#Segundo Ejercico de Probabilidad y Estadistica ISW
#U1A2
#Histogramas y Tablas de Frecuencia
#Asignar carpeta de trabajo
#Activar paquetes
library(fdth)
## 
## Attaching package: 'fdth'
## The following objects are masked from 'package:stats':
## 
##     sd, var
#Set de Datos Climaticos de Alamos, Sonora.
#Declarar variable
ALAMOS <- c(18.5,19.1,20.2,22.1,24.7,28.1,28.0,27.6,26.8,24.8,21.6,19.1)
dist <- fdt(ALAMOS, breaks = "Sturges") #Calcula la Distribucion de Frecuencias segun Sturge

dist#nos brinda una tabla con los calculos de la distribución de frecuencias.
##     Class limits f   rf rf(%) cf  cf(%)
##  [18.315,20.328) 4 0.33 33.33  4  33.33
##  [20.328,22.341) 2 0.17 16.67  6  50.00
##  [22.341,24.355) 0 0.00  0.00  6  50.00
##  [24.355,26.368) 2 0.17 16.67  8  66.67
##  [26.368,28.381) 4 0.33 33.33 12 100.00
#Donde
#f= frecuencia absoluta
#rf= frecuencia relativa
#rf(%) frecuencia relativa porcentual
#cf= frecuencia acumulada
#cf(%)=frecuencia acumulada porcentual
sort(ALAMOS)
##  [1] 18.5 19.1 19.1 20.2 21.6 22.1 24.7 24.8 26.8 27.6 28.0 28.1
#Generacion de Histograma
#Histograma de Frecuencias Absolutas
hist(ALAMOS, breaks = "Sturges")

#Para mostrar multiples graficas en una ventana

par(mfrow=c(3,2)) # particiona mi ventana grafica en 3x2.
hist(ALAMOS, breaks = "Sturges", col= "purple") #histograma utilizando el numero de clases según Sturge

plot(dist, type="cfh", col="pink")   #histograma de frecuencias acumulada
plot(dist, type="cfp", col="orange")   #poligono de frecuencias acumulado
plot(dist, type="fh", col = "cyan")
plot(dist, type="fh", col="blue") #si desea añadir color a su grafico
plot(dist, type="rfh", col= "red") #Histograma de frecuencia relativa

CANANEA <- c(8.5,10,11.4,15.5,19.3,23.9,23.9,23.4,21.4,17.4,12.4,8.7)
dist2 <- fdt(CANANEA, breaks= "Sturges")

dist2
##     Class limits f   rf rf(%) cf  cf(%)
##    [8.415,11.56) 4 0.33 33.33  4  33.33
##   [11.56,14.705) 1 0.08  8.33  5  41.67
##  [14.705,17.849) 2 0.17 16.67  7  58.33
##  [17.849,20.994) 1 0.08  8.33  8  66.67
##  [20.994,24.139) 4 0.33 33.33 12 100.00
sort(CANANEA)
##  [1]  8.5  8.7 10.0 11.4 12.4 15.5 17.4 19.3 21.4 23.4 23.9 23.9
hist(CANANEA, breaks = "Sturges")

par(mfrow=c(3,2))

hist(CANANEA, breaks = "Sturges", col= "grey")

plot(dist2, type="cfh", col="pink")  
plot(dist2, type="cfp", col="green")  
plot(dist2, type="fh", col = "cyan")
plot(dist2, type="fh", col="yellow") 
plot(dist2, type="rfh", col= "purple")

#Cananea cuenta con una media de temperatura mucho mas baja que la de Alamos
#Y las frecuencias de Cananea que se agarran no son exacta como las de Alamos
#Podemos concluir que Cananea es un lugar menos caliente que Alamos.