基本再生産数

いくつかの感染症に対する集団免疫のしきい値(縦軸)と基本生産数(横軸)の関係

基本再生産数(きほんさいせいさんすう、英語: basic reproduction numberR0と表記され、R noughtあるいはR zeroと読まれる[1])とは、疫学において、感染症に感染した1人の感染者が、誰も免疫を持たない集団に加わったとき、直接感染させる人数の期待値である[2]。1人の患者が平均して何人に感染を広げる可能性があるかを示す[3][4]

また、人口学においては女性が生涯に生む期待女児数(純再生産率)を指す[5]。女性が1人だけ娘を産むときには人口は単純再生産であって、1人より多く産めば拡大再生産となる[5]。すなわち母親世代とその娘世代の総数比を指し、R0 が1より大きければ人口は拡大再生産されるが、1より小さければ縮小再生産される[6][7]

基本再生産数は、人口学で生まれた概念であるが、感染症疫学でも基本概念である[6]R0 は、人口学では人口増加の閾値として用いられ、進化生物学においては侵入生物ないし突然変異体適応度と解釈され、感染症疫学では感染症の侵入条件や臨界免疫化割合を定量化するものと見なされる[8]

疫学

[編集]
は、1人の人間から感染された人数の平均値である。たとえば、エボラのは2であるため、平均するとエボラに感染した1人の人間は他の2人の人間にエボラを感染させる。

定義

[編集]

すべての個体が感染症に対する感受性英語版[注釈 1]をもつ集団内で、一感染個体により直接生み出される感染個体数の平均と考えられる[9][要非一次資料]。この定義は他の個体はすべて感染しておらず、(先天性の、あるいはワクチン接種による)免疫をもっていない状況を記述している。オーストラリア保健省英語版による定義などでは「感染症の伝播に対する計画的な介入」のないことが付け加えられてる[10]

定義により基本再生産数 R0ワクチン接種により変更することはできない。また基本再生産数 R0無次元数であり、倍加時間[11] のような時間の単位[12] をもつ割合ではないことに注意が必要である。

基本再生産数 R0 は、環境因子や感染集団の行動による影響も受けるため、病原体に対する生物学的な定数ではない。さらに、基本再生産数 R0 の値は通常、数理モデルから推定されるので、推定値は使用されたモデルや他のパラメータの値に依存する。したがって、文献における値は特定の文脈においてのみ意味があり、古い値を使用したり、異なるモデルに基づく値を比較したりするべきではない[13]。また基本再生産数 R0 自体は集団内における感染症の蔓延する速度を推定するものではない。

感染症モデル

[編集]

基本再生産数 R0 の最も重要な用途は新興感染症が集団内に蔓延するかどうかを決定することと、感染症を撲滅するためには集団のどのくらいの割合にワクチン接種をして免疫化すべきなのかを決定することである。一般に使用される感染症モデル英語版では R0 > 1 のとき、感染症は集団内で蔓延をはじめ、R0 < 1 のときには蔓延しない。一般的に基本再生産数 R0 の値が大きいほど、流行を抑制するのは難しい。単純なモデルでは、感染症の持続的な蔓延を防ぐためには、実質的に免疫化される(=感受性のない)べき集団の割合は集団免疫閾値

1 − 1/R0

を超えなくてはならない[14];一方でエンデミック定常状態における感受性個体数の割合は 1/R0 である。

基本再生産数は、感染個体の感染力の持続時間・微生物の感染性・感染個体が接触している集団における感受性個体数を含むいくつもの要因による影響を受ける。

疫学上、感染症流行を予測・抑制することは公衆衛生上で重要な課題であり、その問題に対して、個体群生態学で使用された数理モデル安定性分析などが感染症疫学でも流用できたことから、この研究分野が発展する事となった[15]

感染個体が単位時間あたり平均 β の感染を生み出す接触をし、感染性期間が平均 τ であると仮定する。このとき基本再生産数は

R0 = βτ

