cocolog:95369157

岩見＆中岡＆岩波『ウイルス感染の数理モデルとシミュレーション』を読んだ。感染症の社会シミュレーションに興味を持って手に取ったが、この本は細胞内のウィルス動態が中心でその点ではアテが外れた。が、参考にできそうな部分も多かった。生物学に関しては基本飛ばして読んだ。 (JRF 7075)

JRF 2025年4月11日 (金) Tweet

『ウイルス感染の数理モデルとシミュレーション - データを定量的に理解する 』(岩見 真吾 ＆ 中岡 慎治 ＆ 岩波 翔也 著, 共立出版, 2024年2月)
https://www.amazon.co.jp/dp/4320115538
https://7net.omni7.jp/detail/1107479533

JRF2025/4/111410

感染症の社会シミュレーションをしたいと計画している。人口の様子から感染を実際に起こさせるエージェント・ベース・シミュレーションをしたいと思っている。そのために以前↓(参: [cocolog:95052633](2024年9月))…みたいな準備はしていた。

《JRF-2018/simple_synthetic_population: 仮想合成人口個票の簡易な実装》
https://github.com/JRF-2018/simple_synthetic_population

JRF2025/4/118366

また近時、大久保 潤『確率的シミュレーション』という本も読んだ([cocolog:95358180](2025年4月))。

今回、この本は生物学寄りの専門書で私には難しすぎるとは思ったのだが、いちおう「目を通す」だけはしてみた。社会シミュレーションとしての側面はほぼなく、細胞内のウィルス動態の常微分連立方程式(たまに偏微分もある)による数理モデルが議論の中心になる。その数学的側面は興味深く読んだが、細胞内の生物学については正直ほとんど飛ばして読むことになった。

JRF2025/4/113940

コロナ禍の最中に書かれた本で、コロナ(COVID-19)にも最後の2章が割かれている。そこは特に興味深く読んだ。ただ、基本は、HIV (AIDS) や C型肝炎(HCV) などが著者らのメインの研究領域のようだ。

ここからは、引用しながら軽くコメントしていく。

JRF2025/4/117160

……。

＞AIC = -2 log f(y|θ*) + 2k＜(p.30)

Kullback-Leibler 情報量については、上の simple_synthetic_population でもちょっと参考にしたのでわかるのだが、AIC (赤池情報量基準) については、毎度いろいろな文献に載ってそのたび確認するのだが、どうも忘れてしまう。

JRF2025/4/114111

今回のこの式は最初どう読めばいいのかピンとこなかった。そこで Grok さんに聞いたところ、log f(y|θ*) は、尤度関数だということで、やっと理解が追いついた。この本では、この直前のところに別の意味の f() が出ていて、それに引っぱられて、意味不明になっていたのだった。ここの文脈で言えば、f() は L(y|θ*) と書いてくれたほうがわかりやすかっただろう。

JRF2025/4/115776

……。

＞FACS 解析＜(p.43)

知らない言葉。ググると…、

「FACS（フローサイトメトリーによる細胞の分選(または、蛍光活性化セルソーティング)）解析は、細胞を1個ずつ測定して、細胞の特性を解析する技術です。細胞の大きさや数、周期、表面マーカー量などを測定し、細胞を特定の集団に選別することができます。」…とのこと。

JRF2025/4/110783

……。

＞ギレスピーアルゴリズム＜(p.51)

直近で大久保 潤『確率的シミュレーション』という本も読んだ([cocolog:95358180](2025年4月))ところ。先に読んどいて良かった。

JRF2025/4/115054

……。

＞ここで分布遅れをもつ微分方程式(4.6)から常微分方程式(4.7)を導出した方法は、"linear-chain-trick" と呼ばれており、数理モデルを用いたさまざまな数理生物学の研究で広く使われている[48]。＜(p.84)

まず偏微分(連立)方程式(4.4)を作り、「特性曲線法」を用いて分布遅れを持つ(常)微分(連立)方程式(4.6)に変形する。そこにガンマ分布を仮定し、それを整数的なアーラン分布に置き換えることで、過程を分解したような常微分(連立)方程式(4.7)に変形する。これが見事だと思った。

JRF2025/4/112690

偏微分方程式というとシュレジンガー方程式みたいな x と t などのものを思い出すが、ここでは、時間遅れが(時間とは別の)もう一つの変数となる偏微分方程式で、そういうのアリなんだ…という発見(思い出し？)があった。

JRF2025/4/115476

……。

＞反応拡散方程式を数値的に解く場合は、偏微分方程式の数値解放に関する知見が必要である。境界条件が有界な矩形や円領域など比較的単純な場合、領域の分割に支障がないため、有限差分法など基本的な数値解法を利用すればよい。一方、境界の形状が複雑な場合は、一般に有限要素法(finite element method)が用いられる。

JRF2025/4/116784

有限要素法では、偏微分方程式を積分形式で定式化された弱形式で表現し、領域を有限要素(2次元の場合は三角形など)に細分化して近似的に離散差分方程式を構成し、領域全体で重ね合わせた離散差分方程式を解くことで、近似解を得る。

近年、有限要素法に基づく数値計算をするうえで、FreeFEM++ といった非常に使いやすいソフトウェアがフリーで利用できる(《コード 5-2: 拡散方程式の数値計算の実装》参照)。
＜(p.134)

JRF2025/4/110625

有限要素法。いつか学ばないといけないと思いつつ、ずっと学ばないで来ている。もう縁はないのかな…。

JRF2025/4/113927

……。

社会シミュレーションでなくウィルス動態的シミュレーションが社会シミュレーション的に役に立つ場面があるようだ。

＞利用できるデータは発症後のウィルス量のみに限定されているが、数理モデルを用いれば、過去の感染時刻を推算することができる。

JRF2025/4/118272

(…)

客観的な指標であるバイオマーカーに基づく提案アプローチでは、"暴露時刻" という不確実性を回避できる。暴露時刻は感染者の想起に依存した主観的なもので、接触者追跡におけるインタビューに内在する1つの問題である。もちろん、COVID-19 に限らず、特に感染イベントを直接観察することが実用的でない感染症の場合、ウィルス量を用いた潜伏期間の推定は非常に有用であると考えられる。
＜(p.246)

JRF2025/4/116160

コロナなどは体内のウィルス量の推移がつりがね状でモデル化でき、例えば二点、調べれば、いつ感染したかが予測できる。それを利用して、例えば、感染したのが帰ってくる前の海外であったかどうかなどを判定できる。

JRF2025/4/116384

また、だいたいの感染パターンが予測できるため、隔離終了を判定する固定期間法(発症後10日間隔離)と検査法(PCR検査で感染がないことを確認して隔離終了とする方法)をすり抜けるだいたいの割合などが予測でき、それに基づいて、隔離期間や検査回数を提案できたりするようだ。そういうことができるのかという驚きがあった。

JRF2025/4/119174