folks,
if i have a matrix of means:
a <- matrix(1:4, 2)
and a matrix of std deviations:
b <- matrix(5:8, 2)
and i want to create a matrix X of random variates such that X[i, j] is a draw from normal distribution with mean = a[i, j] and std dev =  b[i, j], i think i can do this?
X <- matrix(rnorm(4, mean = a, sd = b), ncol = 2)
now if I want to create 10 such matrices, I could possibly do this:
lapply(1:10, function(x) matrix(rnorm(4, mean = a, sd = b), ncol = 2))  
is there a better way to do this with perhaps just one call to the rnorm() function?
cheers,
murali 		 	   		  
	[[alternative HTML version deleted]]

