メイン kwtest.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="xxp.js">document.write("xxp.js ファイルが見つかりません??<br>")</script>
<script src="io.js">document.write("io.js ファイルが見つかりません??<br>")</script>

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

function compareNumbers(a, b) {
//  return a - b
  if (a > b) return 1
  else if (a == b) return 0
  else return -1
}

function kwtest(ng, n, key, g, x)
{
  var r = new Array(ng)
  var num = new Array(ng)
  for (i = 0; i < ng; i++) {
    r[i] = num[i] = 0
  }

  avg = (n+1)/2
  v = t = s = 0
  for (i = 0; i < n; i++) {
    k = g[i]
    num[k]++
    small = equal = 0
    for (j = 0; j < n; j++) {
      if (x[j] < x[i]) {
        small++
      }
      else if (x[j] == x[i]) {
        equal++
      }
    }
    rank = small+(equal+1)/2
    r[k] += rank
    v += Math.pow(rank-avg, 2)
    t += equal*equal-1
  }
  for (i = 0; i < ng; i++) {
    s += r[i]*r[i]/num[i]
    r[i] /= num[i]
  }
  v /= n-1
  s = 12*s/(n*n+n)-3*(n+1)
  c = 1-t/(n*n*n-n)
  s /= c

  printf("         %10s  %5s\n", "コード", "ケース数")
  for (i = 0; i < ng; i++) {
    printf("第%2i 群  %10s  %5i\n", i+1, key[i], num[i])
  }
  printf("\n検定統計量 S = %g%s\n", s, t == 0 ? "" : " (同順位の修正済み)")
  printf("P 値 = %g\n\n", xxp(s, ng-1))

  printf("対比較\n%14s %15s %10s\n", "比較対象", "統計量", "P 値")
  for (i = 0; i < ng-1; i++) {
    for (j = i+1; j < ng; j++) {
      s = Math.pow(r[i]-r[j], 2)/(1/num[i]+1/num[j])/v
      printf("第%2i 群 vs. 第%2i 群 %10.6g %10.6g\n", i+1, j+1, s, xxp(s, ng-1))
    }
  }
}

function calc(data_string)
{
  var data, nc, i, j, ng = 0
  var group = new Array()
  var datax = new Array()
  var key = new Array()
  printf("★ クラスカル・ウォリス検定 ★\n\n")
  if ((data = getdata(data_string, 2)) == false) return
  nc = data.length
  for (i = 0; i < nc; i++) {
    key[ng++] = group[i] = data[i][0]
    for (j = 0; j < ng; j++) {
      if (key[j] == group[i]) {
        if (j < ng-1) {
          ng--
        }
        break
      }
    }
    datax[i] = data[i][1]
  }
  if ((key.length = ng) == 1) {
    printf("群の数が1です。検定できません。\n")
    return
  }
  key.sort(compareNumbers)
  for (i = 0; i < nc; i++) {
    for (j = 0; j < ng; j++) {
      if (group[i] == key[j]) {
        group[i] = j
        break
      }
    }
  }
  kwtest(ng, nc, key, group, datax)
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>生データに基づくクラスカル・ウォリス検定</b></font> <a  href="src/kwtest.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>
<table><tr>
<td><input type="button" name="calcurate" value="計算開始" onClick="calc(this.form.data.value)">  </td>
<td><input type="button" name="clear" value="入力欄クリア" onClick="this.form.data.value=''">  </td>
<td><input type="button" name="clear" value="出力欄クリア" onClick="this.form.result.value=''">  </td>
<td><img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/kwtest.html">使用法</a></td>
</tr></table>

<p>
<table><tr>
<td>入力欄<br><textarea name="data" ROWS=25 COLS=15></textarea></td>
<td>出力欄<br><textarea name="result" ROWS=25 COLS=80></textarea></td>
</tr>
<tr><td colspan=2>入力欄には,群変数と分析対象変数を1行に一対ずつの値を<a href="exa/kugirimoji.html">区切り文字</a>で区切って入力</td>
</tr></table>
</form>

<p><hr noshade>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Average/kwtest.html">手法の解説</a><br>
<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
サブ xxp.js   Last modified: Mar 25, 2004
サブ io.js   Last modified: Mar 25, 2004

・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI

Made with Macintosh