[R-br] Ajuda em desenho de experimentos

[Walmes Zeviani]

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 em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================


2011/10/6 Bernardo Rangel Tura <tura em 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 em 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111007/aa1f4d1f/attachment.html>