#Folder de trabajo
setwd("~/pye1pm")

#Script oara realizar histogramas y polígonos de frecuencia
#de datos climáticos de la estación Álamos, Sonora.

#para activar la biblioteca "fdth": frecuencias, distribuciones
#tablas e histogramas

#Paquete para realizar frecuencias, distribuciones, tablas,
#histogramas
library(fdth)
## 
## Attaching package: 'fdth'
## The following objects are masked from 'package:stats':
## 
##     sd, var
#Datos normales mensuales de la estación de Álamos
#Estos datos los vamos a guardar en una variable

dist <- 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)

#Los datos corresponden en orden secuencial a los 12 meses del año (De enero a febrero)

#Generación de tabla de frecuencias

tabla <- fdt(dist, breaks = "Sturges")
tabla
##     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= frecuecia absoluta
#rf= frecuencia relativa
#rf(%)= frecuencia relativa porcentual
#cf= frecuencia acumulada
#cf(%)= frecuencia acumulada porcentual

#Sort sirve para ordenar los números (datos) de menor a mayor
sort (dist)
##  [1] 18.5 19.1 19.1 20.2 21.6 22.1 24.7 24.8 26.8 27.6 28.0 28.1
#Histograma de frecuencia acumulada

plot(tabla, type="cfh")