🤖 AI Summary
In longitudinal studies, time-varying endogenous covariates—whose evolution depends on their own and outcome history—combined with asynchronous measurements and missing data induce bias in generalized linear mixed model (GLMM) estimation, especially for rare or degenerative diseases. To address this, we propose two novel models: the Joint Mixture Model (JMM) and the Joint Scaling Model (JSM), the first to extend these frameworks to continuous endogenous covariates. Both explicitly model contemporaneous and lagged associations between covariates and non-Gaussian outcomes. By leveraging shared random effects structures and scaling factors, they yield interpretable marginal association coefficients that capture dynamic, time-dependent relationships. Bayesian inference is performed efficiently via integrated nested Laplace approximation (INLA). Simulation studies and analysis of natural history data from Duchenne muscular dystrophy demonstrate unbiased estimation of population-level contemporaneous and lagged effects, substantially enhancing clinical interpretability and practical utility.
📝 Abstract
In longitudinal studies, time-varying covariates are often endogenous, meaning their values depend on both their own history and that of the outcome variable. This violates key assumptions of Generalized Linear Mixed Effects Models (GLMMs), leading to biased and inconsistent estimates. Additionally, missing data and non-concurrent measurements between covariates and outcomes further complicate analysis, especially in rare or degenerative diseases where data is limited. To address these challenges, we propose an alternative use of two well-known multivariate models, each assuming a different form of the association. One induces the association by jointly modeling the random effects, called Joint Mixed Model (JMM); the other quantifies the association using a scaling factor, called Joint Scaled Model (JSM). We extend these models to accommodate continuous endogenous covariates and a wide range of longitudinal outcome types. A limitation in both cases is that the interpretation of the association is neither straightforward nor easy to communicate to scientists. Hence, we have numerically derived an association coefficient that measures the marginal relation between the outcome and the endogenous covariate. The proposed method provides interpretable, population-level estimates of cross-sectional associations (capturing relationships between covariates and outcomes measured at the same time point) and lagged associations (quantifying how past covariate values influence future outcomes), enabling clearer clinical insights. We fitted the JMM and JSM using a flexible Bayesian estimation approach, known as Integrated Nested Laplace Approximation (INLA), to overcome computation burden problems. These models will be presented along with the results of a simulation study and a natural history study on patients with Duchenne Muscular Dystrophy.