Search

Purpose: To describe a procedure for incorporating parametric functions into individual-level simulation models to sample time to event when age-specific rates are available but not the individual data. Methods: Using age-specific event rates, regression analysis was used to parametrize parametric survival distributions (Weibull, Gompertz, log-normal, and log-logistic), select the best fit using the R 2 statistic, and apply the corresponding formula to assign random times to events in simulation models. We used stroke rates in the Spanish population to illustrate our procedure. Results: The 3 selected survival functions (Gompertz, Weibull, and log-normal) had a good fit to the data up to 85 y of age. We selected Gompertz distribution as the best-fitting distribution due to its goodness of fit. Conclusions: Our work provides a simple procedure for incorporating parametric risk functions into simulation models without individual-level data. Highlights: We describe the procedure for sampling times to event for individual-level simulation models as a function of age from parametric survival functions when age-specific rates are available but not the individual data We used linear regression to estimate age-specific hazard functions, obtaining estimates of parameter uncertainty. Our approach allows incorporating parameter (second-order) uncertainty in individual-level simulation models needed for probabilistic sensitivity analysis in the absence of individual-level survival data. [ABSTRACT FROM AUTHOR]

Background: Compartmental infectious disease (ID) models are often used to evaluate nonpharmaceutical interventions (NPIs) and vaccines. Such models rarely separate within-household and community transmission, potentially introducing biases in situations in which multiple transmission routes exist. We formulated an approach that incorporates household structure into ID models, extending the work of House and Keeling. Design: We developed a multicompartment susceptible-exposed-infectious-recovered-susceptible-vaccinated (MC-SEIRSV) modeling framework, allowing nonexponentially distributed duration in exposed and infectious compartments, that tracks within-household and community transmission. We simulated epidemics that varied by community and household transmission rates, waning immunity rate, household size (3 or 5 members), and numbers of exposed and infectious compartments (1–3 each). We calibrated otherwise identical models without household structure to the early phase of each parameter combination's epidemic curve. We compared each model pair in terms of epidemic forecasts and predicted NPI and vaccine impacts on the timing and magnitude of the epidemic peak and its total size. Meta-analytic regressions characterized the relationship between household structure inclusion and the size and direction of biases. Results: Otherwise similar models with and without household structure produced equivalent early epidemic curves. However, forecasts from models without household structure were biased. Without intervention, they were upward biased on peak size and total epidemic size, with biases also depending on the number of exposed and infectious compartments. Model-estimated NPI effects of a 60% reduction in community contacts on peak time and size were systematically overestimated without household structure. Biases were smaller with a 20% reduction NPI. Because vaccination affected both community and household transmission, their biases were smaller. Conclusions: ID models without household structure can produce biased outcomes in settings in which within-household and community transmission differ. Highlights: Infectious disease models rarely separate household transmission from community transmission. The pace of household transmission may differ from community transmission, depends on household size, and can accelerate epidemic growth. Many infectious disease models assume exponential duration distributions for infected states. However, the duration of most infections is not exponentially distributed, and distributional choice alters modeled epidemic dynamics and intervention effectiveness. We propose a mathematical framework for household and community transmission that allows for nonexponential duration times and a suite of interventions and quantified the effect of accounting for household transmission by varying household size and duration distributions of infected states on modeled epidemic dynamics. Failure to include household structure induces biases in the modeled overall course of an epidemic and the effects of interventions delivered differentially in community settings. Epidemic dynamics are faster and more intense in populations with larger household sizes and for diseases with nonexponentially distributed infectious durations. Modelers should consider explicitly incorporating household structure to quantify the effects of non-pharmaceutical interventions (e.g., shelter-in-place). [ABSTRACT FROM AUTHOR]

In an introductory tutorial, we illustrated building cohort state-transition models (cSTMs) in R, where the state transition probabilities were constant over time. However, in practice, many cSTMs require transitions, rewards, or both to vary over time (time dependent). This tutorial illustrates adding 2 types of time dependence using a previously published cost-effectiveness analysis of multiple strategies as an example. The first is simulation-time dependence, which allows for the transition probabilities to vary as a function of time as measured since the start of the simulation (e.g., varying probability of death as the cohort ages). The second is state-residence time dependence, allowing for history by tracking the time spent in any particular health state using tunnel states. We use these time-dependent cSTMs to conduct cost-effectiveness and probabilistic sensitivity analyses. We also obtain various epidemiological outcomes of interest from the outputs generated from the cSTM, such as survival probability and disease prevalence, often used for model calibration and validation. We present the mathematical notation first, followed by the R code to execute the calculations. The full R code is provided in a public code repository for broader implementation. [ABSTRACT FROM AUTHOR]

