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
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@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@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}
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 <dmarcelino@live.com>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 <sergio.edfisica@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@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}
_______________________________________________ R-br mailing list R-br@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.
-- 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
participantes (2)
-
Daniel Marcelino -
Sérgio Henrique almeida da silva ju