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