メイン cross3.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 save_index(x, nc, ix)
{
var nx, i, j
nx = 0
for (i = 0; i < nc; i++) {
ix[nx] = x[i]
for (j = 0; j <= nx; j++) {
if (ix[j] == x[i]) {
break
}
}
if (j == nx) {
nx++
}
}
return nx
}
function sort_index(ix, nx)
{
var i, j, mini, minx, temp
for (i = 0; i < nx-1; i++) {
mini = i
minx = ix[i]
for (j = i+1; j < nx; j++) {
if (ix[j] < minx) {
mini = j
minx = ix[j]
}
}
if (mini != i) {
temp = ix[i]
ix[i] = minx
ix[mini] = temp
}
}
}
function get_index(x, ix, nx)
{
var i
for (i = 0; i < nx; i++) {
if (x == ix[i]) {
return i
}
}
}
function clear_table(table, nx, ny)
{
var i, j
for (i = 0; i <= nx; i++) {
for (j =0; j <= ny; j++) {
table[i][j] = 0
}
}
}
function count(x, y, nc, table, ix, nx, iy, ny)
{
var i, j, k
for (k = 0; k < nc; k++) {
i = get_index(x[k], ix, nx)
j = get_index(y[k], iy, ny)
table[i][j]++
}
}
function sum(table, nx, ny)
{
var i, j
for (i = 0; i < nx; i++) {
for (j = 0; j < ny; j++) {
table[i][ny] += table[i][j]
}
for (j = 0; j <= ny; j++) {
table[nx][j] += table[i][j]
}
}
}
function sep(n)
{
var i
for (i = 0; i < n; i++) {
printf("%10s", "----------")
}
printf("\n")
}
function print_table(table, ix, iy, nx, ny, v1, v2)
{
var i, j
printf("\n==================\n\n")
sep(ny+2)
printf("%"+(ny*5+15)+"s\n", "第 "+v2+" 列")
printf("%13s", "")
sep(ny)
printf("%10s", "第 "+v1+" 列")
for (j = 0; j < ny; j++) {
printf("%10s", iy[j])
}
printf("%10s\n", "合計")
sep(ny+2)
for (i = 0; i < nx; i++) {
printf("%10s", ix[i])
for (j = 0; j <= ny; j++) {
printf("%10i", table[i][j])
}
printf("\n")
}
sep(ny+2)
printf("%10s", "合計")
for (j = 0; j <= ny; j++) {
printf("%10s", table[nx][j])
}
printf("\n")
sep(ny+2)
}
function chi_sq_test(table, nx, ny)
{
var i, j, chi, ex, df, one, five
chi = one = five = 0
for (i = 0; i < nx; i++) {
for (j = 0; j < ny; j++) {
ex = table[i][ny]*table[nx][j]/table[nx][ny]
chi += Math.pow(table[i][j]-ex, 2)/ex
if (ex < 5) {
five++
if (ex < 1) {
one++
}
}
}
}
df = (nx-1)*(ny-1)
printf("カイ二乗値 = %.3f, 自由度 = %i, P 値 = %.3f\n", chi, df, xxp(chi, df))
if (one > 0) {
printf("期待値が1未満のセルがあります。検定結果は慎重に解釈すべきです。\n")
}
if (five >= nx*ny*0.2) {
printf("期待値が5未満のセルが，全体のセルの20%%以上あります。検定結果は慎重に解釈すべきです。\n")
}
}
function getdata_cross(data_string, checkflag)
{
var lines, vars, nv, nc, j0, nv0, i, j
var delimiter = ","
var data_matrix = new Array()
var newline2
if (data_string.indexOf(delimiter) == -1) {
if (data_string.indexOf("\t") != -1) {
delimiter = "\t"
}
else {
delimiter = " "
}
}
newline2 = whichCRLF()
lines = data_string.split(newline2)
nc = 0
for (i = 0; i < lines.length; i++) {
if (lines[i] != "") {
vars = lines[i].split(delimiter)
nv = vars.length
data_matrix[nc] = new Array()
j0 = 0
for (j = 0; j < nv; j++) {
if (vars[j] != "") {
data_matrix[nc][j0++] = vars[j]
}
}
nc++
if (checkflag >= 0) {
if (nc == 1) {
nv0 = j0
}
else if (checkflag == 0 && j0 != nv0) {
printf("入力欄の %i 行目に記述されたデータの個数が前の行と違います。\n入力行：%s\n", i+1, lines[i])
return false
}
else if (checkflag > 0 && j0 != checkflag) {
printf("入力欄の %i 行目に記述されたデータの個数が %i ではありません。\n入力行：%s\n", i+1, checkflag, lines[i])
return false
}
}
}
}
if (nc == 0) {
printf("データが入力されていません。\n")
return false
}
return data_matrix
}
function calcul(data, nc, v1, v2)
{
var i, nx, ny, ix, iy, x, y, table
ix = new Array(100)
iy = new Array(100)
x = new Array(nc)
y = new Array(nc)
for (i = 0; i < nc; i++) {
x[i] = data[i][v1-1]
y[i] = data[i][v2-1]
}
nx = save_index(x, nc, ix)
ny = save_index(y, nc, iy)
sort_index(ix, nx)
sort_index(iy, ny)
table = makeMatrix(nx+1, ny+1)
clear_table(table, nx, ny)
count(x, y, nc, table, ix, nx, iy, ny)
sum(table, nx, ny)
print_table(table, ix, iy, nx, ny, v1, v2)
chi_sq_test(table, nx, ny)
}
function calc(data_string, v10, v20)
{
var data, nc, nv
var v1, v2
if ((data = getdata_cross(data_string, 0)) == false) return
nc = data.length
nv = data[0].length
if (nc < 2) {
printf("ケース数が１以下です\n")
return
}
if (v10 == "" || v20 == "") {
for (v1 = 1; v1 < nv; v1++) {
for (v2 = v1+1; v2 <= nv; v2++) {
calcul(data, nc, v1, v2)
}
}
}
else {
if (isNaN(v1 = parseInt(v10)) || v1 < 0 || v1 > nv || v1 != Math.floor(v1) || isNaN(v2 = parseInt(v20)) || v2 < 0 || v2 > nv || v2 != Math.floor(v2)) {
printf("集計する列の指定がへんです\n")
return
}
calcul(data, nc, v1, v2)
}
}
//-->
</script>
</head>
<body bgcolor="#ffffff">
<font size="+2"><b>クロス集計・独立性の検定</b></font> <a href="src/cross3.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>
カテゴリー化してあるデータ行列を入力欄に入力し，クロス集計する2変数の列を指定して，クロス集計表を行い，独立性の検定を行う。<br>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Cross/cross.html">手法の解説</a> <img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/cross3.html">使用法</a>
<hr>
<table><tr>
<td><input type="button" name="calcurate" value="計算開始" onClick="calc(this.form.data.value, this.form.v1.value, this.form.v2.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 nowrap></td>
</tr></table>
<p>
クロス集計する2変数の指定：<input name="v1" value="" size=3>列目 と <input name="v2" value="" size=3>列目のデータ<br>
入力欄<br><textarea name="data" ROWS=20 COLS=80></textarea><p>
出力欄<br><textarea name="result" ROWS=30 COLS=80></textarea>
</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
