Results of PCA function

Question

Tiago Dias 2018-6-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/407641-results-of-pca-function

评论： Tiago Dias 2018-6-28

X_imputed.mat

Hello, I got some trouble understanding why I get a different result than the PCA function that is incorporated on Matlab. I have uploaded the data that i worked on to get the results

load('X_imputed');
[coeff,score,latent] = pca(X_imputed,'algorithm','eig');

If I got it right, coeff is the Loadings.

I tried to replicate the results on my own, I got the same absolute value that coeff and score, but the signals are different, and I would like to know what I did wrong? or perhaps in my calculation there a sign that is wrong, because i got the absolute valuthe same.

meanX = mean(X_imputed);
Z = (X_imputed-ones(75,1)*mean(X_imputed)); % Centering
covZ = cov(Z); % Covariance
[vec_p,lambda_aux] = eig(covZ);
val_p = diag(lambda_aux);
[lambda,ind_ord] = sort(val_p,'descend');
L1 = vec_p(:,ind_ord); % Loadings
T1 = Z * L1; % Scores

Thanks for your time and help.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

David Goodmanson 2018-6-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/407641-results-of-pca-function#answer_326553

编辑：David Goodmanson 2018-6-27

Hi Tiago,

It's good to verify things as you are doing, and I don't believe that there is any problem. The eigenvalues agree. The columns of the coeff matrix (your L1) are normalized eigenvectors with real coefficients. Assuming they should stay real, eigenvectors are still undetermined within a factor of +-1. In this case, coeff and L1 differ by a factor of -1 in columns 2,5,6,7,8. Nothing wrong with that. Accordingly, score and your T1 differ by a factor of -1 in those same columns.