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核心使用它需要的Linpack，Lapack等例程。

如果您需要不同的東西，您可能需要讓自己其他Linpack等例程，或連接到提供更多的庫。

Doug Bates將Eigen庫包裝在RcppEigen包中，這可能對您有些幫助。 Eigen似乎既強大又相當有特色，同時經過高度優化。