Front dynamics in fractional-order epidemic models

A number of recent studies suggest that human and animal mobility patterns exhibit scale-free, Lévy-flight dynamics. However, current reaction-diffusion epidemics models do not account for the superdiffusive spread of modern epidemics due to Lévy flights. We have developed a SIR model to simulate the spatial spread of a hypothetical epidemic driven by long-range displacements in the infective and susceptible populations. The model has been obtained by replacing the second-order diffusion operator by a fractional-order operator. Theoretical developments and numerical simulations show that fractional-order diffusion leads to an exponential acceleration of the epidemic's front and a power-law decay of the front's leading tail. Our results indicate the potential of fractional-order reaction-diffusion models to represent modern epidemics.     

Inference for epidemic models with time‐varying infection rates: Tracking the dynamics of oak processionary moth in the UK

1. Invasive pests pose a great threat to forest, woodland, and urban tree ecosystems. The oak processionary moth (OPM) is a destructive pest of oak trees, first reported in the UK in 2006. Despite great efforts to contain the outbreak within the original infested area of South‐East England, OPM continues to spread.
2. Here, we analyze data consisting of the numbers of OPM nests removed each year from two parks in London between 2013 and 2020. Using a state‐of‐the‐art Bayesian inference scheme, we estimate the parameters for a stochastic compartmental SIR (susceptible, infested, and removed) model with a time‐varying infestation rate to describe the spread of OPM.
3. We find that the infestation rate and subsequent basic reproduction number have remained constant since 2013 (with R0 between one and two). This shows further controls must be taken to reduce R0 below one and stop the advance of OPM into other areas of England.
4. Synthesis. Our findings demonstrate the applicability of the SIR model to describing OPM spread and show that further controls are needed to reduce the infestation rate. The proposed statistical methodology is a powerful tool to explore the nature of a time‐varying infestation rate, applicable to other partially observed time series epidemic data.

     

Backward bifurcations and multiple equilibria in epidemic models with structured immunity

Many disease pathogens stimulate immunity in their hosts, which then wanes over time. To better understand the impact of this immunity on epidemiological dynamics, we propose an epidemic model structured according to immunity level that can be applied in many different settings. Under biologically realistic hypotheses, we find that immunity alone never creates a backward bifurcation of the disease-free steady state. This does not rule out the possibility of multiple stable equilibria, but we provide two sufficient conditions for the uniqueness of the endemic equilibrium, and show that these conditions ensure uniqueness in several common special cases. Our results indicate that the within-host dynamics of immunity can, in principle, have important consequences for population-level dynamics, but also suggest that this would require strong non-monotone effects in the immune response to infection. Neutralizing antibody titer data for measles are used to demonstrate the biological application of our theory.     

Functional mixed effects models

In this article, a new class of functional models in which smoothing splines are used to model fixed effects as well as random effects is introduced. The linear mixed effects models are extended to nonparametric mixed effects models by introducing functional random effects, which are modeled as realizations of zero-mean stochastic processes. The fixed functional effects and the random functional effects are modeled in the same functional space, which guarantee the population-average and subject-specific curves have the same smoothness property. These models inherit the flexibility of the linear mixed effects models in handling complex designs and correlation structures, can include continuous covariates as well as dummy factors in both the fixed or random design matrices, and include the nested curves models as special cases. Two estimation procedures are proposed. The first estimation procedure exploits the connection between linear mixed effects models and smoothing splines and can be fitted using existing software. The second procedure is a sequential estimation procedure using Kalman filtering. This algorithm avoids inversion of large dimensional matrices and therefore can be applied to large data sets. A generalized maximum likelihood (GML) ratio test is proposed for inference and model selection. An application to comparison of cortisol profiles is used as an illustration.     

EM-REML estimation of covariance parameters in Gaussian mixed models for longitudinal data analysis

