全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
2003 0
2015-06-29
求大神解惑。
例题:生成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;
程序不知道哪里出错了,运行不出结果,求解答,谢谢。

二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群