Random generation with a given cov structure(Cholesky decomposition)

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

Note: some features may need SAS 9.2  IMLplus

matrix A is a given COV.  x1, x2, x3 (col1 col2 col3) are randomly generated having the given COV.

proc iml;
A = {25  0  5,
      0  4  6,
      5  6 59};
U = root(A);
print U;
call randseed(2);
x1=j(500000,1);
x2=x1;
x3=x1;
call randgen(x1, 'normal');
call randgen(x2, 'normal');
call randgen(x3, 'normal');
z=(x1||x2||x3)*U;
para={'col1', 'col2','col3'};
para={'col1' 'col2' 'col3'};
create t from z;
append from z;
close t;
submit para;
proc corr data=t cov;
var &para ;
run;
endsubmit;

quit;

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：composition Generation Structure Cholesky position Structure random Generation Cholesky cov