で与えられる。この単純な式は R0 を減らし、最終的には感染の伝播を減らすいくつかの方法を示唆する。単位時間あたりの接触を減らす(たとえば、伝播が他者との接触を必要とする場合には、家に留まる)、または(防護具などによって)感染を生み出す接触の割合を減らすことで、単位時間あたりの感染を生み出す接触の平均 β を減らすことができる。また感染個体をできるだけ早く発見し、隔離/治療/殺処分(動物の場合によくある)することで、感染性期間の平均 τ を減らすことができる。

区画モデル

[編集]

疫学における区画モデルの代表例としてSIRモデルがある。SIRモデルは感受性 (Susceptible)、感染した (Infected)、隔離された(Removed: 回復または死亡)人の個体数の経時的な疾患動態を説明する。なお、SIRモデルの R(0)と基本再生産数の R0 を混同しないように注意。前者は t = 0 での隔離者数を表わす。

感染症の基本再生産数値

[編集]
介入前の有名な感染症のR0の値と集団免疫のしきい値
感染症 伝染 R0 集団免疫のしきい値[注釈 2]
麻疹(はしか) エアロゾル 12–18[16][17] 92–94%
水痘(水ぼうそう) エアロゾル 10–12[18] 90–92%
流行性耳下腺炎(おたふく風邪) 呼吸器飛沫英語版 10–12[19] 90–92%
COVID-19デルタ株 呼吸器飛沫とエアロゾル (5-9.5)[20] 80–89%
風疹 呼吸器飛沫 6–7[注釈 3] 83–86%
ポリオ 糞口経路 5–7[注釈 3] 80–86%
百日咳 呼吸器飛沫 5.5[25] 82%
天然痘 呼吸器飛沫 3.5–6.0[26] 71–83%
COVID-19アルファ株 呼吸器飛沫とエアロゾル 4–5[27] 75–80%
HIV/AIDS 体液 2–5[28] 50–80%
COVID-19野生型 呼吸器飛沫とエアロゾル[29] 2.9 (2.43.4)[30] 65% (5870%)
SARS 呼吸器飛沫 2–4[31] 50–75%
ジフテリア 唾液 2.6 (1.74.3)[32] 62% (4177%)
インフルエンザスペインかぜ 呼吸器飛沫 2–3[33] 50–67%
風邪 呼吸器飛沫 2–3[34] 50–67%
エボラ2014年のエボラ出血熱の流行 体液 1.8 (1.41.8)[35] 44% (3144%)
インフルエンザ(2009年のパンデミック株 呼吸器飛沫 1.6 (1.32.0)[4] 37% (2551%)
インフルエンザ (季節株) 呼吸器飛沫 1.3 (1.21.4)[36] 23% (1729%)
ニパウイルス 体液 0.48[37] 0%[注釈 4]
MERS 呼吸器飛沫 0.47 (0.290.80)[38] 0%[注釈 4]

限界

[編集]

数理モデル、特に微分方程式から計算をするとき、R0 であると言われる数の多くは、実際には単なる閾値であって、2次感染者数の平均ではない。そのような閾値を数理モデルから導出する多くの手法があるが、常に R0 の真値を与えるものはほとんどない。これはマラリアのように宿主間に媒介動物がいる場合は特に問題が多い[39]

これらの閾値は疾病が絶滅(R0 < 1)するか流行(R0 > 1)するかを決定するが、一般にそれらを異なる疾病で比較することはできない。したがって、上の表にある値は慎重に用いる必要がある。(特に数理モデルから計算された場合。)

手法には生存関数ヤコビ行列の最大固有値の再配置・next-generation method[40]内的自然増加率からの計算[41]エンデミック定常状態の感受性人口からの計算・平均感染年齢からの計算[42]・最終規模方程式からの計算などがある。多くの手法は、たとえ同じ微分方程式から出発したとしても、互いに一致しない[要出典]。実際に2次感染者の平均を計算しているものはさらに少ない。基本再生産数 R0 が実地で観測されることは滅多になく、通常は数理モデルから計算され、これにより有用性は著しく制限される[43]

実効再生産数

[編集]

基本再生産数を、実効再生産数(effective reproduction number)R と混同してはならない。実効再生産数とは、感染個体がすでに存在するかもしれない現在の集団内で、一感染個体により直接生み出される感染個体数の平均である。

広東省疾病予防コントロールセンターによれば、「伝染性を表現するために実効再生産数 (R) がより一般的に使用される。実効再生産数は感染症例ごとによって発生する二次症例数の平均として定義される。管理措置がない場合、R = R0 χ (χは感受人口の割合) となる。」[44][要非一次資料]。例えば2019-nCoVの実効再生産数は2.9で、SARSの実効再生産数は1.7となる[44]

人口学

[編集]

人口学において、基本再生産数 R0 は、1人の女性が生まれてから a 歳まで生き延びる確率(生残率)l(a) と、年齢ごとの年齢別女児出生率 β(a) の積を総和した数字

である。約2.08倍すると男女込みの平均出生児数または、少子化の資料に頻出する合計特殊出生率TFR)となる[6][45]

