<div dir="ltr"><div>Um acho que eu entendi, será que é possível Walmes vc retirar mais algumas dúvidas que me surgiram.<br><br><font color="#000000"><span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt"><span>1)<span style="font:7pt "Times New Roman""> </span></span></span><span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt">Você
modelou os dados da seguinte forma (y = Xb + Zp + ɛ).</span></font><font color="#000000" face="Times New Roman" size="3">
</font><p style="margin:0cm 0cm 0pt 64.35pt;text-indent:0cm"><span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt"><font color="#000000">em que: </font></span></p><font color="#000000" face="Times New Roman" size="3">
</font><p style="margin:0cm 0cm 0pt 64.35pt;text-indent:0cm"><span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt"><font color="#000000">y: vetor de dados da variável resposta.</font></span></p>
<font color="#000000" face="Times New Roman" size="3">
</font><p style="margin:0cm 0cm 0pt 64.35pt;text-indent:0cm"><span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt"><font color="#000000">b: vetor de efeitos fixos (média geral (µ)
+ nitrogênio (N) + potássio (K) + Tempo (T) + N:K + N:T + K:T + K:N:T).</font></span></p><font color="#000000" face="Times New Roman" size="3">
</font><p style="margin:0cm 0cm 0pt 64.35pt;text-indent:0cm"><span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt"><font color="#000000">p: vetor de efeitos das parcelas (N:K:repetições).</font></span></p>
<font color="#000000" face="Times New Roman" size="3">
</font><p style="margin:0cm 0cm 0pt 64.35pt;text-indent:0cm"><span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt"><font color="#000000">X e Z são as matrizes de incidência dos
referidos efeitos. <span> </span></font></span></p><font color="#000000" face="Times New Roman" size="3"></font><font color="#000000"><span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt"><span>2)<span style="font:7pt "Times New Roman""> </span></span></span><span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt">A
variância do erro (σ<sup>2</sup><sub>ɛ </sub>= 5,471) é o meu testador de todos
os efeitos fixos. Fiquei com esta dúvida, pois não conheço as esperanças de
quadrado médio deste modelo.</span></font><font size="3"><font face="Times New Roman"><font color="#000000">
<span style="line-height:150%;font-family:"Times New Roman","serif";font-size:12pt"><br>3) Generalizando para
casos em que eu tenha vários níveis para os meus fatores e queira realizar um
procedimento de comparações múltiplas precisarei calcular algum tipo de erro
combinado. </span></font></font></font><br><br></div><div>Desde já agradeço pela ajuda.<br><br></div><div>Obrigado.<br></div><div><br><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 2 de junho de 2014 11:45, walmes . <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 dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><span style="font-family:courier new,monospace"><span style="font-family:trebuchet ms,sans-serif">Baseado nas discussões, segue um exemplo de análise.</span><br>
</span></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><span style="font-family:courier new,monospace"><br>##-----------------------------------------------------------------------------<br><br>
pkg <- c("lattice", "nlme", "doBy")<br>sapply(pkg, require, character.only=TRUE)<div class=""><br><br>data <- structure(list(N = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 100L, 100L,<br>100L, 100L, 100L, 100L, 100L, 100L), K = c(0L, 0L, 0L, 0L, 50L,<br>
50L, 50L, 50L, 0L, 0L, 0L, 0L, 50L, 50L, 50L, 50L), Tempo = c(10L,<br>10L, 30L, 30L, 10L, 10L, 30L, 30L, 10L, 10L, 30L, 30L, 10L, 10L,<br>30L, 30L), resp = c(24.30206, 32.96457, 36.58875, 42.1957, 43.33543,<br>44.00179, 45.53739, 51.26205, 54.23789, 54.94291, 55.0155, 57.36688,<br>
58.63317, 66.21197, 68.81468, 69.32965)), .Names = c("N", "K",<br>"Tempo", "resp"), class = "data.frame", row.names = c(NA, -16L<br>))<br><br></div>str(data)<br><br>## Gráfico.<br>
xyplot(resp~Tempo|N, groups=K, data=data, type=c("p","a"))<br><br>## Tabela de frequência.<br>ftable(formula=Tempo~N+K, data=data)<br><br>## Os fatores tem 2 níveis, não requer conveter para fator.<br>
<br>data<br><br>## Não informação de unidades experimentais, então assumir que está na<br>## ordem correta e colocar o indicador de repetição.<br><br>data$rept <- 1:2<br>data$parcela <- with(data, interaction(N, K, rept))<br>
str(data)<br><br>## Especificação e ajuste.<br>m0 <- lme(resp~N*K*Tempo, random=~1|parcela, data=data)<br><br>## Dignóstico simples.<br>r <- residuals(m0)<br>f <- fitted(m0)<br>plot(r~f)<br>plot(sqrt(abs(r))~f)<br>
qqnorm(r)<br><br>## Teste de Wald para os termos de efeito fixo.<br>anova(m0)<br><br>## Quadro de estimativas e medidas de ajuste.<br>summary(m0)<br><br>## Predição.<br>pred <- data.frame(N=c(0,100,50,50,50,50),<br> K=c(25,25,0,50,25,25),<br>
Tempo=c(20,20,20,20,10,30))<br><br>pred$y <- predict(m0, newdata=pred, level=0)<br>pred<br><br>##-----------------------------------------------------------------------------</span><br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif;display:inline">
À disposição.<span class="HOEnZb"><font color="#888888"><br>Walmes.</font></span></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">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><br clear="all"><br>-- <br><div style="text-align:center">
<font><b>José Henrique Soler Guilhen</b></font><br><font>Graduando em Agronomia</font><font style="font-size:10pt" face="Verdana"></font><br></div><div style="text-align:center"><font style="font-size:10pt" face="Verdana"><font style="font-size:10pt">Universidade Federal do Espírito Santo - CCAUFES</font></font><br>
</div><font style="font-size:10pt" face="Verdana"><br></font>
</div>