相关系数与协方差间的转换

首先,相关系数和协方差间的数学关系为:

其次,上述数学表达式在sas/iml中的算法表达如下,R为相关系数矩阵,S为对角阵, COV协方差矩阵

其中:

 

最后,上代码:

***** Program 4.10 PROC IML ***********;
***** Converting a correlation matrix to covariance matrix, and vice versa;
***** Example data from Table 4.4 ;
****** Part I: Converting correlation matrix to covariance matrix;
PROC IML;

*** define the correlation matrix;
R={1.00 0.70 0.20,
   0.70 1.00 0.40,
   0.20 0.40 1.00};

*** define the diagonal matrix with standard deviations on the diagonal;
S={15 0 0,
    0 10 0,
    0 0 1};

*** obtain the covariance matrix;	
COV=S*R*S;  

*** print the covariance matrix;
PRINT COV; 
RUN;

***** Part II: Converting covariance matrix to correlation matrix;
PROC IML;

*** define the covariance matrix;
COV={225 105 3,
     105 100 4,
       3   4 1};
	   
*** obtain the matrix with standard deviations on the diagonal;
S=SQRT(DIAG(COV));

*** the inverse of S matrix;
S_INV=INV(S); 

*** obtain correlation matrix;
R=S_INV*COV*S_INV; 

** print out the three matrices;
PRINT COV S R;
RUN;


FROM 《SAS for Monte Carlo Studies》,P96-97

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值