Uncertainty quantification for the random HIV dynamical model driven by drug adherence.

In HIV drug therapy, the high variability of CD4+ T cells and viral loads brings uncertainty to the determination of treatment options and the ultimate treatment efficacy, which may be the result of poor drug adherence. We develop a dynamical HIV model coupled with pharmacokinetics, driven by drug adherence as a random variable, and systematically study the uncertainty quantification, aiming to construct the relationship between drug adherence and therapeutic effect. Using adaptive generalized polynomial chaos, stochastic solutions are approximated as polynomials of input random parameters. Numerical simulations show that results obtained by this method are in good agreement, compared with results obtained through Monte Carlo sampling, which helps to verify the accuracy of approximation. Based on these expansions, we calculate the time-dependent probability density functions of this system theoretically and numerically. To verify the applicability of this model, we fit clinical data of four HIV patients, and the goodness of fit results demonstrate that the proposed random model depicts the dynamics of HIV well. Sensitivity analyses based on the Sobol index indicate that the randomness of drug effect has the greatest impact on both CD4+ T cells and viral loads, compared to random initial values, which further highlights the significance of drug adherence. The proposed models and qualitative analysis results, along with monitoring CD4+ T cells counts and viral loads, evaluate the influence of drug adherence on HIV treatment, which helps to better interpret clinical data with fluctuations and makes several contributions to the design of individual-based optimal antiretroviral strategies.

Combining the dynamic model and deep neural networks to identify the intensity of interventions during COVID-19 pandemic.

During the COVID-19 pandemic, control measures, especially massive contact tracing following prompt quarantine and isolation, play an important role in mitigating the disease spread, and quantifying the dynamic contact rate and quarantine rate and estimate their impacts remain challenging. To precisely quantify the intensity of interventions, we develop the mechanism of physics-informed neural network (PINN) to propose the extended transmission-dynamics-informed neural network (TDINN) algorithm by combining scattered observational data with deep learning and epidemic models. The TDINN algorithm can not only avoid assuming the specific rate functions in advance but also make neural networks follow the rules of epidemic systems in the process of learning. We show that the proposed algorithm can fit the multi-source epidemic data in Xi'an, Guangzhou and Yangzhou cities well, and moreover reconstruct the epidemic development trend in Hainan and Xinjiang with incomplete reported data. We inferred the temporal evolution patterns of contact/quarantine rates, selected the best combination from the family of functions to accurately simulate the contact/quarantine time series learned by TDINN algorithm, and consequently reconstructed the epidemic process. The selected rate functions based on the time series inferred by deep learning have epidemiologically reasonable meanings. In addition, the proposed TDINN algorithm has also been verified by COVID-19 epidemic data with multiple waves in Liaoning province and shows good performance. We find the significant fluctuations in estimated contact/quarantine rates, and a feedback loop between the strengthening/relaxation of intervention strategies and the recurrence of the outbreaks. Moreover, the findings show that there is diversity in the shape of the temporal evolution curves of the inferred contact/quarantine rates in the considered regions, which indicates variation in the intensity of control strategies adopted in various regions.

Mathematical modeling the order of driver gene mutations in colorectal cancer.

Tumor heterogeneity is a large obstacle for cancer study and treatment. Different cancer patients may involve different combinations of gene mutations or the distinct regulatory pathways for inducing the progression of tumor. Investigating the pathways of gene mutations which can cause the formation of tumor can provide a basis for the personalized treatment of cancer. Studies suggested that KRAS, APC and TP53 are the most significant driver genes for colorectal cancer. However, it is still an open issue regarding the detailed mutation order of these genes in the development of colorectal cancer. For this purpose, we analyze the mathematical model considering all orders of mutations in oncogene, KRAS and tumor suppressor genes, APC and TP53, and fit it on data describing the incidence rates of colorectal cancer at different age from the Surveillance Epidemiology and End Results registry in the United States for the year 1973-2013. The specific orders that can induce the development of colorectal cancer are identified by the model fitting. The fitting results indicate that the mutation orders with KRAS → APC → TP53, APC → TP53 → KRAS and APC → KRAS → TP53 explain the age-specific risk of colorectal cancer with very well. Furthermore, eleven pathways of gene mutations can be accepted for the mutation order of genes with KRAS → APC → TP53, APC → TP53 → KRAS and APC → KRAS → TP53, and the alternation of APC acts as the initiating or promoting event in the colorectal cancer. The estimated mutation rates of cells in the different pathways demonstrate that genetic instability must exist in colorectal cancer with alterations of genes, KRAS, APC and TP53.

