[R-br] Ajuda com particionamento de soma de quadrados usando nlme

[Olympio Neto]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120324/681a8e4c/attachment.html>