基本再生産数は、母親世代とその娘世代の総数比を指し、R0 が1より大きければ人口は拡大再生産されるが、1より小さければ縮小再生産される[6]。長期的にみれば、R0 > 1であれば人口は増加し、R0 < 1 であれば人口減少がおきる[6]

日本人口の2005年の R0 は0.61で、これは母親世代の人口の6割程の数の娘しか生まれてこないことを意味する[6]

学説史

[編集]

人口のレベルで感染症がどのように広がっていくのか、また制御するための介入についての数理モデルを使った研究は、18世紀のスイスの数学者ダニエル・ベルヌーイ天然痘死亡率が人間の寿命に与える影響研究に遡り、これは区画モデルの端緒だった[5]。天然痘による死亡リスクというものがなくなった場合に寿命はどう延びるか、つまり、いくつかの競合する死因がある場合に、そのうちの一つが取り除かれたら寿命がどの程度延びるか、という競合リスクモデルの研究の端緒となった[5]

ダニエル・ベルヌーイ1766年フランス王立科学アカデミーで天然痘罹患率と死亡率データ分析「小虫を死因とした死亡例の新解析[疑問点](Essai d’une nouvelle analyse de la mortalite causee par la petite verole)」を発表した[46]。これは予防接種の効果計測や、観察の打ち切り(センサリング censoring)[47][48][49] などの知見を含めた最初期の統計学上の解析となった。

脚注

[編集]

注釈

[編集]
  1. ^ 病原体に対する免疫のないこと[2]
  2. ^ Calculated using p = 1 - 1 / R0.
  3. ^ a b From a module of a training course[21] with data modified from other sources.[22][23][24]
  4. ^ a b When R0 < 1.0, the disease naturally disappears.

出典

