<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Pedro, tem um livro Probability and Statistics With R, que tem o pacote PASWR, veja os exemplos abaixo:   (paginas 506 e 503  respectivamente). <br><br>Espero que te ajude de alguma forma. <br><br><span style="font-family: courier,monaco,monospace,sans-serif;">require(PASWR)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">#ANOVA</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">MEANS<-c(390, 405, 415, 410)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">a<-length(MEANS)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family:
 courier,monaco,monospace,sans-serif;">n<-6</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">N=a*n</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">SD<-20</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">alpha<-0.05</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">Y<-rep(MEANS,rep(n,a))</span><br style="font-family: courier,monaco,monospace,sans-serif;"><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">treat<-factor(rep(1:a, rep(n,a)))</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family:
 courier,monaco,monospace,sans-serif;">SStreat<-summary(aov(Y~treat))[[1]][1,2]</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">lambda<-SStreat/SD^2</span><br style="font-family: courier,monaco,monospace,sans-serif;"><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">cv<-qf(1-alpha,a-1,N-a)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">power<-1-pf(cv,a-1, N-a, ncp=lambda)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">power</span><br style="font-family: courier,monaco,monospace,sans-serif;"><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family:
 courier,monaco,monospace,sans-serif;">power.anova.test(groups=a,n=n, between.var=var(MEANS),within.var=SD^2)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">#TESTE T</span><br style="font-family: courier,monaco,monospace,sans-serif;"><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">MEANSt<-c(405,390)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">at<-length(MEANSt)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">nt<-6</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family:
 courier,monaco,monospace,sans-serif;">Nt<-at*nt</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">dfe<-Nt-at</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">SDt<-10</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">Yt<-rep(MEANSt,rep(nt,at))</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">treatt<-factor(rep(1:at,rep(nt,at)))</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">SStreatt<-summary(aov(Yt~treatt))[[1]][1,2]</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family:
 courier,monaco,monospace,sans-serif;">lambdat<-SStreatt/SDt^2</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">Gamma<-sqrt(lambdat)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">Gamma</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">cvt<-qt(1-alpha,dfe)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">Power<-1-pt(cvt,dfe,ncp=Gamma)</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family: courier,monaco,monospace,sans-serif;">Power</span><br style="font-family: courier,monaco,monospace,sans-serif;"><span style="font-family:
 courier,monaco,monospace,sans-serif;">power.t.test(n=6,delta=15, sd=10, alternative="one.sided")</span><br style="font-family: courier,monaco,monospace,sans-serif;"><br><br><br>Edson Lira<br>
Estatístico<br>
Manaus-Amazonas<br><br>--- Em <b>sex, 1/4/11, Pedro Emmanuel Alvarenga Americano do Brasil <i><emmanuel.brasil@gmail.com></i></b> escreveu:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>De: Pedro Emmanuel Alvarenga Americano do Brasil <emmanuel.brasil@gmail.com><br>Assunto: [R-br] poder a posteriori e regressão multinomial: como responder a um parecerista?<br>Para: "R BR lista" <r-br@listas.c3sl.ufpr.br><br>Data: Sexta-feira, 1 de Abril de 2011, 14:46<br><br><div id="yiv1592725925"><font color="#000066"><font face="courier new,monospace">Amigos de R,<br><br>Eu escrevi um relatorio de uma investigação em que um dos objetivos era tentar explorar potenciais preditores de um subgrupo de uma serie de casos. Inclui na investigação uma serie de pessoas sob suspeita de uma doença. Após a investigação inicial, apenas 81 sujeitos apresentaram investigação inicial inconclusiva
 (este era o subgrupo de interesse). Eu achei que era um numero muito reduzido, e mesmo assim eu tentei rodar uma multinomial com trés desfechos categóricos (doente/nao doente/abandono). CLaro que nao foi possivel rodar a regressão com múltiplos preditores e dos modelos univariados que nao deram erro na função, nenhum foi significativo. Tentei colocar isso nos resultados e na discussão, porque era um objetivo inicial, mas obviamente não havia muito o que discutir sobre esse aspecto. No entanto, um dos parecerista impôs a questão de saber qual seria o beta (poder) para essa análise. No início eu não havia pensado calcular o poder para essa análise, porque era um objetivo secundário, porque os dados eram retrospectivos a partir dos registros médicos, e porque não ha dados similares na literatura que poderiam basear uma estimativa de tamanho amostral de forma razoável.<br>

<br>Eu sou da política de sempre tentar agradar os pareceristas, salvo se for uma coisa descabida e absurda, como por exemplo que modificaria os objetivos iniciais. <br><br>Perguntas: (1) è razoável estimar o poder depois que as análises ja foram conduzidas? Ja li algumas coisas sobre estimativa de poder 'post-hoc' e essas nao foram muito animadoras. (2) procurei mas nao encontrei uma função n R que estimasse poder (ou tamnho amostral ) de regressões multinomiais. Alguém conhece alguma? <br>

<br>Se alguem achar que não é razoável estimar o poder nessa situação, poeria indicar argumentos razoáveis que pudesse utilizar com o parecerista?<br>     <br><br clear="all"></font></font>Abraço forte e que a força esteja com você,<br>

<br>Dr. Pedro Emmanuel A. A. do Brasil<br>Instituto de Pesquisa Clínica Evandro Chagas<br>Fundação Oswaldo Cruz<br>Rio de Janeiro - Brasil<br>Av. Brasil 4365<br>Tel 55 21 3865-9648<br>email: <a rel="nofollow" ymailto="mailto:pedro.brasil@ipec.fiocruz.br" target="_blank" href="/mc/compose?to=pedro.brasil@ipec.fiocruz.br">pedro.brasil@ipec.fiocruz.br</a><br>

email: <a rel="nofollow" ymailto="mailto:emmanuel.brasil@gmail.com" target="_blank" href="/mc/compose?to=emmanuel.brasil@gmail.com">emmanuel.brasil@gmail.com</a><br><br>---Apoio aos softwares livres<br><a rel="nofollow" target="_blank" href="http://www.zotero.org">www.zotero.org</a> - gerenciamento de referências bibliográficas. <br>

<a rel="nofollow" target="_blank" href="http://www.broffice.org">www.broffice.org</a> ou <a rel="nofollow" target="_blank" href="http://www.openoffice.org">www.openoffice.org</a> - textos, planilhas ou apresentações.<br><a rel="nofollow" target="_blank" href="http://www.epidata.dk">www.epidata.dk</a> - entrada de dados.<br>

<a rel="nofollow" target="_blank" href="http://www.r-project.org">www.r-project.org</a> - análise de dados.<br><a rel="nofollow" target="_blank" href="http://www.ubuntu.com">www.ubuntu.com</a> - sistema operacional<br>
</div><br>-----Anexo incorporado-----<br><br><div class="plainMail">_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="/mc/compose?to=R-br@listas.c3sl.ufpr.br">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/cgi-bin/mailman/listinfo/r-br</a><br></div></blockquote></td></tr></table>