R中的全奇异值分解

Question

在大多数应用程序（尤其是统计应用程序）中，thin SVD就足够了。然而，有时需要完整的SVD以获得矩阵（及其共轭）的零空间的正交。似乎R中的svd（）只返回精简版本。是否可以生产完整版？还有其他选择吗？

Answer 1

library(sos)
> findFn("svd NULL space")
found 47 matches;  retrieving 3 pages

这看起来很重要：

MSBVAR null.space Find the null space of a matrix

As does this function in MASS

Answer 2

R Core使用Linpack，Lapack等所需的例程。

如果您需要不同的东西，您可能需要获得其他Linpack等例程，或连接到提供更多的库。

Doug Bates刚刚将Eigen文件夹包裹在RcppEigen包中，这可能会为您提供一些帮助。 Eigen在高度优化的同时似乎既强大又功能强大。