Model analysis and data validation of structured prevention and control interruptions of emerging infectious diseases.

The design of optimized non-pharmaceutical interventions (NPIs) is critical to the effective control of emergent outbreaks of infectious diseases such as SARS, A/H1N1 and COVID-19 and to ensure that numbers of hospitalized cases do not exceed the carrying capacity of medical resources. To address this issue, we formulated a classic SIR model to include a close contact tracing strategy and structured prevention and control interruptions (SPCIs). The impact of the timing of SPCIs on the maximum number of non-isolated infected individuals and on the duration of an infectious disease outside quarantined areas (i.e. implementing a dynamic zero-case policy) were analyzed numerically and theoretically. These analyses revealed that to minimize the maximum number of non-isolated infected individuals, the optimal time to initiate SPCIs is when they can control the peak value of a second rebound of the epidemic to be equal to the first peak value. More individuals may be infected at the peak of the second wave with a stronger intervention during SPCIs. The longer the duration of the intervention and the stronger the contact tracing intensity during SPCIs, the more effective they are in shortening the duration of an infectious disease outside quarantined areas. The dynamic evolution of the number of isolated and non-isolated individuals, including two peaks and long tail patterns, have been confirmed by various real data sets of multiple-wave COVID-19 epidemics in China. Our results provide important theoretical support for the adjustment of NPI strategies in relation to a given carrying capacity of medical resources.

Hormetic and synergistic effects of cancer treatments revealed by modelling combinations of radio - or chemotherapy with immunotherapy.

BACKGROUND: Radio/chemotherapy and immune systems provide examples of hormesis, as tumours can be stimulated (or reduced) at low radio/chemical or antibody doses but inhibited (or stimulated) by high doses. METHODS: Interactions between effector cells, tumour cells and cytokines with pulsed radio/chemo-immunotherapy were modelled using a pulse differential system. RESULTS: Our results show that radio/chemotherapy (dose) response curves (RCRC) and/or immune response curves (IRC) or a combination of both, undergo homeostatic changes or catastrophic shifts revealing hormesis in many parameter regions. Some mixed response curves had multiple humps, posing challenges for interpretation of clinical trials and experimental design, due to a fuzzy region between an hormetic zone and the toxic threshold. Mixed response curves from two parameter bifurcation analyses demonstrated that low-dose radio/chemotherapy and strong immunotherapy counteract side-effects of radio/chemotherapy on effector cells and cytokines and stimulate effects of immunotherapy on tumour growth. The implications for clinical applications were confirmed by good fits to our model of RCRC and IRC data. CONCLUSIONS: The combination of low-dose radio/chemotherapy and high-dose immunotherapy is very effective for many solid tumours. The net benefit and synergistic effect of combined therapy is conducive to the treatment and inhibition of tumour cells.

Exploring COVID-19 transmission patterns and key factors during epidemics caused by three major strains in Asia.

