[R] Monte Carlo Simulation

[Shane Phillips]

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 variables...

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.). 

Please help!

Thanks!

Shane