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

#Nombre: Víctor Hugo Robles Lizárraga
#Carrera: Ingenieria Industrial y de Sistemas
#Matrícula: 229287
#Materia: Probabilidad y estadística

#Descripción del sitio de estudio: Yécora se encuentra a más de 
#1500 metros de altura y se ubica a tan solo 11 kilómetros de la 
#frontera con Chihuahua. Tiene aproximadamente 6089 habitantes.

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

dist <- c(6.3, 6.7, 8.9, 11.9, 14.7, 19.4, 20.2, 19.7, 18.0, 13.0, 8.2, 6.4)

#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(%)
##   [6.237,9.07) 5 0.42 41.67  5  41.67
##    [9.07,11.9) 1 0.08  8.33  6  50.00
##   [11.9,14.74) 2 0.17 16.67  8  66.67
##  [14.74,17.57) 0 0.00  0.00  8  66.67
##   [17.57,20.4) 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]  6.3  6.4  6.7  8.2  8.9 11.9 13.0 14.7 18.0 19.4 19.7 20.2
#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")

#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")