Biochemical Oxygen Demand

La demanda bioquímica de oxígeno es un parámetro que mide la cantidad de dioxígeno consumido al degradar la materia orgánica de una muestra líquida. Es la materia susceptible de ser consumida u oxidada por medios biológicos que contiene una muestra líquida, disuelta o en suspensión.

OBD

Paquetes

library(pacman)
p_load("datasets","DT","fdth")

Datos

data("BOD")
dim(BOD)
## [1] 6 2

Tabla de datos

datatable(BOD)

Demanda de oxigeno

table(BOD)
##     demand
## Time 8.3 10.3 15.6 16 19 19.8
##    1   1    0    0  0  0    0
##    2   0    1    0  0  0    0
##    3   0    0    0  0  1    0
##    4   0    0    0  1  0    0
##    5   0    0    1  0  0    0
##    7   0    0    0  0  0    1

Medidas de tendencia central

summary(BOD)
##       Time           demand     
##  Min.   :1.000   Min.   : 8.30  
##  1st Qu.:2.250   1st Qu.:11.62  
##  Median :3.500   Median :15.80  
##  Mean   :3.667   Mean   :14.83  
##  3rd Qu.:4.750   3rd Qu.:18.25  
##  Max.   :7.000   Max.   :19.80

Ilustrando las MTCcongrafica de caja y bigotes

boxplot(BOD)
#Para dibujar una linea roja de la media del largo del petalo, la linea representa la medida de tendencia central 
abline(h=mean(iris$Petal.Length), col="red")

## Contando la demanda de Oxigeno

table(BOD$demand)
## 
##  8.3 10.3 15.6   16   19 19.8 
##    1    1    1    1    1    1

Varianza

var(BOD$demand)
## [1] 21.44267

Desviación estandar

sd(BOD$demand)
## [1] 4.630623

Tabla de frecuencias, histogramas y poligonos

tabla <- fdt(BOD$demand, breaks = "Sturges")
tabla
##  Class limits f   rf rf(%) cf  cf(%)
##      [8.2,11) 2 0.33 33.33  2  33.33
##       [11,14) 0 0.00  0.00  2  33.33
##       [14,17) 2 0.33 33.33  4  66.67
##       [17,20) 2 0.33 33.33  6 100.00

#Donde #f= frecuencia absoluta #rf= frecuencia relativa #rf(%) frecuencia relativa porcentual #cf= frecuencia acumulada #cf(%)=frecuencia acumulada porcentual

Histogramas

Histograma de frecuencia absluta

plot(tabla,type= "fh")

## Histograma de frecuencia relativa

plot(tabla,type= "rfh")

## Histograma de frecuencia Acomulada

plot(tabla,type= "cfh")

# POLIGONOS Poligonos de frecuencia absoluta

plot(tabla,type= "fp")

## Poligonos de frecuencia relativa

plot(tabla,type= "rfp")

## Poligonos de frecuencia Acomulada

plot(tabla,type= "cfp")