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?