The COVID-19 epidemic has lasted for more than two years since the outbreak in late 2019. An urgent and challenging question is how to systematically evaluate epidemic developments in different countries, during different periods, and to determine which measures that could be implemented are key for successful epidemic prevention. In this study, SBD distance-based K-shape clustering and hierarchical clustering methods were used to analyse epidemics in Asian countries. For the hierarchical clustering, epidemic time series were divided into three periods (epidemics induced by the Original/Alpha, Delta and Omicron variants separately). Standard deviations, the Hurst index, mortality rates, peak value of confirmed cases per capita, average growth rates, and the control efficiency of each period were used to characterize the epidemics. In addition, the total numbers of cases in the different countries were analysed by correlation and regression in relation to 15 variables that could have impacts on COVID-19. Finally, some suggestions on prevention and control measures for each category of country are given. We found that the total numbers of cases per million of a population, total deaths per million and mortality rates were highly correlated with the proportion of people aged over 65 years, the prevalence of multiple diseases, and the national GDP. We also found significant associations between case numbers and vaccination rates, health expenditures, and stringency of control measures. Vaccinations have played a positive role in COVID-19, with a gradual decline in mortality rates in later periods, and are still playing protective roles against the Delta and Omicron strains. The stringency of control measures taken by a government is not an indicator of the appropriateness of a country's response to the outbreak, and a higher index does not necessarily mean more effective measures; a combination of factors such as national vaccination rates, the country's economic foundation and the availability of medical equipment is also needed. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

Coupling an individual adaptive-decision model with a SIRV model of influenza vaccination reveals new insights for epidemic control.

Past seasonal influenza epidemics and vaccination experience may affect individuals' decisions on whether to be vaccinated or not, decisions that may be constantly reassessed in relation to recent influenza related experience. To understand the potentially complex interaction between experience and decisions and whether the vaccination rate is likely to reach a critical coverage level or not, we construct an adaptive-decision model. This model is then coupled with an influenza vaccination dynamics (SIRV) model to explore the interaction between individuals' decision-making and an influenza epidemic. Nonlinear least squares estimation is used to obtain the best-fit parameter values in the SIRV model based on data on new influenza-like illness (ILI) cases in Texas. Uncertainty and sensitivity analyses are then carried out to determine the impact of key parameters of the adaptive decision-making model on the ILI epidemic. The results showed that the necessary critical coverage rate of ILI vaccination could not be reached by voluntary vaccination. However, it could be reached in the fourth year if mass media reports improved individuals' memory of past vaccination experience. Individuals' memory of past vaccination experience, the proportion with histories of past vaccinations and the perceived cost of vaccination are important factors determining whether an ILI epidemic can be effectively controlled or not. Therefore, health authorities should guide people to improve their memory of past vaccination experience through media reports, publish timely data on annual vaccination proportions and adjust relevant measures to appropriately reduce vaccination perceived cost, in order to effectively control an ILI epidemic.

Assessing the transmissibility of epidemics involving epidemic zoning.

BACKGROUND: Epidemic zoning is an important option in a series of measures for the prevention and control of infectious diseases. We aim to accurately assess the disease transmission process by considering the epidemic zoning, and we take two epidemics with distinct outbreak sizes as an example, i.e., the Xi'an epidemic in late 2021 and the Shanghai epidemic in early 2022. METHODS: For the two epidemics, the total cases were clearly distinguished by their reporting zone and the Bernoulli counting process was used to describe whether one infected case in society would be reported in control zones or not. Assuming the imperfect or perfect isolation policy in control zones, the transmission processes are respectively simulated by the adjusted renewal equation with case importation, which can be derived on the basis of the Bellman-Harris branching theory. The likelihood function containing unknown parameters is then constructed by assuming the daily number of new cases reported in control zones follows a Poisson distribution. All the unknown parameters were obtained by the maximum likelihood estimation. RESULTS: For both epidemics, the internal infections characterized by subcritical transmission within the control zones were verified, and the median control reproduction numbers were estimated as 0.403 (95% confidence interval (CI): 0.352, 0.459) in Xi'an epidemic and 0.727 (95% CI: 0.724, 0.730) in Shanghai epidemic, respectively. In addition, although the detection rate of social cases quickly increased to 100% during the decline period of daily new cases until the end of the epidemic, the detection rate in Xi'an was significantly higher than that in Shanghai in the previous period. CONCLUSIONS: The comparative analysis of the two epidemics with different consequences highlights the role of the higher detection rate of social cases since the beginning of the epidemic and the reduced transmission risk in control zones throughout the outbreak. Strengthening the detection of social infection and strictly implementing the isolation policy are of great significance to avoid a larger-scale epidemic.

Lessons drawn from Shanghai for controlling highly transmissible SARS-CoV-2 variants: insights from a modelling study.

