eu já encontrei a solução em breve posto pra vocês Obrigado Daniel
Em 2 de março de 2013 20:51, Daniel Marcelino <[email protected]>escreveu: > Sergio, > 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. > 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. > > p<- 2.5/100 > reps <- 2000 > resultados <- numeric(reps) > for (i in 1:reps) { > participacao <-c(rep(0,600),rep(1,1400)) > x <- 0 > while (participacao==0) { > participacao <- participacao + rbinom(1,1,p) > x <- x + 1 > } > x <- > resultados[i] <- x > } > > mean(resultados) > > > > 2013/3/2 Sérgio Henrique almeida da silva ju <[email protected]> > >> Olá >> >> Como posso simular 2000 amostras aleatórias com 2000 sujeitos cada uma? >> Quero modelar um desfecho binário (óbitos e não óbitos). >> >> A distribuição dos meus dados é: >> ### Meus Dados >> id <- rep(1:2000) >> servidor <- c(rep('Apoio',292),rep('Tecnico',776), rep('Superior',932)) >> participacao <- c(rep(0,600),rep(1,1400)) # 0 = Não participou do estudo >> e 1 = Participou do estudo >> obitos <- ????? >> >> 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? >> >> Obrigado >> >> >> >> -- >> Sérgio Henrique Almeida da Silva Junior >> Doutorando em Epidemiologia em Saúde Pública >> Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ >> http://lattes.cnpq.br/1611345552843383 >> Tel: (21) 68463637 / 94429486 >> >> _______________________________________________ >> R-br mailing list >> [email protected] >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> > > > > -- > "Small steps toward a much better world" > > \begin{signature} > Daniel Marcelino > Land Phone 1+514 343 6111 #3799 > 3200 Jean Brillant, Office C5071 > Montreal, QC; H3T 1N8 > Canada > \end{signature} > > _______________________________________________ > R-br mailing list > [email protected] > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 / 94429486
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
