メイン　corr5.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="inv.js">document.write("inv.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 "000000000000000".substring(0, ww)+i
}

function vname(i)
{
  return "Var"+zfi(2, i+1)
}

function matout(r, nv)
{
  var j, k
  printf("     ")
  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")
  }
}

function calc(data_string)
{
  var epsinv = 1e-8
  var i, j, k, nc, nv, data, mean, variance, sd, r, r2, isw
  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)
  r = makeMatrix(nv, nv)
  r2 = 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
  }
  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])
  }

  if (isw == 1) {
    printf("\nエラー：分散が０になる変数があります。相関係数行列は完全には計算されません。\n")
  }

  for (j = 0; j < nv; j++) {
    for (k = 0; k < nv; k++) {
      r2[j][k] = r[j][k] = (sd[j]*sd[k] == 0) ? "div0" : r[j][k]/((nc-1)*sd[j]*sd[k])
    }
  }

  if (isw == 1) {
    printf("\n★ 相関係数行列 ★\n\n")
    matout(r, nv)
    printf("\n上の出力結果で[NaN.000]が表示されている変数の組では，相関係数が計算できないことを表しています。\n")
    return
  }
  if (inv(r2, nv, epsinv) == 0) {
    printf("\n★ 相関係数・偏相関係数・重相関係数 ★\n\n")
    for (i = 0; i < nv; i++) {
      r[i][i] = Math.sqrt(1-1/r2[i][i])
      for (j = i+1; j < nv; j++) {
        r[j][i] = -r2[i][j]/Math.sqrt(r2[i][i]*r2[j][j])
      }
    }
    matout(r, nv)
    printf("\n下三角行列 r(x,y) は2変数 x, y 以外の全ての変数を制御したときの x, y の間の偏相関係数\n上三角行列は（単純）相関係数\n対角成分 r(x,x) は，変数 x と他の全ての変数との重相関係数\n")
  }
  else {
    printf("\nエラー：相関係数行列の逆行列が求まらないので，偏相関係数と重相関係数は計算できません。\n")
    printf("\n★ 相関係数行列 ★\n\n")
    matout(r, nv)
  }
  printf("\n")
  sep(40)
  printf("\n")
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>相関係数行列，偏相関係数行列，重相関係数</b></font>　<a  href="src/corr5.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>
複数の変数データを入力し，変数間の相関係数，偏相関係数，重相関係数を求める。<br>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Soukan/pearson.html">手法の解説(1)</a>　<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Soukan/partial-corr.html">手法の解説(2)</a>　<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Soukan/multiple-corr.html">手法の解説(3)</a>　<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/corr.html">使用法</a>
<hr>
<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>
</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="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>

サブ　inv.js　　　Last modified: Mar 25, 2004
サブ　io.js　　　Last modified: Mar 25, 2004
・ 直前のページへ戻る　　・ E-mail to Shigenobu AOKI

Made with Macintosh