# 1. #Folder de trabajo 
setwd("~/pye1pm")
#Nombre:Santiago Arredondo Lopez.
#Carrera:Ing. Civil.
#Matricula:00000218082
#Materia:Probabilidad y estadistica.
#Descrippcion del sitio de estudio: es una pueblo mexicano ubicado en la región
#este del estado de Sonora, en la zona de la Sierra Madre Occidental. El pueblo 
#es la cabecera municipal y localidad más habitada del municipio de Arivechi. 
#Según datos del Censo de Población y Vivienda realizado en 2020 por el Instituto
#Nacional de Estadística y Geografía (INEGI), Arivechi tiene una población de 635 
#habitantes, siendo una de las cabeceras municipales más pequeñas y menos
#pobladas del estado.

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

dist <- c(12.8,14.5,17.7,21.5,26.0,29.8,28.8,27.9,26.7,22.8,16.7,13.0,21.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(%)
##  [12.67,16.16) 3 0.23 23.08  3  23.08
##  [16.16,19.64) 2 0.15 15.38  5  38.46
##  [19.64,23.13) 3 0.23 23.08  8  61.54
##  [23.13,26.61) 1 0.08  7.69  9  69.23
##   [26.61,30.1) 4 0.31 30.77 13 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] 12.8 13.0 14.5 16.7 17.7 21.5 21.5 22.8 26.0 26.7 27.9 28.8 29.8
#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")