[R] Monte Carlo Simulation
SPhillips at Lexington1.net
Sat Apr 16 00:23:05 CEST 2011
I apologize for my last post. here is the script I forgot to paste!
eth=runif(1*1000, min=1, max=4)
cogat=rnorm(1*1000, 100, 16)
map=rnorm(1*1000, 200, 9)
simtest=data.frame (subject=subject, treat=treat, gender=gender, eth=round(eth,digits=0), cogat=round(cogat,digits=0),map=round(map,digits=0),growth)
Once again the issues are the correlations, the recoding of the categorical variable and the iterations.
Thanks to all who have helped so far! You are all so smart!
From: Charles Annis, P.E. [charles.annis at statisticalengineering.com]
Sent: Friday, April 15, 2011 3:00 PM
To: Shane Phillips; r-help at r-project.org
Subject: RE: [R] Monte Carlo Simulation
What have you tried so far?
It is often helpful to begin with a much simpler problem, then add
complexity incrementally until you've constructed the desired model.
Charles Annis, P.E.
Charles.Annis at StatisticalEngineering.com
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
Behalf Of Shane Phillips
Sent: Friday, April 15, 2011 2:50 PM
To: r-help at r-project.org
Subject: [R] Monte Carlo Simulation
Hello, R friends...
I am very new to R, and I need some help. I am trying to construct a
simulation for my dissertation.
I need to create 1000 datasets of 1000 subjects with the following
Treatment variable - Drawn from a binomial distribution (1 run, prob=.13)
Covariate 1 - Drawn from a normal distribution (mean=100, sd=16)
Covariate 2 - Drawn from a normal distribution (mean=200, sd=9)
Covariates 1 and 2 need to be correlated (say, r=.80)
Covariate 3 - Drawn from a binomial distribution (1 run, prob=.5)
Covariate 4 - Drawn from a distribution of discrete variables where 1 has an
80% chance of being selected, 2 - 10%, 3 - 5% and 4 - 5%. This variable
would need to be recoded into 4 binary variables.
Covariate 5 - Drawn from a normal distribution (mean=84, sd=2)
Covariate 6 - Drawn from a binomial distribution (1 run, prob=.15)
Covariate 6 needs to correlate with Covariate 2 (r=.70, or so)
I need each dataset saved as a new datafile with an iterative filename (e.g.
sample1, sample2, etc.).
R-help at r-project.org mailing list
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
More information about the R-help