首先,要把你的结果整理成如下矩阵的形式:
x{x11 x12 x13 x14 x15 x16
x21 x22 x23 x24 x25 x26
xn1 xn2 xn3 xn4 xn5 xn6},并且每次x 的矩阵都是不一样的
其次,用sas iml语言计算不同x矩阵的协方差,程序如下:
proc iml;
start corr(x);
x_sum=x[+,];
x_cov=(t(x)*x-t(x_sum)*x_sum/nrow(x))/(nrow(x)-1);print x_cov;
finish corr;
你输入不同的矩阵x,运行
run corr(x);
就得到不同的协方差。
注:如果对数据进行标化,协方差就是-1 到 1 的相关系数 r
[此贴子已经被作者于2009-3-4 14:43:32编辑过]
xmok77 金钱 +100 奖励 2009-3-4 18:57:32
xmok77 金钱 +100 奖励 2009-3-4 19:07:46