Effects of voluntary event cancellation and school closure as countermeasures against COVID-19 outbreak in Japan

Background To control the COVID-19 outbreak in Japan, sports and entertainment events were canceled and schools were closed throughout Japan from February 26 through March 19. That policy has been designated as voluntary event cancellation and school closure (VECSC). Object This study assesses VECSC effectiveness based on predicted outcomes. Methods A simple susceptible–infected–recovered model was applied to data of patients with symptoms in Japan during January 14 through March 26. The respective reproduction numbers for periods before VECSC (R0), during VECSC (Re), and after VECSC (Ra) were estimated. Results Results suggest R0 before VECSC as 2.534 [2.449, 2.598], Re during VECSC as 1.077 [0.948, 1.228], and Ra after VECSC as 4.455 [3.615, 5.255]. Discussion and conclusion Results demonstrated that VECSC can reduce COVID-19 infectiousness considerably, but after VECSC, the value of the reproduction number rose to exceed 4.0.


Introduction
The initial case of COVID-19 in Japan was an airline passenger returning from Wuhan, China on January 3, 2020. The Ministry of Health, Labour and Welfare (MHLW) in Japan would report 3,906 cases in Japan by April 6, 2020, including asymptomatic cases but pointedly excluding those from a large cruise ship: the Diamond Princess [1]. Sports and entertainment events were canceled in Japan for two weeks from February 26 through March 11 in compliance with a government advisory. It had extended until March 19. At that time, it was advised that small business and private meetings be cancelled voluntarily: the measure was not enacted as a law or enforced by authorities. Therefore, people were not arrested or cited even if they did not comply with this government advisory [2]; the effort depended entirely on the voluntary compliance of individuals. Moreover, this was the first measure necessitating voluntary event cancellation. Those characteristics of the measures complicate the ex ante prediction of the proportion of events that were cancelled and the extent to which contact among people was reduced. Moreover, as announced March 3, almost all schools were closed from the middle of March through spring vacation (early April) as a measure to control the spread of COVID-19. Even though young people can be infected and can transmit the virus to adults, school closure effects were questionable ex ante because schoolchildren were not regarded as the most susceptible age class for COVID-19 [3][4][5]. Therefore, these policies must necessarily be evaluated ex post. The policies are known collectively as voluntary event cancellation and school closure (VECSC). After the VECSC period, the VECSC measures ceased, but a tacit norm of cancellation of large events continued.
If the reproduction number under the VECSC environment was less than one, then the outbreak can be expected to have been contained. Alternatively, even if the reproduction number under them were markedly less than the reproduction number before VECSC but greater than one, one could have expected instead that the outbreak would be exacerbated.
This study was conducted to evaluate VECSC based on the epidemic curve because no evidence exists to reflect the number of cancelled events. The obtained result might be expected to contribute to a government decision to implement VECSC as a countermeasure in Japan.

