topic creating replicate/ sample ID in PROC IML in SAS/IML Software and Matrix Computations

creating replicate/ sample ID in PROC IML

SWEETSAS — Sat, 12 Sep 2015 12:11:44 GMT

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;

Re: creating replicate/ sample ID in PROC IML

Rick_SAS — Sat, 12 Sep 2015 23:40:18 GMT

See @IanWakeling and my response at https://communities.sas.com/t5/SAS-IML-Software-and-Matrix/Generating-exponential-and-gamma-distribution/td-p/215439

Re: creating replicate/ sample ID in PROC IML

SWEETSAS — Mon, 14 Sep 2015 01:45:39 GMT

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

Best,