#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.

#A) Análisis de frecuencia, tablas, histrogramas y polígonos


#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
library(modeest)
## 
## Attaching package: 'modeest'
## The following object is masked from 'package:fdth':
## 
##     mfv
#Lo que hacemos cuando escribimos "library" es abrir el "libro" que acabamos
#de comprar

#Datos normales mensuales de la estación de Yécora
#Estos datos los vamos a guardar en una variable

#Se agregó un dato (8.9) para utilizar la moda
dist <- c(6.3, 6.7, 8.9, 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) 6 0.46 46.15  6  46.15
##    [9.07,11.9) 1 0.08  7.69  7  53.85
##   [11.9,14.74) 2 0.15 15.38  9  69.23
##  [14.74,17.57) 0 0.00  0.00  9  69.23
##   [17.57,20.4) 4 0.31 30.77 13 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  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")

#B)Medidas de tendencia central 

#Media
mean(dist)
## [1] 12.48462
#Mediana
median(dist)
## [1] 11.9
#Moda
#Valor más frecuente
mlv(dist, method="mfv") 
## [1] 8.9
#mfv "most frequent value"

#Cuantiles
#Valor mínimo (MIn.), 1st Qu.(Primer cuantil)
# Median (Mediana), Mean (promedio), 3rd Qu. (tercer cuantil)
# Max. (Valor máximo)
summary(dist)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    6.30    8.20   11.90   12.48   18.00   20.20
#Gráfica de caja y bigote
boxplot(dist)

#C)Medidas de de dispersión

#Varianza
var(dist)
## [1] 29.01141
#Si a la varianza le sacamos raíz cuadrada obtenemos la desviación estándar

#Desviación estándar
sd(dist)
## [1] 5.386224
#Grafico de dispersión lineal
plot(dist)