#Conversión a Data Frame de los datos del archivo estatura.txt .
estatura <-read.table("C:/estatura.txt")

#convertir data frame en numeric Array para cálculos.
Estatura<-data.matrix(estatura, rownames.force=NA)

#Frecuencia del array Estatura.
freq<-table(Estatura)

#Rango de los datos.
rango<-(max(Estatura) - min(Estatura))

#Media de los datos.
media<-sapply(Estatura,mean)

#Desviación estándar de los datos.
desvest<-sapply(Estatura,sd)

#Varianza de los datos.
varianza<-sapply(Estatura,var)

#cálculo de los cuartiles
quantile(Estatura)

#generación del Histograma de Estatura
histEstatura<-hist.default(Estatura, freq= NULL, axes=TRUE, plot = TRUE,right= TRUE , col ="deepskyblue", breaks=seq(151,191,by=8),
ylab="Frecuencia", xlab="Estatura en cm", main="Histograma de Estatura")

#generación del Diagrama de cajas de Estatura
boxplot.default(Estatura,range =1.5, width =NULL, vardwidth= FALSE, plot =TRUE, col = "azure",horizontal = TRUE,
main = "Diagrama de Cajas de Estatura", xlab="Estatura en cm")

#Salto de intervalos:
salto<-((max(Estatura)- min(Estatura))/5)

#Generación frecuencia absoluta de los datos.
fAbsolute<-table(Estatura)

#Generación frecuencia relativa de los datos.
fRelative<-(table(Estatura)/115)

#Generación frecuencia absoluta acumulada de los datos.
fACummulative<-cumsum(fAbsolute)

#Generación frecuencia relativa acumulada de los datos.
fRCummulative<-(cumsum(fAbsolute)/115)

#Generación del gráfico de Ojiva.
#Requiere paquete agricolae
ojiva.freq(histEstatura, type = "b", col = "red", xlab="Estatura en cm", ylab= "Frecuencia Acumulada", main = "Gráfico de Ojiva")
grid(col="blue")