#SEGUNDO EJERCICIO DE PYE ISW.
#U1A2.
#HISTOGRAMAS Y TABLAS DE FRECUENCIA.
#Asignar folder de trabajo.
setwd("~/PyEitson")
#activar paquetes
library (fdth)
##
## Attaching package: 'fdth'
## The following objects are masked from 'package:stats':
##
## sd, var
#set climáticos 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, breacks="Sturges") #calcula la distribución de frecuencias según la regla de sturge
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
#nos brinda una tabla con los calculos de la distribución de frecuencias.
#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, breacks="Sturges")
## Warning in plot.window(xlim, ylim, "", ...): "breacks" is not a graphical
## parameter
## Warning in title(main = main, sub = sub, xlab = xlab, ylab = ylab, ...):
## "breacks" is not a graphical parameter
## Warning in axis(1, ...): "breacks" is not a graphical parameter
## Warning in axis(2, ...): "breacks" is not a graphical parameter
#Para mostrar multiples graficas en una ventana
par(mfrow=c(3,2)) # particiona mi ventana grafica en 3x2.
hist(ALAMOS, breaks = "Sturges") #histograma utilizando el numero de clases según Sturge
plot(dist, type="cfh") #histograma de frecuencias acumulada
plot(dist, type="cfp") #poligono de frecuencias acumulado
plot(dist, type="fh")
plot(dist, type="fh", col="blue") #si desea añadir color a su grafico
plot(dist, type="rfh") #histograma de frecuencias relativas
##########################################################################################################
# A C T I V I D A D
#set climáticos de Agua prieta, sonora.
#https://smn.conagua.gob.mx/es/informacion-climatologica-por-estado?estado=son
#declarar variable
APRIETA <- 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, 17.5)
dist <- fdt(APRIETA, breacks="Sturges") #calcula la distribución de frecuencias según la regla de sturge
dist
## Class limits f rf rf(%) cf cf(%)
## [7.92,11.83) 3 0.23 23.08 3 23.08
## [11.83,15.741) 2 0.15 15.38 5 38.46
## [15.741,19.651) 3 0.23 23.08 8 61.54
## [19.651,23.562) 1 0.08 7.69 9 69.23
## [23.562,27.472) 4 0.31 30.77 13 100.00
#nos brinda una tabla con los calculos de la distribución de frecuencias.
#Donde
#f= frecuencia absoluta
#rf= frecuencia relativa
#rf(%) frecuencia relativa porcentual
#cf= frecuencia acumulada
#cf(%)=frecuencia acumulada porcentual
sort(APRIETA)
## [1] 8.0 8.2 9.8 12.2 12.9 16.6 17.5 18.2 21.1 23.6 25.9 26.1 27.2
#GENERACION DE HISTOGRAMA
#histograma de frecuencias absolutas
hist (APRIETA, breacks="Sturges")
## Warning in plot.window(xlim, ylim, "", ...): "breacks" is not a graphical
## parameter
## Warning in title(main = main, sub = sub, xlab = xlab, ylab = ylab, ...):
## "breacks" is not a graphical parameter
## Warning in axis(1, ...): "breacks" is not a graphical parameter
## Warning in axis(2, ...): "breacks" is not a graphical parameter
#Para mostrar multiples graficas en una ventana
par(mfrow=c(3,2)) # particiona mi ventana grafica en 3x2.
hist(APRIETA, breaks = "Sturges") #histograma utilizando el numero de clases según Sturge
plot(dist, type="cfh") #histograma de frecuencias acumulada
plot(dist, type="cfp") #poligono de frecuencias acumulado
plot(dist, type="fh")
plot(dist, type="fh", col="green") #si desea añadir color a su grafico
plot(dist, type="rfh") #histograma de frecuencias relativas
#EN CONCLUSIÓN EN AGUA PRIETA LAS PERSONAS ESTÁN EXPUESTAS
#A BAJAS TEMPERATURAS QUE OCILAN ENTRE LOS 8 A 27 GRADOS, DONDE LA
#MINIMA EN ALAMOS ES DE 18.5 Y LA MÍNIMA EN AGUA PRIETA ES DE 8.0
#AQUÍ PODEMOS NOTAR LA CLARA DIFERENCIA QUE HAY EN LAS TEMPERATURAS DE AMBOS
#LUGARES.