Fast Bayesian parameter estimation for stochastic logistic growth models

Application of logistic differential equation models for early warning of infectious diseases in Jilin Province

Background

There is still a relatively serious disease burden of infectious diseases and the warning time for different infectious diseases before implementation of interventions is important. The logistic differential equation models can be used for predicting early warning of infectious diseases. The aim of this study is to compare the disease fitting effects of the logistic differential equation (LDE) model and the generalized logistic differential equation (GLDE) model for the first time using data on multiple infectious diseases in Jilin Province and to calculate the early warning signals for different types of infectious diseases using these two models in Jilin Province to solve the disease early warning schedule for Jilin Province throughout the year.

Methods

Collecting the incidence of 22 infectious diseases in Jilin Province, China. The LDE and GLDE models were used to calculate the recommended warning week (RWW), the epidemic acceleration week (EAW) and warning removed week (WRW) for acute infectious diseases with seasonality, respectively.

Results

Five diseases were selected for analysis based on screening principles: hemorrhagic fever with renal syndrome (HFRS), shigellosis, mumps, Hand, foot and mouth disease (HFMD), and scarlet fever. The GLDE model fitted the above diseases better (0.80 ≤ R2 ≤ 0.94, P <  0. 005) than the LDE model. The estimated warning durations (per year) of the LDE model for the above diseases were: weeks 12–23 and 40–50; weeks 20–36; weeks 15–24 and 43–52; weeks 26–34; and weeks 16–25 and 41–50. While the durations of early warning (per year) estimated by the GLDE model were: weeks 7–24 and 36–51; weeks 13–37; weeks 11–26 and 39–54; weeks 23–35; and weeks 12–26 and 40–50.

Conclusions

Compared to the LDE model, the GLDE model provides a better fit to the actual disease incidence data. The RWW appeared to be earlier when estimated with the GLDE model than the LDE model. In addition, the WRW estimated with the GLDE model were more lagged and had a longer warning time.

Efficient Bayesian inference of fully stochastic epidemiological models with applications to COVID-19

Epidemiological forecasts are beset by uncertainties about the underlying epidemiological processes, and the surveillance process through which data are acquired. We present a Bayesian inference methodology that quantifies these uncertainties, for epidemics that are modelled by (possibly) non-stationary, continuous-time, Markov population processes. The efficiency of the method derives from a functional central limit theorem approximation of the likelihood, valid for large populations. We demonstrate the methodology by analysing the early stages of the COVID-19 pandemic in the UK, based on age-structured data for the number of deaths. This includes maximum a posteriori estimates, Markov chain Monte Carlo sampling of the posterior, computation of the model evidence, and the determination of parameter sensitivities via the Fisher information matrix. Our methodology is implemented in PyRoss, an open-source platform for analysis of epidemiological compartment models.

Optimal Sampling Regimes for Estimating Population Dynamics

Ecologists are interested in modeling the population growth of species in various ecosystems. Specifically, logistic growth arises as a common model for population growth. Studying such growth can assist environmental managers in making better decisions when collecting data. Traditionally, ecological data is recorded on a regular time frequency and is very well-documented. However, sampling can be an expensive process due to available resources, money and time. Limiting sampling makes it challenging to properly track the growth of a population. Thus, this design study proposes an approach to sampling based on the dynamics associated with logistic growth. The proposed method is demonstrated via a simulation study across various theoretical scenarios to evaluate its performance in identifying optimal designs that best estimate the curves. Markov Chain Monte Carlo sampling techniques are implemented to predict the probability of the model parameters using Bayesian inference. The intention of this study is to demonstrate a method that can minimize the amount of time ecologists spend in the field, while maximizing the information provided by the data.

Early warning of hand, foot, and mouth disease transmission: A modeling study in mainland, China

BACKGROUND

Hand, foot, and mouth disease (HFMD) is a global infectious disease; particularly, it has a high disease burden in China. This study was aimed to explore the temporal and spatial distribution of the disease by analyzing its epidemiological characteristics, and to calculate the early warning signals of HFMD by using a logistic differential equation (LDE) model.

METHODS

