最近在学习多元统计分析的主成分分析时,发现需要经常计算矩阵的特征值和特征向量,自己就找了下用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: