Discussion :

[Daniel.R]

Bonjour à tous,
j'ai un tableau m*n comportant pour chaque case soit un 0 soit un 1.
J'aimerai comparé le profil de mes lignes entre elles, et lorsque certaines lignes ont le même profil, les extraire de mon tableau (en les mettant dans un autre tableau, ainsi j'aurai autant de "sous-tableau" que de profils différents).

Je me suis dit qu'il faudrait prendre la première ligne, la comparer avec les autres, l'extraire elle et ses petites copines, et recommencer l'opération (avec la seconde ligne qui sera passé du coup en tête de tableau) jusqu'à ce que mon tableau initiale soit vide. (s'il y a d'autres "plan" à suivre, je suis aussi preneur )

Soit A mon tableau :
Soit C et row mon nombre de ligne

Code :

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

1
2
3
4
5
6
7
%Pour comparer le profil :
for C = 1:size(A,1)
      Cligne{C} = A(C,:);
   	C_ref{C} = Cligne{C}
   	for row = 2:size(A,1)
   	   Rligne{row} = A(row,:);
   	   corrsp{row} = comp_ref{C} & Rligne{row}

c'est à partir de là que je suis plutôt perdu.

-> De un ma première comparaison/correspondance me donne [], puisque je pars sur

Code :

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

for row = 2:size(A,1)

(en effet je ne vais pas comparer ma première ligne avec ma première ligne. Ou alors je supprime cette comparaison par la suite ?).

-> De deux, je voudrais créer un nouveau tableau avec les lignes identiques. J'ai donc fait :

Code :

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

Nmatrice.reseau{C} = [comp_ref{C},comp_compare]

mais il semble que cela ne marche pas :

Code :

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

1
2
3
4
" ??? Conversion to cell from double is not possible.
 
Error in ==> C:\MATLABR11\toolbox\SPM99\danasupprm.m
On line 19  ==>       Nmatrice.reseau{C} = [comp_ref{C},comp_compare] "

-> De 3, je ne sais pas comment supprimer les lignes "à profil identiques" sans que ca fasse planter le bazar : une fois certaines lignes extraites, ma size(A,2) n'est plus la même. Il faudrait une sorte de mise à jour avant le premier "for"

Si quelqu'un a des conseils pour me faire avancer dans ce truc...
merci d'avance !
Dan

[Jean Dumoncel]

Bonjour,

j'ai un tableau m*n comportant pour chaque case soit un 0 soit un 1.
J'aimerai comparé le profil de mes lignes entre elles, et lorsque certaines lignes ont le même profil, les extraire de mon tableau (en les mettant dans un autre tableau, ainsi j'aurai autant de "sous-tableau" que de profils différents).

Le même profil, cela veut dire que tous les éléments de deux lignes sont égaux? Si oui, ceci peut t'intéresser :

Code :

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

1
2
help unique
[b, m, n] = unique(A, 'rows')

[Daniel.R]

merci.
Ce n'est pas tout à fait cela, mais ce n'est pas grave puisque je souhaite rajouter un truc :
j'aimerai calculer le coefficient de corrélation entre chacune de mes lignes (ce que je n'arrive pas à faire d'ailleurs), et en fonction du résultat extraires ces fameuses lignes de ma matrice...

[Daniel.R]

Pour préciser un peu :

