メイン　armitage.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">
<!--

function calc(data_string)
{
  var data, nc, i, j
  var x = new Array()
  var n = new Array()
  var r = new Array()
  var p = new Array()
  var nt, xt, rt, p_bar, x_bar, b1, b2, b3, b, a, chi_t, chi_q, chi_h

  printf("★ コクラン・アーミテージの検定 ★\n\n")
  printf("%5s%12s%12s%8s%12s\n", "群", "外的基準値", "標本サイズ", "陽性数", "比率")

  if ((data = getdata(data_string, 3)) == false) return
  nc = data.length
  
  nt = rt = xt = 0
  for (i = 0; i < nc; i++) {
    x[i] = data[i][0]
    n[i] = data[i][1]
    r[i] = data[i][2]
    xt += n[i]*x[i]
    nt += n[i]
    rt += r[i]
    for (j = 0; j < i; j++) {
      if (x[i] == x[j]) {
        printf("外的基準値が同じものが複数個あります\n")
        return
      }
    }
    if (Math.floor(n[i]) != n[i] || Math.floor(r[i]) != r[i] || n[i] <= 0 || r[i] < 0 || r[i] > n[i]) {
      printf("標本サイズ，陽性数が変です   標本サイズ = %i   陽性数 = %i\n", n[i], r[i])
      return
    }
    p[i] = r[i]/n[i]
    printf("%5i%12.6g%12i%8i%12.5f\n", i+1, x[i], n[i], r[i], p[i])
  }
  p_bar = rt/nt
  x_bar = xt/nt
  b1 = b2 = b3 = 0
  for (i = 0; i < nc; i++) {
    b1 += n[i]*(p[i]-p_bar)*(x[i]-x_bar)
    b2 += n[i]*Math.pow(x[i]-x_bar, 2)
    b3 += r[i]*r[i]/n[i]
  }
  b = b1/b2
  a = p_bar-b*x_bar
  chi_t = b*b*b2/(p_bar*(1-p_bar))
  chi_h = nt*nt*(b3-rt*rt/nt)/(rt*(nt-rt))
  chi_q = chi_h-chi_t
  printf("\n切片 a = %g　　傾き b = %g\n\n", a, b)
  printf("%16s%12s%8s%12s\n", "要因", "カイ二乗値", "自由度", "有意確率")
  printf("%16s%12.5f%6i%14.5f\n", "傾き", chi_t, 1, xxp(chi_t, 1))
  printf("%16s%12.5f%6i%14.5f\n", "直線からの乖離", chi_q, nc-2, xxp(chi_q, nc-2))
  printf("%16s%12.5f%6i%14.5f\n\n", "合計（非一様性）", chi_h, nc-1, xxp(chi_h, nc-1))
  sep(50)
  printf("\n")
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>コクラン・アーミテージの検定</b></font>　<a  href="src/armitage.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/armitage.html">使用法</a></td>
</tr></table>

<p>
<table><tr>
<td>入力欄<br><textarea name="data" ROWS=25 COLS=25></textarea></td>
<td>出力欄<br><textarea name="result" ROWS=25 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/Hiritu/Armitage.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

Made with Macintosh