Discussion :

[soeursourire]

Bonjour,

On m'a conseillé d'expliquer mon problème dans cette rubrique alors voici...

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 observations par 10000. J'obtiens donc COEFF, LATENT et SCORE. Maintenant j'aurais voulu faire une rotation varimax sur mes données. Donc première question:

1. Il est juste de faire une PCA spatial puis une rotation varimax mais dans ce cas il faut bien faire une rotation sur les loadings (COEFF) ET sur les scores n'est ce pas? puis ensuite la PCA temporelle sur les scores qui ont subis la rotation? Et pour finir une rotation sur les SCORE et les COEFF de cette PCA temporelle?

2. J'ai donc essayé ca:

Code :

Code :

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

1
2
COEFFnew = rotatefactors(COEFF);
SCOREnew = rotatefactors(SCORE);

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

[Jerome Briot]

Cette discussion fait suite à celle-ci : Matlab et varimax

On m'a conseillé d'expliquer mon problème dans cette rubrique alors voici...

Devinez qui ?

[FR119492]

Salut!
Dut t'a dit de passer sur le forum Algo et il a tout à fait raison. Mais maintenant, tu ne sembles pas t'être rendu compte que tu n'étais plus sur le forum Matlab. Tu me sembles présenter un problème de méthodologie que l'on rencontre très souvent sur ce site: tu cherches à savoir comment programmer avant de te demander ce que tu dois programmer. Tu as un problème (de traitement d'images, si j'ai bien compris). Commence par déterminer quels algorithmes tu dois appliquer à tes données, et quand tu seras absolument certain d'avoir achevé cette étape. Ensuite, tu pourras commencer à écrire ton programme, en évitant d'utiliser des "boîtes noires" (c'est peut-être plus rapide, mais on risque de ne pas comprendre ce qu'on fait).
Jean-Marc Blanc

[soeursourire]

Ben pourtant ca me parait clair dans ma tete ce que je veux faire... J'ai tout un tas d'images qui ont été enregistrées sous 3 conditions differentes (50 images par observation). Je veux donc
1. Faire une ACP spatiale sur toutes mes observations avec mes pixels comme variables pour voir OU je trouve des variations dans l'image (ce serait bien d'avoir des zones differentes selon les conditions). -> LOADINGS + SCORES
2. Faire une rotation varimax sur les résultats de ces données
3. Faire une ACP temporelle sur les résultats de l'ACP précédentes pour voir QUAND ces variations se passent dans le temps ->LOADINGS + SCORES
4. De nouveau une rotation varimax sur les resultats pour que ce soit plus clair..

[FR119492]

Salut!

pourtant ca me parait clair dans ma tete

Pour moi, ce n'est pas du tout clair. J'explicite pourquoi:

100x100pixels

une matrice de 3150 observations par 10000

D'abord, ce n'est pas une matrice 3150x10000 mais un tableau, semble-t-il plutôt 3150x100x100.

J'ai utilisé une ACP

Comment, par quel algorithme?

COEFF, LATENT et SCORE

Comment as-tu défini et calculé ces trois variables?

faire une rotation varimax

Quelles sont les formules qui définissent cette rotation?

??? 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);

Tu sembles avoir utilisé la fonction rotatefactors sans trop te demander ce qu'il y avait dedans. Or les messages d'erreur montrent de manière évidente que cette fonction fait appel à la fonction svd (valeurs singulières). C'est donc à toi de décortiquer rotatefactor pour voir pourquoi la matrice à laquelle tu appliques svd est erronée.

Jean-Marc Blanc

[soeursourire]

Effectivement ce n'etait pas clair..
Alors mes images font 100pixels de hauteur par 100pixels de largeur mais je les transforme en vecteur colonne. Et oui j'ai donc un tableau de 3150 par 10000.
J'utilise princomp pour calculer mon ACP et COEFF SCORE et LATENT sont les resultats de cette ACP.
J'utilise bien rotatefactors qui par defaut utilise la rotation varimax. Et svd y'a pas le code avec dans Matlab apparemment.. ce qui est un peu embêtant pour traquer ce qui ce passe dans la fonction non?