Pessoal, recorro novamente a vocês para me ajudarem em um problema que já foi objeto de outra mensagem que postei aqui.

Estou com um modelo de ANOVA com medidas repetidas em que um dos tratamentos é aninhado a outro. O que gostaria de fazer é particionar o tratamento aninhado como se faz com o comando aov (..., split..) (Vi exemplo no site R-idiculas (excelente, por sinal)).
A estrutura dos meus dados é a seguinte:

Id: Unidade Amostral
#Estrutura dos dados:
set.seed(13032012)
d=expand.grid(Trat1=gl(3,18),Tempo=gl(2,1))
d$Trat2=gl(2,9,108)
d$Id=as.factor(rep(1:54,2))
d$Y=abs(rnorm(108,c(10,30,8),10))

#Verificar estrutura para a unidade amostral i=3, por exemplo
d[d$Id==3,]: Tratamento1 =1, Tratamento2=1 Tempo=1, 2.

#Ajuste que fiz usando lme

library(nlme)
 f=lme(Y~(Trat1/Trat2)*Tempo, random=~1|Id,d)

#Gostaria de realizar algo do tipo, só que com o modelo misto usando o pacote nlme
m1=aov(Y~(Trat1/Trat2)*Tempo,d)
names(coef(m1))[5:7]
grep("Trat11", names(coef(m1))[5:7])
grep("Trat12", names(coef(m1))[5:7])
grep("Trat13", names(coef(m1))[5:7])

### As partições, mas sem considerar o modelo misto (queria algo semelhante).

summary(m1)
summary(m1, split=list("Trat1:Trat2"=list(
                         " Trat1-1"=c(1),
                         "Trat2-2"=c(2),
                         "Trat3-3"=c(2)

                         )))

Fico na dúvida, se devo fazer isso na interação tripla, já que nos meus dados reais foi significativa ou se os coeficientes estimados delas também devam entrar no contraste, enfim, eu sempre tenho muitas dúvidas em contrastes no R.

Alguém sabe como posso fazer isso? Já tentei os contrastes com multcomp, mas não é isso que quero... gostaria de ter partição da soma dos quadrados, semelhante ao exemplo com o aov.

Alguém pode me ajudar?

Desde já, muito obrigado.

Olympio