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 <[email protected]> escreveu: De: Pedro Emmanuel Alvarenga Americano do Brasil <[email protected]> Assunto: [R-br] poder a posteriori e regressão multinomial: como responder a um parecerista? Para: "R BR lista" <[email protected]> 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: [email protected] email: [email protected] ---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 [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