BACKGROUND: The continuous emergence of novel SARS-CoV-2 variants with markedly increased transmissibility presents major challenges to the zero-COVID policy in China. It is critical to adjust aspects of the policy about non-pharmaceutical interventions (NPIs) by searching for and implementing more effective ways. We use a mathematical model to mimic the epidemic pattern of the Omicron variant in Shanghai to quantitatively show the control challenges and investigate the feasibility of different control patterns in avoiding other epidemic waves. METHODS: We initially construct a dynamic model with a core step-by-step release strategy to reveal its role in controlling the spread of COVID-19, including the city-based pattern and the district-based pattern. We used the least squares method and real reported case data to fit the model for Shanghai and its 16 districts, respectively. Optimal control theory was utilized to explore the quantitative and optimal solutions of the time-varying control strength (i.e., contact rate) to suppress the highly transmissible SARS-CoV-2 variants. RESULTS: The necessary period for reaching the zero-COVID goal can be nearly 4 months, and the final epidemic size was 629,625 (95%CI: [608,049, 651,201]). By adopting the city-based pattern, 7 out of 16 strategies released the NPIs more or earlier than the baseline and ensured a zero-resurgence risk at the average cost of 10 to 129 more cases in June. By adopting the district-based pattern, a regional linked release can allow resumption of social activity to ~ 100% in the boundary-region group about 14 days earlier and allow people to flow between different districts without causing infection resurgence. Optimal solutions of the contact rate were obtained with various testing intensities, and higher diagnosis rate correlated with higher optimal contact rate while the number of daily reported cases remained almost unchanged. CONCLUSIONS: Shanghai could have been bolder and more flexible in unleashing social activity than they did. The boundary-region group should be relaxed earlier and more attention should be paid to the centre-region group. With a more intensive testing strategy, people could return to normal life as much as possible but still ensure the epidemic was maintained at a relatively low level.

Statistical Analysis of Two-Compartment Pharmacokinetic Models with Drug Non-adherence.

Poor drug adherence is considered one of major barriers to achieving the clinical and public health benefits of many pharmacotherapies. In the current paper, we aim to investigate the impact of dose omission on the plasma concentrations of two-compartment pharmacokinetic models with two typical routes of drug administration, namely the intravenous bolus and extravascular first-order absorption. First, we reformulate the classical two-compartment pharmacokinetic models with a new stochastic feature, where a binomial random model of dose intake is integrated. Then, we formalize the explicit expressions of expectation and variance for trough concentrations and limit concentrations, with the latter proved of the existence and uniqueness for steady-state distribution. Moreover, we mathematically demonstrate the strict stationarity and ergodicity of trough concentrations as a Markov chain. In addition, we numerically simulate the impact of drug non-adherence to different extents on the variability and regularity of drug concentration and compare the drug pharmacokinetic preference between one and two compartment pharmacokinetic models. The results of sensitivity analysis also suggest the drug non-adherence as one of the most sensitive model parameters to the expectation of limit concentration. Our modelling and analytical approach can be integrated into the chronic disease models to estimate or quantitatively predict the therapy efficacy with drug pharmacokinetics presumably affected by random dose omissions.

Threshold conditions for curbing COVID-19 with a dynamic zero-case policy derived from 101 outbreaks in China.

