Bayesian modeling of associations in bivariate piecewise linear mixed-effects models.

Longitudinal processes rarely occur in isolation; often the growth curves of 2 or more variables are interdependent. Moreover, growth curves rarely exhibit a constant pattern of change. Many educational and psychological phenomena are comprised of different developmental phases (segments). Bivariate piecewise linear mixed-effects models (BPLMEM) are a useful and flexible statistical framework that allow simultaneous modeling of 2 processes that portray segmented change and investigates their associations over time. The purpose of the present study was to develop a BPLMEM using a Bayesian inference approach allowing the estimation of the association between the error variances and providing a more robust modeling choice for the joint random-effects of the 2 processes. This study aims to improve upon the limitations of the prior literature on bivariate piecewise mixed-effects models, such as only allowing the modeling of uncorrelated residual errors across the 2 longitudinal processes and restricting modeling choices for the random effects. The performance of the BPLMEM was investigated via a Monte Carlo simulation study. Furthermore, the utility of BPLMEM was illustrated by using a national educational dataset, Early Childhood Longitudinal Study-Kindergarten Cohort (ECLS-K), where we examined the joint development of mathematics and reading achievement scores and the association between their trajectories over 7 measurement occasions. The findings obtained shed new light on the relationship between these 2 prominent educational domains over time. (PsycInfo Database Record (c) 2022 APA, all rights reserved).

Dependence of COVID-19 Policies on End-of-Year Holiday Contacts in Mexico City Metropolitan Area: A Modeling Study.

Background. Mexico City Metropolitan Area (MCMA) has the largest number of COVID-19 (coronavirus disease 2019) cases in Mexico and is at risk of exceeding its hospital capacity in early 2021. Methods. We used the Stanford-CIDE Coronavirus Simulation Model (SC-COSMO), a dynamic transmission model of COVID-19, to evaluate the effect of policies considering increased contacts during the end-of-year holidays, intensification of physical distancing, and school reopening on projected confirmed cases and deaths, hospital demand, and hospital capacity exceedance. Model parameters were derived from primary data, literature, and calibrated. Results. Following high levels of holiday contacts even with no in-person schooling, MCMA will have 0.9 million (95% prediction interval 0.3-1.6) additional COVID-19 cases between December 7, 2020, and March 7, 2021, and hospitalizations will peak at 26,000 (8,300-54,500) on January 25, 2021, with a 97% chance of exceeding COVID-19-specific capacity (9,667 beds). If MCMA were to control holiday contacts, the city could reopen in-person schools, provided they increase physical distancing with 0.5 million (0.2-0.9) additional cases and hospitalizations peaking at 12,000 (3,700-27,000) on January 19, 2021 (60% chance of exceedance). Conclusion. MCMA must increase COVID-19 hospital capacity under all scenarios considered. MCMA's ability to reopen schools in early 2021 depends on sustaining physical distancing and on controlling contacts during the end-of-year holiday.

How do Covid-19 policy options depend on end-of-year holiday contacts in Mexico City Metropolitan Area? A Modeling Study.

