メイン gomperts.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">
<!--
var para_a, para_b, para_c
var epsinv = 1e-8
function mreg(data, nc)
{
var mean, variance, sd, ss, r, nv, i, j, k, b0, b, nv1, temp2
nv = 3
mean = new Array(nv)
variance = new Array(nv)
sd = new Array(nv)
ss = new Array(nv)
b = 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])
}
}
}
for (j = 0; j < nv; j++) {
ss[j] = r[nv-1][j]
variance[j] = r[j][j]/(nc-1)
sd[j] = Math.sqrt(variance[j])
}
for (j = 0; j < nv; j++) {
for (k = 0; k < nv; k++) {
r[j][k] /= (nc-1)*sd[j]*sd[k]
}
}
nv1 = nv-1
if (inv(r, nv1, epsinv) == 1) {
printf("逆行列を求められません\n")
return 1
}
b0 = mean[nv1]
for (i = 0; i < nv1; i++) {
temp2 = 0
for (j = 0; j < nv1; j++) {
temp2 += r[i][j]*r[j][nv1]
}
b[i] = temp2*sd[nv1]/sd[i]
b0 -= b[i]*mean[i]
}
para_a = b0
para_b = b[0]
para_c = -b[1]
return 0
}
function init(data, nc)
{
var m = Math.floor(nc/3)
var i, s1, s2, s3
s1 = s2 = s3 = 0
for (i = 0; i < m; i++) {
s1 += Math.log(data[i])
s2 += Math.log(data[i+m])
s3 += Math.log(data[i+2*m])
}
return Math.log((s1-s2)/(s2-s3))/m
}
function calc(data_string)
{
var data, nc, y, i, j, d, u0, u1, loop
if ((data = getdata(data_string, 0)) == false) return
nc = data.length
if (nc <= 3) {
printf("有効データが %i 個しかないので，分析できません\n", nc)
return
}
y = new Array(nc)
for (i = 0; i < nc; i++) {
if (data[i][0] <= 0) {
printf("従属変数は正の値でなくてはいけません\n入力行： %g\n", data[i][0])
return
}
y[i] = data[i][0]
}
d = new Array(nc)
for (i = 0; i < nc; i++) {
d[i] = new Array(3)
}
printf("★ ゴンペルツ曲線へのあてはめ ★\n\n")
u1 = u0 = init(y, nc)
printf("初期値 b = %.6g\n", u0)
for (loop = 1; ; loop++) {
if (loop % 20 == 0 && confirm("反復回数が "+loop+" 回になりましたが，まだ収束しません。続けますか？") == false) break
for (i = 0; i < nc; i++) {
d[i][0] = Math.exp(-u1*(i+1))
d[i][1] = (i+1)*Math.exp(-u1*(i+1))
d[i][2] = Math.log(y[i])
}
if (mreg(d, nc) == 1) {
return
}
u0 = u1
u1 += para_c/para_b
if (Math.abs((u0-u1)/u0) < 0.000001) break
}
para_a = Math.exp(para_a)
para_b = Math.exp(para_b)
para_c = u1
printf("予測値 = %g・%g ^ exp( - %g x)\n\n", para_a, para_b, para_c)
printf("%12s%12s%12s\n", "X", "Y", "予測値")
for (i = 0; i < nc; i++) {
printf("%12.6g%12.6g%12.6g\n", i+1, y[i], para_a*Math.pow(para_b,Math.exp(-para_c*(i+1))))
}
}
//-->
</script>
</head>
<body bgcolor="#ffffff">
<font size="+2"><b>ゴンペルツ曲線へのあてはめ</b></font> <a href="src/gomperts.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>
入力欄に，従属変数の値を１行に一つずつの値を<a href="exa/kugirimoji.html">区切り文字</a>で区切って入力し，ゴンペルツ曲線 <tt>Y = a*b^exp(-c*X)</tt> にあてはめを行う（本プログラムは線形近似によるものなので，<a href="fit-gomperts.html">非線形回帰</a>を行う方が望ましい）。<br>
。<br>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="../lecture/Regression/growth/gomperts.html">手法の解説</a> <img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/gomperts.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>
<table><tr>
<td>入力欄<br><textarea name="data" ROWS=15 COLS=15></textarea></td>
<td>出力欄<br><textarea name="result" ROWS=15 COLS=80></textarea></td>
</tr></table>
</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
