sas产生多维正态分布数据的问题

shiyiming

大家好！
求助大家，我需要这样一个数据，即5维正态分布随机数
其中第一列数据服从N(33,10^2) ,第二列数据服从N(90,22^2),第三列N(110,56^2),第四列N(30,8^2),第五列N(10,14^2).
在sas中，怎么实现这个目地，恳求高人指教，十分感激！

shiyiming

要看你 5组数据之间是否有相关性，如果没有，简单一点，就单组模拟，合并就行了

shiyiming

proc iml;
mu={33,90,110,30,10} ;
sigma={10,22,56,8,14};
xxx=normal(repeat(0,10,5)) ;
xxx=mu+sigma*xxx;
print xxx;
quit;
多多指教！

shiyiming

If there are not independent, use multivariate transformation, e.g., rotation.

David

shiyiming

to oakland4788:  What's the mean?

shiyiming

各位大侠多谢，
顺便再问一下，如果已知其相关矩阵，怎么解决此问题？

shiyiming

*设(X1,X2,X3)~N((1 2 3)',XX);
proc iml;
mu={1,2,3};
x={10 2 3,4 5 6,7 8 9};
xx=t(x)*x;
r=t(root(xx));
rr=r*normal(repeat(0,3,1)) +mu;
print rr;
quit;

我不清楚SAS中这方面的命令，但茆诗松的《高等数理统计》的第六章最后有这方面的算法，希望有用！