質問編集履歴
5
タイトル変更
test
CHANGED
|
@@ -1 +1 @@
|
|
|
1
|
-
S
|
|
1
|
+
SIRモデルを扱った新型コロナウイルスの新規感染者シミュレーション 最小二乗法
|
test
CHANGED
|
File without changes
|
4
タイトルの変更
test
CHANGED
|
@@ -1 +1 @@
|
|
|
1
|
-
SEIR 最小二乗法による曲線のフィッティング
|
|
1
|
+
SEIR 新型コロナウイルス 最小二乗法による曲線のフィッティング
|
test
CHANGED
|
File without changes
|
3
解決した画像の添付
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -80,6 +80,5 @@
|
|
|
80
80
|
上記にも述べましたが曲線は折れました。
|
|
81
81
|
|
|
82
82
|
### 補足情報（FW/ツールのバージョンなど）
|
|
83
|
-
|
|
84
|
-
|
|
83
|
+
解決できた画像を添付します。
|
|
85
|
-
|
|
84
|
+
|
2
質問2は解決しました
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -16,12 +16,10 @@
|
|
|
16
16
|
SEIRlist = odeint(SEIReq, (S0, E0, I0, R0), t, args = (beta, ipusilon, gamma))
|
|
17
17
|
return SEIRlist[:,2]
|
|
18
18
|
|
|
19
|
-
2. 今回フィッティングさせて際に t = 1 あたりで曲線が緩やかになっていないのが気になります。
|
|
19
|
+
2. (解決済み)今回フィッティングさせて際に t = 1 あたりで曲線が緩やかになっていないのが気になります。
|
|
20
20
|
今までSIRモデルでフィッティングさせた際に緩やかな曲線を描いていましたのでなぜこのようになるのか
|
|
21
21
|
わかる方がいらっしゃったら教えていただきたいです。
|
|
22
22
|
尚、初期条件を色々いじっても曲線は折れました。
|
|
23
|
-
|
|
24
|
-
|
|
25
23
|
|
|
26
24
|
```
|
|
27
25
|
|
1
def Iの箇所でSEIR[:1] → SEIR[:2]に変更しました。
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -14,7 +14,7 @@
|
|
|
14
14
|
自分なりにodeint関数も少し勉強したのですがイマイチわかっていないと思われます。
|
|
15
15
|
def I(t, beta, ipusilon, gamma):
|
|
16
16
|
SEIRlist = odeint(SEIReq, (S0, E0, I0, R0), t, args = (beta, ipusilon, gamma))
|
|
17
|
-
return SEIRlist[:,
|
|
17
|
+
return SEIRlist[:,2]
|
|
18
18
|
|
|
19
19
|
2. 今回フィッティングさせて際に t = 1 あたりで曲線が緩やかになっていないのが気になります。
|
|
20
20
|
今までSIRモデルでフィッティングさせた際に緩やかな曲線を描いていましたのでなぜこのようになるのか
|
|
@@ -63,7 +63,7 @@
|
|
|
63
63
|
|
|
64
64
|
def I(t, beta, ipusilon, gamma):
|
|
65
65
|
SEIRlist = odeint(SEIReq, (S0, E0, I0, R0), t, args = (beta, ipusilon, gamma))
|
|
66
|
-
return SEIRlist[:,
|
|
66
|
+
return SEIRlist[:,2]
|
|
67
67
|
|
|
68
68
|
optparams, cov = optimize.curve_fit(I, t, normalized_cases, )
|
|
69
69
|
print('R0=',optparams[0]/optparams[2])
|