メイン　kappa.html　　　Last modified: Sep 01, 2009
<html>
<head>
<meta http-equiv="Content-Type" content="text/html;CHARSET=EUC-JP">
<link rel="shortcut icon" href="../favicon.ico">
<title>JavaScript</title>
<script src="gxp.js">document.write("gxp.js ファイルが見つかりません??<br>")</script>
<script src="io.js">document.write("io.js ファイルが見つかりません??<br>")</script>

<script language="JavaScript">
<!--

function outmat(r, nr, nc, str)
{
  var j, k
  printf("\n%s\n\n", str)
  for (j = 0; j < nr; j++) {
    for (k = 0; k < nc; k++) {
      printf("%5i ", r[j][k])
    }
    printf("\n")
  }
}

function calc(data_string)
{
  var data, nr, nc, i, j, minrc, maxrc, rsum, csum, n, sumrc, trace, kstat, n03, ase0, z, p, n12, n13, ase1
  var r = new Array()
  if ((data = getdata(data_string, 0)) != false) {
    nr = data.length
    nc = data[0].length
    n = 0
    for (i = 0; i < nr; i++) {
      r[i] = new Array(nc)
      for (j = 0; j < nc; j++) {
        n += r[i][j] = data[i][j]
      }
    }

    outmat(r, nr, nc, "★ 入力された分割表 ★")

    minrc = Math.min(nr, nc)
    maxrc = Math.max(nr, nc)
    rsum = new Array(maxrc)
    csum = new Array(maxrc)
    for (i = 0; i < maxrc; i++) {
      rsum[i]  = csum[i] = 0
    }
    for (i = 0; i < nr; i++) {
      for (j = 0; j < nc; j++) {
        rsum[i] += r[i][j]
        csum[j] += r[i][j]
      }
    }
    sumrc = trace = 0
    for (i = 0; i < minrc; i++) {
      sumrc += rsum[i]*csum[i]
      trace += r[i][i]
    }
    kstat = (n*trace-sumrc)/(n*n-sumrc)
    n03 = n12 = n13 = 0
    for (i = 0; i < minrc; i++) {
      n03 += rsum[i]*(rsum[i]+csum[i])*csum[i]
      n12 += r[i][i]*(rsum[i]+csum[i])
    }
    for (i = 0; i < nr; i++) {
      for (j = 0; j < nc; j++) {
        n13 += r[i][j]*Math.pow(rsum[j]+csum[i], 2)
      }
    }
    ase0 = Math.sqrt((n*n*sumrc+sumrc*sumrc-n03*n)/(n*Math.pow(n*n-sumrc, 2)))
    z = kstat/ase0
    p = gxp(Math.abs(z))*2
    ase1 = Math.sqrt(n*(trace*(n-trace)/Math.pow(n*n-sumrc, 2)+2*(n-trace)*(2*trace*sumrc-n*n12)/Math.pow(n*n-sumrc, 3)+Math.pow(n-trace, 2)*(n*n13-4*sumrc*sumrc)/Math.pow(n*n-sumrc, 4)))
    printf("\nカッパ統計量 κ = %g　標準誤差 = %g\nκ= 0 の検定　z = %g　P 値 = %g\n\n", kstat, ase0, z, p)
    printf("カッパ統計量の 95%% 信頼区間 [%g, %g]\n（信頼区間の計算に使われた標準誤差 = %g）\n", kstat-1.96*ase1, kstat+1.96*ase1, ase1)
  }
  sep(40)
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>κ統計量（正方行列でない場合も含む）</b></font>　<a  href="src/kappa.html"><img src="png/src.png" width=35 height=11 alt="src" align=top></a>　　　　　Last modified: Jun 01, 2006<hr noshade><p>
<font color="#ff0000" size="+2">以下のプログラムのサポートは終了しました。自己責任でお使い下さい。</font>

<form name=Result>
<script language="JavaScript">
<!--
//-->JavaScript がサポートされていないブラウザですか？
</script>
<input type="button" name="calcurate" value="計算開始" onClick="calc(this.form.data.value)">　　
<input type="button" name="clear" value="入力欄クリア" onClick="this.form.data.value=''">　　
<input type="button" name="clear" value="出力欄クリア" onClick="this.form.result.value=''">　　
<img src="../gra/button3.png" width=9 height=9 alt="・"> <A HREF="exa/kappa.html">使用法</a><p>
入力欄（分割表を<a href="exa/kugirimoji.html">区切り文字</a>で区切って，行単位に入力）<br><textarea name="data" ROWS=10 COLS=80></textarea><p>
出力欄<br><textarea name="result" ROWS=30 COLS=80></textarea>
</form>

<p><hr noshade>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <A HREF="javascript:history.go(-1)">直前のページへ戻る</A>　　<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../mail.html">E-mail to Shigenobu AOKI</a>
<p><center><IMG SRC="../gra/ume5.png" width=121 height=37 ALT="Made with Macintosh"></center>
</body>
</html>

サブ　gxp.js　　　Last modified: Mar 25, 2004
サブ　io.js　　　Last modified: Mar 25, 2004
・ 直前のページへ戻る　　・ E-mail to Shigenobu AOKI

Made with Macintosh