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은 고도로 최적화 된 동시에 강력하고 상당히 기능이 풍부합니다.