メイン　kw-test.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 calc(data_string)
{
  var i, j, g, eq, lt, ord, sx = 0, c = 0, p
  var ns, R, tie, data, nc, i, ng = 0
  var group = new Array()
  var datax = new Array()
  var ns = new Array()
  var R = new Array()
  var tie = new Array()
  printf("★ クラスカル・ウォリス検定 ★\n\n")
  if ((data = getdata(data_string, 2)) == false) return
  nc = data.length
  for (i = 0; i < nc; i++) {
    j = group[i] = data[i][0]
    ns[j] = isNaN(ns[j]) ? 1 : ns[j]+1
    datax[i] = data[i][1]
  }
  printf("%8s  %14s  %10s\n", "群", "群変数の値", "データ個数")
  for (i in ns) {
    printf("%8s  %10s  %10i\n", "第"+(++ng)+"群", i, ns[i])
  }
  printf("\n")
  if (ng == 1) {
    printf("群の数が１です。検定できません。\n")
    return
  }
  for (i = 0; i < nc; i++) {
    eq = lt = 0
    for (j = 0; j < nc; j++) {
      if (datax[j] < datax[i]) {
        lt++
      }
      else if (datax[j] == datax[i]) {
        eq++
      }
    }
    ord = lt+(eq+1)/2
    j = group[i]
    R[j] = isNaN(R[j]) ? ord : R[j]+ord
    tie[ord] = isNaN(tie[ord]) ? 1 : tie[ord]+1
  }
  for (i in R) {
    sx += R[i]*R[i]/ns[i]
  }
  for (i in tie) {
    c += Math.pow(tie[i], 3)-tie[i]
  }
  p = xxp(sx = (12*sx/(nc*nc+nc)-3*(nc+1)) / (1-c/(nc*nc*nc-nc)), ng-1)
  printf("検定統計量 = %g　（自由度 %i のカイ二乗分布に従う）\n      P 値 = %.5f\n", sx, ng-1, p)
  sep2(60)
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>生データに基づくクラスカル・ウォリス検定</b></font>　<a  href="src/kw-test.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/kw-test.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>入力欄には，群変数と分析対象変数を１行に一対ずつの値を<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