Which matrix do you need to "output"? Is it X below. Will REPLICATE be set to 1000 when you get this working as desired.

On Wed, Mar 26, 2008 at 10:35 AM, Doaa <d_faik@hotmail.com> wrote: > i made a program on proc iml that generate a data from a multivariate > normal distribution, that i will use in simulation study iwant to > generate about 1000 data sets. i made it by loop statment but my > problem now how to save each data set with a separeate name and how to > reopen these data sets again . i think it is done by a loop but i > don't know how to do it.please replay as sson as possible > my code was > PROC IML; > do Replicate=1 to 10; > seed=12345; > p=3; > n=40; > h=((n-1)**2/n)* betainv(0.05,(p/2),((n-p-1)/2)); > m={1 4 2}; > s={ 1 .4 .16, > .4 1 .4, > .16 .4 1}; > u=ROOT(s); > x=J(n,1,1)*m + NORMAL(J(n,p,seed))*u; > end; > run; > quit; >