# 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
library(modeest)
## 
## Attaching package: 'modeest'
## The following object is masked from 'package:fdth':
## 
##     mfv
#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")

#b)Medias de tendencia central

#1) media
mean (dist)
## [1] 21.51538
#2) Mediana
median(dist)
## [1] 21.5
# 3) Moda
mlv(dist, method="mfv")  #mfv significa "most frequent value"
## [1] 21.5
# 4) Cuantiles
summary(dist)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   12.80   16.70   21.50   21.52   26.70   29.80
# 5) Grafico de caja y bigote
boxplot(dist)

# c) Medidas de dispersion

# 1 ) Varianza
var(dist)
## [1] 37.21141
# Si la varianza le sacamos raiz cuadrada, como resultado
#vamos a optener la:

# 2) Desviacion estandar
sd(dist)
## [1] 6.100116
# Grafico de dispersion lineal
plot(dist)