メイン kw.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">
<!--
//----------------------------------
var nr, nc, um, t, gamma2, half, rank_sum
function zfi(w, i)
{
var ww = (""+i).length
ww = ww > w ? 0 : w-ww
return "000000000000000".substring(0, ww)+i
}
function vname(i, str)
{
return "第"+zfi(2, i+1)+str
}
//----------------------------------
function marginals()
{
var i, j
for (i = 0; i < nr; i++) {
for (j = 0; j < nc; j++) {
t[i][nc] += t[i][j]
}
for (j = 0; j <= nc; j++) {
t[nr][j] += t[i][j]
}
}
}
function kw_h(u)
{
var i, j
var r, retv = 0.0
for (i = 0; i < nr; i++) {
r = -t[i][nc]*half
for (j = 0; j < nc; j++) {
r += rank_sum[j][u[i][j]]
}
retv += r*r/t[i][nc]
}
return retv
}
function asymptotic()
{
var i, j, df, p, stat_val
printf("★ Kruskal-Wallis の検定 ★\n ")
for (j = 0; j < nc; j++) {
printf(" %6s", vname(j, "列"))
}
printf(" %6s\n", "合計")
for (i = 0; i < nr; i++) {
printf("%6s", vname(i, "群"))
for (j = 0; j < nc; j++) {
printf(" %6i", t[i][j])
}
printf(" %6i\n", t[i][j])
}
printf("%6s", "合計 ")
for (j = 0; j < nc; j++) {
printf(" %6i", t[i][j])
}
printf(" %6i\n\n", t[i][j])
stat_val = kw_h(t)
p = xxp(stat_val*gamma2, nr-1)
printf(" カイ二乗値 = %.5f\n", stat_val*gamma2)
printf(" 自由度 = %i\n", nr-1)
printf(" P 値 = %.5f\n\n", p)
}
//----------------------------------
function pow3(x)
{
return x*x*x
}
function calc(data_string)
{
var i, j, k, t0, ntot, tie, q
if ((t0 = getdata(data_string, 0)) != false) {
if ((nr = t0.length) == 0) {
printf("分割表が入力されていません\n")
return
}
nc = t0[0].length
if (nr < 2 || nc < 2) {
printf("２行２列以上の大きさの分割表を対象にします\n")
return
}
t = new Array(nr+1)
for (i = 0; i <= nr; i++) {
t[i] = new Array(nc+1)
for (j = 0; j <= nc; j++) {
t[i][j] = 0
}
}
for (i = 0; i < nr; i++) {
for (j = 0; j < nc; j++) {
t[i][j] = t0[i][j]
}
}
marginals()
q = new Array(nc)
rank_sum = new Array(nc)
gamma2 = tie = 0
ntot = t[nr][nc]
for (j = 0; j < nc; j++) {
q[j] = tie+(t[nr][j]+1)*0.5
gamma2 +=pow3(t[nr][j])-t[nr][j]
tie += t[nr][j]
rank_sum[j] = new Array(t[nr][j]+1)
for (k = 0; k <= t[nr][j]; k++) {
rank_sum[j][k] = k*q[j]
}
}
gamma2 = 12/(ntot*(ntot+1)*(1-gamma2/(pow3(ntot)-ntot)))
half = (ntot+1)*0.5
asymptotic()
}
sep(40)
}
//-->
</script>
</head>
<body bgcolor="#ffffff">
<font size="+2"><b>分割表形式でのクラスカル・ウォリスの検定</b></font> <a href="src/kw.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/kw.html">使用法</a><p>
<table>
<tr>
<td>入力欄<br><textarea name="data" ROWS=20 COLS=15></textarea></td>
<td>出力欄<br><textarea name="result" ROWS=20 COLS=80></textarea></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
