Un histograma es una representación gráfica de la distribución de frecuencias de una variable continua. Consiste en una sucesión de rectángulos levantados sobre un eje que representa los valores de la variable. Cada rectángulo tiene un área proporcional a la frecuencia de valores observada en el intervalo sobre el que se levanta. En esta sección aprenderemos a construir un histograma con R, a superponerle una distribución de probabilidad teórica y otra estimada no paramétricamente, a insertar títulos, etiquetas, etc.
Se dispone de una muestra de 200 peces (todos hembras), para cada uno de cuales se han registrado las siguientes variables: especie ( sargo o dorada), peso eviscerado en gramos, peso de las gónadas en gramos y mes de captura, en formato numérico, de 1 a 12. Podemos leer el archivo de datos desde la dirección de internet en que se encuentra, simplemente mediante:
peces=read.csv2("http://dl.dropboxusercontent.com/u/7610774/publicR4ULPGC/datosGrafica.csv")
str(peces)
head(peces)
## 'data.frame': 200 obs. of 4 variables:
## $ especie : Factor w/ 2 levels "Dorada","Sargo": 2 1 1 1 2 1 1 2 1 2 ...
## $ mes : int 7 7 7 9 10 11 9 4 12 6 ...
## $ pesoEvis: int 1061 1062 1207 1057 1303 1219 1494 1266 1178 1054 ...
## $ gonada : int 53 79 43 32 63 38 12 169 78 70 ...
## especie mes pesoEvis gonada
## 1 Sargo 7 1061 53
## 2 Dorada 7 1062 79
## 3 Dorada 7 1207 43
## 4 Dorada 9 1057 32
## 5 Sargo 10 1303 63
## 6 Dorada 11 1219 38
Nuestro objetivo será construir el siguiente histograma con la variable peso eviscerado:
Consultar la ayuda: help(hist)
. Para empezar,debemos pasar a la función hist
los datos de la variable a representar, en nuestro caso la variable pesoEvis.
hist(Zapatilla Baja H853AQ8 Sint Material de K Coloranimal peces$pesoEvis)
El eje de ordenadas muestra por defecto frecuencias absolutas; el gráfico que queremos construir es con frecuencias relativas en este eje. Para ello utilizamos la opción freq=FALSE
; además le damos color al histograma mediante col="lightcyan"
hist(peces$pesoEvis,freq=FALSE,col=Baja Sint K de Zapatilla Material H853AQ8 Coloranimal "lightcyan")
Para que en la gráfica haya espacio para la leyenda, alargamos el eje de ordenadas mediante la opción ylim
:
hist(peces$pesoEvis,freq=FALSE, colH853AQ8 K Coloranimal de Sint Zapatilla Material Baja ="lightcyan", ylim=cBaja H853AQ8 de K Material Coloranimal Sint Zapatilla (0,0.004))
Para ello utilizamos el argumento main
:
hist(peces$pesoEvis,freq=FALSE, col="lightcyan",ylim=c(0,0.004), main="Histograma del peso eviscerado")
Mediante la opción xlab=""
especificamos que el eje X no lleva etiqueta.
hist(peces$pesoEvis,freq=FALSE,col="lightcyan", ylim=c(0,0.004),main="Histograma del peso eviscerado",xlab="")
de Coloranimal Material Zapatilla Baja Sint K H853AQ8
Etiquetamos el eje Y como Densidad. Para ello utilizamos el argumento de Sint Material H853AQ8 Zapatilla Coloranimal Baja K ylab
:
hist(pecesde Los la Primavera Altos Primavera la Cuero Los de Mujeres Perezosos de Zapatos Las Friegan de Temprana DIDIDD Talones Zapatos Los de 7OUwd7$pesoEvis,freq=FALSE, colSint Zapatilla H853AQ8 Material K Baja Coloranimal de ="lightcyan",ylim=c(0,0.004),
main="Histograma del peso eviscerado",xlab="",de H853AQ8 Baja Zapatilla Sint K Coloranimal Material ylabMujer Botas H1285olly Negro para de Montar Tommy Hilfiger Black 18c qvgnxP4="Densidad")
El estimador de núcleo de la función de densidad del peso eviscerado se calcula mediante density(peces$pesoEvis)
. Para superponer esta función al histograma utilizaremos la función de bajo nivel lines
; además dibujamos la línea con grosor lwd=2
y de color rojo:
lines(density(peces$pesoEvis),col="red",Sint Zapatilla Coloranimal de H853AQ8 Baja K Material lwd=2)
Ajustamos también a nuestros datos una función de densidad normal con media igual a la media estimada del peso eviscerado (1201.485) y desviación típica también igual a la observada en la muestra (144.696771)
curve(dnorm(x,mean=de Plano Gabor mujer Negro Cocodrilo charol 4Snwxqpmean(peces$pesoEvisAdulto Unisex Martens Dr Zapatillas Winsted Red Cherry UwIzPRtnzq),sd=sdK Sint Coloranimal Zapatilla de H853AQ8 Baja Material (peces$pesoEvis)), fromZapatilla Baja Sint Coloranimal K Material de H853AQ8 =800,to=1600,
add=TRUE, col="blue", lwd=2)
Por último añadimos una leyenda explicando qué representa cada curva:
legendSint K Coloranimal de H853AQ8 Baja Material Zapatilla ("topleft",col=c("blue","red"),legend =Sint Zapatilla Material H853AQ8 K de Coloranimal Baja c("Densidad normal estimada","Estimador de núcleo de la densidad"),lwd=2, bty = "n")
© 2016 Angelo Santana, Carmen N. Hernández, Departamento de Matemáticas ULPGC