j'ai un tableau m*n comportant pour chaque case soit un 0 soit une valeur différente de 0 (en fait soit <-1.6 soit >1.6).
Chacune de mes m lignes est une observation de n valeurs (ce n'est pas une mesure répétée, ce sont en l'occurrence des niveaux d'activité pour différentes aires cérébrales).
J'aimerai savoir à quel point le profil de mes lignes se ressemblent en calculant un coefficient de corrélation (des différentes lignes entre elles). (c'est-à-dire savoir à quel point se ressemblent mes différentes carte d'activité (mes observations / lignes)
Malheureusement, j'ai beau avoir chercher sur matlab, ce site, google, etc..., je ne sais comment écrire mes lignes de commandes pour avoir mon résultats souhaité.
Si quelqu'un a une piste...
d'avance merci

[Jean Dumoncel]

Code :

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

help corrcoef

Par contre je ne comprends pas bien pourquoi tu utilises autant de tableau de cellules... Pourrais tu nous montrer le code complet que tu as codé jusqu'à maintenant? (même si il ne fonctionne pas)

[Daniel.R]

Salut,
merci pour la fusion (je pensais qu'à question différente, sujet différent, mais c'est noté ).

Alors, bon, je suis stupide : je cherchais à ce qu'il me calcul les coefficient de corrélation en une fois(de toute ma matrice). Finalement, il suffit de lui donner l'ordre de le faire ligne par ligne... désolé, mais j'y connais rien en programmation et j'apprends en autodidacte (donc par essai/erreur) ^^

[Daniel.R]

bonsoir à tous,
alors j'ai réussi à faire ce que je voulais (enfin je crois), mais il y a une chose que je n'arrive pas :
je voudrais qu'à chaque extraction de lignes (le nombre de lignes extraites pouvant varier à chaque fois que la boucle reprend) de ma matrice initiale, il me sauvegarde ma sous-matrice sous "fichiers1" puis "fichiers2", puis "fichiersn" jusqu'à ce que ma boucle s'arrêtent.
Voici mon code :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
groupe.reseau_s_p : [];
groupe.nom = [];

% corrélation des réseaux seuillés et exprimés en %
	for C = 1:inf
   	if size(matrice.reseau_s_p,1) == 0
   	break
 	   else indy = zeros(size(matrice.reseau_s_p,1),1);
   	comp_ref = matrice.reseau_s_p(1,:)
   	for row = 1:size(matrice.reseau_s_p,1)
   		Rligne{row} = matrice.reseau_s_p(row,:);
      	coeff = corrcoef(comp_ref,Rligne{row})
      	if coeff(1,2) >= 0.75  %0.17334537 %16,515162
      	indy(C) = 1;
      	indy(row) = 1;
      	comp_compare = Rligne{row}
      	elseif coeff(1,2) <= -0.75
      	indy(C) = 1;
      	indy(row) = 1;
      	comp_compare{row} = Rligne{row}           
      	else
      	end
   	end 
   	groupe.reseau_s_p = matrice.reseau_s_p(find(indy > 0),:);
   	groupe.nom = matrice.nom(find(indy > 0),:);
   	matrice.reseau_s_p(find(indy > 0),:) = [];
   	matrice.nom(find(indy > 0,:) = [];
   	end  
   	cd C:\mat
   	if size(matrice.reseau_s_p,1) > 0
   	   indy = zeros(size(matrice.reseau_s_p,1),1);
   	else
   	end
   	save (?????, groupe.reseau_s_p, '-ASCII')
   	save(?????,groupe.nom, '-ASCII')
   	groupe.reseau_s_p = [];
   	groupe.reseau_s_p = [];
   end

plus qu'à sauvegarder et c'est fini, donc si quelqu'un a une idée... merci d'avance et encore merci à ceux qui m'ont aidé jusqu'à présent.

Dan

[Jerome Briot]

Pense à indenter correctement ton code... sous l'éditeur de MATLAB, tu fais Ctrl+A puis Ctrl+I

Voici ce qu ça donne :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
groupe.reseau_s_p : [];
groupe.nom = [];
 
% corrélation des réseaux seuillés et exprimés en %
for C = 1:inf
    if size(matrice.reseau_s_p,1) == 0
        break
    else indy = zeros(size(matrice.reseau_s_p,1),1);
        comp_ref = matrice.reseau_s_p(1,:)
        for row = 1:size(matrice.reseau_s_p,1)
            Rligne{row} = matrice.reseau_s_p(row,:);
            coeff = corrcoef(comp_ref,Rligne{row})
            if coeff(1,2) >= 0.75  %0.17334537 %16,515162
                indy(C) = 1;
                indy(row) = 1;
                comp_compare = Rligne{row}
            elseif coeff(1,2) <= -0.75
                indy(C) = 1;
                indy(row) = 1;
                comp_compare{row} = Rligne{row}
            else
            end
        end
        groupe.reseau_s_p = matrice.reseau_s_p(find(indy > 0),:);
        groupe.nom = matrice.nom(find(indy > 0),:);
        matrice.reseau_s_p(find(indy > 0),:) = [];
        matrice.nom(find(indy > 0,:) = [];
    end
    cd C:\mat
    if size(matrice.reseau_s_p,1) > 0
        indy = zeros(size(matrice.reseau_s_p,1),1);
    else
    end
    save (?????, groupe.reseau_s_p, '-ASCII')
    save(?????,groupe.nom, '-ASCII')
    groupe.reseau_s_p = [];
    groupe.reseau_s_p = [];
end

Maintenant pour répondre à ta question, essaie ceci :

Code :

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

1
2
filename = sprintf('fichier%d.txt',C);
save (filename, 'groupe.reseau_s_p', '-ASCII')

Enfin c'est l'idée...

Et pas besoin d'utiliser CD dans ton code.

Pour enregistrer un fichier dans un dossier, tu fais :

Code :

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

1
2
3
filename = sprintf('fichier%d.txt',C);
pathname = 'C:\mat'
save (fullfile(pathname,filename), 'groupe.reseau_s_p', '-ASCII')

Tu peux aussi te passer de la fonction FIND dans ton code

Code :

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

groupe.reseau_s_p = matrice.reseau_s_p(find(indy > 0),:);

devient simplement :

Code :

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

groupe.reseau_s_p = matrice.reseau_s_p(indy > 0,:);

C'est ce qu'on appelle l'indexage logique