最近在學習多元統計分析的主成分分析時,發現需要經常計算矩陣的特徵值和特徵向量,自己就找了下用R來做計算的函數。
我們可以用sigen()函數來計算特徵對。
#創建一個矩陣
a <- matrix(c(11,sqrt(3)/2,3/2,sqrt(3)/2,21/4,5*sqrt(3)/4,3/2,5*sqrt(3)/4,31/4),
ncol=3,
nrow=3
)
#使用函數eigen()計算特徵對
ev <- eigen(a)
#訪問列表val項,即特徵值
ev$val
#訪問列表vec項,即特徵向量
ev$vec
Out:
a矩陣:
ev: