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]]