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

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

function getmean(data, nc)
{
  var i, sum = 0
  for (i = 0; i < nc; i++) {
    sum += data[i]
  }
  return sum/nc
}

function getvar(data, nc, mean)
{
  var i, sum = 0
  for (i = 0; i < nc; i++) {
    sum += Math.pow(data[i]-mean, 2)
  }
  return sum/(nc-1)
}

function calc(data_string)
{
  var i, nc, data
  var data1 = new Array()
  var data2 = new Array()
  var mean1, mean2, mean2, mean3, var1, var2, var3, t, p
  if ((data = getdata(data_string, 2)) != false) {
    nc = data.length
    if (nc < 2) {
      printf("有効データ組数が1以下です\n")
    }
    else {
      for (i = 0; i < nc; i++) {
        data1[i] = data[i][0]
        data2[i] = data[i][1]
      }
      mean1 = getmean(data1, nc)
      var1 = getvar(data1, nc, mean1)
      mean2 = getmean(data2, nc)
      var2 = getvar(data2, nc, mean2)
      for (i = 0; i < nc; i++) {
        data1[i] -= data2[i]
      }
      mean3 = getmean(data1, nc)
      var3 = getvar(data1, nc, mean3)
      if (mean1 == mean2) {
        printf("平均値の差は0です\n")
      }
      else if (var3 == 0.0) {
        printf("対応のあるデータの差の分散が0だったので,検定できません。\n")
      }
      else {
        t = Math.abs(mean3)/Math.sqrt(var3/nc)
        p = txp(t, nc-1)
        printf("★ 対応のある場合の平均値の差の検定 ★\n\n")
        printf("%4s%13s%13s%13s\n", "変数", "平均値", "不偏分散", "標準偏差")
        printf(" X1 %13.7g%13.7g%13.7g\n", mean1, var1, Math.sqrt(var1))
        printf(" X2 %13.7g%13.7g%13.7g\n", mean2, var2, Math.sqrt(var2))
        printf(" 差 %13.7g%13.7g%13.7g\n\n", mean3, var3, Math.sqrt(var3))
        printf("有効データ組数 = %i\n", nc)
        printf("t0 = %.5f\n", t)
        printf(" p = %.5f\n", p)
      }
    }
  }
  sep(50)
  printf("\n")
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>対応のある場合の平均値の差の検定</b></font> <a  href="src/pttest.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/pttest.html">使用法</a></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=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/paired-t-test.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>

サブ txp.js   Last modified: Mar 25, 2004
サブ io.js   Last modified: Mar 25, 2004

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

Made with Macintosh