SAS中文论坛

标题: sas产生多维正态分布数据的问题 [打印本页]

作者: shiyiming 时间: 2010-5-11 22:45
标题: sas产生多维正态分布数据的问题
大家好！
求助大家，我需要这样一个数据，即5维正态分布随机数
其中第一列数据服从N(33,10^2) ,第二列数据服从N(90,22^2),第三列N(110,56^2),第四列N(30,8^2),第五列N(10,14^2).
在sas中，怎么实现这个目地，恳求高人指教，十分感激！

作者: shiyiming 时间: 2010-5-12 09:08
标题: Re: sas产生多维正态分布数据的问题
要看你 5组数据之间是否有相关性，如果没有，简单一点，就单组模拟，合并就行了

作者: shiyiming 时间: 2010-5-12 23:14
标题: Re: sas产生多维正态分布数据的问题
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 时间: 2010-5-13 06:24
标题: Re: sas产生多维正态分布数据的问题
If there are not independent, use multivariate transformation, e.g., rotation.

David

作者: shiyiming 时间: 2010-5-13 09:06
标题: Re: sas产生多维正态分布数据的问题
to oakland4788: What's the mean?

作者: shiyiming 时间: 2010-5-13 11:57
标题: Re: sas产生多维正态分布数据的问题
各位大侠多谢，
顺便再问一下，如果已知其相关矩阵，怎么解决此问题？

作者: shiyiming 时间: 2010-5-13 19:58
标题: Re: sas产生多维正态分布数据的问题
*设(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中这方面的命令，但茆诗松的《高等数理统计》的第六章最后有这方面的算法，希望有用！

欢迎光临 SAS中文论坛 (http://www.mysas.net/forum/)