1. 程式人生 > >協方差矩陣和相關係數矩陣(R語言)

協方差矩陣和相關係數矩陣(R語言)

一、協方差矩陣

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