<div dir="ltr"><p>Caros,</p><p> alguém sabe como faço para calcular o tamanho de uma 
amostra para um experimento com 3 tratamentos e um controle, rodados em 
paralelo, para um nível de significância de 5% e poder do teste de 80%?</p><p>
 A variável é uma proporção (taxa de usuários que se cadastraram num 
site). Eu sei calcular o tamanho da amostra para um teste de duas 
proporções, mas aqui eu quero testar quatro proporções e escolher a 
maior das quatro.</p><div class=""><p> Eu tenho a 
proporção esperada de cada um dos tratamentos e controle (seja p1 o 
controle, e p2, p3 e p4 o tratamento, p1 = .2, p2= .22, p3=.24 e p4 = 
.3).</p><p> O caminho mais simples é calcular o tamanho da amostra para 
duas proporções, considerando p1 e p4 (maior e menor, sendo que o menor é
 também o controle). Mas parece que assim estou subestimando a 
variabilidade em p2 e p3 e a possibilidade de que eles sejam maiores que
 p4 e p1. Além disso, parece haver o problema de múltiplas comparações.</p><p>
 Li um pouco e eu poderia tentar a correção de Bonferroni para múltiplos
 experimentos, e calcular o tamanho da amostra para duas propoções com 
p1 e p4, e usar alhpa_novo = alpha/m, em que alpha é o nível de 
significância original (5%) e m é o número de tratamentos (4... ou seria
 3?). Porém esse método não leva em conta os valores esperados de p2 e 