BACKGROUND: With more than 20 million residents, Mexico City Metropolitan Area (MCMA) has the largest number of Covid-19 cases in Mexico and is at risk of exceeding its hospital capacity in late December 2020. METHODS: We used SC-COSMO, a dynamic compartmental Covid-19 model, to evaluate scenarios considering combinations of increased contacts during the holiday season, intensification of social distancing, and school reopening. Model parameters were derived from primary data from MCMA, published literature, and calibrated to time-series of incident confirmed cases, deaths, and hospital occupancy. Outcomes included projected confirmed cases and deaths, hospital demand, and magnitude of hospital capacity exceedance. FINDINGS: Following high levels of holiday contacts even with no in-person schooling, we predict that MCMA will have 1·0 million (95% prediction interval 0·5 - 1·7) additional Covid-19 cases between December 7, 2020 and March 7, 2021 and that hospitalizations will peak at 35,000 (14,700 - 67,500) on January 27, 2021, with a >99% chance of exceeding Covid-19-specific capacity (9,667 beds). If holiday contacts can be controlled, MCMA can reopen in-person schools provided social distancing is increased with 0·5 million (0·2 - 1·0) additional cases and hospitalizations peaking at 14,900 (5,600 - 32,000) on January 23, 2021 (77% chance of exceedance). INTERPRETATION: MCMA must substantially increase Covid-19 hospital capacity under all scenarios considered. MCMA's ability to reopen schools in mid-January 2021 depends on sustaining social distancing and that contacts during the end-of-year holiday were well controlled. FUNDING: Society for Medical Decision Making, Gordon and Betty Moore Foundation, and Wadhwani Institute for Artificial Intelligence Foundation. RESEARCH IN CONTEXT: Evidence before this study: As of mid-December 2020, Mexico has the twelfth highest incidence of confirmed cases of Covid-19 worldwide and its epidemic is currently growing. Mexico's case fatality ratio (CFR) - 9·1% - is the second highest in the world. With more than 20 million residents, Mexico City Metropolitan Area (MCMA) has the highest number and incidence rate of Covid-19 confirmed cases in Mexico and a CFR of 8·1%. MCMA is nearing its current hospital capacity even as it faces the prospect of increased social contacts during the 2020 end-of-year holidays. There is limited Mexico-specific evidence available on epidemic, such as parameters governing time-dependent mortality, hospitalization and transmission. Literature searches required supplementation through primary data analysis and model calibration to support the first realistic model-based Covid-19 policy evaluation for Mexico, which makes this analysis relevant and timely.Added value of this study: Study strengths include the use of detailed primary data provided by MCMA; the Bayesian model calibration to enable evaluation of projections and their uncertainty; and consideration of both epidemic and health system outcomes. The model projects that failure to limit social contacts during the end-of-year holidays will substantially accelerate MCMA's epidemic (1·0 million (95% prediction interval 0·5 - 1·7) additional cases by early March 2021). Hospitalization demand could reach 35,000 (14,700 - 67,500), with a >99% chance of exceeding current capacity (9,667 beds). Controlling social contacts during the holidays could enable MCMA to reopen in-person schooling without greatly exacerbating the epidemic provided social distancing in both schools and the community were maintained. Under all scenarios and policies, current hospital capacity appears insufficient, highlighting the need for rapid capacity expansion.Implications of all the available evidence: MCMA officials should prioritize rapid hospital capacity expansion. MCMA's ability to reopen schools in mid-January 2021 depends on sustaining social distancing and that contacts during the end-of-year holiday were well controlled.

Nonidentifiability in Model Calibration and Implications for Medical Decision Making.

BACKGROUND: Calibration is the process of estimating parameters of a mathematical model by matching model outputs to calibration targets. In the presence of nonidentifiability, multiple parameter sets solve the calibration problem, which may have important implications for decision making. We evaluate the implications of nonidentifiability on the optimal strategy and provide methods to check for nonidentifiability. METHODS: We illustrate nonidentifiability by calibrating a 3-state Markov model of cancer relative survival (RS). We performed 2 different calibration exercises: 1) only including RS as a calibration target and 2) adding the ratio between the 2 nondeath states over time as an additional target. We used the Nelder-Mead (NM) algorithm to identify parameter sets that best matched the calibration targets. We used collinearity and likelihood profile analyses to check for nonidentifiability. We then estimated the benefit of a hypothetical treatment in terms of life expectancy gains using different, but equally good-fitting, parameter sets. We also applied collinearity analysis to a realistic model of the natural history of colorectal cancer. RESULTS: When only RS is used as the calibration target, 2 different parameter sets yield similar maximum likelihood values. The high collinearity index and the bimodal likelihood profile on both parameters demonstrated the presence of nonidentifiability. These different, equally good-fitting parameter sets produce different estimates of the treatment effectiveness (0.67 v. 0.31 years), which could influence the optimal decision. By incorporating the additional target, the model becomes identifiable with a collinearity index of 3.5 and a unimodal likelihood profile. CONCLUSIONS: In the presence of nonidentifiability, equally likely parameter estimates might yield different conclusions. Checking for the existence of nonidentifiability and its implications should be incorporated into standard model calibration procedures.