Sergio, <div>Não sei se entendi direito o seu problema. Também não sei o porquê da coluna id. Você terá que ter uma população maior do que a amostra para selecionar amostras de 2k. Rascunhei  alguma coisa. </div><div>Veja que a rbinom é como se atirasse um moeda. Para cada vez que a moeda fosse atirada incrementar-se-ia a probabilidade que guardamos em x.</div>

<div><div><br></div><div><div>p<- 2.5/100</div><div>reps <- 2000</div><div>resultados <- numeric(reps)</div><div>for (i in 1:reps) {</div><div>participacao <-c(rep(0,600),rep(1,1400))</div><div>x <- 0</div>

<div>while (participacao==0) {</div><div>participacao <- participacao + rbinom(1,1,p)</div><div>  x <- x + 1</div><div>}</div><div>x <-</div><div>resultados[i] <- x</div><div>}</div></div></div><div><br></div>

<div>mean(resultados)</div><div><br></div><div><br><br><div class="gmail_quote">2013/3/2 Sérgio Henrique almeida da silva ju <span dir="ltr"><<a href="mailto:sergio.edfisica@gmail.com" target="_blank">sergio.edfisica@gmail.com</a>></span><br>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Olá<div><br></div><div>Como posso simular 2000 amostras aleatórias com 2000 sujeitos cada uma? Quero modelar um desfecho binário (óbitos e não óbitos).</div>


<div><br></div><div>A distribuição dos meus dados é:<br><div><div>
### Meus Dados</div><div>id <- rep(1:2000)</div><div>servidor <- c(rep('Apoio',292),rep('Tecnico',776), rep('Superior',932))</div><div>participacao <- c(rep(0,600),rep(1,1400)) # 0 = Não participou do estudo e 1 = Participou do estudo</div>



</div><div>obitos <- ?????</div><div><br></div><div>Como posso distribuir meus óbitos em relação a variável participacao, sendo uma chande de 2,5 de óbitos de quem não participou em comparação a quem participou?</div>


<div>
<br></div><div>Obrigado</div><span><font color="#888888"><div><br></div><div>  </div><div><br></div>-- <br>Sérgio Henrique Almeida da Silva Junior<br>Doutorando em Epidemiologia em Saúde Pública<br>Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ<br>




<a href="http://lattes.cnpq.br/1611345552843383" target="_blank">http://lattes.cnpq.br/1611345552843383</a> <br>Tel: <a href="tel:%2821%29%2068463637" value="+12168463637" target="_blank">(21) 68463637</a> / 94429486<br>




</font></span></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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>
Leia o guia de postagem (<a 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></blockquote></div><br><br clear="all"><div><br></div>-- <br>


"Small steps toward a much better world"<br><br>\begin{signature}<br>Daniel Marcelino<br>Land Phone 1<a href="tel:%2B514%20343%206111%20%233799" value="+5143436111" target="_blank">+514 343 6111 #3799</a><br>3200 Jean Brillant,  Office C5071<br>

Montreal, QC; H3T 1N8<br>Canada<br>\end{signature}
</div>