# Segundo ejecicio de PYE de ing. software
#U1A2
#Histogramas y tablas de frecuencia
#asignar folder de trabajo
setwd("~/probabilidad y estadistica")
#activar paquetes
library(fdth)
##
## Attaching package: 'fdth'
## The following objects are masked from 'package:stats':
##
## sd, var
#set de datos climaticos de alamos sonora
# declaracion de 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 la regla de 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
#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
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
#generacion de hsitograma
#histograma de frecuencias absolutas
hist(ALAMOS, breaks = "Sturges")
#para mostar 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", col="yellow")# histograma de frecuencias acumulada
plot(dist, type="cfp", col="green")# poligono de frecuencias acumulado
plot(dist, type="fh")
plot(dist, type="fh", col="blue")#si se añade un color al grafico
plot(dist, type="rfh")
#Aplicamos lo mismo pero ahora con los datos de banamichi
banamichi =c(13.8,15.2,17.3,20.8,24.9,29.6,29.5,28.4 ,27.4 ,23.2 ,17.7 ,13.9)
dist1 = fdt(banamichi, breaks = "Sturges")
hist(banamichi, breaks = "Sturges")
par(mfrow=c(3,2))
hist(banamichi, breaks = "Sturges")
#colocaremos los histogramas y el ploligono
plot(dist1, type="cfh", col="yellow")
plot(dist1, type="cfp", col="green")
plot(dist1, type="fh", col="blue")
plot(dist1, type="rfh")
sort(banamichi)
## [1] 13.8 13.9 15.2 17.3 17.7 20.8 23.2 24.9 27.4 28.4 29.5 29.6
#como se puede observar banimichi suele ser mas frio en enero y mas caliente en junio que alamos,
#ademas de tener una media de temperatura mas baja a comparacion de alamos.
