メイン exp2.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="io.js">document.write("io.js ファイルが見つかりません??<br>")</script>

<script language="JavaScript">
<!--

function exp2(x, y, n, data, model)
{
  var i, mx, my, vx, vxy, a, b, b1, yval
  mx = my = vx = vxy = 0
  for (i = 0; i < n; i++) {
    mx += x[i]; my += y[i]
  }
  mx /= n
  my /= n

  for (i = 0; i < n; i++) {
    vx += Math.pow(x[i]-mx, 2)
    vxy += (x[i]-mx)*(y[i]-my)
  }
  vx /= n-1
  vxy /= n-1
  if (vx == 0) {
    printf("あてはめできません。\n")
  }
  else {
    a = vxy/vx
    b = my-a*mx
    b1 = Math.exp(b)
    printf("%s\n\n", model)
    printf("a = %g\n", b1)
    printf("b = %g\n\n", a)
    printf("★ 予測値および残差\n\nケース       独立変数       従属変数        予測値          残差\n")
    for (i = 0; i < n; i++) {
      yval = b1*Math.pow(data[i][0], a)
      printf("%6i%15.7g%15.7g%15.7g%15.7g\n",i+1, x[i], data[i][1], yval, data[i][1]-yval)
    }
  }
}

function calc(data_string, model)
{
  var i, nc, nv, data, x, y
  if ((data = getdata(data_string, 2)) == false) return
  nc = data.length
  nv = data[0].length
  if (nc < 2) {
    printf("ケース数が1以下です\n")
    return
  }

  x = new Array(nc)
  y = new Array(nc)
  for (i = 0; i < nc; i++) {
    if (data[i][0] <= 0 || data[i][1] <= 0) {
      printf("独立変数と従属変数はともに正の値でないといけません\n")
      return
    }
    x[i] = Math.log(data[i][0])
    y[i] = Math.log(data[i][1])
  }

  exp2(x, y, nc, data, model)

  printf("\n")
  sep(40)
  printf("\n")
}
//-->
</script>
</head>

<body bgcolor="#ffffff">
<font size="+2"><b>指数曲線回帰</b></font> <a  href="src/exp2.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>
指数モデル <tt>Y = a・X^b</tt><p>

<table><tr>
<td><input type="button" name="calcurate" value="計算開始" onClick="calc(this.form.data.value, '指数モデル Y = a・X^b')">  </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/corr2.html">使用法</a></td>
</tr></table>
<p>
<table>
<tr>
<td>入力欄<br><textarea name="data" ROWS=20 COLS=15></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="../lecture/Regression/expreg/type3.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>

サブ io.js   Last modified: Mar 25, 2004

・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI

Made with Macintosh