Punto 2 - Modelo de Diseño de Experimentos

load("C:/Users/57318/Downloads/moluscos.RData")
attach(BD_moluscos)
summary(BD_moluscos)
##      c_agua      molusco              cons_o      
##  Min.   : 50   Length:48          Min.   : 1.800  
##  1st Qu.: 50   Class :character   1st Qu.: 6.312  
##  Median : 75   Mode  :character   Median : 9.700  
##  Mean   : 75                      Mean   : 9.305  
##  3rd Qu.:100                      3rd Qu.:11.232  
##  Max.   :100                      Max.   :18.800
library(ggplot2)

BD_moluscos$c_agua=as.character(BD_moluscos$c_agua)

summary(cons_o)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   1.800   6.312   9.700   9.305  11.232  18.800
table(molusco)
## molusco
##  A  B 
## 24 24
table(c_agua)
## c_agua
##  50  75 100 
##  16  16  16
tapply(cons_o, molusco, mean)
##         A         B 
## 10.000417  8.609167
boxplot(cons_o~molusco)

tapply(cons_o, c_agua,mean)
##       50       75      100 
## 12.25062  6.99250  8.67125
tapply(cons_o, c_agua,max)
##   50   75  100 
## 18.8 13.2 14.0
tapply(cons_o, c_agua, min)
##   50   75  100 
## 6.38 1.80 3.68
boxplot(cons_o~c_agua)

boxplot(cons_o~molusco)

mod4=lm(cons_o~molusco+c_agua)
anova(mod4)
## Analysis of Variance Table
## 
## Response: cons_o
##           Df Sum Sq Mean Sq F value   Pr(>F)   
## molusco    1  23.23  23.227  2.0427 0.159846   
## c_agua     1 102.50 102.495  9.0139 0.004362 **
## Residuals 45 511.69  11.371                    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
library(agricolae)

compara6=LSD.test(mod4,"c_agua")
compara6$groups
##       cons_o groups
## 50  12.25062      a
## 100  8.67125      b
## 75   6.99250      b

#Análisis de varianza Valor p= 0 (Asteriscos también) Hay diferencias significativas #Como el analisis tiene asteriscos quiere decir que existe una interaccion