R

时间:2015-06-07 09:29:00

标签: r

当我尝试使用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

1 个答案:

答案 0 :(得分: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