メイン pca.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">
<!--
function zfi(w, i)
{
var ww = (""+i).length
ww = ww > w ? 0 : w-ww
return "0000000000000".substring(0, ww)+i
}
function fname(i)
{
return "PCA"+zfi(2, i+1)
}
function vname(i)
{
return "Var"+zfi(2, i+1)
}
function n_factor(m, ne0)
{
var ne
if (isNaN(ne = parseInt(ne0)) || ne < 1 || ne > m) {
ne = m
}
return ne
}
function calc(data_string, ne0)
{
var data, nc, nv, i, j, k, mean, variance, sd, sdp, r, s, v, e, isw, tm
if ((data = getdata(data_string, 0)) == false) return
nc = data.length
nv = data[0].length
if (nc < 2) {
printf("ケース数が１以下です\n")
return
}
mean = new Array(nv)
variance = new Array(nv)
sd = new Array(nv)
sdp = new Array(nv)
r = makeMatrix(nv, nv)
// Netscape communicator のバグ？配列ベクトルは初期化されていない
for (i = 0; i < nv; i++) {
mean[i] = 0
for (j = 0; j < nv; j++) {
r[i][j] = 0
}
}
for (i = 0; i < nc; i++) {
for (j = 0; j < nv; j++) {
mean[j] += data[i][j]
}
}
for (j = 0; j < nv; j++) {
mean[j] /= nc
}
for (i = 0; i < nc; i++) {
for (j = 0; j < nv; j++) {
for (k = 0; k < nv; k++) {
r[j][k] += (mean[j]-data[i][j])*(mean[k]-data[i][k])
}
}
}
isw = 0
for (j = 0; j < nv; j++) {
variance[j] = r[j][j]/(nc-1)
sd[j] = Math.sqrt(variance[j])
if (variance[j] <= 0) isw = 1
}
if (isw == 1) {
printf("分散が０になる変数があります\n")
return
}
for (j = 0; j < nv; j++) {
for (k = 0; k < nv; k++) {
r[j][k] /= (nc-1)*sd[j]*sd[k]
}
}
printf("★ 主成分分析 ★\n\n")
printf("有効データ組数 = %i\n\n", nc)
printf("%5s%15s%15s%15s\n", "変数", "平均値", "不偏分散", "標準偏差")
for (j = 0; j < nv; j++) {
printf("%5s%15.7g%15.7g%15.7g\n", vname(j), mean[j], variance[j], sd[j])
sdp[j] = Math.sqrt(variance[j]*(nc-1)/nc) // sdp は母標準偏差！
}
printf("\n相関係数行列\n\n ")
for (j = 0; j < nv; j++) {
printf("%8s", vname(j))
}
printf("\n")
for (j = 0; j < nv; j++) {
printf("%5s", vname(j))
for (k = 0; k < nv; k++) {
printf("%8.3f", r[j][k])
}
printf("\n")
}
ne = n_factor(nv, ne0)
v = makeMatrix(ne, nv)
e = new Array(ne)
if (eigen(r, nv, ne, 1e-10, e, v) == 1) {
printf("固有値が求められませんでした\n")
return
}
printf("\n主成分負荷量\n\n ")
for (j = 0; j < ne; j++) {
printf("%10s", fname(j))
}
printf("%10s\n", "寄与率")
for (j = 0; j < nv; j++) {
tm = 0
printf("%6s", vname(j))
for (k = 0; k < ne; k++) {
tm += v[k][j]*v[k][j]*e[k]
printf("%10.5f", v[k][j]*Math.sqrt(e[k]))
}
printf("%10.5f\n", tm)
}
printf("%6s", "固有値")
for (j = 0; j < ne; j++) {
printf("%10.5f", e[j])
}
printf("\n%6s", "寄与率")
s = 0
for (j = 0; j < ne; j++) {
printf("%10.2f", 100*e[j]/nv)
}
printf("\n%6s", "累積％")
s = 0
for (j = 0; j < ne; j++) {
s += e[j]
printf("%10.2f", 100*s/nv)
}
printf("\n\n主成分得点係数\n\n ")
for (j = 0; j < ne; j++) {
printf("%10s", fname(j))
}
printf("\n")
for (j = 0; j < nv; j++) { // v[k][j] を主成分得点係数に変換
printf("%6s", vname(j))
for (k = 0; k < ne; k++) {
// v[k][j] /= Math.sqrt(e[k]) // 主成分得点の分散を1にするときに必要
printf("%10.5f", v[k][j])
}
printf("\n")
}
printf("\n主成分得点\n\nケース")
for (j = 0; j < ne; j++) {
printf("%10s", fname(j))
}
printf("\n")
for (i = 0; i < nc; i++) {
printf("%6i", i+1)
for (k = 0; k < ne; k++) {
s = 0
for (j = 0; j < nv; j++) {
s += v[k][j]*(data[i][j]-mean[j])/sdp[j]
}
printf("%10.5f", s)
}
printf("\n")
}
}
//-->
</script>
</head>
<body bgcolor="#ffffff">
<font size="+2"><b>主成分分析</b></font> <a href="src/pca.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, this.form.ne.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/corr.html">使用法</a></td>
</tr></table>
<p>
抽出する主成分の個数<input name="ne" value="" size=5>（省略時は変数の個数と同じ）<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="../lecture/PCA/index.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
