[R-br] simulação de dados

Daniel Marcelino dmarcelino em live.com
Sábado Março 2 20:51:49 BRT 2013


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 <sergio.edfisica em gmail.com>

> 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
> R-br em listas.c3sl.ufpr.br
> 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}
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130302/0ea56c17/attachment.html>


Mais detalhes sobre a lista de discussão R-br