My codes also demonstrate that the first random number seed determines the random number stream in PROC IML using RANDGEN function.

proc iml;
ntot=j(20,10,.);
fill=j(20,1,.);
seed=1824;
do i=1 to 10;
call randseed(seed);
call randgen(fill,'normal',0,1);
ntot[,i]=fill;
end;
print ntot;
quit;

proc iml;
ntot1=j(20,10,.);
fill=j(20,1,.);
seed1=1822;
do i=1 to 10;
seed1=seed1+2*i;
call randseed(seed1);
call randgen(fill,'normal',0,1);
ntot1[,i]=fill;
end;
print ntot1;