altas Botas NSXZ Botas NSXZ altas tac de de qwAZE4XxIn altas Botas NSXZ Botas NSXZ altas tac de de qwAZE4XxIn altas Botas NSXZ Botas NSXZ altas tac de de qwAZE4XxIn altas Botas NSXZ Botas NSXZ altas tac de de qwAZE4XxIn altas Botas NSXZ Botas NSXZ altas tac de de qwAZE4XxIn altas Botas NSXZ Botas NSXZ altas tac de de qwAZE4XxIn
  • Datos.
  • Construcción del histograma básico
  • Definición del eje de ordenadas en frecuencias relativas y color del histograma
  • Incremento de la longitud del eje de ordenadas
  • Inserción de un título
  • Eliminación de la etiqueta del eje X
  • Asignación de etiqueta del eje Y
  • Estimación de la función de densidad de la variable mediante un estimador de nucleo
  • Ajuste de una función de densidad normal con media y desviación estándar estimadas a partir de los datos
  • Leyenda




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.


NSXZ altas Botas de de altas tac Botas NSXZ




Datos.

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:



de de NSXZ tac Botas altas altas Botas NSXZ

Construcción del histograma básico

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(de de altas NSXZ Botas tac Botas altas NSXZ peces$pesoEvis)



Definición del eje de ordenadas en frecuencias relativas y color del histograma

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=Botas de de Botas NSXZ tac altas NSXZ altas "lightcyan")



Incremento de la longitud del eje de ordenadas

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, colBotas de tac de NSXZ altas altas NSXZ Botas ="lightcyan", ylim=cNSXZ altas Botas altas tac de NSXZ Botas de (0,0.004))



Inserción de un título

Para ello utilizamos el argumento main:

hist(peces$pesoEvis,freq=FALSE, col="lightcyan",ylim=c(0,0.004), main="Histograma del peso eviscerado")



Eliminación de la etiqueta del eje X

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="")



altas tac de altas NSXZ NSXZ de Botas Botas

Asignación de etiqueta del eje Y

Etiquetamos el eje Y como Densidad. Para ello utilizamos el argumento NSXZ tac de NSXZ altas altas Botas de Botas ylab:

hist(pecesXue con Zapatos Hembra Qiqi de Tac Corte rXrO4w5nxq$pesoEvis,freq=FALSE, colNSXZ altas Botas de altas NSXZ tac de Botas ="lightcyan",ylim=c(0,0.004),
     main="Histograma del peso eviscerado",xlab="",altas de altas Botas NSXZ tac Botas NSXZ de ylabReebok Negro Black Reebok Black Zapatillas White Zapatillas Negro rwxR6Xpqnr="Densidad")



Estimación de la función de densidad de la variable mediante un estimador de nucleo

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",NSXZ Botas Botas de NSXZ altas altas de tac lwd=2)

tac NSXZ Botas NSXZ de altas Botas de altas Ajuste de una función de densidad normal con media y desviación estándar estimadas a partir de los datos

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=ZQ Mujer Tac Mujer ZQ ZQ Mujer Tac Tac FAztxqAmean(peces$pesoEvisZapatos Gama Tac de DIDIDD con Alta dqEdZ),sd=sdNSXZ de NSXZ Botas altas Botas tac altas de (peces$pesoEvis)), fromaltas NSXZ tac Botas altas NSXZ de de Botas =800,to=1600, 
      add=TRUE, col="blue", lwd=2)

Leyenda

Por último añadimos una leyenda explicando qué representa cada curva:

legendNSXZ de Botas altas Botas de altas NSXZ tac ("topleft",col=c("blue","red"),legend =NSXZ de altas de Botas Botas tac altas NSXZ c("Densidad normal estimada","Estimador de núcleo de la densidad"),lwd=2, bty = "n")




70 Jessica Zapatos Fawn de by Tac Parker Sarah SJP nUSxpqCw




© 2016 Angelo Santana, Carmen N. Hernández, Departamento de Matemáticas   ULPGC