26-05

set.seed(123)
porc_ger= c(
  rnorm(40,60,6),
  rnorm(40,70,7),
  rnorm(40,80,8))

acido=gl(3,40,120,c("co","c1","c2"))

datos = data.frame(acido, porc_ger)
head (datos)

table(datos$acido)
datos_des=datos [-c(50,111,120),]
table(dato_des$acido)

mod1=aov(porc_ger~acido,datos)
summary(mod1)

boxplot(datos$porc_ger~datos$acido)

##analisis de varianza balanceado(con los datos desbalanceados)

mod2=aov(porc_ger~acido,datos_des)
summary(mod2)

mod3= lm(porc_ger~acido,datos_des)
library(car)
mod3_res=Anova(mod3,type='II')
mod3_res

##ahora con bloques

set.seed(123)
porc_ger= c(
  rnorm(40,60,6),
  rnorm(40,70,7),
  rnorm(40,80,8))


bloq=gl(3,40,120,c("c0","c1","c2"))
acido=gl(4,10,120,c("c0","c1","c2","c3"))
datos=data.frame(acido, bloq, porc_ger)
datos_des= datos[-sample(120,5),]
datos_des

table(datos_des$bloq, datos_des$acido)

mod1=aov(porc_ger~bloq*acido,datos_des)
summary(mod1)

mod2=lm(porc_ger~bloq*acido,datos_des)

mod2_res=Anova(mod2,tyoe="II")


mod2_res

##ajustando desbalanceado con diferente orden

mod3=lm(porc_ger~bloq+acido+bloq:acido,datos_des)
Anova(mod3,tyoe="II")

mod3=lm(porc_ger~bloq+acido+bloq:acido,datos_des)
Anova(mod3,tyoe="II")

#con covariable

diseno factorial simple en bloques completos generalizados y al azar, desbalanceado y con la tecnica analisis de covarianza

set.seed(123)
porc_ger= c(
  rnorm(40,60,6),
  rnorm(40,70,7),
  rnorm(40,80,8))

diam_med=sort(rnorm(120,12,1.3))

bloq=gl(3,40,120,c("c0","c1","c2"))
acido=gl(4,10,120,c("c0","c1","c2","c3"))

datos=data.frame(acido, bloq, porc_ger, diam_med)
datos_des= datos[-sample(120,5),]

table(datos_des$bloq,datos_des$acido)

mod1=lm(porc_ger~bloq+acido+bloq:acido,datos_des)
Anova(mod3,tyoe="II")