向高手请教：如何产生两列不相关的随机数？

10# wkn1986
是啊！不相关不就是相关系数为0吗。

data binormal;
a1=sqrt((1+0)/2);*此处的0可以控制相关系数;
a2=sqrt((1-0)/2);
do i=1 to 100;
rd1=rand('uniform');
rd2=rand('uniform');
e1=a1*rd1+a2*rd2;
e2=a1*rd1-a2*rd2;
output;
end;
keep e1 e2;
run;
proc corr data=binormal;
var e1 e2;
run;
运行一下就知道了

个人观点，假设你所需要产生的两个序列的相关系数矩阵为A，现在可以对A进行Cholesky decomposition成B，即有A=B'B,现在生成两个随机数列X1,X2,然后用X1*B,X2*B就是你所需要求解的最终两个序列。
proc iml;
a={1 0.5,0.5 1};
b=root(a);
c=rannor(repeat(10,100,2));
d=c*b;
print d;
quit;

12# andy162639
该程序得到的相关系数时高时低，不稳定。不知怎么回事？

13# cz851218
我对高手给出的程序修改为产生均匀分布的随机序列：
proc iml;
a={1 0,0 1};
b=root(a);
c=ranuni(repeat(10,100,2));
d=c*b;
print d;
quit;
得到的结果比较满意，相关系数为0.00460，Prob > |r| under H0: Rho=0，由P=0.9638接受H0。
万分感谢！

13# cz851218
有关该程序，我还有一个问题，通过PROC IML得到的两列随机数放在了变量D下，能否把这两列随机数直接输出到一个给定的数据集中？

13# cz851218
问题已解决，多谢！

tj0412ymy 发表于 2010-6-20 17:28
2# wkn1986
由于现要求随机数列是不相关的，所以它们之间的相关系数应该为0！如何才能办到？

What you need is to review your test book because you don't understand the result even it is perfect correct.

18# bobguy
YOU SAID IT. 相关系数为0，不一定意味着不相关，只能说明它们之间没有线性关系；而若相关系数不为0，只能说明它们之间存在着不同程度的线性关系！即相关系数只是针对线性关系而言的。但在很多文献上看到的所谓的“不相关”，也大概多数均指的是没有线性关系！

小白下，问个问题
两列随机数，相关系数为0的概率是多少啊？求高手回答......