r - Rでの完全特異値分解

Question

ほとんどのアプリケーション（特に統計アプリケーション）では、薄いSVDで十分です。ただし、行列（およびその共役）の零空間のオルソベーシスを取得するために、完全なSVDが必要になる場合があります。Rのsvd（）はシンバージョンのみを返すようです。フルバージョンを作成することは可能ですか？代替案はありますか？

Answer 1

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

これはポイントに見えます：

MSBVARnull.space行列の零空間を見つけます

この機能はMASSでも同様です。

Answer 2

R Coreは、Linpack、Lapackなどの必要なルーチンを使用します。

別の何かが必要な場合は、おそらく他のLinpackなどのルーチンを入手するか、より多くを提供するライブラリに接続する必要があります。

Doug Batesは、EigenライブラリをRcppEigenパッケージにラップしました。Eigenは、高度に最適化されている一方で、強力であり、かなり機能的であるように見えます。