質問編集履歴
2
バグフィックス
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -42,10 +42,6 @@
|
|
|
42
42
|
up = 2*sigma**2
|
|
43
43
|
down = (1/2)*sigma**2
|
|
44
44
|
|
|
45
|
-
# カイ二乗分布の方程式は計算に非常に時間がかかるのでテイラー展開(n=20)
|
|
46
|
-
print("expanding")
|
|
47
|
-
f = sum([f.subs(x, x+(up+down)/2).taylor_term(i, x) for i in range(20)])
|
|
48
|
-
|
|
49
45
|
# 評価するパラメータ
|
|
50
46
|
par = s**2
|
|
51
47
|
|
|
@@ -57,6 +53,10 @@
|
|
|
57
53
|
up = up*g
|
|
58
54
|
down = down*g
|
|
59
55
|
|
|
56
|
+
# カイ二乗分布の方程式は計算に非常に時間がかかるのでテイラー展開(n=20)
|
|
57
|
+
print("expanding")
|
|
58
|
+
f = sum([f.subs(x, x+(up+down)/2).taylor_term(i, x) for i in range(20)])
|
|
59
|
+
|
|
60
60
|
# 分布関数を定積分
|
|
61
61
|
print("now integrating")
|
|
62
62
|
F = sym.integrate(f.subs(x, x-(up+down)/2), (x, down, up))
|
1
バグフィックス
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -38,14 +38,14 @@
|
|
|
38
38
|
f = sst.ChiSquared(x, n_2-1)
|
|
39
39
|
f = sst.density(f)(x)
|
|
40
40
|
|
|
41
|
+
# 上限値と下限値の設定
|
|
42
|
+
up = 2*sigma**2
|
|
43
|
+
down = (1/2)*sigma**2
|
|
44
|
+
|
|
41
45
|
# カイ二乗分布の方程式は計算に非常に時間がかかるのでテイラー展開(n=20)
|
|
42
46
|
print("expanding")
|
|
43
47
|
f = sum([f.subs(x, x+(up+down)/2).taylor_term(i, x) for i in range(20)])
|
|
44
48
|
|
|
45
|
-
# 上限値と下限値の設定
|
|
46
|
-
up = 2*sigma**2
|
|
47
|
-
down = (1/2)*sigma**2
|
|
48
|
-
|
|
49
49
|
# 評価するパラメータ
|
|
50
50
|
par = s**2
|
|
51
51
|
|