Significant Excess Mortality probably due to COVID-19 in Tokyo, Japan until March, 2021 ======================================================================================= * Junko Kurita * Tamie Sugawara * Yasushi Ohkusa ## Abstract **Background** No remarkable mortality attributable to COVID-19 confirmed by PCR test has been observed in Japan. **Object** We sought to quantify excess mortality using the National Institute of Infectious Diseases (NIID) model. **Method** We applied the NIID model to deaths of all causes from 1987 up through May, 2021 for the whole of Japan and up through March, 2021 for Tokyo. **Results** Results in Japan show very few excess mortality in August and October, 2020 and May, 2021. It was estimated as 12, 104 and 260, respectively. Conversely, in Tokyo, 1323 excess mortality was detected **Discussion and Conclusion** We detected substantial excess mortality in Tokyo but a few in Japan. It might be important to continue to monitor excess mortality of COVID-19 carefully hereafter. Keywords * excess mortality * COVID-19 * all cause death * stochastic frontier estimation * NIID model * Tokyo * Japan ## 1. Introduction To date, excess mortality has mainly been used to assess the social effects of influenza activity [1–6]. However, since the emergence of COVID-19, excess mortality attributable to COVID-19 has been attracting attention [7] as a measure of the total effects of the disease because it can reflect cases which have not been identified as polymerase chain reaction (PCR) positive. Especially in Japan, PCR tests administered per capita have been few. Therefore, concern has arisen about the possibility that some deaths caused by COVID-19 have not been recognized heretofore. Moreover, excess mortality related to COVID-19 might be expected to contribute to evaluation of vaccine effects. For these evaluations, the estimated excess mortality without the effects of a vaccine should be regarded as a baseline. Nevertheless, no such a trial has been undertaken to date. This study might be the first trial to measure that figure in Japan. As of the end of June, 2021, the COVID-19 outbreak showed about 800 thousand patients and about 15 thousand deaths from the outbreak have been reported in Japan. Although Japan has about one third of the population of the U.S., these figures are vastly different in scale from those of the U.S., which has reported 33.5 million cases of morbidity and 600 thousand cases of mortality (deaths) [8]. In light of the much lower number of patients in Japan, some criticism has arisen that low PCR testing rates might have led to the lower number of documented patients [9]. In this sense, one might regard the number of deaths as reflecting the actual situation in Japan, but with no testing-related bias. Concerning deaths, the case-fatality rate (CFR) is about 5%. In fact, the CFRs in both countries are not much different. The lower PCR testing in Japan might be related to some problems. Therefore, we specifically examined excess mortality attributable to COVID-19 in Japan, irrespective of the cause of death. In Japan, excess mortality was estimated using the National Institute of Infectious Diseases (NIID) model [10], which has been the official procedure for more than ten years. It was applied to two data sources: the national monthly deaths of all causes and the respective weekly pneumonia and influenza deaths in the 21 largest cities and their total. The latter is published regularly in Japanese during the influenza season as [https://www.niid.go.jp/niid/ja/flu-m/2112-idsc/jinsoku/131-flu-jinsoku.html](https://www.niid.go.jp/niid/ja/flu-m/2112-idsc/jinsoku/131-flu-jinsoku.html). Unfortunately, that publication ceased in March 2020 because it is intended for influenza. The first peak in Japan was April 3, 2020[11]: excess mortality cannot be detected until March. Instead, we applied NIID model to the all causes of death in the whole of Japan to evaluate impact of the outbreak of COVID-19. ## 2. Method Excess mortality is defined as the difference between the actual number of deaths and an epidemiological threshold. The epidemiological threshold is defined as the upper bound of the 95% confidence interval (CI) of the baseline. The baseline is defined as the number of deaths that are likely to have occurred if an influenza outbreak had not occurred. Therefore, if the actual deaths are fewer than the epidemiological threshold, then excess mortality is not inferred. The data used for this study were monthly deaths of all causes from 1987 through May 2021[12]. NIID model, the Stochastic Frontier Estimation [13–19], is presented as ![Formula][1] ![Formula][2] where *D**t* represents all causes of death in month/year *t, T**t* denotes the linear time trend, and *M**it* is the dummy variable for a month, which is one if *t* is the *i-*th month and otherwise zero. Moreover, *ν**t* and *ω**t* are stochastic variables as *ν**t* ∼*N*(*0, μ*2) and *ω**t* ∼*N*(*0,ξ*2); they are mutually independent. Although *ν**t* represents stochastic disturbances, *ω**t* denotes non-negative deaths attributable to influenza. These disturbance terms in this model are parameterized by two parameters: *ξ*/*μ* and (*μ**2*+*ξ**2*)**.*5*. If the null hypothesis *ξ*/*μ*=0 is not rejected, then the Stochastic Frontier Estimation model is inappropriate. Study areas were the whole of Japan and its capital, Tokyo. Study period foe estimation was from 1987 to May 2021 for the whole of Japan and up through March, 2021 for Tokyo. We adopted 5% as significant level. ## 3. Results Table 1 summarized the estimation results in the whole of Japan and Table 2 for Tokyo. Figure 1 presents observed deaths, the estimated baseline, and its threshold in Japan. Figure 2 specifically depicts the last year in Japan. We found 12 and 104 excess mortality in August and October, 2020, respectively. and 260 in May, 2021. View this table: [Table 1](http://medrxiv.org/content/early/2021/08/24/2020.07.09.20143164/T1) Table 1 NIID Model estimation results since 1987 until March 2021 in Japan View this table: [Table 2](http://medrxiv.org/content/early/2021/08/24/2020.07.09.20143164/T2) Table 2 NIID Model estimation results since 1987 until December 2020 in Tokyo ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/24/2020.07.09.20143164/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/08/24/2020.07.09.20143164/F1) Figure 1: Observations of the estimated baseline and threshold since 1987 until Mary 2021 in Japan Note: The blue line represents observations. The orange line represents the estimated baseline. The gray line shows its threshold. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/24/2020.07.09.20143164/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/08/24/2020.07.09.20143164/F2) Figure 2: Observation of the estimated baseline and threshold since January 2020 in Japan Note: The blue line represents observations. The orange line represents the estimated baseline. The gray line shows its threshold. Figure 3 and 4 showed the estimated result in Tokyo. We found 595 excess mortality in August and 150 excess mortality in September 76 in October, 458 in December, 44 in January which were 6.3, 1.7, 0.8, 4.1 and 0.4 % of the baseline. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/24/2020.07.09.20143164/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/08/24/2020.07.09.20143164/F3) Figure 3: Observations of the estimated baseline and threshold since 1987 until February 2021 in Tokyo Note: The blue line represents observations. The orange line represents the estimated baseline. The gray line shows its threshold. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/24/2020.07.09.20143164/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/08/24/2020.07.09.20143164/F4) Figure 4: Observation of the estimated baseline and threshold since May 2019 until February 2021 in Tokyo Note: The blue line represents observations. The orange line represents the estimated baseline. The gray line shows its threshold. ## 4. Discussion This study applied the NIID model to all causes of death to detect excess mortality attributable to COVID-19. We found 376 excess mortality in August and October, 2020, and in May, 2021. Moreover, we found 1323 excess mortality in Tokyo. It may not be a few excess mortality. Until January, 2021. Tokyo had 861 deaths due to COVID-19 confirmed PCR test. It means that excess motility was larger than confirmed death due to COVID-19. It may imply that some mortality due to COVID-19 were not captured as COVID-19 death, because all death was not tested by PCR. Or, outbreak of COVID-19 may increase suicide in Japan probably because of too restriction to business such as restaurant or sightseeing [20]. The amount of excess mortality in Tokyo may imply that such an indirect effect of COVID-19 outbreak might be larger than COVID-19 death confirmed by PCR test. It might be increasing after September, especially in winter, even though excess mortality in the whole of Japan were very few at least until May, 2021. On the other hand, about 10200 mortality cases caused by COVID-19 confirmed by PCR testing, were reported in MHLW officially as of May, 2021 throughout Japan [12]. Therefore, even if COVID-19 actually caused external mortality, neither the NIID model nor another statistical model might would detect significant effects attributable to COVID-19. Of course, some deaths from COVID-19 might not have been tested and might have been excluded from these official numbers of mortality cases. However, total deaths might include deaths without diagnosis as COVID-19, but actually those associated with COVID-19. Estimation results suggest that such unrecognized deaths associated with COVID-19 are not significant, even if they exist. Some researchers in Japan have emphasized considerable excess mortality from all causes of death through March, 2021 of around 28.6 thousand at maximum due to COVID-19 [21, 22] using the Farrington algorithm [23] and EuroMOMO [24], which was more than three times larger than the number of death confirmed by PCR test until March, 2021, 9113. Their study measured excess mortalities as the gap between observation and beeline, not threshold as, in prefectures where observation was higher than threshold. Therefore, their estimated too huge excess mortality may seriously mislead the risk participation for COVID-19 among the general population. In particular, in Tokyo, they found 2219 excess mortality until January, 2021. It was more than 2.5 times higher than PCR confirmed death in Tokyo, 861. At the same time, they also found approximately 57.6 thousands negative excess mortality in Japan until March, 2021 and 3679 negative excess mortality in Tokyo until March, 2021. Conversely, we found only 3 persons as negative excess mortality in February, 2020, and 751 negative excess mortality in March, 2020, in the whole of Japan, as shown in Figure 2 and 4. Such a huge negative excess mortality may doubt validity of their procedures. Their estimated baseline might be upward biased. It probably suggested that their adopted procedure have upper biased for excess mortality comparison with NIID model, which was suggested logically [25]. Moreover, they used only five years to estimate and thus volatility in data might be too small to obtain more appropriate threshold. Using pneumonia death data instead of total death data might be better to evaluate excess mortality caused by COVID-19. However, application rule of the International Classification of Diseases was revised on January 2017, after which pneumonia deaths decreased by approximately 25%. April 2020 was the fourth April since that of 2017. However, excess mortality in pneumonia death should be our next challenge. ## 5. Conclusion We found substantial excess mortality since the outbreak of COVID-19 had emerged in Tokyo between August and October and December. It might be increasing rapidly in winter season. Conversely, we found few excess mortality in the whole of Japan. It should be important to continue to monitor excess mortality of COVID-19 carefully hereafter. The present study is based on the authors’ opinions: it does not reflect any stance or policy of their professionally affiliated bodies. ## Data Availability Ministry of Health, Labour and Welfare. Preliminary statistics on demographics [https://www.mhlw.go.jp/toukei/list/81-1a.html](https://www.mhlw.go.jp/toukei/list/81-1a.html) ## 7. Conflict of interest The authors have no conflict of interest to declare. ## 8. Ethical considerations All information used for this study was published on the web site of MHLW [12]. Therefore, no ethical issue is presented. ## 6. Acknowledgement We acknowledge Dr. Nobuhiko Okabe, Kawasaki City Institute for Public Health, Dr.Kiyosu Taniguchi, National Hospital Organization Mie National Hospital, and Dr.Nahoko Shindo, WHO for their helpful support. ## Footnotes * We extend data adding one month. * Received July 9, 2020. * Revision received August 23, 2021. * Accepted August 24, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## 9. References 1. 1).Lin HC, Chiu HF, Ho SC, Yang CY. Association of influenza vaccination and reduced risk of stroke hospitalization among the elderly: a population-based case-control study. Int J Environ Res Public Health 2014; 11: 3639–49. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F24%2F2020.07.09.20143164.atom) 2. 2).Asghar Z, Coupland C, Siriwardena N. Influenza vaccination and risk of stroke: Self-controlled case-series study. Vaccine. 2015; 33: 5458–63. 3. 3).Riedmann EM. Influenza vaccination reduces risk of heart attack and stroke. Hum Vaccin Immunother 2013; 9: 2500. 4. 4).Kwok CS, Aslam S, Kontopantelis E, Myint PK, Zaman MJ, Buchan I, Loke YK, Mamas MA. Influenza, influenza-like symptoms and their association with cardiovascular risks: a systematic review and meta-analysis of observational studies. Int J Clin Pract. 2015;69:928–37. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/ijcp.12646&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25940136&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F24%2F2020.07.09.20143164.atom) 5. 5).Muhammad S, Haasbach E, Kotchourko M, Strigli A, Krenz A, Ridder DA, Vogel AB, Marti HH, Al-Abed Y, Planz O, Schwaninger M. Influenza virus infection aggravates stroke outcome. Stroke. 201142 783–91. 6. 6).Assad F, Cockburn WC, Sundaresan TK. Use of excess mortality from respiratory diseases in the study of influenza. Bull WHO 1973; 49: 219–33. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=4546520&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F24%2F2020.07.09.20143164.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1973T646500001&link_type=ISI) 7. 7).US Center for Disease Control and Prevention. Excess Deaths Associated with COVID-19. [https://www.cdc.gov/nchs/nvss/vsrr/covid19/excess\_deaths.htm](https://www.cdc.gov/nchs/nvss/vsrr/covid19/excess_deaths.htm) [accessed on July 15, 2020] 8. 8).Japan Times. Japan’s daily PCR test capacity tops 20,000. [https://www.japantimes.co.jp/news/2020/05/16/national/japans-daily-pcr-test-capacity-20000/#.XuAwImeP6AA](https://www.japantimes.co.jp/news/2020/05/16/national/japans-daily-pcr-test-capacity-20000/#.XuAwImeP6AA) [accessed on July 10, 2021] 9. 9).World Health Organization. Coronavirus disease (COVID-2019) situation reports. [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports/](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports/) [accessed on July 22,2021] 10. 10).Sugawara T, Ohkusa Y. Comparison of Models for Excess Mortality of Influenza Applied to Japan. Journal of Biosciences and Medicines, 2019, 7, 13–23. doi:10.4236/jbm.2019.76002 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4236/jbm.2019.76002&link_type=DOI) 11. 11).Kurita J, Sugawara T, Ohkusa Y. Mobility data can explain the entire COVID-19 outbreak course in Japan. medRxiv 2020.04.26.20081315; doi: [https://doi.org/10.1101/2020.04.26.20081315](https://doi.org/10.1101/2020.04.26.20081315) 12. 12).Ministry of Health, Labour and Welfare. Preliminary statistics on demographics[https://www.mhlw.go.jp/toukei/list/81-1a.html](https://www.mhlw.go.jp/toukei/list/81-1a.html) (in Japanese) [accessed on July 26, 2020] 13. 13).Aiger AD, Lovell K, Schmitidt P. Formulation and estimation of stochastic frontier production function models. Journal of Econometrics 1977; 21–37. 14. 14).Jondrow J, Lovell K, Materov S, Schmidt P. On the estimation of technical inefficiency in the stochastic frontier production function model. Journal of Econometrics 1982; 233–9. 15. 15).Li T, Rosenman R. Cost inefficiency in Washington Hospitals: A stochastic frontier approach using panel data, Health Care Management Science 2001; 4: 73–81. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1011493209102&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11393744&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F24%2F2020.07.09.20143164.atom) 16. 16).Newhouse JP. Frontier Estimation: How useful a tool for health economics? Journal of Health Economics 1994; 13: 317–22. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0167-6296(94)90030-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10138857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F24%2F2020.07.09.20143164.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994PQ49300004&link_type=ISI) 17. 17).Shelton Brown H. Managed care and technical efficiency. Health Economics. 2003; 12: 149–58. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hec.712&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12563662&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F24%2F2020.07.09.20143164.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000181005400006&link_type=ISI) 18. 18).Jacobs R. Alternative methods to examine hospital efficiency: Data envelopment analysis and stochastic frontier analysis. Health Care Management Science. 2001; 4: 103–15 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1011453526849&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11393739&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F24%2F2020.07.09.20143164.atom) 19. 19).Rosko MD. Cost efficiency of US hospitals: A stochastic frontier approach. Health Economics. 2001; 539–51. 20. 20).Kurita J, Sugawara T, Ohkusa Y. Excess Mortality in Suicide caused by COVID-19 in Japan. doi: [https://doi.org/10.1101/2021.02.13.21251670](https://doi.org/10.1101/2021.02.13.21251670) 21. 21).Excess and Exiguous Deaths Dashboard in Japan. [https://exdeaths-japan.org/](https://exdeaths-japan.org/) [accessed on July 7, 2021] 22. 22).Kawashima T, Nomura S, Tanoue Y, Yoneoka D, Eguchi A, Ng CFS, Matsuura K, Shi S, Makiyama K, Uryu S, Kawamura Y, Takayanagi S, Gilmour S, Miyata H, Sunagawa T, Takahashi T, Tsuchihashi Y, Kobayashi Y, Arima Y, Kanou K, Hashizume M. Excess All-Cause Deaths during Coronavirus Disease Pandemic, Japan, January-May 2020(1). Emerg Infect Dis. 2021 Mar;27(3):789–795. doi: 10.3201/eid2703.203925. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2703.203925&link_type=DOI) 23. 23).Center of Disease Control and Prevention, Excess Deaths Associated with COVID-19 [https://data.cdc.gov/NCHS/Excess-Deaths-Associated-with-COVID-19/xkkf-xrst](https://data.cdc.gov/NCHS/Excess-Deaths-Associated-with-COVID-19/xkkf-xrst) [accessed on June 23, 2021] 24. 24).EUROMOMO,EUROMOMO, [https://www.euromomo.eu/](https://www.euromomo.eu/) [accessed on June 23, 2021] 25. 25).Sugawara T, Ohkusa Y. Comparison of Models for Excess Mortality of Influenza Applied to Japan. Journal of Biosciences and Medicines, 2019, 7, 13–23. doi:10.4236/jbm.2019.76002 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4236/jbm.2019.76002&link_type=DOI) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif