数学 統計に詳しい人が語るコロナウイルス

1 名前：１３２人目の素数さん [2020/02/29(土) 02:18:41 ID:twdO677Q.net] 東大数学科卒の元官僚はこう分析してるが、お前らはどうなると思う？ www.zakzak.co.jp/soc/news/200220/dom2002200003-n2.html 中国国外感染者の中国国内との比率をみると、 １月２０日の数字公表以降は、０・８〜２・６％で比較的安定している。 これは、新型肺炎の感染者のほとんどは中国国内、それも湖北省に集中しているからだ。 ちなみに中国国外での感染者数は、中国国内の１・１％だ（２月１６日現在）。 本コラムで紹介したが、現時点では、最終的な中国国内の感染者数は２０万人超と筆者は推計している。 となると、中国国外の感染者は数千人程度になるだろう。 中国国外のうち日本の比率は１割弱なので、日本の感染者数は数百人程度であろう。 その場合、死者も数人から１０人程度になるだろう。 こうした推計をすると、今の感染者は氷山の一角だと思われるが、今後の増加ペースはどうなるだろうか。 新型コロナウイルスの検査は簡単に行えるので、今後、日本での感染者数は増えていくだろう。 ある時点ではそれがネズミ算的に増えるかのように思える局面もあるだろうが、 筆者の推計が正しければ、現時点ではせいぜい数百人が一つのメドだ。 830 名前：１３２人目の素数さん [2020/07/12(日) 10:18:36.16 ID:JrtRNErT.net] >>780 続き >SIR微分方程式と原子核崩壊 ・SIRモデル。 　感染症数理モデルは� 831 名前：Vンプルで 　集団モデルのSIRが基本モデルだ。 　微分方程式は(ミクロな)個人間の接触　 　が均等より非線型項βSIがあり 　非線型。 ・原子核崩壊。 　放射平衡は 微分方程式d N1/dt=-λ1N1 　が基本で相対論的場の量子論. 　断面積から導かれる。 　親→娘→孫の線型方程式。 ・SIRモデルの局所線型化。 　SIRモデルで初期t≒0の時 　S(t)≒S(0)=1 で線型方程式になる 　SIRモデル(t≒0)と放射平衡は 　共に線型微分方程式として 　対応できる ・マスク効果は個人間接触の相互作用. 　断面積パラメータで扱え 　集団のSIRモデルで感染率パラメータ 　となる。 　ミクロ相互作用の接触は密度に依存 　する 以上は数理モデルの概論です 　 [] [ここ壊れてます] 832 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 10:41:01.76 ID:9O0Ntz/e.net] >>780 レスありがとうございます。 K値とは合羽値の頭文字出会ったとｗ >694のκ値とは別物。 833 名前：１３２人目の素数さん [2020/07/12(日) 10:57:43 ID:JrtRNErT.net] >>782 K値は根拠不明 834 名前：１３２人目の素数さん [2020/07/12(日) 14:27:41.83 ID:yl7MdC0U.net] 工学の人が好きそうな、いきあたりばったり的なパラメータだよね＞K値 理学系なら演繹的に捻り出すモデルを考えろよ、と言いたい。 835 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 14:54:44 ID:qAF2CiMg.net] モデルを考えている間に感染が広がったり経済的損失が拡大してしまったら役に立たない 病気の原因はわかっても患者が死んだら意味がない 学問的には意味があるだろうけど実践的には意味がない モデルを考えるのは後からでもいい 836 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 15:17:13 ID:Wf1trDX+.net] 日本に必要なのは、都道府県別実効再生産数を1未満に保つ、 仮に1を上回ったら2週間以内に1以下に下げる っていう数値目標を加えること 実効再生産数が1を超えると2週間後に感染者が目に見えて増大、4週間後には急拡大する 感染者数とかだけで見てると先手先手をうった対応が出来ず後手後手にまわるのは いまの東京を見れば明白 東京は5月後半に実効再生産数が継続的に1を超えて、 感染拡大するのは明白だったのにろくな対策しなかった 837 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 15:19:19 ID:Wf1trDX+.net] 実効再生産数の目標には外部要因も考慮する必要がある たとえば、自分ところの都道府県は実質的には実効再生産数が1以下で、 自分ところの都道府県内での感染は抑制されてるが、 他県、たとえば東京からどんどん感染者が流入してくるので数値的には1を超えた、みたいな場合ね 838 名前：１３２人目の素数さん [2020/07/12(日) 17:30:49.38 ID:PYIegNFg.net] >SIRモデルと再生産数 潜伏期間が0や住民間で均等な接触 を仮定し集団を対象にしたシンプルな 非線型微分方程式。 生物の基本的な数理モデル ・防疫戦略の基本的な道具の一つ ・インドボンベイのペスト流行でも 　流行曲線を再現した ・体内のウィルス感染モデルにも 　使用されている。 ・初期に非感染者S≒1の場合 　d I/dt=(β-γ) I =(Ro-1)γI 　基本再生産数Ro≡β/γ 　Ro>1の時　感染者 Iが増加 ・実効再生産数R 　感染原因が住民間の接触だから　 　接触を減らせばRが減少する。 　8割接触が減少すれば 　R=Ro×(1-0.8)=0.2R だからマスクをして三密はさける。 (物理公式と同じく意味もわからず 公式へ数値を入れ結果を振り回さない、 Rは指標の一つでデータの取り方に 依存する>マスコミ素人学者政治家) ・都道府県単位で平均された 　データ指標には注意が必要。 　ウィルス拡散は放射能拡散と 同様に県境がないし業種交通網 　年齢他にも依存する ・検査対象データは検査数対象 　特性に依存するし東京以外は 　戦略的検査数が少なく過ぎると 　思われる 　抗体検査では大阪>東京>宮城 　だった ・SIRモデルから潜伏期間EのSEIR モデル、進行波解がある(SIR＋拡散 項)モデル、再帰性があるインフル 　エンザモデルなどがあり個人行動 　をシュミレーションしたibmモデル 　も以前からある。 　新型インフルエンザ対策として 　政府へ提言されていた 　マスク備蓄もね 839 名前：１３２人目の素数さん [2020/07/12(日) 17:31:39.31 ID:PYIegNFg.net] >>787 検査していますか？ 840 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 17:43:16.95 ID:PLqlj++l.net] >>788 海外から流入するパラメータの違うIを加えたモデルってあるかな。 841 名前：１３２人目の素数さん [2020/07/12(日) 18:02:36 ID:PYIegNFg.net] >>790 移流項＋SIRモデルやibmモデルかな。 SIRモデルは人口一定の閉鎖地域だから 感染者Iの流入流出項を加える。 (物理化学モデルの類似なら 閉じ込められた分子が衝突して化学 反応をしている。 更に外部から分子が流入するとか) 西浦教授が海外からの水際対策を説明 するため玉を流してぶつけていたが 西浦教授に聞いてください。 842 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 18:06:06 ID:4t109+lN.net] >>787 その数が分かっているなら、Rtの分子から流入分は除くべき。 西浦先生の計算では、輸入分の患者はそう取り扱っていたと思う。 843 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 18:18:56.02 ID:PLqlj++l.net] >>791 レスありがとうございます。 SEIRに人口増加の項を加えたモデルはみたことがあるけど IBMって今はレノボではと思ったら別のIBMだったｗ こういうのがみつかりました。 https://bradduthie.github.io/blog/individual-based-models-in-r/ 理解できるかわからないけど読んでみます。 844 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 21:07:07.01 ID:PLqlj++l.net] RにIBMのパッケージがあった。 https://cran.r-project.org/web/packages/ibm/index.html 845 名前：１３２人目の素数さん [2020/07/13(月) 09:41:22 ID:dpBf9py+.net] >>785 だから、実践なんか二の次で真理を求めるのが理学。 846 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 19:17:11.63 ID:ux2hxf89.net] K値による予測を使うべきでない理由 katukawa.com/?p=6228 847 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 19:21:00.33 ID:FAV+lxAc.net] >>795 真理って言っても確率的に大雑把にしか言えないと思う 個々の人から排出される飛沫にウイルスがどの程度含まれて それかが別の人の体内にどの程度入って その体内でウイルスがどう動いて免疫がどう働いて 結果的に症状や他の人に感染させるかとか 具体的な事は複雑で判らないと思う 真理にはたどり着けない 848 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 19:33:06.81 ID:tR49RdlH.net] まぁ何をどう言い繕ってみても再生産数を1未満、最低でも1以下にしない限り社会は崩壊する。 再生産数が1より大きい状態では何をどうやっても社会を維持することはできない。 849 名前：１３２人目の素数さん [2020/07/13(月) 20:39:03.48 ID:loo9MSbm.net] >>796 >K値というものが物理学者によって 提唱されて、大阪などでは利用されて いるようです 阪大核物理 Nが提唱し大阪府がK値を 利用していると明記してください。　 京大RIMSのMグループはIUTが数学 のabc予想から天気予報の改善まで 有効と主張しています。 関西の数理科学 850 名前：はアホくさく悲惨です がコロナ対策は命を守る問題です。 数学と物理は感染症分野は素人だから 生物数理モデルから個人的に学ぶべき であくまで素人の範囲でしょう。 生物数理モデルは故山口昌哉京大教授 がカオスを含め解析し弟子も活躍し 今では世界的に正統派です。 西浦北大教授の感染症対策の仕事や話も 山口先生の仕事の発展があると思って います。 [] [ここ壊れてます] 851 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 20:46:24.24 ID:DPFA8Q9h.net] All analyses were performed using R version 3.6.3 (R Foundation). なんて書いてあるんで読んでみたけど 割り算に使っただけみたいだ。 https://jamanetwork.com/journals/jama/fullarticle/2768351 対照群もなんにもない。小学生に夏休みの自由研究みたいだ。 852 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 20:48:12.18 ID:DPFA8Q9h.net] >>800 平均値と標準偏差の計算に使っただけだろうな。 853 名前：１３２人目の素数さん mailto:sage [2020/07/14(火) 14:27:26.92 ID:cML8ukUt.net] K値は基本役に立たない 一番役に立つのは実効再生産数 実効再生産数をモニタリング指標、および数値目標にして管理するのが最適 K値予測が当たるのは、一時的に実効再生産数が1を超えて、また1以下に戻る、みたいな場合 もちろん、実効再生産数が1を超えたら行政がすぐい対応して1以下に減らす対策をとることが前提 行政が1か月以上実効再生産数を減らす対策をとらない東京には当てはまらない 854 名前：１３２人目の素数さん [2020/07/14(火) 22:48:00.75 ID:MNivzfqp.net] SIRモデルによるCov19対策 ・法令システム.医療システム 　行政システムによる。 　SIRモデルから状況分析と予測より 感染者集団Ｙを減少させる つまりd Ｙ/dt<0にする。 SIRモデルで( S.I.R)→(X. Y. Z) とした。 ・SIRモデルは集団免疫モデルで 　対策がない場合を予測する。 　対策を考える場合も使える 　第1式　新規感染者の増加　 　第2式　感染対策の効果 　　　　=新規感染者-隔離感染者 　第3式　医療システム ・対策のための実効再生産数R 　定性的な式 　個人間は均等な接触から R→Ro×(個人間の接触回数割合) =(感染率)/(回復率)]× (個人間の均等な　　接触回数割合) SIRモデルのβは広義の(感染率)× (個人間の均等な接触回数割合) ・第2式で dＹ/dt<0のためには 　β感染率→小　接触回数→小 　γ隔離率→大　 　結果　Rは減少する γ→大の対策について ・第2式　 d Ｙ/dt=βS Y-γＹ 戦略的PCR検査対策。 隔離率には抗体獲得.ワクチンと 検査隔離がある。 抗体集団検査が約大阪015>東京0.1> 宮城0.01パーセントだから集団免疫 なしかつワクチンがない。 クラスターが予想される地域店などへ 特にPCR検査を増やし感染者を 病院.ホテルへ隔離することは隔離者が 増加し Yが減少する。 ・第3式医療システム d Z/dt=γY隔離感染者 検査隔離(γ→大)で 隔離者増大より病床数が満杯になり 結果医療システムが崩壊するから 軽症ならホテル隔離へ。 東京方針はRを小さくする。 ・大阪K値について 　データへ指数関数を当てはめた 現象論に基づくから対策に繋がらない 　戦略的PCR検査なし 本音は 経済優先でコロナはただの風邪路線 855 名前：１３２人目の素数さん [2020/07/15(水) 11:05:27.21 ID:+c3dKonF.net] >>797 単純な原理から、いかに自然現象をより高い精度で近似できるかって ことが、自然現象を理解するってことでしょ。科学における「真理」 という言葉は、そういうただし書きつきだと理解しないと。 そこらへんは数学とは事情が異なる。 自然科学において、絶対的な真理など空疎な言葉。そんなものは、 宗教屋が唱えるおまじないでしかない。 856 名前：１３２人目の素数さん mailto:sage [2020/07/15(水) 14:49:16 ID:nLcbkn+C.net] >>804 同じようにある条件下で成り立つ法則を数式などで数量も含めて表現したものと思う ニュートン力学は量子領域や宇宙領域では成り立たないだろうし 量子領域の物理で人間サイズを扱おうとするとシュレーディンガーの猫みたいになると思う ウイルスと人間の免疫の活動の結果として健康に大きな影響なければ普通の風邪などと同様に無視していいからミクロな視点も無視はできない ウイルスが体内でどのように増殖するのかも知るのも真理を知ることにはなると思う マクロはミクロの積み重ねだから 個人の行動変容が感染拡大に影響するのも確かなこと 同じウイルスと人間の活動をどんな観点で記述するかの違いの様に思う 同じ現象だから互いに整合するはずだと思う 857 名前：１３２人目の素数さん mailto:sage [2020/07/15(水) 17:41:50.43 ID:Y5eCgaZb.net] K値予測 katukawa.com/wp-content/uploads/fig6-1024x574.png 858 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 05:27:26 ID:f7C3/cjh.net] PCR検査、アメリカでは200人だかの検体をまとめて検査してるそうだ。 つまり検査結果が陰性なら全員がセーフ、陽性なら誰か陽性の人が1人以上いるぞ！ってこと 複数の検査を 検査1,2,3,4,5,6,7,8,9,10 検査a,b,c,d,e,f,g,h,i,j のように受けることでシロの人の洗い出しが素早くできたりはしないでしょうか またはその逆にクロのひとの洗い出し 859 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 07:13:17.08 ID:C2pOqxoK.net] >>807 サンプルを混ぜる手法はバランスパズル・偽コイン問題やiPS細胞の実験手法で有名よね www.jscf.org/jscf/katudou/Walk_Again_2008/08_10_5-1.pdf 860 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 08:22:41.73 ID:muripj9V.net] >>807 >>633の式を当てはめると、アメリカは陽性率の事前確率を25/100万ぐらいに見積もってるのだな。 861 名前：１３２人目の素数さん [2020/07/16(木) 09:19:06 ID:dzMbDhwI.net] >>805 マクロなシステムの振る舞いをミクロな構造から説き起こそうとする 演繹的な方法論には限界があるが、それでもそれしかないからね。 モデル化ってのまさにそういうことでしょ。 マクロなシステムの予測に役立つk値のようなパラメータを経験的に 見出すような、帰納的なやり方もありなのは確かにそうなんだろう。 だけど、熱力学が統計力学で置き換えられたように、そこでとどまっ ていては現象を理解した気にはなれんのよね。 やはり普遍的な基本法則から（たとえ数値シミュレーションの結果 であっても）演繹的に導かれるような説明ができないと理解した気 になれない。 862 名前：１３２人目の素数さん [2020/07/16(木) 09:22:20 ID:dzMbDhwI.net] >>807 へー、実際にやってるんだ。 863 名前：１３２人目の素数さん [2020/07/16(木) 14:37:18 ID:tJZRomem.net] >>803 訂正 ばつ.大阪015>東京0.1> 宮城0.01 まる.大阪0.17>東京0.1>宮城0.03 864 名前：１３２人目の素数さん [2020/07/16(木) 18:04:59.93 ID:d2bliECv.net] K値モデルによる予想。 >8府県(大阪兵庫京都愛知神奈川 埼玉千葉福岡) 7月9日ごろにピークアウト →大阪府認定のk値だが7月15日 　大阪61人では、大風呂敷を広げて 　根拠に欠ける曲線を入れても予測 　能力に全く欠ける。 　現象論的にも過去に検証がなく 　思いつきのレベルだが大阪府が　 　K値を全国へばらまいた 　第1波の厚労省集団抗体検査で 　大阪府が約0.17パーセントから 　推測する。 　府民882万人より15000人の潜在 　感染者がいて隔離感染者2126 人　　 (7.14)だから大阪府に潜在感染者が　　 　12.8千人もいた。 　医療システムが脆弱では戦略的PCR 　検査が不可能でk値やメディア出演 　のパフォーマンスよりその場を 　しのいできた。 　大阪夜の街ホストクラブで感染は 　徳島の指摘から発覚した。 　大阪府は関空の水際の規制緩和を 　要望など経済優先に前のめり 　よく知らんが都構想の都合がある 　のだろうか 865 名前：１３２人目の素数さん [2020/07/16(木) 18:11:15.37 ID:d2bliECv.net] >>813 つづき 移流項があるSIRモデル。 ・移流項は個人が一様性等方の運動を 　して個人間の感染,相互作用が仮定 される、 　ブラウン運動も一様等方の運動。 　移流項d∇^2(S、 I、R)より 　∂S/∂ t=-βSI＋d(∇^2)S 　 ∂ I/∂t=βSI-γY＋d(∇^2 )I 　∂R/∂t=γY ＋d(∇^2)R まず移流項なし閉鎖地域のSIR方程式を 解くこと。 解き方は力学系.近似解.数値計算などが あり力学系ベクトル場から閾値の存在が すぐわかる。 Murray数理生物学は感染症の章に ペスト狂犬病などがあり 空間.時空間の発展モデルがある。 更にロジスティック曲線もある。 念のためゴンペルツ曲線は全くない。 866 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 22:14:21 ID:lhXKpfMQ.net] >>807 検体をまとめて調査する手法は、検査能力が足りないときに検査数拡大したいときによく使われる手法 中国では、武漢の集団PCR検査で使ってた 韓国も、3月ころの検査数足りないときにやってた しかし200人ってのはちょっと信憑性かけるのでは？ アメリカの陽性率だと、200人もまぜたらほぼ全部陽性になって検査にならない 867 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 22:19:18.21 ID:lhXKpfMQ.net] 日本は学校なんかの全生徒検査とかがあれば、 検体混ぜるのを使えばいいとおもうよ 感染者と同じクラスの児童や担任、クラス外で親密な人→検体まとめずに個別検査、 それ以外→検体混ぜて検査 みたいな感じで ただし、紙とFAXで検査してる日本じゃ、検体ミックスなんてやったら保健所がパンクしたりミスがたくさん出そう ちゃんとコンピュータ化してれば、検体ミックスが有効に使える 868 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 05:10:32.66 ID:l8XG6+rW.net] 水中でも感染力を失わないならやばそう。 https://www.medrxiv.org/content/10.1101/2020.06.17.20133504v1 869 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 13:35:13 ID:nAA2dWKD.net] K値モデル完全失敗、西浦予想に近い形になる K値モデルは、感染者が増えたらあわてて政府・行政が対策するなら、K値モデルに近くなる 現実は、感染者増えても自治体も政府も無対策どころかさらに感染者ふやす対策をとる 870 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 15:50:29.96 ID:3kM9UPoB.net] >>815 >>816 毒ワイン問題の帰結から類推すると、陽性と出るか陰性と出るか1/2の確率になる� 871 名前：ｭらいの 人数分を混ぜて検査するのが、最適な方法ではないだろうか。 [] [ここ壊れてます] 872 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 16:10:22.26 ID:l/H5WeF7.net] >>819 >>419 873 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 16:30:29.34 ID:PvhkH04x.net] 誰かk値理論の論文読んで解説してくれん？ なんかアホらしくて探す気にも読む気にもならん。 1) 感染者数の推移についてどんなモデル曲線を採用してるのか、どんなパラメータを持ってるのか 2) そのパラメータをどんなモニタリング指標を使ってどう統計処理して推定するのか？ もちろん2)はk値とかいうのを使うんやろうけど 874 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 16:58:10.33 ID:l/H5WeF7.net] ツイッターかnoteかで解読してるヤツおらんかった？ 875 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 17:04:19.18 ID:PvhkH04x.net] なんかどうせデタラメやろうから調べる気にすらならん。 検索して読む時間全部無駄になる予感しかせん 876 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 17:08:29.39 ID:l/H5WeF7.net] 実際に大阪とかはK値モデル採用して意思決定してるんだから、その予測が出来るから無駄ではないと思うけど。 877 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 17:17:55.08 ID:PvhkH04x.net] つまりは大阪がどんなデタラメ理論に基づいて動いてるかわかる事に意味があると？ イヤ、ちょっとご勘弁w 878 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 17:48:59 ID:l/H5WeF7.net] >>825 発言権を得て大阪以外も追随する可能性があるし、トンデモ理論の成り立ちとかわかるから面白いと思うけどね。（大阪専門家委員にはウイルス学の宮坂いるんだけどこの人もちょっとアホだからK値に感心している。あと山中伸弥教授もいたようだから、K値と接触して宣伝する可能性があって危険。） 879 名前：１３２人目の素数さん [2020/07/17(金) 18:13:37 ID:WeNCWXbE.net] >>802 >行政が1か月以上実効再生産数を 減らす対策をとらない東京には 当てはまらない デマですね。 あなたはなぞのK値でも信奉してください 東京都は戦略的にPCR検査を増やし 隔離感染者を増大していますね。 これは基本的なSIRモデルで有効な 対策で実効再生産数を小さくします 880 名前：１３２人目の素数さん [2020/07/17(金) 19:43:00.71 ID:Ew4KkMJN.net] >>821 >>796にK値のポイントが整理されて いますね。 「論文」も検索できるし 881 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 21:58:21.84 ID:4GCjkvh6.net] >> 882 名前：828 なんかアホらしすぎて読む気にならんのよねぇ。 ただ例え便所の落書きとはいえ、k値理論に口出しするなら、やっぱ読まんといかんかなぁと言う気もするし。 なんでこんな怪しさ満載の理論採用されてるんかねぇ？ そもそも実績としてもことごとく予測値外してるのに [] [ここ壊れてます] 883 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 22:53:16.49 ID:l8XG6+rW.net] >>821 >796にでてくるゴンペルツ曲線 https://en.wikipedia.org/wiki/Gompertz_function ここに利用法として Examining disease spread と書いてあるな。 884 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 23:19:19.44 ID:CMxFmiKu.net] >>830 thx 未知パラメータ３つか。 コレをK値で推定すると。 しかしコレ増加しっぱなし、減少しっぱなしやん。 それでみんなどこから起点日に取るかの恣意性があるって言ってんのか。 それとも起点日もパラメータにしてちゃんと推計値出すのかな？ 読んだ？と思われる人のレスだと無さそうだなww 885 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 00:58:41 ID:nqMg34Lo.net] >>827 東京はクラスター対策がいい加減すぎて、PCR検査で実効再生産数を十分低くすることは難しい 中国並みに厳格なクラスター対策しない限り、東京のようなザルクラスター対策では 大して実効再生産数下げられない 886 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 01:04:31 ID:nqMg34Lo.net] ちなみに北京は、6月中旬からはじまった感染クラスターで、 普段の感染対策と、大量のPCR検査を組み合わせて抑え込みに成功した 高リスクグループから� 887 名前：�に毎日数十万件とかいう検査をおこなった結果抑え込めたわけで、 いまは低リスクグループの大量検査を実施中で、すでに検査数が1000万件を超えてる [] [ここ壊れてます] 888 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 07:39:06 ID:kzNofEnA.net] >>831 ゴンペルツ曲線って新規発生数じゃなくて累積患者数を表しているんじゃないかな？ ゴンペルツ分布というのがあるらしいが 新規発生数=pdf 累積患者数=cdfとすると パラメータによってpdfは単調にも一峰性にもなるみたい。 https://en.m.wikipedia.org/wiki/Gompertz_distribution 889 名前：１３２人目の素数さん [2020/07/18(土) 07:49:49 ID:b9ODbCJm.net] >>832 中国は共産党独裁国家だし 大阪は検査もろくにしなかった。 法令システムでは 措置法は規制がなくザルで厚労省は 食中毒並みのPCR検査義務と休業命令 も通達しない。 国、政府、国会の法改正責任が重大なのにGO TO キャンペーンへ、 まあコロナはただの風邪路線なのだろう 890 名前：１３２人目の素数さん [2020/07/18(土) 07:55:01 ID:b9ODbCJm.net] >>829 論文を読むとインチキ手口がわかるよw 891 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 08:49:38.59 ID:bK/VOAbE.net] >>834 そうかもね。よんだ？ >>836 自信が確信になってもなぁw 892 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 09:17:12 ID:bK/VOAbE.net] >>830 しかしなぁ sir方程式の解はピーク時点は左右対称になんかならんからな。 もちろん現実世界では再生産数は社会の対策で変動してるからsirの解通りには変動しないけど。 なのでどこでピークアウトするかは理論的に予測できるハズはほとんどないんだけどな。 893 名前：１３２人目の素数さん [2020/07/18(土) 10:06:46 ID:b9ODbCJm.net] >>837 K値論文の仮定が根拠不明だ。 SIRは数値計算を行なっている。 実際の測定値は検査などによる隔離感染者数だな 感染者数Iは推定で対策よりγが 変わればSIR予測も変わる。 894 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 11:49:40.36 ID:ILbvGgBu.net] >>839 まぁまだ論文読んでないけど(そして読む気も消え失せたけどw)大体インチキなのは確定したな。 モデルの曲線に何とるにせよ、そのモデルの曲線はなんかの意味もって選定しないと話にならん。 sirのモデルではsir方程式 d/dt I = β SI - γ I (感染者の増加分は感受性宿主数と感染者数に比例し、減少分(死亡、回復する数)は感染者数に比例する をとる。 「そやろな」と納得いく式だし、実績もある。 この方程式に従えば感受性宿主数がγ/βに到達するまでピークアウトしない。 にもかかわらず実際の第一波でははるかに早い段階でピークアウトしてる。 それは感染率βが定数ではないからだ。 行動制限などによって接触率が下がれば感染率βも下がる。 第一波が早期にピークアウトしたのは緊急事態宣言下で市民が接触率を下げた効果。 つまりどこでピークアウトするかは市民の生活様式によるものでとても現在の新規感染者数から予測できるようなものじゃない。 もうそれを現在の新規感染者すうから計算できるK値なるもので予測できると言ってる時点でインチキ確定と言っていいやろ。 もう全く読む気が失せたww 895 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 11:51:05.18 ID:ILbvGgBu.net] あ、達するまでというか減るまでね。 896 名前：１３２人目の素数さん [2020/07/18(土) 14:36:31 ID:+nayeGtY.net] K値論文 Indicator K to analyze the transition of COVID-19 spread In this study, we introduce a new indicator called the K value defined by K(d) = 1−N(d − 7)/N(d), where d is the number of days from the reference d 897 名前：ate, and N (d) and N (d − 7) are the total number of infected people on days d and (d − 7), respectively. Since N (d) is greater than N (d − 7) during the period from the initiation of spread to convergence, K takes a value between 0 and 1. ( infected people 感染者) k = 1 + 2.88K′ (Fig. S1). この公式を図へあてはめたのか、 図から帰納的に公式を導いたのか、 特に2.88の由来はなんだろう？？ [] [ここ壊れてます] 898 名前：１３２人目の素数さん [2020/07/18(土) 23:28:48.51 ID:b9ODbCJm.net] K値について調べてみました https://tatsuharug.com/kvalue 899 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 01:35:49 ID:m6+VaLKl.net] >>842 Indicator K の K はなんの K ? 900 名前：１３２人目の素数さん [2020/07/19(日) 06:37:20 ID:GNq9W41i.net] KはあまカッパのKからともいわれている が、よくしらない。 阪大か大阪府へきいてみたら 901 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 07:29:35 ID:w6BaQv4H.net] >>842 R0 = 2.5であとは適当にパラメータ設定したSEIRモデルにK値を重ねてグラフにしてみた。 https://i.imgur.com/5AsiQxW.png どうやって予想に使うんだろう？ 902 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 08:00:38 ID:w6BaQv4H.net] >>846 変数名が別のデータと被っていたので rm(list=ls())でメモリクリアしてから 作図し直した。 https://i.imgur.com/KpwhZ0G.png 903 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 08:24:34 ID:w6BaQv4H.net] >>847 K値の基準となる日を１週間ずつずらしてK値の変動をグラフにしたみた。 https://i.imgur.com/fkN9JLg.png 904 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 21:02:46 ID:iNIKliXH.net] K値モデルが正しく、なにもしなくて収束するなら、 アメリカやブラジルはとっくに収束してるわ K値モデルがダメな点は、ただの一学者の主張じゃなくて、 大阪府知事が信じて取り入れてしまった点だな やっぱ実効再生産数がいちばん大事 905 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 21:20:04.76 ID:iNIKliXH.net] なぜいまだに政府も東京都も大阪府も、 実効再生産数をモニタリング指標および政策目標につかわないのかまったく不明 東京と大阪以外は実効再生産数0.7、東京と大阪は実効再生産数0.8 これくらいを政策目標にすればいいのに 一時的に実効再生産数が0.5以下になるまで規制して感染者を大きく減らし、 その後は実効再生産数0.8目標でいいとおもう 906 名前：１３２人目の素数さん [2020/07/19(日) 21:56:54.03 ID:GNq9W41i.net] >実効再生産数 山中京大教授は実効再生産数のエクセル用公式を見つけて数値いれて威張って いますが、 実効再生産数は感染症数理モデルに基づきますのでマネすることはおやめください 907 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 00:09:27.09 ID:Ze2tCmmK.net] 疫学者が頑張って計算する正確な実効再生産数じゃなくても、 日々の確認陽性者数をエクセルに入れて処理する簡単な実効再生産数でも十分役に立つんだよ ネットで公表されてる簡易モデルは、実際の確認感染者数の増減とだいたい一致してる 908 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 00:13:48.19 ID:Ze2tCmmK.net] 日本の都道府県別実効再生産数 ttps://rt-live-japan.com/ ttps://covid.gutas.net/ アメリカの州ごとの実効再生産数 ttps://rt.live/ 世界の実効再生産数 ttps://covid.gutas.net/world 疫学者ががんばって計算する実効再生産数じゃなくても、 このレベルの簡易モデルで、十分通用するのが今回の新型コロナだよ 日本の場合、感染者がかなり少ない都道府県はまともに算出されないが、 それ以外はだいたい妥当な値になってる 909 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 00:30:32.49 ID:Ze2tCmmK.net] ドイツやイタリアは実効再生産数ベースに対策を考えてるのでは？ 実効再生産数が上がらないように対策してる感じがする 日本は、日本全国および東京都の実効再生産数が継続的に1を超えてるのに、 感染対策をやらずに自粛解除だけ進めていったのがこの事態を引き起こした 本来、自粛解除しても実効再生産数が１を超えないように感染対策をした上で 自粛解除をする必要があった 910 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 01:56:42.64 ID:L7jTb9sR.net] そう、なのに日本だとわけのわからんモニタリング指標に従って対策立ててる アホじゃないかと 911 名前：１３２人目の素数さん [2020/07/20(月) 08:24:49.04 ID:IvKiJrmj.net] 感染者数　成長曲線の3パターン ・指数関数型 マルサス人口論　原子核崩壊　 　官邸　アベノミクスカンフル剤 　コロナ対策より 　株価が指数関数型の上昇希望 ・ロジスティック曲線 　dX/dt=N(1-X/N)X SIモデル=SIRモデルの第1式.第2式 (隔離率γ=0)の感染者増加曲線。 PCR検査より感染者の隔離対策は 過去に実証された有効対策であり 　SIRモデルに基づく ・ゴンペルツ曲線 　高齢者の人口減少 　腫瘍ガン細胞の体積増加 　大阪府地方で流行っているK値 　=感染者数の新興占いの式。 　神奈川も流行り始めた模様 912 名前：１３２人目の素数さん [2020/07/20(月) 08:44:21 ID:IvKiJrmj.net] 東京都のコロナ対策はSIRモデルが 重要な役割をしている。 913 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 08:55:29 ID:mgzJmLOO.net] SIRって言ってもRは当分の間ほぼ0のままだから、事実上SIモデルだよな。 914 名前：１３２人目の素数さん [2020/07/20(月) 09:31:31 ID:IvKiJrmj.net] 積極的にPCR検査しない道県は SIモデル。 府はガン細胞増殖のK値モデル 915 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 10:24:10.81 ID:jF/Tht4x.net] γは回復率やろ 平均回復日数の逆数 916 名前：１３２人目の素数さん [2020/07/20(月) 10:47:58 ID:AYBv/aE8.net] γは隔離率とも回復率ともいう。 抗体ができる、死亡、隔離より感染しない者を表わす係数 917 名前：１３２人目の素数さん [2020/07/20(月) 11:45:14 ID:AYBv/aE8.net] ・SIRモデル 閉鎖区域の総人口Nが一定。 S＋ I＋R= N 一定。 感染可能者S→潜期→ウィルス排出 →発病I→一部回復(抗体).隔離.死亡 などR、 ウィルス排出→発病の期間0とする。 第1式　dX/dt=-βXY 第2式dY/dt=βXY-γY 第3式d Z/dt=γY. ・発展SIRSモデル 　一部R→Sがある場合 918 名前：１３２人目の素数さん mailto:sage [2020/07/23(木) 02:43:39 ID:iTEunCvL.net] 大阪モデル完全崩壊だな K値唱える人がいるのはいいが、それを政治が採用するとかやばすぎる 実効再生産数を見るのが最適 ただし、実効再生産数はクラスター対策が不可能になったら突然急騰するので、 クラスター対策できる上限には注意が必要 919 名前：１３２人目の素数さん [2020/07/23(木) 18:41:45.80 ID:bfbrr4ZK.net] 実効再生産数は数理モデルによる。 数学として SIモデルからロジスティック方程式 となり感染率βを小さくする対策 新しい日常不要不急不急の外出禁止 が必要だ。 更にオイラー差分を使い離散ロジスティック方程式が得られる。 シャルコフスキーの定理、リ.ヨーク の定理のカオスの意味で β Nよりカオスになる場合があるが、 感染症数理モデルとしては感染症 の基礎知識が必要だろう。 920 名前：１３２人目の素数さん [2020/07/24(金) 18:33:22.00 ID:g3DPzG/z.net] K値って、大阪だけじゃなく神奈川も使ってるのか 知事がアホだとこういったオカルトに走り出す 921 名前：１３２人目の素数さん mailto:sage [2020/07/24(金) 21:41:46.70 ID:E42oSbW+.net] 世間は無駄に大騒ぎだが以下の事実がある 感染者増と全く比例してない おまえら全く割出せなかっただろ？ つまり今後の予測もギブアップだよな 死者数 7月*1日　*1人 7月*2日　*1人 7月*3日　*1人 7月*4日　*0人 7月*5日　*0人 7月*6日　*1人 7月*7日　*2人 7月*8日　*1人 7月*9日　*1人 7月10日　*0人 7月11日　*0人 7月12日　*0人 7月13日　*0人 7月14日　*2人 7月15日　*1人 7月16日　*0人 7月17日　*0人 7月18日　*0人 7月19日　*0人 7月20日　*3人 7月21日　*1人 7月22日　*1人 922 名前：１３２人目の素数さん mailto:sage [2020/07/24(金) 22:44:07.67 ID:leZfEHQl.net] おいおまえら また新しいの出てきたぞ K値　改め　Rw値 wwwww 923 名前： [] [ここ壊れてます] 924 名前：１３２人目の素数さん mailto:sage [2020/07/24(金) 23:11:05 ID:kQJgVFx5.net] 医療提供体勢が整っているから死者数が少なくなるんじゃね？ 925 名前：１３２人目の素数さん [2020/07/24(金) 23:11:27 ID:2uKiwAiC.net] 死亡者は重症感染者数による。 第1波の重症者は回復か死亡より 激減したが、今後第2波が GO TOキャンペーンで拡大し 若年軽症から感染した重症者 (特に高齢者)が増加するだろう。 東洋経済の実効再生産数公式 西浦教授監修とのこと 7月23日 全国　1.4 ・東京1.20神奈川1.19埼玉1.19千葉1.1 ・大阪2.46兵庫2.23京都1.5 ・愛知9.51 岐阜3.88 ・北海道1.29 福岡 K値予想の大阪 新規感染者149人 陽性率13パーセント。 大阪と官邸はコロナはただの風邪路線「GO TOキャンペーン」 全国へコロナウィルスばらまき 926 名前：１３２人目の素数さん mailto:sage [2020/07/24(金) 23:13:52 ID:kQJgVFx5.net] 薬やその使い方も少しわかってきただろうし 血栓に注意が必要ということもわかってきただろう サイトカインストームにステロイドを使うとか 回復した人の血漿を治療に使うとか 対処法が変われば結果も変わるだろう 927 名前：１３２人目の素数さん [2020/07/24(金) 23:14:10 ID:2uKiwAiC.net] >>869 追加 福岡3.63 928 名前：１３２人目の素数さん [2020/07/25(土) 01:04:16.54 ID:3rBi8tDq.net] >>870 確かに。 対症療法とはいえ、それなりに対策が立てられるようになってんのかもね。 アメリカもUKも新規死者数の減少幅のほうが、感染者の減少幅より大きい （致死率が半分くらいに下がってる？） 929 名前：１３２人目の素数さん [2020/07/25(土) 22:57:57.28 ID:M+D5uJMO.net] 発展形も含め基本のSIRモデルが コロナ対策の基盤で ・第3式から隔離者増は医療システム内。 医療崩壊がおきるとコロナに限らず、 死亡者増加を招く。 野党は日本でも独並みのPCR検査7万件をと喚いているが、現実の医療システム 無視で机上の妄想。 与論島などベッド数が少ないし 医療崩壊危機がすぐおきる。 また個人の任意PCR検査は高額だから 政府か共産党、 まあ国会議員が寄付すれば良い 930 名前：１３２人目の素数さん [2020/07/25(土) 23:04:39.05 ID:M+D5uJMO.net] ・移流項つきSIRモデルは感染経路 対策。 ・GOTOキャンペーン推進の、官邸や 大阪関西愛知などは東京問題 東京由来などと非科学的な発言を している。 東京発新幹線は名古屋関西方面の 感染爆発地域も東北長野新幹線方面では 落ち着いていて実効再生産数もこれを 示している。 第1波後厚労省の集団抗体検査は 宮城0.03<東京0.1<大阪0.17パーセント。 大阪関西愛知などPCR検査も過少で 第1波が地域で燻っていても見逃して いただけだ ・海外からの感染者入国管理は 政府の責任だが成田空港も毎日感染者 が発見されている。 水際対策は行動アプリ義務もなく ズブズブの上更に緩めようとしている。 西浦教授の試算で毎日10人の感染者が 入国すれば3ヶ月で感染爆発だ。 ・野党国会議員は緊急事態宣言中に 新宿歌舞伎町で歓楽するバカ丸出し だが、横田基地の兵隊さんも歓楽し ている？ まず防衛省は米軍が軍人従業員の感染 を管理しているか情報を開示すべきだ。 野党は沖縄の米軍基地だけに興味が あるらしい、、

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef