5.3.- Código R de la entrada

Incluyo aquí el código R de esta entrada para poder practicar con él :

#    Cargo el paquete para el calculo de la moda
#
library(modeest)
#    Cargamos los datos del nivel del lago Huron entre 1875-1972
#
data ("LakeHuron")
#    Las medidas de centralización de los datos
#
paste (' La media es   : ', format(mean(LakeHuron), digits=5))
paste (' La mediana es : ', format(median(LakeHuron), digits=5))
paste (' La moda es    : ', format(mfv(LakeHuron), digits=5))
#    La medida de dispersión de los datos
#
paste(' La desviación standard es : ',format(sd(LakeHuron), digits = 4) )
paste(' La varianza es : ', format(sd(LakeHuron)^2, digits = 4) )
#    Salidas gráficas
#
hist(LakeHuron, main='Nivel de Lago Huron', xlab = 'Nivel',
     ylab = 'DENSIDAD',
     xlim = c(574, 585),
     prob=TRUE)
#    presentamos la densidad de los datos para comprobar la normalidad
lines(density(LakeHuron), col='blue', lwd=2)
#    Comprobación de normalidad
shapiro.test(LakeHuron)
#    añadimos la distribución de datos
rug (LakeHuron, col='red')
#    las medidas de centralización
abline(v=mean(LakeHuron), col='blue')
abline(v=median(LakeHuron), col='green')
abline(v=mfv(LakeHuron), col='red')
#    La media +/- 3 veces la desviación estandard
abline(v=mean(LakeHuron)+(3*sd(LakeHuron)), col='red', lty=2, lwd=2)
abline(v=mean(LakeHuron)-(3*sd(LakeHuron)), col='red', lty=2, lwd=2)
#    Salida grafico cartesiano
plot(x=seq(from=1, to=length(LakeHuron)), y=LakeHuron,
     main = 'Nivel del Lago Huron 1.875 a 1.972',
     xlab = 'Medida',
     ylab =' Valor',
     type='b', ylim = c(574, 585))
#    La media +/- 3 veces la desviación estandard
abline(h=mean(LakeHuron)+(3*sd(LakeHuron)), col='red', lty=2, lwd=2)
abline(h=mean(LakeHuron)-(3*sd(LakeHuron)), col='red', lty=2, lwd=2)
#    La media
abline(h=mean(LakeHuron), col='blue')