BookmarkSubscribeRSS Feed
SWEETSAS
Obsidian | Level 7

The following program works perfectly well; it gives me what I want for just one sample. But I want multiple samples. However, in data step, I can easily create sample ID  using a do "rep=1 to m" statement and analyze by that sample ID. By sample ID, I mean replicate. How do I generate the sample ID (replicate) in IML so that I can analyze my data by sample ID?

 

Thanks in advance:

 

proc iml;
call randseed(1);

sigma = 0.5; /* scale parameter */
E = j(1000, 1);
call randgen(E, "Exponential", sigma); /* E ~ Expo(0.5) */
G = j(1000, 1);
call randgen(G, "Gamma", 2, sigma); /* G ~ Gamma(2, 0.5) */

 

create test var {E G}; append; close;

quit;

2 REPLIES 2
SWEETSAS
Obsidian | Level 7

Great! Thanks everyone. I see it's in page 69, section 4.5.2 (Reshaping Matrices) of "Simulating Data with SAS"

 

Best,

 

J