Walmes, bom dia<br><br>Pois bem, li seu material e consegui multiplicar! Se você puder ajudar, aqui vai uma questão:<br><br>Ajustei um modelo que consta no livro Ramalho et al., 2005, cap 8 pág 126, igual ao que eu transcrevi no 1º email, que eu repito aqui:<br>

<br>\[ y_{ijk} = \mu + t_i + r_j + (tr)_{ij} + s_k + (rc)_{jk} + (ts)_{ik} + e_{ijk} \]<br><br>Estou interassado realmente nos comp da variância... que independe do problema do satterthwaite (já que não envolve a fonte do problema, a meu ver!). Mas uma coisa interessante para se estimar é se a fonte de variação colheita é significativa, assim pretendia estimar a significância dessa fonte. Pela esperança dos QM para o teste F, eu preciso combinar QMs da seguinte forma:<br>

<br>\[ F_{colheitas} = \frac{QM_{colheita} + QM_{erro}} {QM_{bloco x colheita} + QM_{trat x colheita}} \]<br><br>No mesmo livro, Ramalho, et al., 2005, cap 5, pág 76, tem as expressões do Satterthwaite para os dois graus de liberdade (numerador e denominador), que são iguais!.<br>
<br>Aí fui olhar seu código, com os respectivos exemplos...<br><br>Pelo que pude compreender sua função da expressão de sattertwaite estima só o gl do numerador (ou denominador), estou errado? <i>Se eu quiser além desse, estimar o grau de liberdade do denominador, que também é uma combinação de QMs posso usar a mesma expressão?</i><br>

<br>Comparativamente ao seu exemplo no tópico 15, BL é bloco, obviamente, ES é o equivalente ao
meu tratamento, AD é seria o equivalente a minha fonte colheita. Pude perceber que é análogo ao modelo do tópico 15, mas acrescentei mais uma interação, a bloco x colheita (BL:ES). <br>
<br>Agora o que exatamente pega é... o teste F para colheitas eu consigo estimar pela expressão, mas os devidos graus de liberdade (do numerador e do denominador) eu usei a função "satter" e <i>no número de níveis do fator, como em alguns QMs são interações eu pus o produto dos níveis de cada fator envolvido, isso está certo?</i><br>
<br>A propósito olhei o material, muito bom hein... E mais, pelo menos no meu entendimento essa expressão do satterthwaite valia um CMR lá no ridiculas!<br><br>obrigado abraço,<br>FH<br><br><br><div class="gmail_quote">
2011/4/30 Walmes Zeviani <span dir="ltr"><<a href="mailto:walmeszeviani@gmail.com" target="_blank">walmeszeviani@gmail.com</a>></span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

<font face="trebuchet ms,sans-serif">Fernando,<br><br>Veja o tópico 15 e 16 desse material <a href="http://www.leg.ufpr.br/%7Ewalmes/cursoR/cursoR4.pdf" target="_blank">http://www.leg.ufpr.br/~walmes/cursoR/cursoR4.pdf</a>. Fiz uma função para obter o valor do grau de liberdade e quadrado médio combinado para o desdobramento da interação.<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><br><div class="gmail_quote"><div><div></div><div>2011/4/29 Fernando Henrique Toledo <span dir="ltr"><<a href="mailto:fernandohtoledo@gmail.com" target="_blank">fernandohtoledo@gmail.com</a>></span><br>
</div></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div>
Prezados colegas de lista,<br><br>Realizei uma análise em DBC (4 repetições, 10 tratamentos, mais de uma colheita), conforme o modelo de parcelas subdividas no tempo, ou seja: y_{ijk} = \mu + t_i + r_j + (tr)_{ij} + s_k + (rc)_{jk} + (ts)_{ik} + e_{ijk}, assim nas os testes F (e o p-valor), para a fonte de variação colheita (s_k) deve ser feitas de acordo com as esperanças e para tal faz-se necessário a combinação de alguns quadrados médios, portanto tem que se realizar a aproximação por Welch-Satterthwaite, como posso realizar isso no R? Aberto a qualquer sugestão... Grato desde já!<br>



<br>att,<br><font color="#888888">FH<br>
</font><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>
<br></blockquote></div><br>
<br>_______________________________________________<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>
<br></blockquote></div><br>