协方差矩阵和相关系数矩阵(R语言)

版权声明:本文为博主原创文章,未经博主允许不得转载哦!!! https://blog.csdn.net/lph188/article/details/84501481

一、协方差矩阵

1.协方差定义

                                          S = \frac{1}{n - 1}\sum_{i=1}^{n}\left ( x_{(i)} -\bar{x}\right )\left ( x_{(i)} -\bar{x}\right )^{T}

 2.R语言实现

#协方差矩阵
#等价于cov(data)

data<-as.matrix(data)
n<-nrow(data)
mx<-diag(1,n)-matrix(1,n,n)/n
covA<-t(data)%*%mx%*%data/(n-1);covA 

一、相关系数矩阵

1.相关系数定义

                                  r_{jk} = \frac{s_{jk}}{\sqrt{s_{jj}}\sqrt{s_{kk}}}, j,k = 1,2,...,p

2.R语言实现

#相关系数矩阵
#等价于cor(data)

corA<-covA
for(j in 1:n )
{
  for(k in 1:n)
  {
    corA[j,k]<-covA[j,k]/(sqrt(covA[j,j])*sqrt(covA[k,k]))
  }
}
corA

相关资料:数据分析与R软件第二版(李素兰著)每章例题代码和数据,下载地址:https://download.csdn.net/download/lph188/10802159

猜你喜欢

转载自blog.csdn.net/lph188/article/details/84501481
今日推荐