WALMES, muito obrigado!!
Agora a saída do desdobramento está correta.
#### adaptado de Walmes ---##### Fatiar solo/nc.
m1 <- aov(mspa~bloco+np*solo + solo/nc + np:solo:nc, data=a1)
anova(m1)
m1$assign
names(coef(m1))[m1$assign==5] # ir mudando o nº ate achar a inter solo:nc
# buscando as posiçoes por expressões regulares,
grep("ES", names(coef(m1))[m1$assign==5])
grep("EU_0-10", names(coef(m1))[m1$assign==5])
grep("EU_40-60", names(coef(m1))[m1$assign==5])
grep("SE_0-10", names(coef(m1))[m1$assign==5])
grep("SE_40-60", names(coef(m1))[m1$assign==5])
grep("TM", names(coef(m1))[m1$assign==5])
# de posse das posiçoes dos solos informar no summary
summary(m1, split=list("solo:nc"=list(
"ES"= c(1),
"EU_0-10"=c(2),
"EU_40-60"=c(3),
"SE_0-10"=c(4),
"SE_40-60"= c(5),
"TM"= c(6)
)))
Odirley R. Campos
Engenheiro Agrônomo UFV/MG
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
