[R] simulation

[Peter Dalgaard]

Dimitris Rizopoulos wrote:
> maybe you're looking for something like this:
>
> x <- rpois(999, 2000)
> y <- numeric(length(x))
> for (i in seq_along(x))
>      y[i] <- sum(exp(rgamma(x[i], scale = 2, shape = 0.5)))
>
>   
Or use sapply,
sapply(x, function(x) sum(exp(rgamma(x[i], scale = 2, shape = 0.5)) )

or even

replicate(999, sum(exp(rgamma(rpois(1,2000), scale = 2, shape = 0.5)) )