Discussion :

[soeursourire]

Salut!

J'ai utilisé une ACP pour rechercher ce qui varie le plus dans mes images de 100x100pixels. Donc j'ai utilisé princomp sur une matrice de 3150 par 10000. J'obtiens donc COEFF, LATENT et SCORE. Maintenant j'aurais voulu faire une rotation varimax sur mes données. J'ai donc essayé ca:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

COEFFnew = rotatefactors(COEFF);

Malheureusement j'ai ce message d'erreur:

??? Error using ==> svd
Input to SVD must not contain NaN or Inf.

Error in ==> rotatefactors>orthomax at 205
[L, D, M] = svd(A' * (d*B.^3 - gamma*B * diag(sum(B.^2))));

Error in ==> rotatefactors at 152
[B, T] = orthomax(A, 1, normalize, reltol, maxit);

Est-ce que je fais quelque chose de faux? Comment puis-je résoudre ce problème?
Merci d'avance pour votre aide

[Caro-Line]

Le message indique juste qu'il y a des NaN ou des Inf dans le calcul fait par SVD.
C'est donc que ta matrice d'entrée COEFF ne doit pas être "normale".
A quoi ressemble-t-elle ?
Y aurait-il des données inconnues (NaN ou Inf) dedans ?

[soeursourire]

Ben pourtant non a priori:
COEFF est une matrice de 10000x3149 (la sortie de princomp)
et si je fais:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
X=find(isnan(COEFF(:))==1);
Y=find(isinf(COEFF(:))==1);

X et Y sont vides donc il ne devrait pas y avoir de NaN ni de Inf...

[Caro-Line]

Comme c'est l'entrée de SVD qui pose problème il faudrait que tu fasses les calculs en pas à pas pour voir à quel moment ça dégénère.
Soit tu le fais en mode debug, soit tu écris en ligne de commande chaque ligne de la fonction ROTATEFACTORS pour voir où ça coince.

[Jerome Briot]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
X=find(isnan(COEFF(:))==1);
Y=find(isinf(COEFF(:))==1);

A noter que le test s'effectue plus proprement comme ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
X = any(isnan(COEFF(:)));
Y = any(isinf(COEFF(:)));

[soeursourire]

Si j'essaye sur une matrice aléatoire de la meme taille j'ai un problème de mémoire vous croyez que mon problème pourrait venir de la taille de ma matrice?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
C=rand(10000,3149);
Cb=rotatefactors(C);

??? Out of memory. Type HELP MEMORY for your options.

Error in ==> rotatefactors>orthomax at 205
[L, D, M] = svd(A' * (d*B.^3 - gamma*B * diag(sum(B.^2))));

Error in ==> rotatefactors at 152
[B, T] = orthomax(A, 1, normalize, reltol, maxit);

[soeursourire]

En faisant:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
X = any(isnan(COEFFs(:)));
Y = any(isinf(COEFFs(:)));

X=0 et Y=0