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

MATLABpca svd nan inf

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

Best Answer

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!
Related Question