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

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

var which = "不定"
function dispersion(sw)
{
  which = sw
}

function convert(n, u)
{
  if    (which == "分散") return u*n/(n-1)
  else if (which == "不偏分散") return u
  else if (which == "標準偏差") return u*u
  else        return  n*u*u
}

function revconvert(n, u)
{
  if    (which == "分散") return u*(n-1)/n
  else if (which == "不偏分散") return u
  else if (which == "標準偏差") return Math.sqrt(u)
  else        return  Math.sqrt(u/n)
}

function calc(data_string)
{
  var i, j, sb, sw, st, vb, vw, dfb, dfw, f, p, chi, term1, term2, data, ng, ns, means, vars, gmean, gvar, nc
  var datax = new Array()
  printf("★ 一元配置分散分析 ★\n\n")
  if ((data = getdata(data_string, 3)) == false) return
  ng = data.length
  ns = new Array(ng)
  means = new Array(ng)
  vars = new Array(ng)
  gmean = nc = sb = sw = 0
  if (which == "不定") {
    printf("散布度として何を入力するかを選択してから\n「計算開始」ボタンを押してください\n")
    return
  }
  for (j = 0; j < ng; j++) {
    nc += (ns[j] = data[j][0])
    gmean += (means[j] = data[j][1])*ns[j]
    vars[j] = convert(ns[j], data[j][2])
    sw += vars[j]*(ns[j]-1)
  }
  gmean /= nc
  for (j = 0; j < ng; j++) {
    sb += ns[j]*Math.pow(means[j]-gmean, 2)
  }
  st = sb+sw
  gvar = st/(nc-1)
  dfb = ng-1
  dfw = nc-ng
  vb = sb/dfb
  vw = sw/dfw
  f = vb/vw
  p = fxp(f, dfb, dfw)
  printf("%10s %13s %12s %15s\n", "群", "データの個数", "平均値", which)
  for (j = 0; j < ng; j++) {
    printf("%10s %10i %15.6g %15.6g\n", "第"+(j+1)+"群", ns[j], means[j], data[j][2])
  }
  printf("%10s %10i %15.6g %15.6g\n\n", "全体", nc, gmean, revconvert(nc, gvar))
  printf("相関比 = %g\n", sb/st)
  sep(50)
  printf("分散分析表\n\n")
  printf("%4s %10s %7s %10s %10s %9s\n", "要因", "平方和", "自由度", "平均平方", "F 値", "P 値")
  printf("%4s %10.6g %7i %10.6g %10.5f %9.5f\n", "群間", sb, dfb, vb, f, p)
  printf("%4s %10.6g %7i %10.6g\n", "群内", sw, dfw, vw)
  printf("%4s %10.6g %7i %10.6g\n", "全体", st, nc-1, gvar)
  sep(50)
  printf("★ 等分散性の検定 ★\n\n")
  term1 = term2 = 0
  for (j = 0; j < ng; j++) {
    term2 += (ns[j]-1)*Math.log(vars[j])
    term1 += 1/(ns[j]-1)
  }
  chi = (dfw*Math.log(vw)-term2)/(1+(term1-1/dfw)/(3*(ng-1)))
  printf("Bartlett test\n  カイ二乗値 = %g   d.f. = %i   P 値 = %.5f\n\n", chi, ng-1, xxp(chi, ng-1))
  sep2(60)
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>二次データに基づく一元配置分散分析と等分散性の検定</b></font> <a  href="src/onewayANOVA2.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/onewayANOVA2.html">使用法</a></td>
</tr><tr>
<td colspan=4>散布度として入力するのは <input type="radio" name="kind" onClick="dispersion('分散')">分散  <input type="radio" name="kind" onClick="dispersion('不偏分散')">不偏分散  <input type="radio" name="kind" onClick="dispersion('標準偏差')">標準偏差  <input type="radio" name="kind" onClick="dispersion('標準誤差')">標準誤差</td>
</tr></table>

<p>
<table><tr>
<td>入力欄<br><textarea name="data" ROWS=20 COLS=15></textarea></td>
<td>出力欄<br><textarea name="result" ROWS=20 COLS=65></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/oneway-ANOVA.html">手法の解説(1)</a>  <img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Average/Bartlett.html">手法の解説(2)</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
サブ fxp.js   Last modified: Mar 25, 2004
サブ io.js   Last modified: Mar 25, 2004

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

Made with Macintosh