<html><head></head><body><div style="color:#000; background-color:#fff; font-family:bookman old style, new york, times, serif;font-size:16px"><div id="yui_3_16_0_ym19_1_1461439790803_3473">O meu e-mail foi respondido, porém, não sei o motivo de não ter recebido de volta a respostas em meu e-mail! Estou reenviando com a respostas e minha outra colocação. <br></div><div id="yui_3_16_0_ym19_1_1461439790803_3799"><br></div><div id="yui_3_16_0_ym19_1_1461439790803_3801">E-mail original</div><div id="yui_3_16_0_ym19_1_1461439790803_3825"></div><div id="yui_3_16_0_ym19_1_1461439790803_3483">Prezad@s;</div><div id="yui_3_16_0_ym19_1_1461439790803_3490" dir="ltr">
<a id="yui_3_16_0_ym19_1_1461439790803_3484" href="http://r-br.2285057.n4.nabble.com/R-br-Otimizacao-Fatorial-td4665997.html" name="yui_3_16_0_ym19_1_1461280577114_49221" target="_top" rel="nofollow"></a><a id="yui_3_16_0_ym19_1_1461439790803_3485" href="http://r-br.2285057.n4.nabble.com/R-br-Otimizacao-Fatorial-td4665997.html" name="yui_3_16_0_ym19_1_1461280577114_4924" target="_top" rel="nofollow"></a>
estou iniciando um projeto em química onde será montado um
experimento no esquema fatorial onde tenho alguns fatores definidos a
princípio
e gostaria de otimizar este processo. A questão é que tenho 4
respostas que serão medidas para cada tratamento. Gostaria de
saber se a biblioteca rsm, lm e aov seriam a mais indicada para que eu
possa analisar depois de montar o experimento e buscar o deslocamentos
necessários a fim de obter a otimização do processo nos
próximos experimentos. Estou buscando material nesta linha e usando o
R de preferência. Minha dúvida também se refere ao fato que tenho
encontrado materiais com modelos com apenas uma resposta. Se
alguém puder ajudar ficarei grato.</div><div id="yui_3_16_0_ym19_1_1461439790803_3511"><div id="yui_3_16_0_ym19_1_1461439790803_3514"> </div><div id="yui_3_16_0_ym19_1_1461439790803_3828"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1461439790803_3697">O Ivan Bezerra Allaman respondeu da seguinte forma <br></div><div id="yui_3_16_0_ym19_1_1461439790803_3708" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1461439790803_3707" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1461439790803_4011" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1461439790803_4010" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1461439790803_3654">O é que eu entendi sua dúvida...</div><div id="yui_3_16_0_ym19_1_1461439790803_3579"><br id="yui_3_16_0_ym19_1_1461439790803_3580"></div><div id="yui_3_16_0_ym19_1_1461439790803_3581"><div id="yui_3_16_0_ym19_1_1461439790803_3582">vars=list()</div><div id="yui_3_16_0_ym19_1_1461439790803_3583">for(i in 1:30) vars[[i]] <- rnorm(48)</div><div id="yui_3_16_0_ym19_1_1461439790803_3584">names(vars) <- paste('var',1:30,sep='')</div><div id="yui_3_16_0_ym19_1_1461439790803_3585">vars</div><div id="yui_3_16_0_ym19_1_1461439790803_3586"><br id="yui_3_16_0_ym19_1_1461439790803_3587"></div><div id="yui_3_16_0_ym19_1_1461439790803_3588">varis <- do.call('cbind',vars)</div><div id="yui_3_16_0_ym19_1_1461439790803_3589"><br id="yui_3_16_0_ym19_1_1461439790803_3590"></div><div id="yui_3_16_0_ym19_1_1461439790803_3591">da <- expand.grid(FatA=rep(LETTERS[1:6],4),</div><div id="yui_3_16_0_ym19_1_1461439790803_3592"> FatB=letters[1:2])</div><div id="yui_3_16_0_ym19_1_1461439790803_3593"><br id="yui_3_16_0_ym19_1_1461439790803_3594"></div><div id="yui_3_16_0_ym19_1_1461439790803_3595">dad <- data.frame(da,varis)</div><div id="yui_3_16_0_ym19_1_1461439790803_3596"><br id="yui_3_16_0_ym19_1_1461439790803_3597"></div><div id="yui_3_16_0_ym19_1_1461439790803_3598">tabs <- apply(dad[,-c(1:2)],</div><div id="yui_3_16_0_ym19_1_1461439790803_3599"> 2,</div><div id="yui_3_16_0_ym19_1_1461439790803_3600"> function(x) aov(x ~ FatA*FatB,dad))</div><div id="yui_3_16_0_ym19_1_1461439790803_3601"><div id="yui_3_16_0_ym19_1_1461439790803_3846" dir="ltr">lapply(tabs,summary) <br></div><div id="yui_3_16_0_ym19_1_1461439790803_3845" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1461439790803_3909" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1461439790803_3830" dir="ltr">Oi Ivan,</div><div id="yui_3_16_0_ym19_1_1461439790803_3856" dir="ltr">no caso o que procuro é uma função que faça o que você escreveu mesmo, mas no lugar de x e aov seria algo multivariado como abaixo.</div><div id="yui_3_16_0_ym19_1_1461439790803_3999" dir="ltr">fit=manova(varis~da, data= dad)<br></div><div id="yui_3_16_0_ym19_1_1461439790803_4023" dir="ltr"> Neste caso a manova seria pouco informativo sobre o otimização do processo. Queria algum modelo com a estrutura do fit que escrevi, mas com mais informações sobre a otimização individual para cada y os fatores em questão. <br></div><div id="yui_3_16_0_ym19_1_1461439790803_4075" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1461439790803_4076" dir="ltr">obrigado <br></div><div id="yui_3_16_0_ym19_1_1461439790803_3831" dir="ltr"><br></div></div></div></div><div id="yui_3_16_0_ym19_1_1461439790803_3857" class="signature">André Oliveira Souza.
Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES</div></div></body></html>