メイン density.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>density function</TITLE>
<script src="gxp.js">document.write("gxp.js ファイルが見つかりません??<br>")</script>
<script src="xxp.js">document.write("xxp.js ファイルが見つかりません??<br>")</script>
<script src="txp.js">document.write("txp.js ファイルが見つかりません??<br>")</script>
<script src="fxp.js">document.write("fxp.js ファイルが見つかりません??<br>")</script>
<script src="gamma.js">document.write("gamma.js ファイルが見つかりません??<br>")</script>
<script src="beta.js">document.write("beta.js ファイルが見つかりません??<br>")</script>
<script src="io.js">document.write("io.js ファイルが見つかりません??<br>")</script>
<SCRIPT LANGUAGE="JavaScript">
<!--
function pow(a, b)
{
return a == 0 ? 0 : Math.pow(a, b)
}
function calc1(lower0, upper0, ndiv0, mean0, sd0)
{
var lower = parseFloat(lower0)
var upper = parseFloat(upper0)
var ndiv = parseFloat(ndiv0)
var mean = parseFloat(mean0)
var sd = parseFloat(sd0)
var x, temp, temp2
if (isNaN(lower) || isNaN(upper) || isNaN(ndiv) || ndiv <= 0 || lower > upper || isNaN(mean) || isNaN(sd) || sd <= 0) {
printf("指定に誤りがあります。\n")
return
}
temp = 2*sd*sd
temp2 = Math.sqrt(2*Math.PI)*sd
printf("正規分布の密度関数 N(%g, %g)\n", mean, sd)
printf("%15s %15s %15s\n", "x", "f(x)", "F(x)")
for (x = lower; x <= upper; x += ndiv) {
printf("%15.7f %15.7f %15.7f\n", x, Math.exp(-pow(x-mean, 2)/temp)/temp2, 1-gxp((x-mean)/sd))
}
printf("\n")
sep(60)
printf("\n")
}
function calc2(lower0, upper0, ndiv0, df0)
{
var lower = Math.max(0, parseFloat(lower0))
var upper = Math.max(0, parseFloat(upper0))
var ndiv = parseFloat(ndiv0)
var df = parseFloat(df0)
var x, temp
if (isNaN(lower) || isNaN(upper) || isNaN(ndiv) || ndiv <= 0 || lower > upper || isNaN(df) || Math.floor(df) != df) {
printf("指定に誤りがあります。\n")
return
}
temp = pow(2,df/2)*gamma(df/2)
printf("カイ二乗分布の密度関数 自由度 = %i\n", df)
printf("%15s %15s %15s\n", "x", "f(x)", "F(x)")
for (x = lower; x <= upper; x += ndiv) {
printf("%15.7f %15.7f %15.7f\n", x, pow(x, df/2-1)*Math.exp(-x/2)/temp, 1-xxp(x, df))
}
printf("\n")
sep(60)
printf("\n")
}
function calc3(lower0, upper0, ndiv0, df0)
{
var lower = parseFloat(lower0)
var upper = parseFloat(upper0)
var ndiv = parseFloat(ndiv0)
var df = parseFloat(df0)
var x, temp, temp2
if (isNaN(lower) || isNaN(upper) || isNaN(ndiv) || ndiv <= 0 || lower > upper || isNaN(df) || Math.floor(df) != df) {
printf("指定に誤りがあります。\n")
return
}
temp = -(df+1)/2
temp2 = Math.sqrt(df)*Beta(df/2,0.5)
printf("t 分布の密度関数 自由度 = %i\n", df)
printf("%15s %15s %15s\n", "x", "f(x)", "F(x)")
for (x = lower; x <= upper; x += ndiv) {
printf("%15.7f %15.7f %15.7f\n", x, pow(1+x*x/df, temp)/temp2, x > 0 ? 1-0.5*txp(x, df) : 0.5*txp(x, df))
}
printf("\n")
sep(60)
printf("\n")
}
function calc4(lower0, upper0, ndiv0, df10, df20)
{
var lower = Math.max(0, parseFloat(lower0))
var upper = Math.max(0, parseFloat(upper0))
var ndiv = parseFloat(ndiv0)
var df1 = parseFloat(df10)
var df2 = parseFloat(df20)
var x, temp, temp2, temp3, temp4
if (isNaN(lower) || isNaN(upper) || isNaN(ndiv) || ndiv <= 0 || lower > upper || isNaN(df1) || Math.floor(df1) != df1 || isNaN(df2) || Math.floor(df2) != df2) {
printf("指定に誤りがあります。\n")
return
}
temp = pow(df1, df1/2)*pow(df2, df2/2)
temp2 = Beta(df1/2, df2/2)
temp3 = (df1+df2)/2
temp4 = df1/2-1
printf("F 分布の密度関数 自由度 = ( %i , %i )\n", df1, df2)
printf("%15s %15s %15s\n", "x", "f(x)", "F(x)")
for (x = lower; x <= upper; x += ndiv) {
printf("%15.7f %15.7f %15.7f\n", x, temp*pow(x, temp4)/(temp2*pow(df1*x+df2, temp3)), 1-fxp(x, df1, df2))
}
printf("\n")
sep(60)
printf("\n")
}
// --> JavaScript対応ブラウザでないと動作しません。
</SCRIPT>
</HEAD>
<body bgcolor="#ffffff">
<font size="+2"><b>密度関数</b></font> <a href="src/density.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>
正規分布，カイ二乗分布，t 分布，F 分布において，パラメータを指定することにより，密度関数の値を計算する。<br>
<img src="../gra/button3.png" width=9 height=9 alt="・"> <a href="exa/density.html">使用法</a>
<hr>
<table>
<tr>
<td valign=top bgcolor="#ff5555" align=right>正規分布<br>
平均値 <input name="mean" value="0" size=10><br>
標準偏差 <input name="sd" value="1" size=10><br>
<hr>
開始値 <input NAME="lower1" VALUE="-3" SIZE=10><br>
終了値 <input NAME="upper1" VALUE="3" SIZE=10><br>
刻み幅 <input NAME="ndiv1" VALUE="0.5" SIZE=10><br>
<input type="button" name="cl1" value="クリア" onClick="this.form.lower1.value=this.form.upper1.value=this.form.ndiv1.value=this.form.mean.value=this.form.sd.value=''"> <input type="button" name="calcurate1" value="計算開始" onClick="calc1(this.form.lower1.value, this.form.upper1.value, this.form.ndiv1.value, this.form.mean.value, this.form.sd.value)"><br>
</td>
<td valign=top bgcolor="#ef5555" align=right>カイ二乗分布<br>
自由度 <input name="df_c" value="1" size=10><br>
<hr>
開始値 <input NAME="lower2" VALUE="0" SIZE=10><br>
終了値 <input NAME="upper2" VALUE="10" SIZE=10><br>
刻み幅 <input NAME="ndiv2" VALUE="0.5" SIZE=10><br>
<input type="button" name="l2" value="クリア" onClick="this.form.lower2.value=this.form.upper2.value=this.form.ndiv2.value=this.form.df_c.value=''"> <input type="button" name="calcurate2" value="計算開始" onClick="calc2(this.form.lower2.value, this.form.upper2.value, this.form.ndiv2.value, this.form.df_c.value)"><br>
</td>
<td valign=top bgcolor="#df5555" align=right>t 分布<br>
自由度 <input name="df_t" value="1" size=10><br>
<hr>
開始値 <input NAME="lower3" VALUE="-3" SIZE=10><br>
終了値 <input NAME="upper3" VALUE="3" SIZE=10><br>
刻み幅 <input NAME="ndiv3" VALUE="0.5" SIZE=10><br>
<input type="button" name="cl3" value="クリア" onClick="this.form.lower3.value=this.form.upper3.value=this.form.ndiv3.value=this.form.df_t.value=''"> <input type="button" name="calcurate3" value="計算開始" onClick="calc3(this.form.lower3.value, this.form.upper3.value, this.form.ndiv3.value, this.form.df_t.value)"><br>
</td>
<td valign=top bgcolor="#cf5555" align=right>
F 分布<br>
第１自由度 <input name="df1" value="1" size=10><br>
第２自由度 <input name="df2" value="1" size=10><br>
<hr>
開始値 <input NAME="lower4" VALUE="0" SIZE=10><br>
終了値 <input NAME="upper4" VALUE="10" SIZE=10><br>
刻み幅 <input NAME="ndiv4" VALUE="0.5" SIZE=10><br>
<input type="button" name="cl4" value="クリア" onClick="this.form.lower4.value=this.form.upper4.value=this.form.ndiv4.value=this.form.df1.value=this.form.df2.value=''"> <input type="button" name="calcurate4" value="計算開始" onClick="calc4(this.form.lower4.value, this.form.upper4.value, this.form.ndiv4.value, this.form.df1.value, this.form.df2.value)"><br>
</td>
</tr>
<tr>
<td nowrap colspan=4><br>出力欄 <input type="button" name="clear" value="出力欄クリア" onClick="this.form.result.value=''"><br><textarea name="result" rows=30 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>
サブ gxp.js Last modified: Mar 25, 2004
サブ xxp.js Last modified: Mar 25, 2004
サブ txp.js Last modified: Mar 25, 2004
サブ fxp.js Last modified: Mar 25, 2004
サブ gamma.js Last modified: Mar 25, 2004
サブ beta.js Last modified: Mar 25, 2004
サブ io.js Last modified: Mar 25, 2004
直前のページへ戻る E-mail to Shigenobu AOKI
