Gráficas y funciones

En el siguiente caso analizamos gráficas y sentencias para facilitar el uso de R

reporte <-c(
  rep("descortes",12), 
  rep("lento",42), 
  rep("fria",5),
  rep("estrecho",20),
  rep("atmosfera",10))


reporte
##  [1] "descortes" "descortes" "descortes" "descortes" "descortes" "descortes"
##  [7] "descortes" "descortes" "descortes" "descortes" "descortes" "descortes"
## [13] "lento"     "lento"     "lento"     "lento"     "lento"     "lento"    
## [19] "lento"     "lento"     "lento"     "lento"     "lento"     "lento"    
## [25] "lento"     "lento"     "lento"     "lento"     "lento"     "lento"    
## [31] "lento"     "lento"     "lento"     "lento"     "lento"     "lento"    
## [37] "lento"     "lento"     "lento"     "lento"     "lento"     "lento"    
## [43] "lento"     "lento"     "lento"     "lento"     "lento"     "lento"    
## [49] "lento"     "lento"     "lento"     "lento"     "lento"     "lento"    
## [55] "fria"      "fria"      "fria"      "fria"      "fria"      "estrecho" 
## [61] "estrecho"  "estrecho"  "estrecho"  "estrecho"  "estrecho"  "estrecho" 
## [67] "estrecho"  "estrecho"  "estrecho"  "estrecho"  "estrecho"  "estrecho" 
## [73] "estrecho"  "estrecho"  "estrecho"  "estrecho"  "estrecho"  "estrecho" 
## [79] "estrecho"  "atmosfera" "atmosfera" "atmosfera" "atmosfera" "atmosfera"
## [85] "atmosfera" "atmosfera" "atmosfera" "atmosfera" "atmosfera"
barplot(table(reporte))

data_resumen <- table(reporte)
data_resumen
## reporte
## atmosfera descortes  estrecho      fria     lento 
##        10        12        20         5        42
contaminante <- c(0.367,3.067,0.767,4.133,2.000,1.900,1.467,
                  3.567,4.200,4.900,0.800,4.633,4.233,4.567,3.333,1.933,2.667,
                  3.333,0.600,4.000,1.133,3.200,2.767,3.033,1.633)
contaminante <- c(0.367,3.067,0.767,4.133,2.000,1.900,1.467,
                  3.567,4.200,4.900,0.800,4.633,4.233,4.567,3.333,1.933,2.667,
                  3.333,0.600,4.000,1.133,3.200,2.767,3.033,1.633)


hist(contaminante,
     col = 2, 
     main = "histograma",
     xlab = "rango",
     ylab= "frecuencias")

d <- length(contaminante)
intervalos <- trunc(d^0.5+.999)
R <- diff(range(contaminante))
A = R/intervalos
bb = min(contaminante) + (0:intervalos)*A



hist(contaminante,
     col = 3, 
     main = "histograma",
     xlab = "rango",
     ylab= "frecuencias",
     breaks = bb)