Question

我有一些带有3个坐标的数据点并且使用PCA函数我通过这样做将其转换为具有2个坐标的点

import numpy as np
from sklearn.decomposition import PCA
X = np.array([[-1, -1, -3], [-2, -1, -1], [-3, -2, -2], [1, 1, 1], [2, 1, 5], [3, 2, 6]]) #data
pca = PCA(n_components=2)
pca.fit(X)
PCA(copy=True, n_components=2, whiten=False)
XT = pca.fit_transform(X)
print XT
#output obtained
#[[-4.04510516 -1.24556106]
#[-2.92607624  0.61239898]
#[-4.55000611  1.13825234]
#[ 0.81687144 -1.11632484]
#[ 4.5401931   0.56854397]
#[ 6.16412297  0.04269061]]

我在特征空间中得到主轴，使用'components_'属性表示数据中的最大方差方向

W = (pca.components_)
print W
# output obtained
#[[ 0.49508794  0.3217835   0.80705843]
# [-0.67701709 -0.43930775  0.59047148]]

现在我想通过这样做使用'components_'属性将第一个点[-1，-1，-3]（这是X中的第一个点）投影到2D子空间

projectedXT_0 = np.dot(W,X[0])
print projectedXT_0
#output obtained
#[-3.23804673 -0.65508959]

#expected output
#[-4.04510516 -1.24556106]

我没有达到我的预期，显然我在使用'components_'属性计算projectionPoint时做错了。请证明使用'components_'属性来获得一个点的投影。

注意：我知道'transform'功能会这样做，但我想使用'components_'属性。

Answer 1

你忘了减去平均值。

查看pca transform的来源：

if self.mean_ is not None:
    X = X - self.mean_
X_transformed = fast_dot(X, self.components_.T)
if self.whiten:
    X_transformed /= np.sqrt(self.explained_variance_)
return X_transformed

如何使用sklearn.decomposition包中的PCA的“components_”属性将新点投射到新的基础？

1 个答案: