# 1. #Folder de trabajo 
setwd("~/Documents/pye1pm")
#Nombre:
#Carrera:
#Matrícula:
#Materia:

#Descripción del sitio de estudio: Los Hornos es un ejido del municipio de 
#Cajeme ubicado en el sur del estado mexicano de Sonora, en la zona del 
#valle del Yaqui. Según los datos del Censo de Población y Vivienda realizado 
#en 2020 por el Instituto Nacional de Estadística y Geografía (INEGI), 
#Los Hornos tiene un total de 707 habitantes.

# Script para realizar histogramas y polígonos de frecuencia 
# de datos climáticos, de la estación de Hornos, 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 estación de Hornos, Sonora.
# Estos datos los vamos a guardar en una variable 

dist <- c(17.1, 18.0, 19.3, 22.1, 25.2, 29.3, 31.3, 31.0, 30.3, 26.5, 21.3, 17.5)

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

# 4. Generación de tabla de frecuencias

tabla <- fdt(dist, breaks = "Sturges") 
tabla
##     Class limits f   rf rf(%) cf  cf(%)
##  [16.929,19.866) 4 0.33 33.33  4  33.33
##  [19.866,22.803) 2 0.17 16.67  6  50.00
##  [22.803,25.739) 1 0.08  8.33  7  58.33
##  [25.739,28.676) 1 0.08  8.33  8  66.67
##  [28.676,31.613) 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.1 17.5 18.0 19.3 21.3 22.1 25.2 26.5 29.3 30.3 31.0 31.3
#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. Polígonos de frecuencia 

# Polígono de frecuencia absoluta 
plot(tabla, type="fp")

# Polígono de frecuencia acumulada 
plot(tabla, type="cfp")

# Polígono de frecuencia relativa 
plot(tabla, type="rfp")