#Segundo ejercicio de la materia PyE
#U1A2
#felipe de jesus escarrega valdez 215312

#folder de trabajo
setwd("~/PYE1112ITSON")
#Paquetes
library(fdth)
## 
## Attaching package: 'fdth'
## The following objects are masked from 'package:stats':
## 
##     sd, var
#Datos
#Temperatura media normal (1951-2010) de la estacion 
#SMN EMA CONAGUA de ALAMOS, SONORA
#https://smn.conagua.gob.mx/tools/RESOURCES/Normales5110/NORMAL26002.TXT

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

#Tabla de frecuencias absolutas 
dist <- fdt(ALAMOS,breaks="Sturges")
dist
##     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
#ordenar datos de menor a mayor 
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
#HISTOGRAMA de frecuencias absolutas 

hist(ALAMOS)

#____________________________________________________
#Elavoracion de histograms y poligonos de frecuencia de alamos
#____________________________________________________


#Histogramde frecuencias absolutas
plot(dist, type = "fh")

#Poligono frecuencias absolutas
plot(dist, type = "fp")

#Histogramde frecuencias acumuladas 
plot(dist, type = "cfh")

#Polignono frecuencias acumuladas 
plot(dist, type = "cfp")

#Histogramde frecuencias relativas  
plot(dist, type = "rfh")

#poligono frecuencias relativas  
plot(dist, type = "rfp")

#_______________________________________________________________________________________________________________
                                           #TAREA
#_______________________________________________________________________________________________________________
#guardar las temperaturas en la variable
AGUA_PRIETA <- c(8.0, 9.8, 12.9, 16.6, 21.1, 25.9, 27.2, 26.1, 23.6, 18.2, 12.2, 8.2)

#hacer clases en la forma sturges 
dist2 <- fdt(AGUA_PRIETA,breaks="Sturges")
dist2
##     Class limits f   rf rf(%) cf  cf(%)
##     [7.92,11.83) 3 0.25 25.00  3  25.00
##   [11.83,15.741) 2 0.17 16.67  5  41.67
##  [15.741,19.651) 2 0.17 16.67  7  58.33
##  [19.651,23.562) 1 0.08  8.33  8  66.67
##  [23.562,27.472) 4 0.33 33.33 12 100.00
#ordenar los dator de menor a mayor
sort(AGUA_PRIETA)
##  [1]  8.0  8.2  9.8 12.2 12.9 16.6 18.2 21.1 23.6 25.9 26.1 27.2
#Histograma de agua prieta
hist(AGUA_PRIETA)

#__________________________________________________________________
#Elavoracion de histograms y poligonos de frecuencia de agua prieta
#__________________________________________________________________


#Histogramde frecuencias absolutas
plot(dist2, type = "fh")

#Poligono frecuencias absolutas
plot(dist2, type = "fp")

#Histogramde frecuencias acumuladas 
plot(dist2, type = "cfh")

#Polignono frecuencias acumuladas 
plot(dist2, type = "cfp")

#Histogramde frecuencias relativas  
plot(dist2, type = "rfh")

#poligono frecuencias relativas  
plot(dist2, type = "rfp")

#_____________________________________________________
#                    comparacion
#_____________________________________________________

# podemos observar que en agua prieta el clima es mas frio que en alamos
#ya que los intervalos de clases de agua prieta son mas chicos.

#podemos observar que agua prieta tiene 7 meses mas frios y 5 meses de calor 
#pero los meses de calor no superan a alamos.

#osesa que en agua prieta tiene un año  25% muy frio 42% fresco y un 33% de calor 
#mientras que en alamos 50% fresco 50% calor, lo que quierede decir que en alamos hace frio o hace calor 

#en alamos es mas bruco el cambio de climas  mientras que en agua prieta es mas escalonado el cambio de climas