シンタックスのいじる場所

In Reply to: Re:weighted kappaの算出方法
投稿者:オロナイン - 投稿日時:2002年11月20日 17時04分56秒

　こんにちは。cooさんの目的に役に立ちましたでしょうか？

　シンタックスについては僕も詳しくないのですが、知ってる範囲での
追加情報です。

　以下のシンタックスをいじるとすると、以下の２点だけだと思います。

①　Cooさんのデータをスペース区切りで作り
「BEGIN DATA 」「END DATA.」の間に入れます。

　例では、５×５でしたが、任意のサイズでOKです。

②列（変数）の数に合わせて、
「DATA LIST LIST / x1 to x5.」
「GET x /var=x1 to x5. 」
の２カ所を変える

　７つの変数なら、 x1 to x7. です。

③あとは、すべて実行でOK

　
　ポイントとしては、中間処理が見えた方が間違いが少な
いので、ところどころに「print」コマンドで、途中結果を
出力します。

　以下の例は、先のシンタックスに、途中経過の出力を加え
たものです。

　ご参考までに(^○^)

* 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.

print x.

compute wt1=make(nrow(x),ncol(x),0).
compute wt2=make(nrow(x),ncol(x),0).
compute prop=x/msum(x).

print wt1.
print wt2.
print prop.

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.

print wt1.
print wt2.

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.

print wk1num.
print wk1den.
print wk1.

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 wk2num.
print wk2den.
print wk2.

print wk1.
print wk2.
end matrix.

この投稿へのコメント