Methods
We applied a simple deterministic susceptible-infected-recovered (SIR) model [6] to the epidemic curve in Japan for its population of 120 million. We assume an incubation period following the empirical distribution from the early stage of the outbreak. A person on the i-th day of incubation will move to a symptomatic or an asymptomatic state with probability of p(i)/S j = i L p(j), where p(i) represents the proportion of the i-th day on the empirical distribution of the incubation period. Also, L denotes the maximum length of the incubation period. Alternatively, a person moves to the i+1 day in the incubation state with one minus p(i)/S j = i L p(j). For the model, symptomatic and asymptomatic states are assumed to continue for seven days, followed by a transition to a recovery state with probability of one [7]. This model ignores outcomes and necessary medical resources: states of death and hospitalization are not incorporated into the model.
Asymptomatic cases are not observable unless complete laboratory-based surveillance is performed. One exceptional study of a sample of elderly people indicated them as 3 cases among 23 [7]. We checked its robustness for the proportion of asymptomatic cases assuming the proportion of 4 cases among 19 from pre-outbreak experiences of Japanese residents of Wuhan [8].
Infectivity by severely infected patients and by mildly infected patients was assumed to be equal. Moreover, we assumed for simplicity that asymptomatic cases have the same capability of infecting others as symptomatic cases have [7]. Therefore, we verified the model robustness through sensitivity analysis of the infectiousness of asymptomatic cases, such as 50% assumed from simulation studies of influenza infections [9][10][11][12][13]. The distribution of infectiousness by symptomatic and asymptomatic cases was assumed to be 30% on the onset day, 20% on the following day, and 10% for the subsequent five days [7]. Regarding its robustness, we also applied a uniform distribution for a week.
Because VECSC policies were administered during February 26-March 19, we divided the full data period into three sub-periods: before VECSC, during VECSC, and after VECSC. Their corresponding reproduction numbers are represented as R 0 , R e , and R a . We modeled the change in the reproduction number based on the prior day's reported number of persons susceptible, incubating, symptomatic, asymptomatic, and in a recovered state.
For this study, R 0 was defined as the number of newly infected persons from one person infected over the infectious period when all persons were susceptible with no countermeasures. Also, R e is defined by the numbers of newly infected persons from one person infected over the infectious period when all persons were susceptible under VECSC. Finally, R a is defined by the numbers of newly infected persons from one person infected during a time of countermeasures after the VECSC period. For estimation, the number of newly infected persons was calculated using R 0 (R e or R a ) × the proportion of susceptible persons among the population × the number of persons with infectiousness weighted by their degree of infectiousness.
The values of R 0 , R e , and R a were sought to fit the data to minimize the sum of absolute values of discrepancies among the bootstrapped epidemic curve and the fitted values. The estimated relative values of the three reproduction numbers were calculated using 10,000 iterations of bootstrapping for empirical distribution of the data for symptomatic patients. In this sense, although this model is deterministic at each bootstrapping iteration, its range can be estimated through 10,000 iterations.
The bootstrapping procedure we used has fully replicated bootstrapping for a fixed number of initial cases. There were N patients in the data, with numbering of the patients from the initial case to the last case. Initially, no patient was on the bootstrapped epidemic curve. If a random variable drawn from a uniform distribution of (0,1) was t included in the internal [i/(N-1), (i+1)/(N-1)], then we added one to the onset date of the i+1th patient to the bootstrapped epidemic curve. We replicated this procedure N-1 times. Thereby, we obtained the bootstrapped epidemic curve with N-1 patients. Finally, we added the initial patient, for whom the onset date was January 14, to the bootstrapped epidemic curve. Consequently, we obtained a bootstrapped epidemic curve with N patients starting from January 14.
We estimated the curve sequentially as described hereinafter. First, we estimated R 0 as the best fit to bootstrapped data for the pre-VECSC period. Then, based on the obtained R 0 and on the course of the outbreak before the VECSC period, we estimated R e as the best fit to bootstrapped data in the VECSC period. Finally, based on the obtained values of R 0 and R e , we estimated R a as the best fit to bootstrapped data for the post-VECSC period. At each step, reproduction numbers were grid searched in the interval of (0,10) by 0.001.
We applied some one-way sensitivity analyses in addition to the base case as explained above: the infectious power of asymptomatic cases was 50% of the power of symptomatic cases; the proportion of asymptomatic cases was 4 /19 of symptomatic cases; infectiousness and infectious patterns showed a uniform distribution for one week. We also estimated the parameters by minimizing the sum of the squared residuals or maximizing Poisson likelihood instead of the difference in absolute values and using the fitted incubation period as a log normal distribution instead of the empirical distribution.
Moreover, we examined the possibility of infection during two days before onset. Its expected number can be represented as , where x(i,t) denotes the number of persons at the i-th incubation period on day t. The first term, indicates the expected number of patients who do not show symptoms day t but who will show symptoms on day t+1, except for the third term. p(i)/S j = i L p(j) indicates probability of showing symptoms among patients who did not show symptoms until i-1 th period but will show symptom on day t +1. Similarly, the second term of t) indicates the expected number of patients who do not show symptoms on day t or t +1 but who will show symptoms on day t +2. The third term, x(L-1,t), represents patients who showed no symptom until the last day of the incubation period. Therefore, they will certainly show symptoms on the following day of the last day of incubation period. Their infectiousness was assumed to be the average of the first two days of the symptomatic period. To meet the definition of the reproduction number, infectiousness in the incubation and (a) symptomatic period were adjusted such that the sum of infectiousness from two days before onset to one week following onset would be one.

