メイン 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("類似度行列の対角成分に0でないものがありましたので,0に変更して以下の解析を行いました。\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

Made with Macintosh