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).
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:
En temperaturas altas y bajas los machos consumen más oxígeno en relación a las hembras.
En temperatura media, las hembras sobrepasan el consumo de Oxígeno en comparación a los machos.
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.
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.