PLOS ONE
We also examined incorporation of the imported cases into the model. They were entered into the model as newly symptomatic cases with no history of incubation period when the imported cases were reported. The parameters in this model were estimated as minimizing the sum of the absolute difference among the number of patients including the imported cases and prediction. We adopted 5% as the level for which we inferred significance. We used software (Matlab 2014a; The MathWorks Inc.) to code the model as explained above.

Data source
The numbers of symptomatic patients during January 14-March 26were published by the MHLW [1]. The data were reported publicly from 46 prefectures and some cities. They were summarized in S1 File as of April 6. During this period, 1516 cases were recorded with onset dates. We excluded imported cases and cases representing infected persons from the Diamond Princess because they were presumed not to have been community-acquired in Japan.

Ethical considerations
All information used for this study has been published elsewhere [1]. This study therefore poses no ethical issues.

Results
Fig 1 depicts the empirical distribution of the incubation period among 125 cases for which the exposed date and onset date were published by MHLW. Its mode was six days. The average was 6.1 days. The maximum incubation period length was 19 days. Because the ranges of R e and R a were overlapped in all cases, these factors might not greatly affect both parameters, at least in the region with parameters and distribution that are considered reasonable. The estimated R 0 was the most sensitive. When using Poisson likelihood for evaluation, R 0 was significantly lower than in the base case. Moreover, infectiousness before onset raised R 0 significantly. In other scenarios, the ranges of R 0 were overlapped.

Discussion
We applied a simple deterministic SIR model including asymptomatic cases. A deterministic model might be simpler but more appropriate than a stochastic model for the study including changing parameters and policy evaluation [14].

