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

# Script para realizar histogramas y poligonos de frecuencia 
# de datos climáticos, de la estación de CIUDAD OBREGON (DGE), Sonora.

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

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

dist <- c( 17.4, 18.6, 20.3, 23.5, 26.9, 30.6, 31.9, 31.6, 31.0, 27.5, 22.1, 17.9)

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

# 4. Generaciacion de tabla de frecuencias

tabla <- fdt(dist, breaks = "Sturges") 
tabla
##     Class limits f   rf rf(%) cf  cf(%)
##  [17.226,20.225) 3 0.25 25.00  3  25.00
##  [20.225,23.223) 2 0.17 16.67  5  41.67
##  [23.223,26.222) 1 0.08  8.33  6  50.00
##   [26.222,29.22) 2 0.17 16.67  8  66.67
##   [29.22,32.219) 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

# 5. Ordenar los numeros (datos) de menor a mayor
sort(dist)
##  [1] 17.4 17.9 18.6 20.3 22.1 23.5 26.9 27.5 30.6 31.0 31.6 31.9
#6. Histogramas de frecuencia 

# Histograma de frecuencia absoluta 
plot(tabla, type="fh")

# Histograma de frecuencia acumulada 
plot(tabla, type="cfh")

# Histograma de frecuencia relativa 
plot(tabla, type="rfh")

#7. Poligonos de frecuencia 

# Poligono de frecuencia absoluta 
plot(tabla, type="fp")

# Poligono de frecuencia acumulada 
plot(tabla, type="cfp")

# Poligono de frecuencia relativa 
plot(tabla, type="rfp")