Desdobrando interações

Fazendo o desdobramento dos graus de liberdade para estudar os níveis de B dentro de A, pode-se utilizar os seguintes procedimentos no R.

 ex.b.d.a<-aov(resp~bloco+ a/b +bloco*a,data=ex09)
 summary(ex.b.d.a)
summary(ex.b.d.a,split=list("a:b"=list(B.d.A1=c(1,5,9),B.d.A2=c(2,6,10),B.d.A3=c(3,7,11),B.d.A4=c(4,8,12))))

Lembre-se que aqui interessa apenas concluir sobre o desdobramento.

Para estudar os níveis de A dentro de B, usa-se:

 ex.a.d.b<-aov(resp~bloco+b/a +bloco*a -a,data=ex09)
 summary(ex.a.d.b)
 ex.ab<-summary(ex.a.d.b,split=list("b:a"=list(A.d.B1=c(1,5,9),A.d.B2=c(2,6,10),A.d.B3=c(3,7,11),A.d.B4=c(4,8,12))))

Observe que nesta saída, o teste F não está correto, pois é necessário utilizar o resíduo médio e os graus de liberdade médio para testar os efeitos de A dentro de B.

 
 ex09.av1 <- aov(resp ~ bloco + a*b + bloco*a, data=ex09)
 summary(ex09.av1)

 ex09.anova<-anova(ex09.av1)
 res.medio<-(ex09.anova$Mean[5]+3*ex09.anova$Mean[6])/4
 res.medio

 gl.medio<-27

O teste F pode ser calculado da seguinte maneira:

 FAB1<-468.1/res.medio;FAB1
 FAB2<-137.7/res.medio;FAB2
 FAB3<-108.3/res.medio;FAB3
 FAB4<-430.9/res.medio;FAB4

E os respectivos p-valores associados a cada F calculado:

 pvalorAB1 <- 1 - pf(FAB1, 3, gl.medio);  pvalorAB1
 pvalorAB2 <- 1 - pf(FAB2, 3, gl.medio);  pvalorAB2
 pvalorAB3 <- 1 - pf(FAB3, 3, gl.medio);  pvalorAB3
 pvalorAB4 <- 1 - pf(FAB4, 3, gl.medio);  pvalorAB4