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