Fig 1. Empirical distribution of incubation period published by Ministry of Health, Labour and
Welfare, Japan. Bars shows distribution of incubation period among 91 cases whose exposed date and onset date were published by Ministry of Health, Labour and Welfare, Japan. Because the patients whose incubation was longer than 14 days were just 4%, they were included as bars on 14. https://doi.org/10.1371/journal.pone.0239455.g001 An earlier study [15][16][17] estimated R 0 for COVID-19 as 2.24-3.58 in Wuhan. Our R 0 obtained for the period before VECSC was similar. However, an earlier other study [18] found the average number of secondary cases infected by a single primary case found in contact tracing in Japan as 0.6 before the VECSC period, with 80% of patients found to have infected no one. If the average number of secondarily infected persons was less than one, then the outbreak might have stalled on its own. Given such an eventuality, no countermeasures would have been required. However, based on results of a study uploaded to the site of Japanese Public Health Association, health authorities in the public health sector or local government were urged to conduct contact tracing to detect clusters [19].
Infectiousness two days before symptom onset has been described in the literature [7,20]. We applied it as sensitivity analysis. It showed slightly higher reproduction numbers than the base model, although the estimated R e and R a were found to have a non-significant difference from the base case because a longer period of infectiousness implies fewer infected people on a day when the reproduction number does not change. However, the epidemic curve was given. Therefore, the estimated reproduction number probably increased by extending the period of infectiousness from seven days to nine days. Moreover, when using Poisson likelihood for evaluation, R 0 alone was found to be significantly lower than in the base case. Although the difference was not large, it might engender some difference attributable to selection of the evaluation function.
Among children and younger adults, the proportion of asymptomatic infected people might be larger than among elderly people, but that point remains uncertain. Lower susceptibility of children [3][4][5] might imply a higher proportion of asymptomatic cases in children. In this sense, the finding of 3 cases divided by 23 among elderly people might be the lower bound of the proportion for the general population. Overall, sensitivity analysis led to similar estimation with the base case. Even though significant difference was found through comparison with the base case only in R 0 , differences among them were not large, quantitatively speaking. Therefore, one can infer that the obtained results are robust.
The average incubation period applied in an earlier study was 5.1 days [21], but it was 6.1 days for the present study. However, the latter number based on the epidemic curve considered in the present study is certainly consistent with the data. It might reflect lower severity among patients in Japan than in other countries such as the US. The estimated value of the reproduction number probably would be smaller if we were to adopt a slightly shorter incubation period, which is a similar result to that found with sensitivity analysis for infectiousness before onset.
Many weeks have passed since the end of the study period. Therefore, the reporting lag almost disappeared. If timely estimation, meaning estimation using data from a day prior, was necessary, then we used an adjusted delay for reported data [8,22,23].
Instead of using the minimized sum of squared residuals, such as maximum likelihood estimation based on a normal distribution, we used the minimized sum of the absolute values of discrepancies for the bootstrapped epidemic curve and the fitted values. In general, minimization of the sum of the absolute values is more robust than minimization of the sum of squares because the absolute value is less sensitive than squared values to the effects of outliers [24][25][26]. The daily epidemic curve sometimes exhibits spikes according to the day of the week, and for other reasons. Especially, few patients per day were reported during the early stage of the outbreak. Therefore, spikes might be quite large. These spikes are probably outliers. They might overly affect the estimator. For this reason, we prefer minimization of the sum of the absolute values to minimization of the sum of squared values when analyzing daily data in the earlier stage of the outbreak.
Moreover, we chose non-parametric approaches using actual data in preference to parametric approaches assuming a particular distribution. Such an assumption might affect results through miss-specification. Therefore, we prefer to use the actual distribution of incubation periods or epidemic curves with no unnecessary or restrictive assumption. Sensitivity analysis using minimum least squares of errors or using Poisson likelihood for maximization showed very similar results to those obtained with the base case. Although R 0 and R a were slightly higher, R e were slightly lower. That result might indicate that some outlier in the pre-VECSC and post-VECSC period pulls the estimator in sensitivity analysis.
The present study has some limitations. The first is that even though we evaluated VECSC, the respective effects of voluntary event cancellation and school closure cannot be distinguished. To distinguish their respective effects, one would have to develop a model including several age classes. School closure mainly affects contact patterns among schoolchildren; voluntary event cancellation mainly affects patterns among adults. Therefore, studies of those respective age groups might elucidate the separate effects of these policies. Such studies stand as challenges for our future study.
The second point is under ascertainment. Although the epidemic curves of COVID-19 in all countries are subject to under ascertainment, it might be very difficult to evaluate the degree to which they are affected. It might bias the estimation result.
A third point is a lack of estimation of the outcomes, such as dead or severe cases or necessary medical resources for the care of COVID-19 patients. We particularly examined how policies affect the reproduction number under the effects of imposed countermeasures. Therefore, we ignored prediction of the entire course of the outbreak and its outcomes such as the number of deaths. Nevertheless, outcomes are expected to be a primary concern for modelling. Moreover, the collapse of medical services can be expected to engender worse outcomes even if the reproduction number remains unchanged. Prediction of the effects of severe policies including lockdowns is anticipated as a challenge to be addressed in our future research.

Conclusion
Results have demonstrated that VECSC can reduce the infectiousness of COVID-19 considerably: approximately to one. However, the figure is probably greater than one. Outbreaks might continue for a long time. Therefore, lockdown policies are expected to be as effective as VECSC if they are executed carefully. After VECSC, the reproduction number escalated again beyond that before VECSC. Similar phenomena might be observed after a lockdown. It is our earnest hope that results of the present study can contribute to governmental policy-making related to lockdown measures and other countermeasures used to combat the spread and destruction of COVID-19.
The present study was based on the authors' opinions. Neither the results nor implications reflect any stance or policy of professionally affiliated bodies.