ABCDEFGHIJKLMNOPQRSTUVWX
1
スプレッドシートを使った新型コロナ感染シミュレーション(その2)
2
「42万死亡推計」の再現
3
4
(1-1) 目的など
5
西浦氏が 2020-04-15 に発表した、「42万死亡推計」をスプレッドシートの計算で再現します。
6
7
また、西浦氏が「42万死亡推計」を説明した Newsweek 記事
8
2020-06-11 「8割おじさん」の数理モデルとその根拠──西浦博・北大教授
9
https://www.newsweekjapan.jp/stories/world/2020/06/8-39.php
10
に掲載のグラフ（↓）を再現します。（雑誌では 2020-06-09 号）
11
12
（グラフが表示されない場合は、「再読み込み」をお試しください。）
13
(figure 1) Newsweek 記事のグラフ
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
縦軸に「10万人当たり」とありますが、これは「人口10万人当たり」であり、
33
「各世代10万人当たり」ではありません。「人口10万人当たり」だと西浦氏のグラフと
34
計算結果がよく一致するからです。「(3)計算とグラフ(dt=0.1)」の計算結果をご参照下さい。
35
36
(1-2) 西浦氏が示しているパラメータが誤っていると考えている事情
37
上の Newsweek 記事で、西浦氏は「42万死亡推計」のアルゴリズムやパラメータの所在を示して
38
います。
39
https://github.com/contactmodel/COVID19-Japan-Reff/blob/master/BerkleyMadonna_May2020.txt
40
どうやらこれは、"Berkeley Madonna" という計算ソフトのスクリプトらしいのですが、私はこの
41
ソフトが手元にないのでスクリプトを直接動作させた確認はしていません。
42
ただしこのファイルの意味するところは、およそ「見れば分かる」と思います。
43
この(その2)では、上のファイルで示されているパラメータと計算式を使っています。
44
ただし私は、西浦氏が別に示している計算の途中経過などを詳しく調べた結果、
45
上のファイルのパラメータの一部に微妙な誤りがあると判断しています。
46
検討の結果、計算の途中経過に適合すると判断できるパラメータを得ることが出来ました。
47
その値と根拠などについては、以下の連ツイで示しています。
48
https://twitter.com/sarkov28/status/1413048029720973312
49
誤ったパラメータの開示は不適切であり、私は不満です。
50
51
(1-3) パラメータの変更方法、コピーや保存の方法
52
このスプレッドシートを「自分の google drive」にコピーすれば、変更したり保存したりできます。
53
gmail のアカウントを持っている人は既に「自分の googlr drive」を持っているのですぐに使えます。
54
コピーするには、画面上方のこのスプレッドシートの名前
55
「スプレッドシートを使った新型コロナ感染シミュレーション(その2)」の下のメニューの
56
「ファイル」の「コピーを作成」を選択して下さい。
57
コピーは、「自分の google drive」の指定したフォルダに作成されると思います。
58
「自分の google drive」は、以下から開くことができます。
59
https://drive.google.com/drive/my-drive?ths=true
60
https://sarkov28.hatenablog.com/entry/2022/04/01/201503 の
61
「開いているグーグル スプレッドシートのコピーを、自分の google drive に作成する方法」
62
にもう少し説明を書きました。関連情報はネットで調べて頂いてもすぐに出てくると思いますが、
63
分からない場合はお問い合わせください。
64
65
(1-4) このスプレッドシートについて
66
3世代SIR 方程式をスプレッドシートで計算し、グラフを描画します。
67
これは、米山隆一 @RyuichiYoneyama さんが 2020-03-29 に Web で示し、
68
twitter で言及した方法シンプルな SIR 方程式を解く方法を、私 @sarkov28 が改変したものです。
69
関連する米山さんの tweet は以下です。
70
https://twitter.com/RyuichiYoneyama/status/1243940595782967302
71
https://twitter.com/RyuichiYoneyama/status/1255291796801847297
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100