Question

我使用sklearn.decomposition.PCA

执行SVD

从SVD的等式

A = U x S x V_t

V_t = V的转置矩阵（对不起，我无法粘贴原始等式）

如果我想要矩阵U，S和V，如果我使用sklearn.decomposition.PCA怎么能得到它？

Answer 1

首先，根据矩阵的大小，PCA的sklearn实现并不总是计算完整的SVD分解。以下摘自PCA's GitHub reciprocity：

svd_solver : string {'auto', 'full', 'arpack', 'randomized'}
        auto :
            the solver is selected by a default policy based on `X.shape` and
            `n_components`: if the input data is larger than 500x500 and the
            number of components to extract is lower than 80% of the smallest
            dimension of the data, then the more efficient 'randomized'
            method is enabled. Otherwise the exact full SVD is computed and
            optionally truncated afterwards.
        full :
            run exact full SVD calling the standard LAPACK solver via
            `scipy.linalg.svd` and select the components by postprocessing
        arpack :
            run SVD truncated to n_components calling ARPACK solver via
            `scipy.sparse.linalg.svds`. It requires strictly
            0 < n_components < X.shape[1]
        randomized :
            run randomized SVD by the method of Halko et al.

此外，它还会对数据执行一些操作（请参阅here）。

现在，如果您想获得U, S, V中使用的sklearn.decomposition.PCA，可以使用pca._fit(X)。例如：

from sklearn.decomposition import PCA
X = np.array([[1, 2], [3,5], [8,10], [-1, 1], [5,6]])
pca = PCA(n_components=2)
pca._fit(X)

打印

(array([[ -3.55731195e-01,   5.05615563e-01],
        [  2.88830295e-04,  -3.68261259e-01],
        [  7.10884729e-01,  -2.74708608e-01],
        [ -5.68187889e-01,  -4.43103380e-01],
        [  2.12745524e-01,   5.80457684e-01]]),
 array([ 9.950385  ,  0.76800941]),
 array([[ 0.69988535,  0.71425521],
        [ 0.71425521, -0.69988535]]))

但是，如果您只想对原始数据进行SVD分解，我建议使用scipy.linalg.svd

通过sklearn.decomposition.PCA执行svd，如何从中获取U S V？

1 个答案: