求大神解惑。
例题:生成1000组相关系数均为0.3的服从多元正态分布的随机向量(X,Y,Z),其中随机变量X~N(0,1),Y~N(10,22),Z~N(20,32)。程序如下: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};         /*定义均值向量*/
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;
程序不知道哪里出错了,运行不出结果,求解答,谢谢。