p3, apenas p1 e p4, e isso nåo parece certo. </p><p> Por fim, vale 
mebrar que o que quero é descobrir o melhor tratamento (conversão), e 
nåo testar se as propoções são iguais (eu já sei que não são exatamente 
iguais).</p><p> Eu fiz um script simulando para o meu caso particular 
vários tamanhos de amostra, e qual a probabilidade de eu escolher o 
melhor tratamento [ ou seja, prob(escolher p4| dados e que p4 é o 
melhor). Mas não estou seguro nem que meu script está certo e, mais do 
que isso, qual a relação entre o número de tratamtnos, e as diferenças 
esperadas no efeito de cada tratamento sobre o tamanho da amostra (no 
caso de duas propoções, é fåcil estabelecer quanto minha amostra aumenta
 pela diferença esperada entre p1 e p2, para um dado alpha e beta).</p><p> Se alguém tiver também referência na litertatura, já ajuda.<br> Obrigado.</p><p>Manoel<br></p><p>segue o script:</p><p>n.sim <-2000 # número de simulações<br>mat <- matrix( nrow=n.sim, ncol=4)<br>n.sample <- seq(500, 5500, by=1000) # tamanho de amostras para ver taxa de acerto para cada tamanho</p><p>## proporção de cada tratamento</p><p>## vou criar 10 combinações de proporção para 4 tratamentos</p><p>## a primeira todos os tratamentos têm a mesma proporção, igual a 20%</p><p>## me ajuda a ter uma base do resultado da simulação.</p><p>set.seed(2) # fixa a semente, para resultados serem reproduíveis<br></p><p>p1 <- c(.2, runif(9, .15, .3))<br>p2 <- c(.2, runif(9, .15, .3))<br>p3 <- c(.2, runif(9, .15, .3))<br>p4 <- c(.2, runif(9, .15, .3))</p><p>mat.p <- matrix(c(p1, p2, p3, p4), nrow=10, ncol=4, byrow=F)<br></p><p># para armazenar qual é o melhor tratamento, entre 10 combinações de efeitos diferentes dos tratamentos<br></p><p>maior <- numeric()</p><p># simulação<br></p><p>for ( i in 1:10) maior[i] <- which.max(mat.p[i,])<br><br>mat.true.max <- matrix(nrow=length(n.sample), ncol=10)<br>num.true.max  <- 0<br><br>for (h in 1:6) {<br>  for ( k in 1:10) {<br>    for ( i in 1:n.sim) {<br>      for ( j in 1:4) {<br>        mat[i,j] <- sum(rbinom(n.sample[h], 1, mat.p[k,j]))<br>      }<br>      num.true.max <- num.true.max + as.numeric(which.max(mat[i,]) == maior[k])<br>    } <br>    mat.true.max[h,k] <- num.true.max/n.sim<br>    num.true.max <- 0<br>  }<br>}</p><p># verifito taxa de acerto para cada um dos tamanhos da amostra. Acerto é escolher a maior proporção.<br></p><p>mat.true.max<br></p></div><div class="gmail_extra"><br><div class="gmail_quote">2015-08-24 19:04 GMT-03:00 Andre Oliveira <span dir="ltr"><<a href="mailto:andreolsouza@yahoo.com.br" target="_blank">andreolsouza@yahoo.com.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="color:rgb(0,0,0);background-color:rgb(255,255,255);font-family:bookman old style,new york,times,serif;font-size:16px"><div>Pessoal boa noite,</div><div>seria de grande valia divulgar as próximas transmissões no grupo rbr. <br></div><div><span></span></div><div> </div><div>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>  <br><div><br><br></div><div style="display:block"> <div style="font-family:bookman old style,new york,times,serif;font-size:16px"> <div style="font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:16px"> <div dir="ltr"> <font face="Arial" size="2"> Em Sábado, 22 de Agosto de 2015 21:26, David Feitosa <<a href="mailto:contato@davidfeitosa.com" target="_blank">contato@davidfeitosa.com</a>> escreveu:<br> </font> </div>  <br><br> <div><div><div><div dir="ltr">Eu tb achei muito boa.<div><br clear="none"><div>Sabia nem que o ShareLatex estava com integração com o R.</div><div>Fico usando o knitr local no RStudio mesmo.</div><div><br clear="none"></div><div>Aqui em Fortaleza não há muitos usuários de R.</div><div>Só conheço duas pessoas por aqui e uma delas sou eu. Hahaha</div><div><br clear="none"></div><div>Python tem recebido mais destaque por aqui do que o R.</div></div><div>Java ainda domina desenvolvimento e pesquisa na região.</div></div><div><br clear="all"><div><div><div dir="ltr"><div><div dir="ltr"><div><div><br clear="none"></div><div>Atenciosamente,</div><div><br clear="none"></div>David F.</div></div></div></div></div></div>
<br clear="none"><div><div>Em 22 de agosto de 2015 17:45, Luis Benites Sánchez <span dir="ltr"><<a rel="nofollow" shape="rect" href="mailto:lbenitesanchez@gmail.com" target="_blank">lbenitesanchez@gmail.com</a>></span> escreveu:<br clear="none"><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Prezado Jobenil,<div><br clear="none"></div><div>Muito obrigado, temos ainda mais palestras no Meetup de usuários de R em São Paulo, ainda temos 13 palestras mais. Em todos os meetup (se os palestrantes concordam) teremos as filmações das reuniões e estarão disponíveis no canal de YouTube do Meetup. </div><div><br clear="none"></div><div>Não esqueçam de junstar-se ao Meetup de usuários de R em São Paulo em nosso site <a rel="nofollow" shape="rect" href="http://www.meetup.com/pt/useR-SP/" target="_blank">www.meetup.com/pt/useR-SP/</a> para que possa conhecer mais detalhes sobre as futuras reuniões  que estão sendo realizadas no IME-USP (Auditório Jacy Monteiro e que são transmitidas ao vivo).</div><div><br clear="none"></div><div>A palestra do 4 e 18 de setembro que minha colega e esposa Rocío Maehara comentou em seu e-mail serão de 14 a 16 horas (sextas).</div><div><br clear="none"></div><div>Um abraço</div><div><br clear="none"></div><div>Luis Benites Sánchez</div><div>Organizador del Meetup de usuários de R em São Paulo</div><div><br clear="none"></div></div><div><br clear="none"><div><span>2015-08-22 17:26 GMT-03:00 Jobenil - Gmail <span dir="ltr"><<a rel="nofollow" shape="rect" href="mailto:pjobenil@gmail.com" target="_blank">pjobenil@gmail.com</a>></span>:<br clear="none"></span><div><div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div lang="PT-BR"><div><div><span style="font-size:11pt">Prezados,<u></u><u></u></span></div><div><span style="font-size:11pt"><u></u> <u></u></span></div><div><span style="font-size:11pt">Acabei de assistir à apresentação dos Professores Paulo, Walmes e Fernando. Parabéns pelo evento e aos organizadores do Meetup pela iniciativa. <u></u><u></u></span></div><div><span style="font-size:11pt">Este é uma grande oportunidade para quem trabalha isolado em alguns rincões do país a fora.<u></u><u></u></span></div><div><span style="font-size:11pt">Abraços<u></u><u></u></span></div><div><span style="font-size:11pt">Jobenil Júnior<u></u><u></u></span></div><div><span style="font-size:11pt"><u></u> <u></u></span></div><div><span style="font-size:11pt"><u></u> <u></u></span></div><div><b><span style="font-size:11pt">De:</span></b><span style="font-size:11pt"> R-br [mailto:<a rel="nofollow" shape="rect" href="mailto:r-br-bounces@listas.c3sl.ufpr.br" target="_blank">r-br-bounces@listas.c3sl.ufpr.br</a>] <b>Em nome de </b>Luiz Roberto Martins Pinto<br clear="none"><b>Enviada em:</b> sexta-feira, 21 de agosto de 2015 19:27<br clear="none"><b>Para:</b> <a rel="nofollow" shape="rect" href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br clear="none"><b>Assunto:</b> Re: [R-br] (Transmissão ao vivo) Meetup R São Paulo Users Group (IME-USP) Sexta 21 de agosto do 2015<u></u><u></u></span></div><div><div><div><u></u> <u></u></div><div><div>Walmes,<u></u><u></u></div><div><div><u></u> <u></u></div></div><div><div>Agradeço o link.<u></u><u></u></div></div><div><div><u></u> <u></u></div></div><div><div>Tive dificuldades para acessar. E ainda não entendi o que aconteceu. <u></u><u></u></div></div><div><div><u></u> <u></u></div></div><div><div>Luiz Roberto.<u></u><u></u></div></div></div><div><div><br clear="all"><u></u><u></u></div><div><div><div><div><div><div>Luiz Roberto Martins Pinto<br clear="none">Prof. Pleno/DCET/UESC<u></u><u></u></div></div><div><div>Laboratório de Estatística Computacional<u></u><u></u></div></div><div><div>Universidade Estadual de Santa Cruz<u></u><u></u></div></div><div><div>Ilhéus-Bahia-Brasil<br clear="none"><br clear="none"><a rel="nofollow" shape="rect" href="mailto:luizroberto.uesc@gmail.com" target="_blank">luizroberto.uesc@gmail.com</a><br clear="none">skype: lrmpinto<u></u><u></u></div></div><div><div style="margin-bottom:12pt"><a rel="nofollow" shape="rect" href="http://lattes.cnpq.br/2732314327604831" target="_blank">http://lattes.cnpq.br/2732314327604831</a> <br clear="none"><br clear="none"><u></u><u></u></div></div></div></div></div></div><div><u></u> <u></u></div><div><div>Em 21 de agosto de 2015 16:54, Walmes Zeviani <<a rel="nofollow" shape="rect" href="mailto:walmeszeviani@gmail.com" target="_blank">walmeszeviani@gmail.com</a>> escreveu:<u></u><u></u></div><blockquote style="border-width:medium medium medium 1pt;border-style:none none none solid;border-color:-moz-use-text-color -moz-use-text-color -moz-use-text-color rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm"><div><div><div style="margin-bottom:12pt"><span>Para os que não estiveram presentes, passaram por dificuldades ou querer rever, o vídeo também está disponível em <a rel="nofollow" shape="rect" href="https://www.youtube.com/watch?v=UgPX49gkby4" target="_blank">https://www.youtube.com/watch?v=UgPX49gkby4</a>. Caso interessar, o endereço memorizável do canal é <a rel="nofollow" shape="rect" href="http://www.leg.ufpr.br/youtube" target="_blank">http://www.leg.ufpr.br/youtube</a>, e o perfil na Google+ é <a rel="nofollow" shape="rect" href="http://www.leg.ufpr.br/g+" target="_blank">http://www.leg.ufpr.br/g+</a>.<u></u><u></u></span></div></div><div><div><span>À disposição.<u></u><u></u></span></div></div><div><div><span>Walmes.<u></u><u></u></span></div></div><div>​<u></u><u></u></div></div><div><br clear="none">_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a rel="nofollow" shape="rect" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a rel="nofollow" shape="rect" 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 clear="none">Leia o guia de postagem (<a rel="nofollow" shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<u></u><u></u></div></blockquote></div><div><u></u> <u></u></div></div></div></div></div></div><br clear="none">_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br clear="none">
<a rel="nofollow" shape="rect" 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 clear="none">
Leia o guia de postagem (<a rel="nofollow" shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br clear="none"></blockquote></div></div></div><br clear="none"><br clear="all"><div><br clear="none"></div>-- <br clear="none"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><pre></pre><div style="color:rgb(136,136,136);font-family:arial;white-space:normal">Luis B. Sánchez</div><span><span style="color:rgb(136,136,136);font-family:arial;white-space:normal">Ph.D Candidate</span></span><div style="color:rgb(136,136,136);font-family:arial;white-space:normal">Dept. of Statistics-IME, University at São Paulo (USP)</div><div style="color:rgb(136,136,136);font-family:arial;white-space:normal">São Paulo, BRAZIL</div><div style="color:rgb(136,136,136);font-family:arial;white-space:normal"><a rel="nofollow" shape="rect" href="mailto:lsanchez@ime.usp.br" style="font-size:12.8px" target="_blank">lsanchez@ime.usp.br</a><br clear="none"></div><div style="color:rgb(136,136,136);font-family:arial;white-space:normal">Skype: lbenitesanchez<br clear="none"></div><div><span style="color:rgb(136,136,136);font-family:arial;white-space:normal;font-size:13px">Twitter: @lbenitesanchez</span></div><div><a rel="nofollow" shape="rect" href="http://www.luis-benites.com/" style="font-family:arial;font-size:13px;white-space:normal" target="_blank">www.luis-benites.com/</a><span style="color:rgb(136,136,136);font-family:arial;white-space:normal;font-size:13px"><br clear="none"></span></div><div style="color:rgb(136,136,136);font-family:arial;white-space:normal"><br clear="none"></div>
</div></div></div></div></div></div></div></div>
</div>
<br clear="none">_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br clear="none">
<a rel="nofollow" shape="rect" 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 clear="none">
Leia o guia de postagem (<a rel="nofollow" shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br clear="none"></blockquote></div><br clear="none"></div></div></div></div><br><div>_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a shape="rect" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a shape="rect" 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 clear="none">Leia o guia de postagem (<a shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</div><br><br></div>  </div> </div>  </div></div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto: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" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature">Manoel Galdino<br><a href="https://sites.google.com/site/galdinomcz/" target="_blank">https://sites.google.com/site/galdinomcz/</a><br></div>
</div></div>