Shane Phillips
2011-Apr-11 20:16 UTC
[R] Help on calculating a variable using random numbers
I'm new to R, but I'm trying to write a program for a dissertation that generates a dataset as follows... subject=1:1000 treat=rbinom(1*1000,1,.13) gender=rbinom(1*1000,1,.5) 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)) simtest I need to add a variable named "growth". If the "treat" variable for an observation is 0 then "growth" needs to be a randomly generated a number from a normal distribution with a mean of .1 and a sd of .03. If the "treat" variable is 1 then "growth" needs to be a randomly generated a number from a normal distribution with a mean of .5 and a sd of .03. Please help!
Dennis Murphy
2011-Apr-11 21:09 UTC
[R] Help on calculating a variable using random numbers
Hi: Try simtest <- transform(simtest, growth = rnorm(1000, m = ifelse(treat == 0, 0.1, 0.5), s = 0.03)) HTH, Dennis On Mon, Apr 11, 2011 at 1:16 PM, Shane Phillips <SPhillips@lexington1.net>wrote:> I'm new to R, but I'm trying to write a program for a dissertation that > generates a dataset as follows... > > subject=1:1000 > treat=rbinom(1*1000,1,.13) > gender=rbinom(1*1000,1,.5) > 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)) > simtest > > I need to add a variable named "growth". If the "treat" variable for an > observation is 0 then "growth" needs to be a randomly generated a number > from a normal distribution with a mean of .1 and a sd of .03. If the > "treat" variable is 1 then "growth" needs to be a randomly generated a > number from a normal distribution with a mean of .5 and a sd of .03. > > Please help! > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]