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