Discussion :

[gc001]

Bonjour à tous,

Je devrais calculer le nombre de bonnes réponse selon le stimulus présenté.

J'ai donc pour simplifier 2 colonnes : l'une (A) reprenant les bonnes ou mauvaises réponses (codées 1 ou 0) et une autre (B) correspondant au stimuli présenté. Cette colonne est un peu plus complexe car elle ne varie pas de 1 à nom nombre de présentation mais elle comprend des répétition du même chiffre (car présentation plusieurs fois du même stimulus). En effet, cette colonne est donc 4, 4, 4, 5, 5, 5, 6, 6,...

Le but serait donc de dire à Matlab : Pour le chiffre 4, tu comptes les bonnes réponses dans la colonne A, puis tu fais le chiffre 5, jusqu'à la fin. De ce fait, je voudrai faire une sommation partitionnée de la colonne A selon le type de stimulus considéré dans la colonne B.

Voici à quoi j'ai pensé :

Code :

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

1
2
3
4
5
 
rep_corr = [];    
    for ii = 4:size(B);    
        rep_corr = sum(A)      
    end

Je suis conscient des lacunes de ce code car lorsque je dis 4:size(B) Matlab risque de prendre en considération toutes les lignes et pas uniquement tout les lignes avec le même nom. Mon problème est que je parviens pas à lui faire prendre en considération uniquement les chiffres identiques dans la colonne B. De plus lorsque je regarde rep_corr, je n'ai qu'une seule ligne alors que je devrais en avoir autant qu'il y a de stimuli présenté.

Quelqu'un aurait un tuyau pour m'aider à avancer?
Je vous remercie d'avance

rq/ Matlab 7.5.0 (R2007B)

[duf42]

Bonjour,

Un exemple pour une valeur:

Code :

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

1
2
 
sum(A(B==3))

Pour définir la liste des valeurs de B tu peux utiliser la fonction UNIQUE et ainsi construire ta boucle.

Au passage je te conseille de lire ceci: FAQ: Indexage linéaire

Bonne apm,
Duf

[Jerome Briot]

C'est un travail pour la fonction HISTC :

Code :

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

1
2
3
4
5
6
7
8
A = [rand(1,8)>.5 ; 4 4 4 5 5 5 6 6].'
 
pA = prod(A,2);
uA = unique(A(:,2));
 
n = histc(pA,uA);
 
[uA n]

[gc001]

Ok, j'y vois un peu plus claire dans l'utilisation de la fonction UNIQUE. Vraiment impressionnant le raisonnement. J'étais en fait parti sur une décomposition de B par le nombre de stimuli en utilisant une boucle mais cela ne tournait toujours pas.

Merci bcp en tout cas! Je vais voir ce que cela donne.

[gc001]

Oui c'est bien cela, ça marche correctement.

Cependant, ce raisonnement n'est valable que si une de mes deux colonnes est composée de 0 et de 1. Dans le cas où on aurait des temps de réaction (colonne C par exemple), y a-t-il quelques choses de plus simple que :

Code :

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

1
2
3
4
 
TR4 = C(1)+C(2)+C(3);
TR5 = C(4)+C(5)+C(6);
TR6 = C(7)+C(8)+C(9);

En gros, je dis que pour connaitre le temps de réaction pour le stimulus 4, il faut ajouter les lignes 1, 2 et 3 de ma colonne C. Ainsi de suite pour les autres stimuli.

En tout cas, je vous remercie pour l'excellent tuyau d'au-dessus!

[Jerome Briot]

Alors il faut utiliser SPARSE :

Code :

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

1
2
3
4
5
A = [rand(1,8) ; 4 4 4 5 5 5 6 6].'
 
[a,b,c] = find(sparse(A(:,2),ones(size(A,1),1),A(:,1)));
 
[a c]

J'admet que ces solutions sont très MATLABiennes et pas vraiment très didactique pour des débutants... mais quand on a appris à se passer des boucles, c'est difficile de revenir en arrière.

[gc001]

Merci pour ta réponse.

Je parviens plus ou moins à comprendre la logique. Et c'est vrai que ça rend les choses plus courtes comparées aux boucles. Pour m'aider à comprendre, j'avais fait la décomposition du find en 2 :

Code :

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

1
2
3
 
X = sparse(A(:,2)
Y = ones(size(A,1),1),A(:,1))

Cependant, cela ne m'aide pas à voir où se fait la sommation (quelle fonction l'induit) car a aucun moment on n'indique sum. Cela rend donc la formule difficile à utiliser si on s'intéresse à un écart-type ou une moyenne.

Se sont des questions superflues mais c'est juste pour comprendre.
Merci d'avance pour les éclaircissements.

[Jerome Briot]

Cependant, cela ne m'aide pas à voir où se fait la sommation (quelle fonction l'induit) car a aucun moment on n'indique sum.

C'est dans la construction de la matrice creuse avec SPARSE

On peut lire dans la documentation de la fonction SPARSE :

S = sparse(i,j,s,m,n,nzmax) uses vectors i, j, and s to generate an m-by-n sparse matrix such that S(i(k),j(k)) = s(k), with space allocated for nzmax nonzeros. Vectors i, j, and s are all the same length. Any elements of s that are zero are ignored, along with the corresponding values of i and j. Any elements of s that have duplicate values of i and j are added together.

C'est très utile pour faire des assemblage de matrice en une seule fois

[gc001]

Oui effectivement! J'ai mal lu.

Y a-t-il d'autres fonctions qui permettent de faire des soustraction, faire une moyenne ou un écart-type? Ou bien le raisonnement dois totalement changer si on veut faire de telles choses?

Je voulais par ailleurs vous demander pourquoi lorsqu'on l'introduit la formule que vous avez proposé dans une boucle (je sais que vous n'aimez pas trop les boucles mais c'est la seule manière que j'ai de faire je pense), cela ne fonctionne pas. Il fait le calcul comme si le if n'existait pas.

Code :

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

1
2
3
4
if res_sorted(:,4) == 1;
[a,b,c] = find(sparse(matrix_RT(:,2),ones(size(matrix_RT,1),1),matrix_RT(:,1)));
rep = [a c]
end

L'idée ici, est d'additionner uniquement les réponses correctes codée par un 1. Pour seulement ensuite faire la moyenne.

Merci beaucoup.