By 31 May 2022, original/Alpha, Delta and Omicron strains induced 101 outbreaks of COVID-19 in mainland China. Most outbreaks were cleared by combining non-pharmaceutical interventions (NPIs) with vaccines, but continuous virus variations challenged the dynamic zero-case policy (DZCP), posing questions of what are the prerequisites and threshold levels for success? And what are the independent effects of vaccination in each outbreak? Using a modified classic infectious disease dynamic model and an iterative relationship for new infections per day, the effectiveness of vaccines and NPIs was deduced, from which the independent effectiveness of vaccines was derived. There was a negative correlation between vaccination coverage rates and virus transmission. For the Delta strain, a 61.8% increase in the vaccination rate (VR) reduced the control reproduction number (CRN) by about 27%. For the Omicron strain, a 20.43% increase in VR, including booster shots, reduced the CRN by 42.16%. The implementation speed of NPIs against the original/Alpha strain was faster than the virus's transmission speed, and vaccines significantly accelerated the DZCP against the Delta strain. The CRN ([Formula: see text]) during the exponential growth phase and the peak time and intensity of NPIs were key factors affecting a comprehensive theoretical threshold condition for DZCP success, illustrated by contour diagrams for the CRN under different conditions. The DZCP maintained the [Formula: see text] of 101 outbreaks below the safe threshold level, but the strength of NPIs was close to saturation especially for Omicron, and there was little room for improvement. Only by curbing the rise in the early stage and shortening the exponential growth period could clearing be achieved quickly. Strengthening China's vaccine immune barrier can improve China's ability to prevent and control epidemics and provide greater scope for the selection and adjustment of NPIs. Otherwise, there will be rapid rises in infection rates and an extremely high peak and huge pressure on the healthcare system, and a potential increase in excess mortality.

The effectiveness of various control strategies: An insight from a comparison modelling study.

Several local outbreaks have occurred and been suppressed with the dynamic zero-COVID policy and widely promoted vaccination program implemented in China. The epidemic duration and final size vary significantly in different cities, which may be attributed to different implementation patterns and intensities of non-pharmaceutical interventions (NPIs). It's important to capture the underlying mechanism to explore more efficient implementation patterns of NPIs in order to prevent the future resurgence. In this study, outbreaks caused by Delta variant in Xi'an, Yangzhou and Guangzhou in 2021 are chosen as the examples. A novel model dividing the population into three groups is proposed to describe the heterogeneity of control interventions. The model is calibrated and key parameters related to NPIs are estimated by using multi-source epidemic data. The estimation results show a lower transmission probability but a higher initial reproduction number in Xi'an. Sensitivity analysis are conducted to investigate the impact of various control measures in different epidemic phases. The results identify the vital role of enhancing closed-off management, strengthening tracing and testing intensities, on shortening the epidemic durations and reducing the final size. Further, we find that sufficiently implemented closed-off management would prevent the city from lockdown. Strengthening the tracing other than the testing strategy in the initial stage is more effective on containing the epidemic in a shorter duration with less infections.

A quantitative method to project the probability of the end of an epidemic: Application to the COVID-19 outbreak in Wuhan, 2020.

The end-of-outbreak declaration is an important part of epidemic control, marking the relaxation or cancellation of prevention and control measures. We propose a probability model to retrospectively quantify the confidence of giving the end-of-outbreak declaration during the COVID-19 epidemic in early 2020 in Wuhan. By using the linear spline, we firstly estimates the time-varying proportion of cases who miss the nonpharmaceutical interventions (NPIs) among all reported cases. Assuming the reproduction numbers being 1.5, 2.0, 3.0, 4.0, 5.0 and 6.0, the respective probability of the end of the COVID-19 outbreak with time after the last reported case can be iteratively computed. Consequently, the varying reproduction numbers produce slightly different increasing patterns of NPI effectiveness, and the end-of-outbreak declarations with 95% confidence are projected consistently earlier than the day when the lockdown was actually lifted. The reason for the timing discrepancy is discussed as well.

Assessing Age-Specific Vaccination Strategies and Post-vaccination Reopening Policies for COVID-19 Control Using SEIR Modeling Approach.

As the availability of COVID-19 vaccines, it is badly needed to develop vaccination guidelines to prioritize the vaccination delivery in order to effectively stop COVID-19 epidemic and minimize the loss. We evaluated the effect of age-specific vaccination strategies on the number of infections and deaths using an SEIR model, considering the age structure and social contact patterns for different age groups for each of different countries. In general, the vaccination priority should be given to those younger people who are active in social contacts to minimize the number of infections, while the vaccination priority should be given to the elderly to minimize the number of deaths. But this principle may not always apply when the interaction of age structure and age-specific social contact patterns is complicated. Partially reopening schools, workplaces or households, the vaccination priority may need to be adjusted accordingly. Prematurely reopening social contacts could initiate a new outbreak or even a new pandemic out of control if the vaccination rate and the detection rate are not high enough. Our result suggests that it requires at least nine months of vaccination (with a high vaccination rate > 0.1%) for Italy and India before fully reopening social contacts in order to avoid a new pandemic.

