メイン cross2.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">
<!--
//----------------------------------
var nr, nc, um, perm_table, t, expectation, criterion, nntrue, nntrue2, log_expmax, ntables, stat_val, nc1, nr1, Fisher = true
var EXPMAX = 1e100
var EPSILON = 1e-14
function zfi(w, i)
{
var ww = (""+i).length
ww = ww > w ? 0 : w-ww
return "000000000000000".substring(0, ww)+i
}
function vname(i, str)
{
return "第"+zfi(2, i+1)+str
}
//----------------------------------
function permutation(n)
{
var i
perm_table = new Array(n+1)
perm_table[0] = perm_table[1] = 0.0
for (i = 2; i <= n; i++) {
perm_table[i] = perm_table[i-1]+Math.log(i)
}
}
function calc_prob(um)
{
var nprod = 0
for (i = 0; i < nr; i++) {
nprod += perm_table[t[i][nc]]
for (j = 0; j < nc; j++) {
nprod -= perm_table[um[i][j]]
}
}
return nprod
}
function found()
{
var i, j
var hh, nprod
if (Fisher == true) {
nprod = calc_prob(um)
if (nprod <= criterion+EPSILON) {
nntrue += Math.exp(nprod-nntrue2*log_expmax)
while (nntrue >= EXPMAX) {
nntrue /= EXPMAX
nntrue2++
}
}
}
else {
hh = chi_sq(um)
if (hh >= stat_val-EPSILON) {
nprod = calc_prob(um)
nntrue += Math.exp(nprod-nntrue2*log_expmax)
while (nntrue >= EXPMAX) {
nntrue /= EXPMAX
nntrue2++
}
}
}
ntables++
}
function search(x, y)
{
var t
if (y == 0) {
found()
}
else if (x == 0) {
if ((t = um[0][0]-um[y][0]) < 0) {
return
}
um[0][0] = t
search(nc1, y-1)
um[0][0] += um[y][0]
}
else {
search(x-1, y)
while (um[y][0] && um[0][x]) {
++um[y][x]
--um[y][0]
--um[0][x]
search(x-1, y)
}
um[y][0] += um[y][x]
um[0][x] +=um[y][x]
um[y][x] = 0
}
}
function exact()
{
var i, j, ntot, denom, denom2
ntables = nntrue = nntrue2 = 0
ntot = t[nr][nc]
permutation(ntot)
denom2 = 0
log_expmax = Math.log(EXPMAX)
denom = perm_table[ntot]
for (j = 0; j < nc; j++) {
denom -= perm_table[t[nr][j]]
}
while (denom > log_expmax) {
denom -= log_expmax
denom2++
}
denom = Math.exp(denom)
criterion = calc_prob(t)
printf("criterion=%g\n", criterion);
um = new Array(nr)
for (i = 0; i < nr; i++) {
um[i] = new Array(nc)
for (j = 0; j < nc; j++) {
um[i][j] = 0
}
}
// 列和
for (i = 0; i < nc; i++) {
um[0][i] = t[nr][i]
}
// 行和 ２．．．ｋ
for (i = 1; i < nr; i++) {
um[i][0] =t[i][nc]
}
nc1 = nc-1
nr1 = nr-1
search(nc1, nr1)
printf("%s の方法による正確な P 値 = %.10g\n", Fisher ? "Fisher" : "Pearson", nntrue / denom * Math.pow(EXPMAX, nntrue2 - denom2))
printf("査察した分割表の個数は %s 個\n", ntables)
}
//----------------------------------
function marginals()
{
var i, j
for (i = 0; i < nr; i++) {
for (j = 0; j < nc; j++) {
t[i][nc] += t[i][j]
}
for (j = 0; j <= nc; j++) {
t[nr][j] += t[i][j]
}
}
}
function make_exp()
{
var i, j
expectation = new Array(nr)
for (i = 0; i < nr; i++) {
expectation[i] = new Array(nc)
for (j = 0; j < nc; j++) {
expectation[i][j] = t[i][nc]*t[nr][j]/t[nr][nc]
}
}
}
function chi_sq(u)
{
var i, j
var temp, retv = 0.0
for (i = 0; i < nr; i++) {
for (j = 0; j < nc; j++) {
temp = expectation[i][j] - u[i][j]
retv += temp*temp/expectation[i][j]
}
}
return retv
}
function asymptotic()
{
var i, j, df, p
printf("★ 独立性の検定（Exact test） ★\n ")
for (j = 0; j < nc; j++) {
printf(" %6s", vname(j, "列"))
}
printf(" %6s\n", "合計")
for (i = 0; i < nr; i++) {
printf("%6s", vname(i, "行"))
for (j = 0; j < nc; j++) {
printf(" %6i", t[i][j])
}
printf(" %6i\n", t[i][j])
}
printf("%6s", "合計 ")
for (j = 0; j < nc; j++) {
printf(" %6i", t[i][j])
}
printf(" %6i\n\n", t[i][j])
make_exp(t)
stat_val = chi_sq(t)
df = (nr-1)*(nc-1)
p = xxp(stat_val, df)
printf(" カイ二乗値 = %.5f\n", stat_val)
printf(" 自由度 = %i\n", df)
printf(" P 値 = %.5f\n\n", p)
}
//----------------------------------
function calc(data_string)
{
var i, j, t0
if ((t0 = getdata(data_string, 0)) != false) {
if ((nr = t0.length) == 0) {
printf("分割表が入力されていません\n")
return
}
nc = t0[0].length
if (nr < 2 || nc < 2) {
printf("２行２列以上の大きさの分割表を対象にします\n")
return
}
t = new Array(nr+1)
for (i = 0; i <= nr; i++) {
t[i] = new Array(nc+1)
for (j = 0; j <= nc; j++) {
t[i][j] = 0
}
}
for (i = 0; i < nr; i++) {
for (j = 0; j < nc; j++) {
t[i][j] = t0[i][j]
}
}
marginals()
asymptotic()
exact()
}
sep(40)
}
//-->
</script>
</head>
<body bgcolor="#ffffff">
<font size="+2"><b>独立性の検定（Exact test）</b></font> <a href="src/cross2.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>
集計済みのクロス集計表を入力し，独立性の検定を行う。漸近近似によるカイ二乗検定を行った上で，正確な検定（Exact test）を行う。<br>
Fisher の方法とは，観察された分割表の生起確率よりも小さい生起確率を持つ分割表の生起確率を合計したものを P 値とする。<br>
Pearson のカイ二乗法とは，観察された分割表のカイ二乗値以上の値を持つ分割表の生起確率を合計したものを P 値とする。<br>
参考文献：<a href="http://www.cytel.com/new.pages/SX.tp.2.html">Mehta の Exact Inference for Categorical Data</a><br>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../exact/exact.html">関連する解説ページ</a> <img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/cross2.html">使用法</a><br>
注意：計算時間が長くなる場合がある（データの個数などからは事前に予測できない）。
<hr>
<input type="radio" name="method" onClick="Fisher=true" checked> Fisher の方法
<input type="radio" name="method" onClick="Fisher=false"> Pearson の方法 <p>
<input type="button" name="calcurate" value="計算開始" onClick="calc(this.form.data.value)">
<input type="button" name="clear" value="入力欄クリア" onClick="this.form.data.value=''">
<input type="button" name="clear" value="出力欄クリア" onClick="this.form.result.value=''">
<p>
<table><tr>
<td>入力欄<br><textarea name="data" ROWS=20 COLS=20></textarea></td>
<td>出力欄<br><textarea name="result" ROWS=20 COLS=80></textarea></td>
</tr></table>
</form>
<p><hr noshade>
<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><IMG SRC="../gra/ume5.png" width=121 height=37 ALT="Made with Macintosh">
</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
