zhangyanrong
下面多元正态分布的程序 ,其中,
xi=C`*(z1//z2//Z3)+u;
m=m//xi`;[/b]不知道其具体含义,尤其是在IML中//的含义,请大侠指教
程序如下:
proc IML;
R={1.00 0.30 0.30,
0.30 1.00 0.30,
0.30 0.30 1.00}; /*定义相关系数矩阵*/
S={1 0 0,
0 2 0,
0 0 3};
E=S*R*S; /*计算协方差矩阵*/
u={0,10,20}; /*定义均值向量*/
do i=1 to 1000;
z1=RANNOR(0);
z2=RANNOR(0);
z3=RANNOR(0);
C=root(E);
xi=C`*(z1//z2//Z3)+u;
m=m//xi`;
end; /*产生1000组多元正态分布的数据*/
create example4 var{x y z}; /*存储到数据集example4中*/
append from m;
run;