メイン dual.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="eigen.js">document.write("eigen.js ファイルが見つかりません??<br>")</script>
<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 lpw = 80
function zfi(w, i)
{
var ww = (""+i).length
ww = ww > w ? 0 : w-ww
return "000000000000000".substring(0, ww)+i
}
function axis_no(i)
{
return "第"+zfi(2, i+1)+"軸"
}
function name(str, i)
{
return str+zfi(2, i+1)
}
function sep3(ndash)
{
printf("-------")
while (ndash--) {
printf("------")
}
printf("\n")
}
function print_table(table, nxr, nxc)
{
var i, j, nxc1
nxc1 = nxc+1
printf("★ 分割表\n")
sep3(nxc1)
printf(" ")
for (i = 0; i < nxc; i++) {
printf("%6s", name("C", i))
}
printf(" 合計\n")
sep3(nxc1)
for (i = 0; i < nxr; i++) {
printf(" %4s ", name("R", i))
for (j = 0; j < nxc1; j++) {
printf("%6i", table[i][j])
}
printf("\n")
sep3(nxc1)
}
printf(" 合計")
for (j = 0; j < nxc1; j++) {
printf("%6i", table[nxr][j])
}
printf("\n")
sep3(nxc1)
}
function AxisName(is, ie, str)
{
var i
printf("%13s", str)
for (i = is; i < ie; i++) {
printf("%10s", axis_no(i))
}
printf("\n")
}
function out_sub(sw, is, ie, nr, nc, v, w, e)
{
var i, j, str
str = new Array("相関比による重み付けをしないもの", "相関比により重み付けをしたもの")
printf("\n★ 最適重みベクトル --- %s\n\n", str[sw])
AxisName(is, ie, "行")
for (j = 0; j < nr; j++) {
printf("%13s", name("R", j))
for (i = is; i < ie; i++) {
printf("%10.5f", sw ? w[j][i]*Math.sqrt(e[i]) : w[j][i])
}
printf("\n")
}
printf("\n")
AxisName(is, ie, "列")
for (j = 0; j < nc; j++) {
printf("%13s", name("C", j))
for (i = is; i < ie; i++) {
printf("%10.5f", sw ? v[i][j]*Math.sqrt(e[i]) : v[i][j])
}
printf("\n")
}
}
function dual(table, nr, nc)
{
var i, j, k, is, ie, naxis, w, e, v, fr, fc, ft, trace, cum1, cum2, work, max_nr_nc
max_nr_nc = Math.max(nr, nc)
w = makeMatrix(max_nr_nc, max_nr_nc)
v = makeMatrix(nc-1, nc)
e = makeVector(nc-1)
fr = makeVector(nr)
fc = makeVector(nc)
naxis = Math.min(nr, nc)-1
ft = table[nr][nc]
for (i = 0; i < nr; i++) {
fr[i] = table[i][nc]
}
for (j = 0; j < nc; j++) {
fc[j] = table[nr][j]
}
trace = 0.0
for (i = 0; i < nc; i++) {
for (j = 0; j < nc; j++) {
work = 0.0
for (k = 0; k < nr; k++) {
work += table[k][i]*table[k][j]/fr[k]
}
work /= Math.sqrt(fc[i]*fc[j])
work -= Math.sqrt(fc[i]*fc[j])/ft
w[i][j] = work
if (i == j) {
trace += work
}
}
}
eigen(w, nc, naxis, 1e-10, e, v)
for (i = 0; i < naxis; i++) {
for (j = 0; j < nc; j++) {
v[i][j] *= Math.sqrt(ft/fc[j])
}
}
for (i = 0; i < naxis; i++) {
for (j = 0; j < nr; j++) {
work = 0.0
for (k = 0; k < nc; k++) {
work += table[j][k]*v[i][k]
}
w[j][i] = work/(fr[j]*Math.sqrt(e[i]))
}
}
print_table(table, nr, nc)
cum2 = 0.0
for (is = ie = 0; ie != naxis; is =ie) {
ie = Math.min(ie+(lpw-13)/10, naxis)
cum1 = cum2
printf("\n")
AxisName(is, ie, " ")
printf(" 相関比の二乗")
for (i = is; i < ie; i++) {
printf("%10.5f", e[i])
}
printf("\n 相関係数")
for (i = is; i < ie; i++) {
printf("%10.5f", Math.sqrt(e[i]))
}
printf("\n 寄与率 ％")
for (i = is; i < ie; i++) {
printf("%10.1f", e[i]/trace*100)
}
cum2 = cum1
printf("\n累積寄与率 ％")
for (i = is; i < ie; i++) {
cum2 += e[i]/trace*100.0
printf("%10.1f", cum2)
}
printf("\n カイ二乗値")
for (i = is; i < ie; i++) {
if (e[i] < 1.0) {
printf("%10.3f", -(ft-1-(nr+nc-1)*0.5)*Math.log(1-e[i]))
}
}
printf("\n 自由度")
for (i = is; i < ie; i++) {
printf("%10i", nr+nc-1-2*(i+1))
}
printf("\n P 値")
for (i = is; i < ie; i++) {
if (e[i] < 1.0) {
printf("%10.5f", xxp(-(ft-1-(nr+nc-1)*0.5)*Math.log(1-e[i]), nr+nc-1-2*(i+1)))
}
}
printf("\n")
out_sub(0, is, ie, nr, nc, v, w, e)
out_sub(1, is, ie, nr, nc, v, w, e)
}
}
// このプログラムは、双対尺度法（西里）のうちの分割表、度数表の形式をとるデータを分析します。
// データ形式を変更すれば数量化III類（によって、同じ結果が得られます。
// 多肢選択データまたは反応型データについては、数量化III類を使用してください。
// その他のデータ型についてはサポートしていません。
// 参考文献
// 西里静彦：質的データの数量化−双対尺度法とその応用−、朝倉書店、１９８２
function calc(data_string)
{
var i, j, t0, t, nr, nc
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]
}
}
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]
}
}
dual(t, nr, nc)
}
sep(40)
}
//-->
</script>
</head>
<body bgcolor="#ffffff">
<font size="+2"><b>双対尺度法（データ形式が分割表・度数表の場合）</b></font> <a href="src/dual.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>
入力欄には分割表の本体部分（合計欄を除く）を，<a href="exa/kugirimoji.html">区切り文字</a>で区切って行単位に入力する。<br>
参考文献： 西里静彦「質的データの数量化−双対尺度法とその応用−」，朝倉書店，1982
<hr>
<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>
サブ eigen.js Last modified: Mar 25, 2004
サブ 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
