求助，如何用SAS求样本数据的协方差阵？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

如题。。困惑中

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：协方差阵 样本数据 协方差 如何用 样本数 求助 SAS 计算 协方差阵

我想到的有2个方法：
方法一：应用corr过程，程序如下
DATA raw;
input x1 x2 x3 @@;
CARDS;
171.0      58.5      81.0            154.8      44.5      74.0
175.0      65.0      87.0            164.0      51.0      72.0
159.0      38.0      71.0            165.2      55.0      79.0
155.3      45.0      74.0            164.5      46.0      71.0
152.0      35.0      63.0            159.1      48.0      72.5
158.3      44.5      75.0            164.2      46.5      73.0
;
PROC CORR COV ;
RUN;

方法二：用iml直接计算，程序如下
proc iml;
start cov(x);
n=nrow(x);
x_sum=x[+,];
x_cov=(t(x)*x-t(x_sum)*x_sum/n)/(n-1);
print x_cov;
finish cov;
x={171.0      58.5      81.0,  
175.0      65.0      87.0,  
159.0      38.0      71.0,  
155.3      45.0      74.0,  
152.0      35.0      63.0,  
158.3      44.5      75.0,
154.8      44.5      74.0,
164.0      51.0      72.0,
165.2      55.0      79.0,
164.5      46.0      71.0,
159.1      48.0      72.5,
164.2      46.5      73.0};
run cov(x);
quit;
在以上的iml，我已经建立里cov模块，以后直接调用就行了。

仅供参考啊！~

楼上高手