<font face="trebuchet ms,sans-serif">Tura,<br><br>A função agricolae::design.bib() faz a matriz de deliamento para experimentos desse tipo. Pórem, para o seu caso sendo 156 níveis de tratamento, com cada bloco (pessoa/entrevistado) contendo 6 níveis, estora a capacidade de cálculo da função. Com 2 e 3 níveis por bloco, os resultados são esses<br>
<br><br>> require(agricolae)<br>> trt <- 1:156<br>> k <- 2<br>> bib <- design.bib(trt, k, number=1, seed=41, kinds="Super-Duper") # seed = 41<br><br>Parameters BIB<br>==============<br>Lambda : 1<br>
treatmeans : 156<br>Block size : 2<br>Blocks : 12090<br>Replication: 155 <br><br>Efficiency factor 0.5032258 <br><br><<< Book >>><br>> str(bib)<br>'data.frame': 24180 obs. of 3 variables:<br>
$ plots: num 1 2 3 4 5 6 7 8 9 10 ...<br> $ block: Factor w/ 12090 levels "1","2","3","4",..: 1 1 2 2 3 3 4 4 5 5 ...<br> $ trt : Factor w/ 156 levels "1","2","3","4",..: 111 72 102 125 88 34 105 19 84 97 ...<br>
> k <- 3<br>> bib <- design.bib(trt, k, number=1, seed=41, kinds="Super-Duper") # seed = 41<br><br>Parameters BIB<br>==============<br>Lambda : 154<br>treatmeans : 156<br>Block size : 3<br>Blocks : 620620<br>
Replication: 11935 <br><br>Efficiency factor 0.6709677 <br><br><<< Book >>><br>> str(bib)<br>'data.frame': 1861860 obs. of 3 variables:<br> $ plots: num 1 2 3 4 5 6 7 8 9 10 ...<br> $ block: Factor w/ 620620 levels "1","2","3","4",..: 1 1 1 2 2 2 3 3 3 4 ...<br>
$ trt : Factor w/ 156 levels "1","2","3","4",..: 53 109 62 79 155 94 24 72 32 91 ...<br>> <br><br>Penso que se você usar modelo misto, você não precisa ter um delineamento tão balanceado assim. Tudo depende dos objetivos da sua pesquisa.<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">2011/10/6 Bernardo Rangel Tura <span dir="ltr"><<a href="mailto:tura@centroin.com.br">tura@centroin.com.br</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Mestres em R e DOE!<br>
<br>
Preciso de uma ajuda para resolver um projeto de pesquisa.<br>
<br>
No meu experimento temos 156 níveis (ou tratamentos) que preciso analisar e só posso pesquisar 6 em cada pessoa entrevistada. Pensei em trabalhar com um modelo de blocos incompletos balanceados (BIB).<br>
<br>
Desta forma seriam necessárias 12.480 entrevistas, pois cada nível seria avaliado por 80 pessoas, mas eu precisaria agora gerar os blocos para o experimento. Como posso fazer isto no R?<br>
<br>
Uma segunda questão, esta mais difícil na minha opinião, estes 156 níveis são divididas em 3 estratos: A, B e C. Algumas pessoas sugerem que estas entrevistas teriam de ter blocos do tipo: 2 níveis A, 2 níveis B e 2 níveis C. Existe algum modelo que permita isto, qual seria? É possível implementar este plano no R como no caso anterior?<br>
<br>
Antecipadamente agradeço<br>
<br>
[]s<br>
Tura<br>
<br>
______________________________<u></u>_________________<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/<u></u>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-<u></u>guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br>