[R-br] multcomp lmer contraste interação..

[Eder David Borges da Silva]

Pessoal da lista:
Minha duvida:
require(lme4)
require(multcomp)
## Exemplo no lme4 moodificado
fm2 <- lmer(Reaction ~ Days/Subject + (0+Days|Subject), sleepstudy)
summary(fm2)
# Tenho interesse neste contraste por exemplo:
summary(glht(fm2,c('Days:Subject3-Days:Subject1 = 0')) )
## Funcionou...
## Agora com meus dados querro coisa semelhante

### Agora com meu conjunto de dados:
DADOSS <-
expand.grid(sitio=1,ponto=as.character(c('SS')),serie=2,camara=as.character(LETTERS[1:5]),tempo=c(0,15,30,45),armazen=c('VIDRO','HEAD','PP'))
DADOSS <-
rbind(DADOSS,expand.grid(sitio=1,ponto=as.character(c('SL')),serie=2,camara=as.character(LETTERS[6:10]),tempo=c(0,15,30,45),armazen=c('VIDRO','HEAD','PP')))
DADOSS <-
rbind(DADOSS,expand.grid(sitio=1,ponto=as.character(c('SL','SS')),serie=2,camara='Z',tempo=1:2,armazen=c('VIDRO','HEAD','PP')))
DADOSS$ID <-
paste(DADOSS[,1],DADOSS[,2],DADOSS[,3],DADOSS[,4],DADOSS[,5],DADOSS[,6],sep='',collapse
= NULL)
DADOSS$camaraM <- paste(DADOSS$sitio,DADOSS$ponto,DADOSS$camara,sep='')
DADOSS$resp <- rnorm(nrow(DADOSS))

model <- lmer(resp~armazen/tempo+(1|camaraM)+(1|ponto),DADOSS)
summary(model)
anova(model)
### Eu tenho interesse apenas em este contraste (e nos outros que envolvem a
interação ):
summary(glht(model,c('armazenVIDRO:tempo-armazenPP:tempo = 0')))
## PQ não funciona????????????(Apesar de ter pouco cabelo, já arranquei
alguns para tentar descobrir....)
# OBS: tempo é númerico mesmo....
## Alguma dica?
## Tem formas mais inteligente de comparar interação com multcomp?
Att
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110628/8a07ca48/attachment.html>