Punto 4

Se colectaron 24 cangrejos para este experimento, con el fin de determinar el efecto que tiene dos factores (temperatura y sexo) sobre la tasa de consumo de oxígeno (mgO2). Se realizará un análisis que permita conocer como es el consumo de O2 de los cangrejos en distintas temperaturas.

require(ggplot2)
require(agricolae)
library(readxl)
Cangrejos <- read_excel("C:/Users/beatriz/beatriz/bioestadistica/parcial 3/Cangrejos.xlsx")
View(Cangrejos)
Cangrejos
## # A tibble: 24 x 3
##    Sexo   Temperatura Datos
##    <chr>  <chr>       <dbl>
##  1 macho  Baja          1.9
##  2 macho  Baja          1.8
##  3 macho  Baja          1.6
##  4 macho  Baja          1.4
##  5 hembra Baja          1.8
##  6 hembra Baja          1.7
##  7 hembra Baja          1.4
##  8 hembra Baja          1.5
##  9 macho  Media         2.3
## 10 macho  Media         2.1
## # ... with 14 more rows

Para el experimento se evaluaron tres variables: temperatura (baja, media y alta), sexo (hembra y macho) y consumo de O2 en mg (Datos).

Exploración de datos

ggplot(Cangrejos, aes(y=Datos, x=Temperatura, fill=Sexo))+ facet_grid(~Temperatura)+ geom_boxplot()+theme_bw()+ ylab("Consumo de O2 (mg)")+ggtitle("Diagrama de cajas")

De acuerdo a la gráfica:

Anova

Modelo4= lm(Datos~Temperatura+Sexo, data = Cangrejos)
anova(Modelo4)
## Analysis of Variance Table
## 
## Response: Datos
##             Df Sum Sq Mean Sq F value    Pr(>F)    
## Temperatura  2 7.5833  3.7917 73.4463 6.108e-10 ***
## Sexo         1 0.0037  0.0037  0.0726    0.7903    
## Residuals   20 1.0325  0.0516                      
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

La tabla anterior indica que el valor P para la temperatura es significativa, por lo tanto, está relacionada con el consumo de oxígeno. Mientras que, el sexo no influye siginificativamente con respecto al consumo de oxígeno dentro del experimento.

Postanova

Postanova4=LSD.test(Modelo4, "Temperatura")
Postanova4
## $statistics
##    MSerror Df     Mean       CV  t.value       LSD
##   0.051625 20 2.345833 9.685742 2.085963 0.2369773
## 
## $parameters
##         test p.ajusted      name.t ntr alpha
##   Fisher-LSD      none Temperatura   3  0.05
## 
## $means
##        Datos       std r      LCL      UCL Min Max   Q25  Q50   Q75
## Alta  3.0125 0.2232071 8 2.844932 3.180068 2.7 3.4 2.875 3.00 3.125
## Baja  1.6375 0.1922610 8 1.469932 1.805068 1.4 1.9 1.475 1.65 1.800
## Media 2.3875 0.2474874 8 2.219932 2.555068 2.0 2.7 2.250 2.40 2.600
## 
## $comparison
## NULL
## 
## $groups
##        Datos groups
## Alta  3.0125      a
## Media 2.3875      b
## Baja  1.6375      c
## 
## attr(,"class")
## [1] "group"

La temperatura alta presenta mayor consumo de Oxígeno con respecto a las demas temperaturas, ya que su promedio es de 3.01 y se clasifica en el grupo a. Las temperaturas media y baja presentan diferencias de medias con respecto a la temperatura alta.