Discussion :

[eviasra]

bonjour;

je voudrai calculer l'indice kappa d'une matrice, mais j'ai une erreur qui me bloque 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
m=[50 0 0;0 49 1;0 1 49]
 
%indice kappa
%calculer P0
N=sum(sum(m));
PO=sum(diag(m))/N;
%------------------
%calculer k
%k=po-pe/1-pe
%pe=somme de pi*p.i/n²
for i=1:3
    for j=1:3
pe= (sum(m(i,:))*sum (m(:,j)))/(1-sum(m(:,j)));
 
    end
end

[phryte]

Bonsoir,
D'après :

https://fr.mathworks.com/matlabcentr...ontent/kappa.m

J'ai essayé cela :

Code :

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

1
2
3
4
5
6
7
8
9
n=sum(m(:)); %Sum of Matrix elements
m=m./n; %proportion
r=sum(m,2); %rows sum
s=sum(m); %columns sum
Ex=r*s; %expected proportion for random agree
f=eye(3,3);
po=sum(sum(m.*f));
pe=sum(sum(Ex.*f));
k=(po-pe)/(1-pe)

[eviasra]

merci phryte
me suis inspiré de ce code, je poste ma solution :

Code :

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

1
2
3
4
5
6
7
8
N=sum(sum(m));
po=sum(diag(m))/N;
pi=sum(m,2);% addition des lignes
p.i=sum(m);% addition des colonnes 
pii=p.i';
p1=sum(pi.*pii);
pe=p1/N^2;
k=(po-pe)/(1-pe);