Coupling the Within-Host Process and Between-Host Transmission of COVID-19 Suggests Vaccination and School Closures are Critical.

Most models of COVID-19 are implemented at a single micro or macro scale, ignoring the interplay between immune response, viral dynamics, individual infectiousness and epidemiological contact networks. Here we develop a data-driven model linking the within-host viral dynamics to the between-host transmission dynamics on a multilayer contact network to investigate the potential factors driving transmission dynamics and to inform how school closures and antiviral treatment can influence the epidemic. Using multi-source data, we initially determine the viral dynamics and estimate the relationship between viral load and infectiousness. Then, we embed the viral dynamics model into a four-layer contact network and formulate an agent-based model to simulate between-host transmission. The results illustrate that the heterogeneity of immune response between children and adults and between vaccinated and unvaccinated infections can produce different transmission patterns. We find that school closures play a significant effect on mitigating the pandemic as more adults get vaccinated and the virus mutates. If enough infected individuals are diagnosed by testing before symptom onset and then treated quickly, the transmission can be effectively curbed. Our multiscale model reveals the critical role played by younger individuals and antiviral treatment with testing in controlling the epidemic.

Statistical analysis of one-compartment pharmacokinetic models with drug adherence.

Pharmacokinetics is a scientific branch of pharmacology that describes the time course of drug concentration within a living organism and helps the scientific decision-making of potential drug candidates. However, the classical pharmacokinetic models with the eliminations of zero-order, first-order and saturated Michaelis-Menten processes, assume that patients perfectly follow drug regimens during drug treatment, and the significant factor of patients' drug adherence is not taken into account. In this study, therefore, considering the random change of dosage at the fixed dosing time interval, we reformulate the classical deterministic one-compartment pharmacokinetic models to the framework of stochastic, and analyze their qualitative properties including the expectation and variance of the drug concentration, existence of limit drug distribution, and the stochastic properties such as transience and recurrence. In addition, we carry out sensitivity analysis of drug adherence-related parameters to the key values like expectation and variance, especially for the impact on the lowest and highest steady state drug concentrations (i.e. the therapeutic window). Our findings can provide an important theoretical guidance for the variability of drug concentration and help the optimal design of medication regimens. Moreover, The developed models in this paper can support for the potential study of the impact of drug adherence on long-term treatment for chronic diseases like HIV, by integrating disease models and the stochastic PK models.

Joint impacts of media, vaccination and treatment on an epidemic Filippov model with application to COVID-19.

A non-smooth SIR Filippov system is proposed to investigate the impacts of three control strategies (media coverage, vaccination and treatment) on the spread of an infectious disease. We synthetically consider both the number of infected population and its changing rate as the switching condition to implement the curing measures. By using the properties of the Lambert W function, we convert the proposed switching condition to a threshold value related to the susceptible population. The classical epidemic model involving media coverage, linear functions describing injecting vaccine and treatment strategies is examined when the susceptible population exceeds the threshold value. In addition, we consider another SIR model accompanied with the vaccination and treatment strategies represented by saturation functions when the susceptible population is smaller than the threshold value. The dynamics of these two subsystems and the sliding domain are discussed in detail. Four types of local sliding bifurcation are investigated, including boundary focus, boundary node, boundary saddle and boundary saddle-node bifurcations. In the meantime, the global bifurcation involving the appearance of limit cycles is examined, including touching bifurcation, homoclinic bifurcation to the pseudo-saddle and crossing bifurcation. Furthermore, the influence of some key parameters related to the three treatment strategies is explored. We also validate our model by the epidemic data sets of A/H1N1 and COVID-19, which can be employed to reveal the effects of media report and existing strategy related to the control of emerging infectious diseases on the variations of confirmed cases.

Effects of medical resource capacities and intensities of public mitigation measures on outcomes of COVID-19 outbreaks.

