¡¡°Ê²¼¤ò¥·¥ó¥¿¥Ã¥¯¥¹¥¨¥Ç¥£¥¿¡¼¤Ë¤Ï¤Ã¤Æ¡¢¡Ö¤¹¤Ù¤Æ¡×¤Ç¼Â¹Ô¤·¤Æ¤ß¤Æ
¤¤¤¿¤À¤±¤Þ¤¹¤«¡©HP¤ÈƱ¤¸¤â¤Î¤Ç¤¹¤¬¡¢¤¦¤Á¤Ç¤Ï·ë²Ì¤Þ¤Ç½Ð¤Þ¤·¤¿¡£
¡¡ÄÉ¿¡§¤À¤ó¤Ü¤µ¤ó¤Î´Ä¶¤À¤È¤É¤¦¤Ç¤¹¤«¡©
* Read in the IxI matrix of counts for I-level scale.
DATA LIST LIST / x1 to x5.
BEGIN DATA
25 5 3 2 0
12 18 8 3 1
0 3 14 4 2
1 3 3 19 4
2 1 5 9 20
END DATA.
* Commands to compute weighted Kappa .
matrix.
GET x /var=x1 to x5.
compute wt1=make(nrow(x),ncol(x),0).
compute wt2=make(nrow(x),ncol(x),0).
compute prop=x/msum(x).
loop i=1 to nrow(x).
loop j=1 to ncol(x).
compute wt1(i,j)=1-(abs(i-j)/(nrow(x)-1)).
compute wt2(i,j)=1-((i-j)/(nrow(x)-1))**2.
end loop.
end loop.
compute wk1num=msum(wt1&*prop)-msum(mdiag(rsum(prop))*wt1*mdiag(csum(prop))).
compute wk1den=1-msum(mdiag(rsum(prop))*wt1*mdiag(csum(prop))).
compute wk1=wk1num/wk1den.
compute wk2num=msum(wt2&*prop)-msum(mdiag(rsum(prop))*wt2*mdiag(csum(prop))).
compute wk2den=1-msum(mdiag(rsum(prop))*wt2*mdiag(csum(prop))).
compute wk2=wk2num/wk2den.
print wk1.
print wk2.
end matrix.