## 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)
)))