This study included datasets of HFMD cases reported in seven regions in Mainland China. The early warning time (week) was calculated using the LDE model with the key parameters estimated by fitting with the data. Two key time points, "epidemic acceleration week (EAW)" and "recommended warning week (RWW)", were calculated to show the early warning time.

RESULTS

The mean annual incidence of HFMD cases per 100,000 per year was 218, 360, 223, 124, and 359 in Hunan Province, Shenzhen City, Xiamen City, Chuxiong Prefecture, Yunxiao County across the southern regions, respectively and 60 and 34 in Jilin Province and Longde County across the northern regions, respectively. The LDE model fitted well with the reported data (R2 > 0.65, P < 0.001). Distinct temporal patterns were found across geographical regions: two early warning signals emerged in spring and autumn every year across southern regions while one early warning signals in summer every year across northern regions.

CONCLUSIONS

The disease burden of HFMD in China is still high, with more cases occurring in the southern regions. The early warning of HFMD across the seven regions is heterogeneous. In the northern regions, it has a high incidence during summer and peaks in June every year; in the southern regions, it has two waves every year with the first wave during spring spreading faster than the second wave during autumn. Our findings can help predict and prepare for active periods of HFMD.

Predicting the Future Course of Opioid Overdose Mortality: An Example From Two US States

BACKGROUND

The rapid growth of opioid abuse and the related mortality across the United States has spurred the development of predictive models for the allocation of public health resources. These models should characterize heterogeneous growth across states using a drug epidemic framework that enables assessments of epidemic onset, rates of growth, and limited capacities for epidemic growth.

METHODS

We used opioid overdose mortality data for 146 North and South Carolina counties from 2001 through 2014 to compare the retrodictive and predictive performance of a logistic growth model that parameterizes onsets, growth, and carrying capacity within a traditional Bayesian Poisson space-time model.

RESULTS

In fitting the models to past data, the performance of the logistic growth model was superior to the standard Bayesian Poisson space-time model (deviance information criterion: 8,088 vs. 8,256), with reduced spatial and independent errors. Predictively, the logistic model more accurately estimated fatality rates 1, 2, and 3 years in the future (root mean squared error medians were lower for 95.7% of counties from 2012 to 2014). Capacity limits were higher in counties with greater population size, percent population age 45-64, and percent white population. Epidemic onset was associated with greater same-year and past-year incidence of overdose hospitalizations.

CONCLUSION

Growth in annual rates of opioid fatalities was capacity limited, heterogeneous across counties, and spatially correlated, requiring spatial epidemic models for the accurate and reliable prediction of future outcomes related to opioid abuse. Indicators of risk are identifiable and can be used to predict future mortality outcomes.

Inferring time-dependent population growth rates in cell cultures undergoing adaptation

Background

The population growth rate is an important characteristic of any cell culture. During sustained experiments, the growth rate may vary due to competition or adaptation. For instance, in presence of a toxin or a drug, an increasing growth rate indicates that the cells adapt and become resistant. Consequently, time-dependent growth rates are fundamental to follow on the adaptation of cells to a changing evolutionary landscape. However, as there are no tools to calculate the time-dependent growth rate directly by cell counting, it is common to use only end point measurements of growth rather than tracking the growth rate continuously.

Results

We present a computer program for inferring the growth rate over time in suspension cells using nothing but cell counts, which can be measured non-destructively. The program was tested on simulated and experimental data. Changes were observed in the initial and absolute growth rates, betraying resistance and adaptation.

Conclusions

For experiments where adaptation is expected to occur over a longer time, our method provides a means of tracking growth rates using data that is normally collected anyhow for monitoring purposes. The program and its documentation are freely available at https://github.com/Sandalmoth/ratrack under the permissive zlib license.

Mean exit time and escape probability for the stochastic logistic growth model with multiplicative α-stable Lévy noise

