Pessoal da lista:<div>Minha duvida:<div><div>require(lme4)</div><div>require(multcomp)</div><div>## Exemplo no lme4 moodificado</div><div>fm2 <- lmer(Reaction ~ Days/Subject + (0+Days|Subject), sleepstudy)</div><div>summary(fm2)</div>
<div># Tenho interesse neste contraste por exemplo:</div><div>summary(glht(fm2,c('Days:Subject3-Days:Subject1 = 0')) )</div><div>## Funcionou...</div><div>## Agora com meus dados querro coisa semelhante</div><div>
<br></div><div>### Agora com meu conjunto de dados:  </div><div>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'))</div>
<div>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')))</div><div>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')))</div>
<div>DADOSS$ID <- paste(DADOSS[,1],DADOSS[,2],DADOSS[,3],DADOSS[,4],DADOSS[,5],DADOSS[,6],sep='',collapse = NULL)</div><div>DADOSS$camaraM <- paste(DADOSS$sitio,DADOSS$ponto,DADOSS$camara,sep='')</div>
<div>DADOSS$resp <- rnorm(nrow(DADOSS))</div><div><br></div><div>model <- lmer(resp~armazen/tempo+(1|camaraM)+(1|ponto),DADOSS)</div><div>summary(model)</div><div>anova(model)</div><div>### Eu tenho interesse apenas em este contraste (e nos outros que envolvem a interação ):</div>
<div>summary(glht(model,c('armazenVIDRO:tempo-armazenPP:tempo = 0')))</div><div>## PQ não funciona????????????(Apesar de ter pouco cabelo, já arranquei alguns para tentar descobrir....)</div><div># OBS: tempo é númerico mesmo....</div>
<div>## Alguma dica?</div></div></div><div>## Tem formas mais inteligente de comparar interação com multcomp?</div><div>Att</div><div><br></div>