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
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.