Background: Historically, correctional facilities have had large outbreaks of respiratory infectious diseases like COVID-19. Hence, importation and exportation of such diseases from correctional facilities raises substantial concern. Methods: We developed a stochastic simulation model of transmission of respiratory infectious diseases within and between correctional facilities and the community. We investigated the infection dynamics, key governing factors, and relative importance of different infection routes (e.g., incarcerations and releases versus correctional staff). We also developed machine-learning meta-models of the simulation model, which allowed us to examine how our findings depended on different disease, correctional facility, and community characteristics. Results: We find a magnification-reflection dynamic: a small outbreak in the community can cause a larger outbreak in the correction facility, which can then cause a second, larger outbreak in the community. This dynamic is strongest when community size is relatively small as compared with the size of the correctional population, the initial community R-effective is near 1, and initial prevalence of immunity in the correctional population is low. The timing of the correctional magnification and community reflection peaks in infection prevalence are primarily governed by the initial R-effective for each setting. Because the release rates from prisons are low, our model suggests correctional staff may be a more important infection entry route into prisons than incarcerations and releases; in jails, where incarceration and release rates are much higher, our model suggests the opposite. Conclusions: We find that across many combinations of respiratory pathogens, correctional settings, and communities, there can be substantial magnification-reflection dynamics, which are governed by several key factors. Our goal was to derive theoretical insights relevant to many contexts; our findings should be interpreted accordingly. Highlights: We find a magnification-reflection dynamic: a small outbreak in a community can cause a larger outbreak in a correctional facility, which can then cause a second, larger outbreak in the community. For public health decision makers considering contexts most susceptible to this dynamic, we find that the dynamic is strongest when the community size is relatively small, initial community R-effective is near 1, and the initial prevalence of immunity in the correctional population is low; the timing of the correctional magnification and community reflection peaks in infection prevalence are primarily governed by the initial R-effective for each setting. We find that correctional staff may be a more important infection entry route into prisons than incarcerations and releases; however, for jails, the relative importance of the entry routes may be reversed. For modelers, we combine simulation modeling, machine-learning meta-modeling, and interpretable machine learning to examine how our findings depend on different disease, correctional facility, and community characteristics; we find they are generally robust. [ABSTRACT FROM AUTHOR]

Decision models can combine information from different sources to simulate the long-term consequences of alternative strategies in the presence of uncertainty. A cohort state-transition model (cSTM) is a decision model commonly used in medical decision making to simulate the transitions of a hypothetical cohort among various health states over time. This tutorial focuses on time-independent cSTM, in which transition probabilities among health states remain constant over time. We implement time-independent cSTM in R, an open-source mathematical and statistical programming language. We illustrate time-independent cSTMs using a previously published decision model, calculate costs and effectiveness outcomes, and conduct a cost-effectiveness analysis of multiple strategies, including a probabilistic sensitivity analysis. We provide open-source code in R to facilitate wider adoption. In a second, more advanced tutorial, we illustrate time-dependent cSTMs. [ABSTRACT FROM AUTHOR]

Most decisions are associated with uncertainty. Value of information (VOI) analysis quantifies the opportunity loss associated with choosing a suboptimal intervention based on current imperfect information. VOI can inform the value of collecting additional information, resource allocation, research prioritization, and future research designs. However, in practice, VOI remains underused due to many conceptual and computational challenges associated with its application. Expected value of sample information (EVSI) is rooted in Bayesian statistical decision theory and measures the value of information from a finite sample. The past few years have witnessed a dramatic growth in computationally efficient methods to calculate EVSI, including metamodeling. However, little research has been done to simplify the experimental data collection step inherent to all EVSI computations, especially for correlated model parameters. This article proposes a general Gaussian approximation (GA) of the traditional Bayesian updating approach based on the original work by Raiffa and Schlaifer to compute EVSI. The proposed approach uses a single probabilistic sensitivity analysis (PSA) data set and involves 2 steps: 1) a linear metamodel step to compute the EVSI on the preposterior distributions and 2) a GA step to compute the preposterior distribution of the parameters of interest. The proposed approach is efficient and can be applied for a wide range of data collection designs involving multiple non-Gaussian parameters and unbalanced study designs. Our approach is particularly useful when the parameters of an economic evaluation are correlated or interact. [ABSTRACT FROM AUTHOR]