This paper presents procedures for implementing the EM algorithm to compute REML estimates of variance covariance components in Gaussian mixed models for longitudinal data analysis. The class of models considered includes random coefficient factors, stationary time processes and measurement errors. The EM algorithm allows separation of the computations pertaining to parameters involved in the random coefficient factors from those pertaining to the time processes and errors. The procedures are illustrated with Pothoff and Roy''s data example on growth measurements taken on 11 girls and 16 boys at four ages. Several variants and extensions are discussed.     

Coexistence of multiple pathogen strains in stochastic epidemic models with density-dependent mortality

Stochastic differential equations that model an SIS epidemic with multiple pathogen strains are derived from a system of ordinary differential equations. The stochastic model assumes there is demographic variability. The dynamics of the deterministic model are summarized. Then the dynamics of the stochastic model are compared to the deterministic model. In the deterministic model, there can be either disease extinction, competitive exclusion, where only one strain persists, or coexistence, where more than one strain persists. In the stochastic model, all strains are eventually eliminated because the disease-free state is an absorbing state. However, if the population size and the initial number of infected individuals are sufficiently large, it may take a long time until all strains are eliminated. Numerical simulations of the stochastic model show that coexistence cases predicted by the deterministic model are an unlikely occurrence in the stochastic model even for short time periods. In the stochastic model, either disease extinction or competitive exclusion occur. The initial number of infected individuals, the basic reproduction numbers, and other epidemiological parameters are important determinants of the dominant strain in the stochastic epidemic model.     

Permutation tests for random effects in linear mixed models

Inference regarding the inclusion or exclusion of random effects in linear mixed models is challenging because the variance components are located on the boundary of their parameter space under the usual null hypothesis. As a result, the asymptotic null distribution of the Wald, score, and likelihood ratio tests will not have the typical χ(2) distribution. Although it has been proved that the correct asymptotic distribution is a mixture of χ(2) distributions, the appropriate mixture distribution is rather cumbersome and nonintuitive when the null and alternative hypotheses differ by more than one random effect. As alternatives, we present two permutation tests, one that is based on the best linear unbiased predictors and one that is based on the restricted likelihood ratio test statistic. Both methods involve weighted residuals, with the weights determined by the among- and within-subject variance components. The null permutation distributions of our statistics are computed by permuting the residuals both within and among subjects and are valid both asymptotically and in small samples. We examine the size and power of our tests via simulation under a variety of settings and apply our test to a published data set of chronic myelogenous leukemia patients.     

Inference for an epidemic when susceptibility varies

A stochastic epidemic model featuring fixed-length latent periods, gamma-distributed infectious periods and randomly varying heterogeneity among susceptibles is considered. A Markov chain Monte Carlo algorithm is developed for performing Bayesian inference for the parameters governing the infectious-period length and the hyper-parameters governing the heterogeneity of susceptibility. This method of analysis applies to a wider class of diseases than methods proposed previously. An application to smallpox data confirms results about heterogeneity suggested by an earlier analysis that relied on less realistic assumptions.     

Inference in spline-based models for multiple time-to-event data, with applications to a breast cancer prevention trial

As part of the National Surgical Adjuvant Breast and Bowel Project, a controlled clinical trial known as the Breast Cancer Prevention Trial (BCPT) was conducted to assess the effectiveness of tamoxifen as a preventive agent for breast cancer. In addition to the incidence of breast cancer, data were collected on several other, possibly adverse, outcomes, such as invasive endometrial cancer, ischemic heart disease, transient ischemic attack, deep vein thrombosis and/or pulmonary embolism. In this article, we present results from an illustrative analysis of the BCPT data, based on a new modeling technique, to assess the effectiveness of the drug tamoxifen as a preventive agent for breast cancer. We extended the flexible model of Gray (1994, Spline-based test in survival analysis, Biometrics 50, 640-652) to allow inference on multiple time-to-event outcomes in the style of the marginal modeling setup of Wei, Lin, and Weissfeld (1989, Regression analysis of multivariate incomplete failure time data by modeling marginal distributions, Journal of the American Statistical Association 84, 1065-1073). This proposed model makes inference possible for multiple time-to-event data while allowing for greater flexibility in modeling the effects of prognostic factors with nonlinear exposure-response relationships. Results from simulation studies on the small-sample properties of the asymptotic tests will also be presented.     

