[R-br] multcomp lmer contraste interação..
Eder David Borges da Silva
eder em leg.ufpr.br
Quarta Junho 29 13:36:50 BRT 2011
Ivan e Walmes,
Vou dar uma olhada nas opções sugeridas, mas o armazenVIDRO não dar certo é
uma coisa meio estranha mesmo assim.
Walmes: Este não é o geral o qual conversamos, é uma outra coisa comparando
tipos de armazenamento, mas conversamos sobre esses dados com calma depois.
Obrigado pelas dicas
Att
Em 29 de junho de 2011 10:35, Walmes Zeviani <walmeszeviani em gmail.com>escreveu:
> É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: walmes em ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110629/963af5bc/attachment.html>
Mais detalhes sobre a lista de discussão R-br