Éder,
Confirmando as suspeitas do Ivan, que sugeriu remover o intercepto, também
acredito que tem algo relacionado ao nível de referência. Veja que o
contraste abaixo dá resultado
> summary(glht(model,c('armazenHEAD:tempo-armazenPP:tempo = 0')))
Simultaneous Tests for General Linear Hypotheses
Fit: lmer(formula = resp ~ armazen/tempo + (1 | camaraM) + (1 | ponto),
data = DADOSS)
Linear Hypotheses:
Estimate Std. Error z value
Pr(>|z|)
armazenHEAD:tempo - armazenPP:tempo == 0 0.08059 0.33029 0.244
0.807
(Adjusted p values reported -- single-step method)
O que envolve o nível armazenVIDRO (que é o nível de referência) não dá
certo. Então a opção seria mudar os níveis ajustando vários lmer() ou
remover o intercepto como Ivan sugeriu. Veja se a remoção do intercepto não
muda o valor dos contrastes (não deveria, mas vale a pena conferir). Outra
saída seria passar matrizes de contrastes para glht().
Fora do objetivo da mensagem, não seria aleatório o efeito de ponto dentro
de câmara conforme conversamos?
À disposição.
Walmes.
==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: [email protected]
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
