メイン　ttest2.html　　　Last modified: Jul 08, 2010
<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="io.js">document.write("io.js ファイルが見つかりません??<br>")</script>
<script src="fxp.js">document.write("fxp.js ファイルが見つかりません??<br>")</script>
<script src="betai.js">document.write("betai.js ファイルが見つかりません??<br>")</script>

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

function txp(t, df)
{
  return betai(df/(df+t*t), df/2, 0.5)
}

function decomp(data, a)
{
  var i, j, n = 0
  for (i = 0; i < data.length; i++) {
    for (j = 0; j < data[i].length; j++) {
      a[n++] = data[i][j]
    }
  }
  return n
}

function getMean(a, n)
{
  var mean = 0
  var i
  for (i = 0; i < n; i++) {
    mean += a[i]
  }
  return mean/n
}

function getVar(a, n, m)
{
  var v = 0
  var i
  if (n <= 1) {
    return 0
  }
  else {
    for (i = 0; i < n; i++) {
      v += Math.pow(a[i]-m, 2)
    }
    return v/(n-1)
  }
}

function v_test2(n1, n2, u1, u2)
{
  var f, df1, df2, p
  df1 = n1-1
  df2 = n2-1
  printf("二群の等分散性の検定\n")
  if (u2 == 0) {
    f = Number.POSITIVE_INFINITY
    p = 0
    printf("　  F 値 = ∞\n")
    printf("　自由度 = ( %i, %i )\n", df1, df2)
    printf("　  P 値 = 0　（両側確率）\n\n")
  }
  else {
    f = u1/u2
    p = fxp(f, df1, df2)
    p = Math.min(p, 1-p)*2
    printf("　  F 値 = %.6g\n", f)
    printf("　自由度 = ( %i, %i )\n", df1, df2)
    printf("　  P 値 = %.6f　（両側確率）\n\n", p)
  }
}

function confirm(n1, m1, u1, n2, m2, u2)
{
  printf("第一群：\n")
  printf("　標本サイズ = %i\n", n1)
  printf("　    平均値 = %.6g\n", m1)
  printf("　  不偏分散 = %.6g\n\n", u1)
  printf("第二群：\n")
  printf("　標本サイズ = %i\n", n2)
  printf("　    平均値 = %.6g\n", m2)
  printf("　  不偏分散 = %.6g\n\n", u2)
}

function effect_size(n1, m1, u1, n2, m2, u2)
{
  var g
  g = (m2-m1)/Math.sqrt(((n1-1)*u1+(n2-1)*u2)/(n1+n2-2))
  printf("effect size g = %.6g\neffect size d = %.6g\n\n", g, (1-3/(4*(n1+n2)-9))*g)
}

function t_test(n1, m1, u1, n2, m2, u2)
{
  var df, ue, t, p, p1, p2
  df = n1+n2-2
  ue = ((n1-1)*u1+(n2-1)*u2)/df
  t = Math.abs(m1-m2)/Math.sqrt(ue*(1/n1+1/n2))
  p = txp(t, df)
  printf("通常の t 検定（等分散性が仮定できるとき）\n")
  printf("　  t 値 = %.5f\n", t)
  printf("　自由度 = %i\n", df)
  printf("　  P 値 = %.5f\n\n", p)
  t = Math.abs(m1-m2)/Math.sqrt(u1/n1+u2/n2)
  df = Math.pow(u1/n1+u2/n2,2)/(Math.pow(u1/n1,2)/(n1-1)+Math.pow(u2/n2,2)/(n2-1))
  p = txp(t, df)
  printf("等分散性が仮定できないとき（Welch の方法）\n")
  printf("　  t 値 = %.5f\n", t)
  printf("　自由度 = %.5f\n", df)
  printf("　  P 値 = %.5f　（小数自由度に対応した正確な値）\n\n", p)
}

function calc(data1_string, data2_string)
{
  var data1, data2, n1, m1, u1, n2, m2, u2
  var array1 = new Array()
  var array2 = new Array()
  if ((data1 = getdata(data1_string, -1)) != false && (data2 = getdata(data2_string, -1)) != false) {
    n1 = decomp(data1, array1)
    n2 = decomp(data2, array2)
    m1 = getMean(array1, n1)
    m2 = getMean(array2, n2)
    u1 = getVar(array1, n1, m1)
    u2 = getVar(array2, n2, m2)
    confirm(n1, m1, u1, n2, m2, u2)
    v_test2(n1, n2, u1, u2)
    t_test(n1, m1, u1, n2, m2, u2)
    effect_size(n1, m1, u1, n2, m2, u2)
  }
  sep(40)
  printf("\n")
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>原データからの t 検定（二群の平均値の差の検定）</b></font>　<a  href="src/ttest2.html"><img src="png/src.png" width=35 height=11 alt="src" align=top></a>　　　　　　　　　　　　　　　　　　　　　　　　　Last modified: Jul 08, 2010<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.data1.value,this.form.data2.value)">　　</td>
<td><input type="button" name="clear" value="入力欄クリア" onClick="this.form.data1.value=this.form.data2.value=''">　　</td>
<td><input type="button" name="clear" value="出力欄クリア" onClick="this.form.result.value=''">　　</td>
<td nowrap><img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/ttest2.html">使用法</a></td>
</tr></table>
<p>

<table>
<tr>
<td valign=top>第一群のデータ<br><textarea name="data1" rows=10 cols=20></textarea><p>
第二群のデータ<br><textarea name="data2" rows=10 cols=20></textarea></td>
<td valign=top rowspan=2>出力欄<br><textarea name="result" rows=24 cols=80></textarea></td>
</tr>
</table>
</form>
<p><hr noshade>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Average/t-test.html">手法の解説(1)</a>　　<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Average/bunsan1.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>

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

Made with Macintosh