メイン princo.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="io.js">document.write("io.js ファイルが見つかりません??<br>")</script>
<script language="JavaScript">
<!--
var lpw = 80
var EPSEIG = 1e-10
function zfi(w, i)
{
var ww = (""+i).length
ww = ww > w ? 0 : w-ww
return "0000000000000".substring(0, ww)+i
}
function vname(str, i)
{
return str+zfi(2, i+1)
}
function matout_sim(a, n)
{
var i, ie, is, j
for (is = ie = 0; ie != n; is = ie) {
ie = Math.min(is+Math.floor((lpw-12)/12), n)
printf("\n")
for (i = is; i < n; i++) {
printf("%-12s", vname("case", i))
for (j = is; j < ie; j++) {
printf("%12.4f", a[j][i])
}
printf("\n")
}
printf(" ")
for (i = is; i < ie; i++) {
printf("%12s", vname("case", i))
}
printf("\n")
}
}
function calc(data_string)
{
var i, ie, ok, ok2, is, j, n, naxis, eigen_value, cum, cum0, h, mean, norm, r, sum, temp, temp2, temp3
printf("★ 主座標分析 ★\n\n")
if ((r = getdata(data_string, 0)) == false) return
n = r.length
cum = makeVector(n)
mean = makeVector(n)
eigen_value = makeVector(n)
h = makeMatrix(n, n)
ok = ok2 = true
for (i = 0; i < n; i++) {
if (r[i][i] != 0) {
r[i][i] = 0
ok = false
}
for (j = 0; j < i; j++) {
if (r[i][j] != r[j][i]) {
ok2 = false
}
}
}
if (ok == false) {
printf("類似度行列の対角成分に０でないものがありましたので，０に変更して以下の解析を行いました。\n")
}
printf("★ 類似度行列\n")
matout_sim(r, n)
if (ok2 == false) {
printf("類似度行列が対称行列ではありませんでした。この分析手法では許されないことです。\n")
return
}
temp2 = 0.0
for (i = 0; i < n; i++) {
temp = 0.0
for (j = 0; j < n; j++) {
temp += r[j][i]
}
mean[i] = temp/n
temp2 += temp
}
temp2 /= n*n
for (i = 0; i < n; i++) {
for (j = 0; j <= i; j++) {
temp3 = r[j][i]
temp3 -= mean[i]+mean[j]
h[i][j] = h[j][i] = temp3+temp2
}
}
eigen(h, n, n, EPSEIG, eigen_value, r)
naxis = 0
sum = 0.0
for (i = 0; i < n; i++) {
if (eigen_value[i] > 1e-5) {
naxis++
sum += eigen_value[i]
}
}
if (naxis <= 0) {
printf("解は求まりません。")
return
}
for (i = 0; i < naxis; i++) {
norm = 0.0
for (j = 0; j < n; j++) {
norm +=r[i][j]*r[i][j]
}
temp = Math.sqrt(eigen_value[i]/norm)
for (j = 0; j < n; j++) {
r[i][j] *= temp
}
}
cum0 = 0.0
sum = 100.0/sum
for (i = 0; i < naxis; i++) {
cum0 += eigen_value[i]
cum[i] = cum0*sum
}
for (is = ie = 0; ie != naxis; is = ie) {
ie = Math.min(is+Math.floor((lpw-12)/12), naxis)
printf("\n★ 解\n\n ")
for (i = is; i < ie; i++) {
printf("%12s", vname("axis", i))
}
printf("\n 固有値")
for (i = is; i < ie; i++) {
printf("%12.3f", eigen_value[i])
}
printf("\n 寄与率％")
for (i = is; i < ie; i++) {
printf("%12.3f", eigen_value[i]*sum)
}
printf("\n累積寄与率％")
for (i = is; i < ie; i++) {
printf("%12.3f", cum[i])
}
printf("\n\n ")
for (i = is; i < ie; i++) {
printf("%12s", vname("axis", i))
}
printf("\n")
for (i = 0; i < n; i++) {
printf("%-12s", vname("case", i))
for (j = is; j < ie; j++) {
printf("%12.5f", r[j][i])
}
printf("\n")
}
}
}
//-->
</script>
</head>
<body bgcolor="#ffffff">
<font size="+2"><b>主座標分析</b></font> <a href="src/princo.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 nowrap><img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/princo.html">使用法</a></td>
</tr></table>
<p>
入力欄<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="../lecture/misc/princo.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>
サブ eigen.js Last modified: Mar 25, 2004
サブ io.js Last modified: Mar 25, 2004
直前のページへ戻る E-mail to Shigenobu AOKI
