Weighted PCA throws error: Input to SVD must not contain NaN or Inf.

14 次查看（过去 30 天）
Daniel Ringel
Daniel Ringel 2017-8-9
评论： Daniel Ringel 2017-8-17
Hello, I started to use Matlab's built in PCA on a large data set in terms of variables (around 70 Observations but 100,000 variables in matrix M).
Using standard PCA I get results (N-1 PC's).
[coeff,scores,latent,~,explained] = pca(M);
When I run weighted PCA on the same data
[wcoeff,scores,latent,~,explained] = pca(M,'VariableWeights','variance');
I get the following error message:
Error using svd
Input to SVD must not contain NaN or Inf.
Error in pca>localSVD (line 477)
[U,sigma,coeff] = svd(x,'econ');
Error in pca (line 347)
[U,sigma, coeff, wasNaN] = localSVD(x, n,...
Error in ScriptA (line 62)
[wcoeff,scores,latent,~,explained] = pca(M,'VariableWeights','variance');
What am I doing wrong? Or what is the problem with matrix M that I need to resolve? Thank you, Daniel

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

采纳的回答

Faiz Gouri
Faiz Gouri 2017-8-17
Your input matrix probably contains infinity or not a number(NaN). Set a breakpoint in your file then run it again. When MATLAB reaches the breakpoint, look at the input matrix values. It will contain at least one element that is Inf, NaN, or -Inf. The SVD function doesn't know how to compute the singular value decomposition of a matrix with a nonfinite element. If you are updating matrix in a loop, the first iteration may not be the one that's throwing the error; in that case, or once you've determined that input matrix does actually contain a nonfinite value, set an error breakpoint to stop as soon as an Inf or NaN value is created. This will indicate where the nonfinite element of X is being introduced. Then figure out, based on your knowledge of your algorithm, how to prevent the nonfinite value from being introduced.
Hope this helps!
  1 个评论
Daniel Ringel
Daniel Ringel 2017-8-17
Thank your for your help. My initial input matrix contains no NaN or Inf values. So they have to be generated during the PCA processing. I wonder why...?
Are you suggesting I manipulate the built in Matlab PCA function by adding break points? How would I do that?

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Dimensionality Reduction and Feature Extraction 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by