[R-br] poder a posteriori e regressão multinomial: como responder a um parecerista?

[Edson Lira]

Pedro, tem um livro Probability and Statistics With R, que tem o pacote PASWR, veja os exemplos abaixo:   (paginas 506 e 503  respectivamente). 

Espero que te ajude de alguma forma. 

require(PASWR)
#ANOVA
MEANS<-c(390, 405, 415, 410)
a<-length(MEANS)
n<-6
N=a*n
SD<-20
alpha<-0.05
Y<-rep(MEANS,rep(n,a))

treat<-factor(rep(1:a, rep(n,a)))
SStreat<-summary(aov(Y~treat))[[1]][1,2]
lambda<-SStreat/SD^2

cv<-qf(1-alpha,a-1,N-a)
power<-1-pf(cv,a-1, N-a, ncp=lambda)
power

power.anova.test(groups=a,n=n, between.var=var(MEANS),within.var=SD^2)

#TESTE T

MEANSt<-c(405,390)
at<-length(MEANSt)
nt<-6
Nt<-at*nt
dfe<-Nt-at
SDt<-10
Yt<-rep(MEANSt,rep(nt,at))
treatt<-factor(rep(1:at,rep(nt,at)))
SStreatt<-summary(aov(Yt~treatt))[[1]][1,2]
lambdat<-SStreatt/SDt^2
Gamma<-sqrt(lambdat)
Gamma
cvt<-qt(1-alpha,dfe)
Power<-1-pt(cvt,dfe,ncp=Gamma)
Power
power.t.test(n=6,delta=15, sd=10, alternative="one.sided")



Edson Lira

Estatístico

Manaus-Amazonas

--- Em sex, 1/4/11, Pedro Emmanuel Alvarenga Americano do Brasil <emmanuel.brasil em gmail.com> escreveu:

De: Pedro Emmanuel Alvarenga Americano do Brasil <emmanuel.brasil em gmail.com>
Assunto: [R-br] poder a posteriori e regressão multinomial: como responder a um parecerista?
Para: "R BR lista" <r-br em listas.c3sl.ufpr.br>
Data: Sexta-feira, 1 de Abril de 2011, 14:46

Amigos de R,

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.



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. 

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? 



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?
     

Abraço forte e que a força esteja com você,



Dr. Pedro Emmanuel A. A. do Brasil
Instituto de Pesquisa Clínica Evandro Chagas
Fundação Oswaldo Cruz
Rio de Janeiro - Brasil
Av. Brasil 4365
Tel 55 21 3865-9648
email: pedro.brasil em ipec.fiocruz.br


email: emmanuel.brasil em gmail.com

---Apoio aos softwares livres
www.zotero.org - gerenciamento de referências bibliográficas. 


www.broffice.org ou www.openoffice.org - textos, planilhas ou apresentações.
www.epidata.dk - entrada de dados.


www.r-project.org - análise de dados.
www.ubuntu.com - sistema operacional


-----Anexo incorporado-----

_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110404/387ce3d3/attachment.html>