R中的特征值和特征向量

Question

当我尝试使用eigen（）函数为以下矩阵计算特征值和特征向量时，与我的预期相反，最小特征值不为零，为什么？

 > lapMatrix
      [,1] [,2] [,3] [,4]
 [1,]  0.6 -0.1 -0.5  0.0
 [2,] -0.1  1.0  0.0 -0.9
 [3,] -0.5  0.0  1.0 -0.5
 [4,]  0.0 -0.9 -0.5  1.4

 > eigen(lapMatrix)
  $values
 [1] 2.277671e+00 1.241300e+00 4.810288e-01 1.297990e-16

 $vectors
            [,1]       [,2]       [,3] [,4]
 [1,]  0.1359786  0.5205445  0.6786333  0.5
 [2,] -0.5412623  0.3853145 -0.5554889  0.5
 [3,] -0.3480021 -0.7447136  0.2725734  0.5
 [4,]  0.7532858 -0.1611455 -0.3957179  0.5

Answer 1

您可以简单地将结果取整 ：

round(eigen(mm)$values,2)
## [1] 2.28 1.24 0.48 0.00

您可以将其包装在一个小函数中：

reigen <- 
  function(...,digits =2) round(eigen(...)$values,digits)

或使用zapsmall如注释中建议的那样）（它将自动确定位数）：

reigen <- 
  function(...,digits =2) zapsmall(eigen(...)$values)

reigen(mm)
[1] 2.278 1.241 0.481 0.000