Pessoal,<div>Rodei agora tirando o intercepto, como esperado os coeficientes estão os mesmo e funcionou os contrastes:</div><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~0+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('armazenHEAD:tempo-armazenPP:tempo = 0',</div><div> 'armazenHEAD:tempo-armazenVIDRO:tempo = 0',</div><div> 'armazenPP:tempo-armazenVIDRO:tempo = 0')))</div>
<div><br></div><div>Valeu pelas dicas</div><div>Faloww</div><div><br><br><div class="gmail_quote">Em 29 de junho de 2011 13:36, Eder David Borges da Silva <span dir="ltr"><<a href="mailto:eder@leg.ufpr.br">eder@leg.ufpr.br</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Ivan e Walmes,<div><br></div><div>Vou dar uma olhada nas opções sugeridas, mas o <span style="font-family:'trebuchet ms', sans-serif">armazenVIDRO não dar certo é uma coisa meio estranha mesmo assim.</span></div>
<div><font face="'trebuchet ms', sans-serif"><br></font></div><div><font face="'trebuchet ms', sans-serif">Walmes: Este não é o geral o qual conversamos, é uma outra coisa comparando tipos de armazenamento, mas conversamos sobre esses dados com calma depois.<br>
</font></div><div><font face="'trebuchet ms', sans-serif">Obrigado pelas dicas<br></font></div><div><div class="gmail_quote">Att</div><div class="gmail_quote"><br></div><div class="gmail_quote">
Em 29 de junho de 2011 10:35, Walmes Zeviani <span dir="ltr"><<a href="mailto:walmeszeviani@gmail.com" target="_blank">walmeszeviani@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div></div><div class="h5">
<font face="trebuchet ms,sans-serif">Éder,<br><br>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<br>
<br><span style="font-family:courier new,monospace">> summary(glht(model,c('armazenHEAD:tempo-armazenPP:tempo = 0')))</span><br style="font-family:courier new,monospace"><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace"> Simultaneous Tests for General Linear Hypotheses</span><br style="font-family:courier new,monospace"><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">Fit: lmer(formula = resp ~ armazen/tempo + (1 | camaraM) + (1 | ponto), </span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace"> data = DADOSS)</span><br style="font-family:courier new,monospace"><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">Linear Hypotheses:</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace"> Estimate Std. Error z value Pr(>|z|)</span><br style="font-family:courier new,monospace"><span style="font-family:courier new,monospace">armazenHEAD:tempo - armazenPP:tempo == 0 0.08059 0.33029 0.244 0.807</span><br style="font-family:courier new,monospace">
<span style="font-family:courier new,monospace">(Adjusted p values reported -- single-step method)</span><br style="font-family:courier new,monospace"><br>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().<br>
<br>Fora do objetivo da mensagem, não seria aleatório o efeito de ponto dentro de câmara conforme conversamos?<br><br>À disposição.<br>Walmes.<br><br clear="all"></font><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>
<br></div></div>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>
</blockquote></div><br></div></div>