Tura, 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

require(agricolae) trt <- 1:156 k <- 2 bib <- design.bib(trt, k, number=1, seed=41, kinds="Super-Duper") # seed = 41

Parameters BIB ============== Lambda : 1 treatmeans : 156 Block size : 2 Blocks : 12090 Replication: 155 Efficiency factor 0.5032258 <<< Book >>>

str(bib) 'data.frame': 24180 obs. of 3 variables: $ plots: num 1 2 3 4 5 6 7 8 9 10 ... $ block: Factor w/ 12090 levels "1","2","3","4",..: 1 1 2 2 3 3 4 4 5 5 ... $ trt : Factor w/ 156 levels "1","2","3","4",..: 111 72 102 125 88 34 105 19 84 97 ... k <- 3 bib <- design.bib(trt, k, number=1, seed=41, kinds="Super-Duper") # seed = 41

Parameters BIB ============== Lambda : 154 treatmeans : 156 Block size : 3 Blocks : 620620 Replication: 11935 Efficiency factor 0.6709677 <<< Book >>>

str(bib) 'data.frame': 1861860 obs. of 3 variables: $ plots: num 1 2 3 4 5 6 7 8 9 10 ... $ block: Factor w/ 620620 levels "1","2","3","4",..: 1 1 1 2 2 2 3 3 3 4 ... $ trt : Factor w/ 156 levels "1","2","3","4",..: 53 109 62 79 155 94 24 72 32 91 ...

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. À 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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2011/10/6 Bernardo Rangel Tura <tura@centroin.com.br>

Mestres em R e DOE!

Preciso de uma ajuda para resolver um projeto de pesquisa.

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).

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?

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?

Antecipadamente agradeço

[]s Tura

______________________________**_________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<http://www.leg.ufpr.br/r-br-guia>) e forneça código mínimo reproduzível.