El tamaño de la muestra es de 24 Cangrejos, los cuales se sometieron al experimineto para evidenciar el efecto de dos variables (temperatura y sexo) sobre la tasa de consumo de oxigeno. Posteriormente se analiza la relacion entre el consumo de oxigeno en los cangrejos respecto a diferentes temperaturas.

require(ggplot2)
## Loading required package: ggplot2
require(agricolae)
## Loading required package: agricolae
## Warning: package 'agricolae' was built under R version 4.0.3
library(readxl)
## Warning: package 'readxl' was built under R version 4.0.3
Crustaceos <- read_excel("C:/Bioestadistica/Crustaceos.xlsx")
View(Crustaceos)
Crustaceos
## # A tibble: 24 x 3
##    Sexo   Temperatura   CO2
##    <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

Segun la tabla para este experimento se evaluaron tres variables: temperatura (baja, media y alta), sexo (hembra y macho) y consumo de O2 en mg (Consumo de O2).

ggplot(Crustaceos, aes(y=CO2, 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 los machos consumen más oxígeno en relación a las hembras. En temperaturas bajas los machos y las hembras consumen poco oxigeno, sin embargo los machos tienen una tasa de consumo de oxigeno mayor. En temperatura media, las hembras sobrepasan la tasa de consumo de Oxígeno en comparación a los machos.

#Anova

Modelo4= lm(CO2~Temperatura+Sexo, data = Crustaceos)
anova(Modelo4)
## Analysis of Variance Table
## 
## Response: CO2
##             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 nos deja evidenciar que el valor P para la temperatura es significativo, es decir, que esta variable se relaciona directamente con el consumo de oxigeno. Por otro lado la variable sexo no se relaciona de manera directa con respecto al consumo de oxigeno.

#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
##          CO2       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
##          CO2 groups
## Alta  3.0125      a
## Media 2.3875      b
## Baja  1.6375      c
## 
## attr(,"class")
## [1] "group"

La anterior tabla nos deja evidenciar, que en temperaturas altas existe un mayor consumo de oxigeno por parte de los crustaceos en relacion a las otras temperaturas(media y baja), ya que su promedio es de 3.01 clasificandose en el grupo a.

R Markdown

This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.

When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

summary(cars)
##      speed           dist       
##  Min.   : 4.0   Min.   :  2.00  
##  1st Qu.:12.0   1st Qu.: 26.00  
##  Median :15.0   Median : 36.00  
##  Mean   :15.4   Mean   : 42.98  
##  3rd Qu.:19.0   3rd Qu.: 56.00  
##  Max.   :25.0   Max.   :120.00

Including Plots

You can also embed plots, for example:

Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.