BACKGROUND: The COVID-19 pandemic is complex and is developing in different ways according to the country involved. METHODS: To identify the key parameters or processes that have the greatest effects on the pandemic and reveal the different progressions of epidemics in different countries, we quantified enhanced control measures and the dynamics of the production and provision of medical resources. We then nested these within a COVID-19 epidemic transmission model, which is parameterized by multi-source data. We obtained rate functions related to the intensity of mitigation measures, the effective reproduction numbers and the timings and durations of runs on medical resources, given differing control measures implemented in various countries. RESULTS: Increased detection rates may induce runs on medical resources and prolong their durations, depending on resource availability. Nevertheless, improving the detection rate can effectively and rapidly reduce the mortality rate, even after runs on medical resources. Combinations of multiple prevention and control strategies and timely improvement of abilities to supplement medical resources are key to effective control of the COVID-19 epidemic. A 50% reduction in comprehensive control measures would have led to the cumulative numbers of confirmed cases and deaths exceeding 590,000 and 60,000, respectively, by 27 March 2020 in mainland China. CONCLUSIONS: Multiple data sources and cross validation of a COVID-19 epidemic model, coupled with a medical resource logistic model, revealed the key factors that affect epidemic progressions and their outbreak patterns in different countries. These key factors are the type of emergency medical response to avoid runs on medical resources, especially improved detection rates, the ability to promote public health measures, and the synergistic effects of combinations of multiple prevention and control strategies. The proposed model can assist health authorities to predict when they will be most in need of hospital beds and equipment such as ventilators, personal protection equipment, drugs, and staff.

Impacts of varying strengths of intervention measures on secondary outbreaks of COVID-19 in two different regions.

By March 2020, China and Singapore had achieved remarkable results in the prevention and control of COVID-19, but in April Singapore's outbreak began to deteriorate, while China's remained controlled. Using detailed data from Tianjin, China, and Singapore, a stochastic discrete COVID-19 epidemic model was constructed to depict the impact of the epidemic. Parameter estimation and sensitivity analysis were developed to study the probability of imported cases inducing an outbreak in relation to different prevention and control efforts. Results show that the resumption of work and the re-opening of schools will not lead to an outbreak if the effective reproduction number is lower than 1 and approaches 0 and tracking quarantine measures are strengthened. Once an outbreak occurs, if close contacts can be tracked and quarantined in time, the outbreak will be contained. If work is resumed and schools are re-opened with the effective reproduction number greater than 1, then it is more likely that a secondary outbreak will be generated. Also, the greater the number of undetected foreign imported cases and the weaker the prevention and control measures, the more serious the epidemic. Therefore, the key to prevention of a second outbreak is to return to work and to re-open schools only after the effective reproduction number is less than 1 for a period, and when tracking quarantine measures have been strengthened. Our model provides a qualitative and quantitative basis for decision-making for the prevention and control of COVID-19 epidemics and the prediction, early warning and risk assessment of secondary outbreaks.

Impact of media reports on the early spread of COVID-19 epidemic.

Media reports can modify people's knowledge of emerging infectious diseases, and thus changing the public attitudes and behaviors. However, how the media reports affect the development of COVID-19 epidemic is a key public health issue. Here the Pearson correlation and cross-correlation analyses are conducted to find the statistically significant correlations between the number of new hospital notifications for COVID-19 and the number of daily news items for twelve major websites in China from January 11th to February 6th 2020. To examine the implication for transmission dynamics of these correlations, we proposed a novel model, which embeds the function of individual behaviour change (media impact) into the intensity of infection. The nonlinear least squares estimation is used to identify the best-fit parameter values in the model from the observed data. To determine impact of key parameters with media impact and control measures for the later outcome of the outbreak, we also carried out the uncertainty and sensitivity analyses. These findings confirm the importance of the responses of individuals to the media reports, and the crucial role of experts and governments in promoting the public under self-quarantine. Therefore, for mitigating epidemic COVID-19, the media publicity should be focused on how to guide people's behavioral changes by experts, and the management departments and designated hospitals of the COVID-19 should take effective quarantined measures, which are critical for the control of the disease.