著者：畝見達夫，作成：令和2年9月1日，編集：12月6日
この文書は個体ベース感染シミュレータ SimEpidemic に使われているモデルについて記述したものである。
シミュレーションの場となる世界は，連続値直行座標で位置を表現する正方形の2次元平面である。 初期状態では人口に対応する数の個体がランダムに配置される。 初期状態における個体の健康状態は，ランダムに選ばれた数人の未発症感染者を除き， 未感染であるものとする。 辺の長さと人口を調整することで人口密度を調整する。
シミュレーション過程は，まず，現在の状態での各個体の相互作用を計算し， その結果と個体自身の状態に応じて次の状態を算出するというステップを繰り返すことで実現する。
各個体は各ステップにおいて，周囲から受ける力ベクトルに従い移動する。 力には，周囲の個体との距離の２乗に反比例する斥力と周囲のもっとも好む個体への引力を合成して力のベクトルを計算し， 質量と摩擦係数を想定した積分の近似計算1により次の位置を算出する。 世界の境界に隣接する場合は境界線からの斥力も受ける。 個体 \(i\) が受ける力 \(F_i\) は， \[F_i = \alpha\cdot\sum_j\frac{p_i - p_j}{|p_i - p_j|^3}+\beta\cdot(p_b-p_i) +\gamma\cdot\sum_{j=1}^4\frac{p_i-w_j}{|p_i-w_j|^3}+\zeta\cdot G_i\] \(p_i\) は個体 \(i\) の位置，\(b\) は周囲で最も好みの個体，\(w_j\) は境界線 \(j\) 上の最も近い点，\(\alpha,\:\beta,\:\gamma,\:\zeta\) は係数である2。 \(G_i\) は後に述べる集会の中心へ向かうベクトルである。
発症していない個体は設定された確率で世界の領域内の別の位置に移動する。 対策としての移動制限は，この確率と距離を調整することで実現する。
学校，職場，イベント，集会，さらに家庭など人が集まる現象を再現するモデルとして， 集会 を確率的に発生される。個々の集会は円形の領域で．大きさ，持続時間，および，強さが割り当てらる。 集会の円周付近の発症していない個体は，集会の中心へ向かう力を受ける。力の大きさは円周上で最も強く， そこからの距離に反比例して小さくなる。 右図に個体が受ける力の分布を示す。 上の式における \(G_i\) は以下のように定義する。
\[G_i=\left\{\begin{array}{ll} 0 & \mbox{if } |p_i - G_p|-G_r \notin [-C_r,C_r] \\ \displaystyle 1 - k & \mbox{else if }|p_i - G_p| > G_r \\ \displaystyle 1 + k & \mbox{else if }G_r > C_r \\ \displaystyle |p_i - G_p|\:/\:G_r& \mbox{otherwise}\\ \end{array}\right.\] ただし \(G_p\) はその個体から最も近い集会の中心，\(G_r\) はその半径，\(C_r\) は円周から影響する領域の境界までの距離を表す定数である。 \(k\) は \((|p_i - G_p|-G_r)\:/\:C_r\) である。
さらに，集会に留まるよう中心から \(G_r\) 以内の領域では摩擦係数が増すものとし，中心で通常の倍の摩擦になるよう距離に比例した補正を適用する。
各個体の健康状態を 1. 未感染，2. 未発症，3. 発症，4. 快復，5. 死亡の５種類に分類する。 未感染の個体と未発症または発症の個体がある距離以内に近づくと，設定された確率で感染する。 感染した個体は個体に割り振られた速度で症状が進み，ある閾値を越えると発症し， さらに次の閾値を越えると死亡する。 感染後，個体に割り振られた快復開始日数に達した個体は，症状が一定速度で快方に向かい， 感染前の健康状態に戻った時点で完治したとみなす。 快復した個体は免疫を獲得したものとみなされるが， 設定された期間を過ぎると免疫を失ったとみなし，未感染の状態にもどる。
感染者から他個体に感染する感染性の強さは，感染から完成性の遅れ期間を経るまでは 0 で， その後ピーク到達期間と潜伏期間の短い方までは，時間経過に比例して強まる。
社会的距離，移動の制限，隔離，接触者追跡を次節に述べる検査とともに設定する。
行動モデルの節で述べた他の個体との衝突回避の斥力の係数を変化させることで社会的距離を表現する。 パラメータとして，斥力の係数の増加率と要請への協力率を設定する。
行動モデルの節で述べた他の場所への移動の距離と頻度を調整することで，移動制限を表現する。
次節で述べる検査の結果が陽性であった個体を隔離する。 隔離された個体は世界とは別に設けられた病院領域へ移動する。 隔離された個体が快復すると元の位置へ戻る。
接触者追跡は，各ステップで感染可能な距離以内に近づいた別の個体を接触者として， 個体が持つリストに追記することで実現する。既にリストにある個体と再度接近した場合は， 過去の記録を消去し，新しい記録で書き換える。ただし，接触から14日を過ぎた記録は消去する。 接近したすべての個体を記録するのではなく，設定された確率により記録するかどうかが確率的に決まる。
検査の対象となる個体の選択方法と確率，検査の感度と特異度を設定する。 対象とする個体は，症状の出た個体，疑いのある個体，陽性の検査結果が判明した個体との接触者 の中から選ぶ。接触者以外の対象は設定された確率により選択する。 接触者の検査確率は，接触者の記録の確率に吸収できるので，別にパラメータを設けない。
潜伏期間や移動距離など分布の形式をとるパラメータは， 最小値，最大値，最頻値の３つの数値の組み合わせにより表現し， その条件を満足するベル型の確率密度関数に従った擬似乱数発生により個々の場合の値を決める。3
© Tatsuo Unemi, 2020. All rights reserved.