#Distribucion de frecuencia de datos de alamos 
ALAMOS <- c(18.5,19.1,20.2,22.1,24.7,28.1,28.0,27.6,26.8,24.8,21.6,19.1)
#Debera instalar el paquete que vamos a llamar
library(fdth) # carga la extensión fdth (para el calculo de distribución de frecuencias.
## 
## Attaching package: 'fdth'
## The following objects are masked from 'package:stats':
## 
##     sd, var
dist <- fdt(ALAMOS,breaks="Sturges") # calcula la distribución de frecuencias utilizando la regla Sturge

dist #nos brinda una tabla con los calculos de la distribución de frecuencias.
##     Class limits f   rf rf(%) cf  cf(%)
##  [18.315,20.328) 4 0.33 33.33  4  33.33
##  [20.328,22.341) 2 0.17 16.67  6  50.00
##  [22.341,24.355) 0 0.00  0.00  6  50.00
##  [24.355,26.368) 2 0.17 16.67  8  66.67
##  [26.368,28.381) 4 0.33 33.33 12 100.00
#Donde
#f= frecuencia absoluta
#rf= frecuencia relativa
#rf(%) frecuencia relativa porcentual
#cf= frecuencia acumulada
#cf(%)=frecuencia acumulada porcentual

par(mfrow=c(3,2)) # particiona mi ventana grafica en 3x2.
hist(ALAMOS, breaks = "Sturges") #histograma utilizando el numero de clases según Sturge

plot(dist, type="cfh")   #histograma de frecuencias acumulada
plot(dist, type="cfp")   #poligono de frecuencias acumulado
plot(dist, type="fh")
plot(dist, type="fh", col="blue") #si desea añadir color a su grafico


#Distribucion de frecuencia de datos de Yecora
YECORA <- c(7.2,8.4,10.2,13.4,17.3,21.5,21.6,21.0,19.4,14.9,9.9,7.2,14.3)
library(fdth)
dist <- fdt(YECORA,breaks="Sturges") 
dist
##     Class limits f   rf rf(%) cf  cf(%)
##   [7.128,10.066) 4 0.31 30.77  4  30.77
##  [10.066,13.003) 1 0.08  7.69  5  38.46
##  [13.003,15.941) 3 0.23 23.08  8  61.54
##  [15.941,18.878) 1 0.08  7.69  9  69.23
##  [18.878,21.816) 4 0.31 30.77 13 100.00
par(mfrow=c(3,2))

hist(YECORA, breaks = "Sturges")

plot(dist, type="cfh")
plot(dist, type="cfp")
plot(dist, type="fh")
plot(dist, type="fh", col="red")