Maximin D-optimal designs for longitudinal mixed effects models

In this article, the optimal selection and allocation of time points in repeated measures experiments is considered. D-optimal cohort designs are computed numerically for the first- and second-degree polynomial models with random intercept, random slope, and first-order autoregressive serial correlations. Because the optimal designs are locally optimal, it is proposed to use a maximin criterion. It is shown that, for a large class of symmetric designs, the smallest relative efficiency over the model parameter space is substantial.     

Nonlinear mixed effects models for repeated measures data

We propose a general, nonlinear mixed effects model for repeated measures data and define estimators for its parameters. The proposed estimators are a natural combination of least squares estimators for nonlinear fixed effects models and maximum likelihood (or restricted maximum likelihood) estimators for linear mixed effects models. We implement Newton-Raphson estimation using previously developed computational methods for nonlinear fixed effects models and for linear mixed effects models. Two examples are presented and the connections between this work and recent work on generalized linear mixed effects models are discussed.     

Evaluating somatic cell scores with a Bayesian Gaussian linear state-space model

Because accurate characterization of health state is important for managing dairy herds, we propose to validate the use of a linear state-space model (LSSM) for evaluating monthly somatic cell scores (SCSs). To do so, we retrieved SCS from a dairy database and collected reports on clinical mastitis collected in 20 farms, during the period from January 2008 to December 2011 in the Walloon region of Belgium. The dependent variable was the SCS, and the independent variables were the number of days from calving, year of calving and parity. The LSSM also incorporated an error-free underlying variable that described the trend across time as a function of previous clinical and subclinical status. We computed the mean sum of squared differences between observed SCS and median values of the posterior SCS distribution and constructed the receiver operating characteristic (ROC) curve for SCS thresholds going from 0 to 6. Our results show SCS estimates are close to observed SCS and area under the ROC curve is higher than 90%. We discuss the meaning of the parameters in light of our current knowledge of the disease and propose methods to incorporate, in LSSM, this knowledge often expressed in the form of ordinary differential equations.     

Generalized linear mixed models for mapping multiple quantitative trait loci

Many biological traits are discretely distributed in phenotype but continuously distributed in genetics because they are controlled by multiple genes and environmental variants. Due to the quantitative nature of the genetic background, these multiple genes are called quantitative trait loci (QTL). When the QTL effects are treated as random, they can be estimated in a single generalized linear mixed model (GLMM), even if the number of QTL may be larger than the sample size. The GLMM in its original form cannot be applied to QTL mapping for discrete traits if there are missing genotypes. We examined two alternative missing genotype-handling methods: the expectation method and the overdispersion method. Simulation studies show that the two methods are efficient for multiple QTL mapping (MQM) under the GLMM framework. The overdispersion method showed slight advantages over the expectation method in terms of smaller mean-squared errors of the estimated QTL effects. The two methods of GLMM were applied to MQM for the female fertility trait of wheat. Multiple QTL were detected to control the variation of the number of seeded spikelets.     

Random effects selection in linear mixed models

We address the important practical problem of how to select the random effects component in a linear mixed model. A hierarchical Bayesian model is used to identify any random effect with zero variance. The proposed approach reparameterizes the mixed model so that functions of the covariance parameters of the random effects distribution are incorporated as regression coefficients on standard normal latent variables. We allow random effects to effectively drop out of the model by choosing mixture priors with point mass at zero for the random effects variances. Due to the reparameterization, the model enjoys a conditionally linear structure that facilitates the use of normal conjugate priors. We demonstrate that posterior computation can proceed via a simple and efficient Markov chain Monte Carlo algorithm. The methods are illustrated using simulated data and real data from a study relating prenatal exposure to polychlorinated biphenyls and psychomotor development of children.