In this paper, we formulate a stochastic logistic fish growth model driven by both white noise and non-Gaussian noise. We focus our study on the mean time to extinction, escape probability to measure the noise-induced extinction probability and the Fokker–Planck equation for fish population [Formula: see text]. In the Gaussian case, these quantities satisfy local partial differential equations while in the non-Gaussian case, they satisfy nonlocal partial differential equations. Following a discussion of existence, uniqueness and stability, we calculate numerical approximations of the solutions of those equations. For each noise model we then compare the behaviors of the mean time to extinction and the solution of the Fokker–Planck equation as growth rate [Formula: see text], carrying capacity [Formula: see text], intensity of Gaussian noise [Formula: see text], noise intensity [Formula: see text] and stability index [Formula: see text] vary. The MET from the interval [Formula: see text] at the right boundary is finite if [Formula: see text]. For [Formula: see text], the MET from [Formula: see text] at this boundary is infinite. A larger stability index [Formula: see text] is less likely leading to the extinction of the fish population.

The study on the early warning period of varicella outbreaks based on logistic differential equation model

Chickenpox is a common acute and highly contagious disease in childhood; moreover, there is currently no targeted treatment. Carrying out an early warning on chickenpox plays an important role in taking targeted measures in advance as well as preventing the outbreak of the disease. In recent years, the infectious disease dynamic model has been widely used in the research of various infectious diseases. The logistic differential equation model can well demonstrate the epidemic characteristics of epidemic outbreaks, gives the point at which the early epidemic rate changes from slow to fast. Therefore, our study aims to use the logistic differential equation model to explore the epidemic characteristics and early-warning time of varicella. Meanwhile, the data of varicella cases were collected from first week of 2008 to 52nd week of 2017 in Changsha. Finally, our study found that the logistic model can be well fitted with varicella data, besides the model illustrated that there are two peaks of varicella at each year in Changsha City. One is the peak in summer–autumn corresponding to the 8th–38th week; the other is in winter–spring corresponding to the time from the 38th to the seventh week next year. The ‘epidemic acceleration week’ average value of summer–autumn and winter–spring are about the 16th week (ranging from the 15th to 17th week) and 45th week (ranging from the 44th to 47th week), respectively. What is more, taking warning measures during the acceleration week, the preventive effect will be delayed; thus, we recommend intervene during recommended warning weeks which are the 15th and 44th weeks instead.

A Bayesian estimation approach for the mortality in a stage-structured demographic model

Control interventions in sustainable pest management schemes are set according to the phenology and the population abundance of the pests. This information can be obtained using suitable mathematical models that describe the population dynamics based on individual life history responses to environmental conditions and resource availability. These responses are described by development, fecundity and survival rate functions, which can be estimated from laboratory experiments. If experimental data are not available, data on field population dynamics can be used for their estimation. This is the case of the extrinsic mortality term that appears in the mortality rate function due to biotic factors. We propose a Bayesian approach to estimate the probability density functions of the parameters in the extrinsic mortality rate function, starting from data on population abundance. The method investigates the time variability in the mortality parameters by comparing simulated and observed trajectories. The grape berry moth, a pest of great importance in European vineyards, has been considered as a case study. Simulated data have been considered to evaluate the convergence of the algorithm, while field data have been used to obtain estimates of the mortality for the grape berry moth.

Subspace identification for a stochastic model of plague

In this paper, a stochastic model of plague is first studied by subspace identification. First, the discrete model of plague is obtained based on the classical model. The corresponding stochastic model is proposed for the existence of stochastic disturbances. Second, for the model, the parameter matrices and noise intensity are obtained. Finally, the simulations of the model show that the subspace identification is more precise than least square method.

Automated inference procedure for the determination of cell growth parameters

The growth rate and carrying capacity of a cell population are key to the characterization of the population's viability and to the quantification of its responses to perturbations such as drug treatments. Accurate estimation of these parameters necessitates careful analysis. Here, we present a rigorous mathematical approach for the robust analysis of cell count data, in which all the experimental stages of the cell counting process are investigated in detail with the machinery of Bayesian probability theory. We advance a flexible theoretical framework that permits accurate estimates of the growth parameters of cell populations and of the logical correlations between them. Moreover, our approach naturally produces an objective metric of avoidable experimental error, which may be tracked over time in a laboratory to detect instrumentation failures or lapses in protocol. We apply our method to the analysis of cell count data in the context of a logistic growth model by means of a user-friendly computer program that automates this analysis, and present some samples of its output. Finally, we note that a traditional least squares fit can provide misleading estimates of parameter values, because it ignores available information with regard to the way in which the data have actually been collected.