U1A2.R

# Segundo ejercicio de PyE ISW
# U1A2
# Isaac Castelo Valenzuela 
# 00000216823

# Histogramas y tablas de frecuencia

# Asignamos el folder de trabajo
setwd("~/PyE Programas")

# 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)
BACANORA <- c(12.6,14.7,16.9,21.4,24.7,30.5,29.7,28.7,27.1,22.6,16.5,12.7,21.5)
BATACOSA <- c(16.0,17.3,19.2,22.1,25.4,29.6,30.2,29.2,28.5,25.2,20.1,16.4,23.3)

#Calcula la distribucion de frecuencias segun la regla de Sturges
dist <- fdt(ALAMOS, breaks="Sturges")
dist <- fdt(BACANORA, breaks="Sturges")
dist <- fdt(BATACOSA, breaks="Sturges")

dist

##   Class limits f   rf rf(%) cf  cf(%)
##  [15.84,18.77) 3 0.23 23.08  3  23.08
##   [18.77,21.7) 2 0.15 15.38  5  38.46
##   [21.7,24.64) 2 0.15 15.38  7  53.85
##  [24.64,27.57) 2 0.15 15.38  9  69.23
##   [27.57,30.5) 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(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

#sort(BACANORA)
#sort(BATACOSA)

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

#####################  BACANORA  #########################
#hist(BACANORA, breaks = "Sturges")


#par(mfrow=c(3,2))
#hist(BACANORA, breaks = "Sturges")

#plot(dist, type="cfh")  
#plot(dist, type="cfp")  
#plot(dist, type="fh")
#plot(dist, type="fh", col="red") 
#plot(dist, type="rfh") 

#####################  BATACOSA  #########################

#hist(BATACOSA, breaks = "Sturges")


#par(mfrow=c(3,2))
#hist(BATACOSA, breaks = "Sturges")

#plot(dist, type="cfh")  
#plot(dist, type="cfp")  
#plot(dist, type="fh")
#plot(dist, type="fh", col="purple") 
#plot(dist, type="rfh")