メイン t-test.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">
<!--
var which = "不偏分散"
function dispersion(sw)
{
which = sw
document.Result.which_disp.value = which
}
function txp(t, df)
{
return betai(df/(df+t*t), df/2, 0.5)
}
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 v_test2(n1, n2, u1, u2)
{
var f, df1, df2, p
df1 = n1-1
df2 = n2-1
printf("二群の等分散性の検定\n")
if (u2 == 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(" %s = %.6g\n\n", which, u1)
printf("第二群：\n")
printf(" 標本サイズ = %i\n", n2)
printf(" 平均値 = %.6g\n", m2)
printf(" %s = %.6g\n\n", which, 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(frm)
{
var n1, m1, u1, n2, m2, u2
if (frm.n1.value == "" || frm.m1.value == "" || frm.u1.value == "" || frm.n2.value == "" || frm.m2.value == "" || frm.u2.value == "") {
printf("入力が不完全ですね\n")
return
}
n1 = parseInt(frm.n1.value)
m1 = parseFloat(frm.m1.value)
u1 = parseFloat(frm.u1.value)
n2 = parseInt(frm.n2.value)
m2 = parseFloat(frm.m2.value)
u2 = parseFloat(frm.u2.value)
if (isNaN(n1) || isNaN(m1) || isNaN(u1) || isNaN(n2) || isNaN(m2) || isNaN(u2)) {
printf("入力が不完全か，数値でないものがあります\n")
return
}
else if (n1 != parseFloat(frm.n1.value) || n2 != parseFloat(frm.n2.value)) {
printf("標本サイズは整数値でなくてはいけません\n")
return
}
else if (n1 <= 0 || n2 <= 0) {
printf("標本サイズが０以下というのは！ ヘンですね\n")
return
}
if (u1 < 0 || u2 < 0) {
printf(which + "が負の数というのは！ ヘンですね\n")
return
}
if (u1 == 0 && u2 == 0) {
printf("両群の" + which + "が共に０というのは！ ヘンですね\n")
return
}
confirm(n1, m1, u1, n2, m2, u2)
u1 = convert(n1, u1)
u2 = convert(n2, u2)
v_test2(n1, n2, u1, u2)
t_test(n1, m1, u1, n2, m2, u2)
effect_size(n1, m1, u1, n2, m2, u2)
printf("--------------------\n")
}
//-->
</script>
</head>
<body bgcolor="#ffffff">
<font size="+2"><b>二群の平均値の差の検定</b></font> <a href="src/ttest.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 valign="top">
<table border=1 cellpadding=3>
<caption>入力欄</caption>
<tr>
<td></td> <th nowrap><input name="dummy1" value="標本サイズ" size=10></th> <th nowrap><input name="dummy2" value="平均値" size=6></th> <th nowrap><input name="which_disp" value="不偏分散" size=8></th>
</tr>
<tr>
<th nowrap>第一群</th>
<td><input name="n1" value="10" size=10></td>
<td><input name="m1" value="2.56" size=10></td>
<td><input name="u1" value="3.45" size=10></td>
</tr>
<tr>
<th nowrap>第二群</th>
<td><input name="n2" value="19" size=10></td>
<td><input name="m2" value="3.16" size=10></td>
<td><input name="u2" value="2.95" size=10></td>
</tr>
</table>
<table cellpadding=3>
<th colspan=2 nowrap align=left><br>散布度としては何を入力しますか</th>
<tr>
<td nowrap><input type="radio" name="which" onClick="dispersion('分散')">分散</td>
<td nowrap><input type="radio" name="which" checked onClick="dispersion('不偏分散')">不偏分散</td>
</tr>
<tr>
<td nowrap><input type="radio" name="which" onClick="dispersion('標準偏差')">標準偏差</td>
<td nowrap><input type="radio" name="which" onClick="dispersion('標準誤差')">標準誤差</td>
</tr>
<tr>
<td nowrap><br><input type="button" name="calcurate" value="計算開始" onClick="calc(this.form)"></td>
<td nowrap><br><input type="button" name="clear" value="入力欄クリア" onClick="this.form.n1.value=this.form.m1.value=this.form.u1.value=this.form.n2.value=this.form.m2.value=this.form.u2.value=''"></td>
<td nowrap><br><input type="button" name="clear" value="出力欄クリア" onClick="this.form.result.value=''"></td>
</tr>
</table>
</td>
<td valign="top" align="center">
出力欄<br><textarea name="result" ROWS=20 COLS=60></textarea>
</td>
</tr>
</table>
</form>
<ol>
<li> まず，「入力欄クリア」と「出力欄クリア」ボタンをクリックし，入力欄と出力欄をクリアします。
<li> 入力する散布度に対応するボタンをクリックします（標準偏差は不偏分散の平方根をとったものです）。
<li> 入力欄の６つの桝目に数値を入力します。
<li> 「計算開始」ボタンをクリックします。
<li> 入力欄をクリアして，別の数値を入力し計算します（出力欄はクリアしなくても，結果が追加されます。また，いつでもクリアして構いません）。
<li> 出力欄の内容は，コピー＆ペーストできます。
</ol>
<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