[編集]
  1. ^ Milligan, Gregg N.; Barrett, Alan D. T. (2015). Vaccinology : an essential guide. Chichester, West Sussex: Wiley Blackwell. pp. 310. ISBN 978-1-118-63652-7. OCLC 881386962 
  2. ^ a b 用語集 2018.
  3. ^ 【図解】R0(基本再生産数)とは? 感染症疫学用語の基礎知識 AFPBB,2020年3月31日 13:33
  4. ^ a b “Pandemic potential of a strain of influenza A (H1N1): early findings”. Science 324 (5934): 1557–61. (June 2009). Bibcode2009Sci...324.1557F. doi:10.1126/science.1176062. PMC 3735127. PMID 19433588. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3735127/. 
  5. ^ a b c d 稲葉寿「感染症の数理東京大学大学院数理科学研究科平成 20 年度年第 7 回例会「医療とアクチュアリー」講演 2009年3月18日。
  6. ^ a b c d e f 稲葉寿. “基本再生産数”. 理学のキーワード. 東京大学大学院理学系研究科・理学部. 2020年4月15日閲覧。
  7. ^ de Boer; Rob J. Theoretical Biology. http://theory.bio.uu.nl/rdb/books/tb.pdf 2007年11月13日閲覧。 
  8. ^ 稲葉寿「基本再生産数 R0 の数学」日本数理生物学会ニュースレターJSMB Newsletter No. 64, pp. 1–7, 2011
  9. ^ Christophe Fraser; Christl A. Donnelly; Simon Cauchemez (19 June 2009). “Pandemic Potential of a Strain of Influenza A (H1N1): Early Findings”. Science 324 (5934): 1557–1561. Bibcode2009Sci...324.1557F. doi:10.1126/science.1176062. PMC 3735127. PMID 19433588. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3735127/. Free text
  10. ^ The reproduction number, “Using Mathematical Models to Assess Responses to an Outbreak of an Emerged Viral Respiratory Disease”, Australian Goverment Department of Health, (2006), ISBN 1-74186-357-0, https://www1.health.gov.au/internet/publications/publishing.nsf/Content/mathematical-models~mathematical-models-models.htm~mathematical-models-2.2.htm 
  11. ^ Siegel, Ethan. “Why 'Exponential Growth' Is So Scary For The COVID-19 Coronavirus” (英語). Forbes. 2020年3月19日閲覧。
  12. ^ Jones, James. “Notes On R0”. Stanford University. 2020年4月13日閲覧。
  13. ^ Delamater, Paul L.; Street, Erica J.; Leslie, Timothy F.; Yang, Y. Tony; Jacobsen, Kathryn H. (January 2019). “Complexity of the Basic Reproduction Number (R0 )”. Emerging Infectious Diseases 25 (1): 1–4. doi:10.3201/eid2501.171901. ISSN 1080-6040. PMC 6302597. PMID 30560777. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6302597/. 
  14. ^ Fine, Paul; Eames, Ken; Heymann, David L. (2011-04-01). “"Herd Immunity": A Rough Guide” (英語). Clinical Infectious Diseases 52 (7): 911–916. doi:10.1093/cid/cir007. ISSN 1058-4838. PMID 21427399. https://academic.oup.com/cid/article/52/7/911/299077. 
  15. ^ 西浦 & 稲葉 2006.
  16. ^ “The basic reproduction number (R0) of measles: a systematic review”. The Lancet. Infectious Diseases 17 (12): e420-e428. (December 2017). doi:10.1016/S1473-3099(17)30307-9. PMID 28757186. 
  17. ^ “Complexity of the Basic Reproduction Number (R0)”. Emerging Infectious Diseases 25 (1): 1–4. (January 2019). doi:10.3201/eid2501.171901. PMC 6302597. PMID 30560777. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6302597/. 
  18. ^ Ireland's Health Services. Health Care Worker Information. https://www.hse.ie/eng/health/immunisation/hcpinfo/guidelines/chapter23.pdf 2020年3月27日閲覧。 
  19. ^ Australian government Department of Health Mumps Laboratory Case Definition (LCD)
  20. ^ McMorrow, Meredith (July 29, 2021). “Improving communications around vaccine breakthrough and vaccine effectiveness”. p. 15. July 30, 2021閲覧。
  21. ^ Centers for Disease Control and Prevention; World Health Organization (2001). "History and epidemiology of global smallpox eradication". Smallpox: disease, prevention, and intervention (training course) (Presentation). Atlanta: Centers for Disease Control and Prevention (published 25 August 2014). cdc:27929. 2017年3月17日時点のオリジナルよりアーカイブ (PDF)。2021年6月17日閲覧
  22. ^ “Herd Immunity: History, Theory, Practice”. Epidemiologic Reviews 15 (2): 265–302. (1993-07-01). doi:10.1093/oxfordjournals.epirev.a036121. ISSN 1478-6729. PMID 8174658. https://academic.oup.com/epirev/article-abstract/15/2/265/440430. 
  23. ^ “National, state, and urban-area vaccination-coverage levels among children aged 19–35 months, United States, 1999”. American Journal of Preventive Medicine 20 (4): 88–153. (2001-05-01). doi:10.1016/S0749-3797(01)00274-4. PMID 12174806. https://www.ajpmonline.org/article/S0749-3797(01)00274-4/pdf. 
  24. ^ “Vaccination coverage among children enrolled in Head Start programs or day care facilities or entering school”. MMWR. CDC surveillance summaries: Morbidity and mortality weekly report. CDC surveillance summaries 49 (9): 27–38. (2000-09-22). ISSN 1545-8636. PMID 11016876. https://www.cdc.gov/mmwr/preview/mmwrhtml/ss4909a2.htm. 
  25. ^ “Incidence and reproduction numbers of pertussis: estimates from serological and social contact data in five European countries”. PLOS Medicine 7 (6): e1000291. (June 2010). doi:10.1371/journal.pmed.1000291. PMC 2889930. PMID 20585374. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2889930/. 
  26. ^ “Transmission potential of smallpox in contemporary populations”. Nature 414 (6865): 748–51. (December 2001). Bibcode2001Natur.414..748G. doi:10.1038/414748a. PMID 11742399. https://www.nature.com/articles/414748a 18 March 2020閲覧。. 
  27. ^ Gallagher, James (12 June 2021). “Covid: Is there a limit to how much worse variants can get?”. BBC News. https://www.bbc.com/news/health-57431420 21 July 2021閲覧。 
  28. ^ Playing the Numbers Game: R0”. National Emerging Special Pathogen Training and Education Center. 12 May 202時点のオリジナルよりアーカイブ。27 December 2020閲覧。 “[...] while infections that require sexual contact like HIV have a lower R0 (2-5).”
  29. ^ “Airborne transmission of SARS-CoV-2”. Science 370 (6514): 303–304. (October 2020). Bibcode2020Sci...370..303P. doi:10.1126/science.abf0521. PMID 33020250. オリジナルの5 October 2020時点におけるアーカイブ。. https://science.sciencemag.org/content/370/6514/303.2 30 October 2020閲覧。. 
  30. ^ “Reproductive number of coronavirus: A systematic review and meta-analysis based on global level evidence”. PLOS ONE 15 (11): e0242128. (2020-11-11). Bibcode2020PLoSO..1542128B. doi:10.1371/journal.pone.0242128. PMC 7657547. PMID 33175914. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7657547/. 
  31. ^ Consensus document on the epidemiology of severe acute respiratory syndrome (SARS). Department of Communicable Disease Surveillance and Response (Technical report). World Health Organization. p. 26. hdl:10665/70863. WHO/CDS/CSR/GAR/2003.11. A number of researchers have estimated the basic reproduction number by fitting models to the initial growth of epidemics in a number of countries. Their observations indicate that the SARS-CoV is less transmissible than initially thought with estimates of Ro in the range of 2-4.
  32. ^ “Clinical and Epidemiological Aspects of Diphtheria: A Systematic Review and Pooled Analysis”. Clinical Infectious Diseases 71 (1): 89–97. (June 2020). doi:10.1093/cid/ciz808. PMC 7312233. PMID 31425581. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7312233/. 
  33. ^ Mills CE, Robins JM, Lipsitch M (2004). “Transmissibility of 1918 pandemic influenza”. Nature 432 (7019): 904–6. doi:10.1038/nature03063. PMID 15602562. http://www.hsph.harvard.edu/disasters/articles/LoreeLipsitch.pdf. 
  34. ^ Magic formula that will determine whether Ebola is beaten”. The Telegraph. Telegraph.Co.Uk. 30 March 2020閲覧。
  35. ^ “A systematic review of early modelling studies of Ebola virus disease in West Africa”. Epidemiology and Infection 145 (6): 1069–1094. (April 2017). doi:10.1017/S0950268817000164. PMID 28166851. "The median of the R0 mean estimate for the ongoing epidemic (overall) is 1.78 (interquartile range: 1.44, 1.80)" 
  36. ^ “Seasonal influenza in the United States, France, and Australia: transmission and prospects for control”. Epidemiology and Infection (Cambridge University Press) 136 (6): 852–64. (June 2008). doi:10.1017/S0950268807009144. PMC 2680121. PMID 17634159. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2680121/. "The reproduction number across influenza seasons and countries lied in the range 0.9–2.0 with an overall mean of 1.3, and 95% confidence interval (CI) 1.2–1.4." 
  37. ^ “The pandemic potential of Nipah virus”. Antiviral Research 100 (1): 38–43. (October 2013). doi:10.1016/j.antiviral.2013.07.011. PMID 23911335. 
  38. ^ “The role of superspreading in Middle East respiratory syndrome coronavirus (MERS-CoV) transmission”. Euro Surveillance 20 (25): 14–8. (June 2015). doi:10.2807/1560-7917.ES2015.20.25.21167. PMID 26132768. 
  39. ^ “The Failure of R0. Computational and Mathematical Methods in Medicine 2011 (527610): 1–17. (2011). doi:10.1155/2011/527610. PMC 3157160. PMID 21860658. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3157160/. 
  40. ^ “The Basic Reproduction Ratio”. Mathematical Epidemiology of Infectious Diseases : Model Building, Analysis and Interpretation. New York: Wiley. (2000). pp. 73–98. ISBN 0-471-49241-8. https://books.google.com/books?id=5VjSaAf35pMC&pg=PA73 
  41. ^ “The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda”. Journal of Theoretical Biology 229 (1): 119–126. (2004). arXiv:q-bio/0503006. doi:10.1016/j.jtbi.2004.03.006. PMID 15178190. 
  42. ^ Ajelli M; Iannelli M; Manfredi P; Ciofi degli Atti, ML (2008). “Basic mathematical models for the temporal dynamics of HAV in medium-endemicity Italian areas”. Vaccine 26 (13): 1697–1707. doi:10.1016/j.vaccine.2007.12.058. PMID 18314231. 
  43. ^ “Perspectives on the Basic Reproductive Ratio”. Journal of the Royal Society Interface 2 (4): 281–93. (2005). doi:10.1098/rsif.2005.0042. PMC 1578275. PMID 16849186. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1578275/. 
  44. ^ a b Liu, Tao; Hu, Jianxiong; Kang, Min; Lin, Lifeng; Zhong, Haojie; Xiao, Jianpeng; He, Guanhao; Song, Tie et al. (2020-01-25). “Transmission dynamics of 2019 novel coronavirus (2019-nCoV)” (英語). bioRxiv: 2020.01.25.919787. doi:10.1101/2020.01.25.919787. https://www.biorxiv.org/content/10.1101/2020.01.25.919787v1. 
  45. ^ 稲葉寿「基本再生産数理論の最近の進歩 (第10回生物数学の理論とその応用)」『数理解析研究所講究録』第1917巻、京都大学数理解析研究所、2014年9月、143-152頁、CRID 1050001335843791104hdl:2433/223338ISSN 1880-2818 
  46. ^ Mem Math Phy Acad Roy Sci Paris(Mémoires de mathématique et de physique, presentés à l'Académie royale des sciences, par divers sçavans & lûs dans ses assemblées.)Blower, S; Bernoulli, D (2004). “An attempt at a new analysis of the mortality caused by smallpox and of the advantages of inoculation to prevent it” (PDF). Reviews in medical virology 14 (5): 275–88. doi:10.1002/rmv.443. PMID 15334536. オリジナルの27 September 2007時点におけるアーカイブ。. https://web.archive.org/web/20070927032605/http://www.semel.ucla.edu/biomedicalmodeling/pdf/Bernoulli%26Blower.pdf. 
  47. ^ 宮田敏 医学統計勉強会
  48. ^ 西川正子「生存時間分析とは
  49. ^ 折笠秀樹「講座 リハビリテーション医学研究で使われる多変量解析・3 ロジスティック回帰とCox回帰」『総合リハビリテーション』第37巻第7号、株式会社医学書院、2009年7月、631-641頁、CRID 1390565159911685504doi:10.11477/mf.1552101548ISSN 03869822 

参考文献

[編集]

関連項目

[編集]