Abstract:Many statistical models have been proposed in the literature for the analysis of longitudinal data. One may propose to model two or more correlated longitudinal processes simultaneously, with a goal of understanding their association over time. Joint modeling is then required to carefully study the association structure among the outcomes as well as drawing joint inferences about the different outcomes. In this study, we sought to model the associations among six nutrition outcomes while circumventing the computational challenge posed by their clustered and high‐dimensional nature. We analyzed data from a 2 × 2 randomized crossover trial conducted in Kenya, to compare the effect of high‐dose and low‐dose iodine in household salt on systolic blood pressure (SBP) and diastolic blood pressure (DBP) in women of reproductive age and their household matching pair of school‐aged children. Two additional outcomes, namely, urinary iodine concentration (UIC) in women and children were measured repeatedly to monitor the amount of iodine excreted through urine. We extended the model proposed by Mwangi et al. (2021, Communications in Statistics: Case Studies, Data Analysis and Applications, 7(3), 413–431) allowing flexible piecewise joint models for six outcomes to depend on separate random effects, which are themselves correlated. This entailed fitting 15 bivariate general linear mixed models and deriving inference for the joint model using pseudo‐likelihood theory. We analyzed the outcomes separately and jointly using piecewise linear mixed‐effects (PLME) model and further validated the results using current state‐of‐the‐art Jones and Kenward methodology (JKME model) used for analyzing randomized crossover trials. The results indicate that high‐dose iodine in salt significantly reduced blood pressure (BP) compared to low‐dose iodine in salt. Estimates for the random effects and residual error components showed that SBP and DBP had strong positive correlation, with effect of the random slope indicating that significantly related outcomes are strongly associated in their evolution. There was a moderately strong inverse relationship between evolutions of UIC and BP both in women and children. These findings confirmed the original hypothesis that high‐dose iodine salt has significant lowering effect on BP. We further sought to evaluate the performance of our proposed PLME model against the widely used JKME model, within the multivariate joint modeling framework through a simulation study mimicking a 2×2 crossover design. From our findings, the multivariate joint PLME model performed exceptionally well both in estimation of random‐effects matrix (G) and Hessian matrix (H), allowing satisfactory model convergence during estimation. It allowed a more complex fit to the data with both random intercepts and slopes effects compared to the multivariate joint JKME model that allowed for random intercepts only. When a hierarchical viewpoint is adopted, in the sense that outcomes are specified conditionally upon random effects, the variance–covariance matrix of the random effects must be positive definite. In some cases, additional random effects could explain much variability in the data, thus improving precision in estimation of the estimands (effect size) parameters. The key highlight in this evaluation shows that multivariate joint JKME model is a powerful tool especially while fitting mixed models with random intercepts only, in crossover design settings. Addition of random slopes may lead to model complexities in most cases, resulting in unsatisfactory model convergence during estimation. To circumvent convergence pitfalls, extention of JKME model to PLME model allows a more flexible fit to the data (generated from crossover design settings), especially in the multivariate joint modeling framework.

A Cross-Validated Ensemble Approach to Robust Hypothesis Testing of Continuous Nonlinear Interactions: Application to Nutrition-Environment Studies

Application of Flexible Matching Strategy to Detect Gene-Environment Interactions for Increasing the Study Power

A robust and adaptive framework for interaction testing in quantitative traits between multiple genetic loci and exposure variables

A multivariate approach to joint testing of main genetic and gene-environment interaction effects

Collapsible Kernel Machine Regression for Exposomic Analyses

Robust semiparametric gene-environment interaction analysis using sparse boosting.

Robust Identification of Gene-Environment Interactions under High-Dimensional Accelerated Failure Time Models

A Robust Approach for Identifying Gene-Environment Interactions for Prognosis

A Varying Coefficient Model to Jointly Test Genetic and Gene–Environment Interaction Effects

Robust Bayesian variable selection for gene-environment interactions

Kernel Machine and Distributed Lag Models for Assessing Windows of Susceptibility to Environmental Mixtures in Children's Health Studies

Bayesian kernel machine regression-causal mediation analysis

A penalized robust method for identifying gene-environment interactions.

Identifying Gene-Environment Interactions for Prognosis Using a Robust Approach

A penalized robust semiparametric approach for gene-environment interactions

A unified method for rare variant analysis of gene-environment interactions

Pairwise fitting of piecewise mixed models for the joint modeling of multivariate longitudinal outcomes, in a randomized crossover trial

Functional varying index coefficient model for dynamic gene-environment interactions

Evaluating the robustness of targeted maximum likelihood estimators via realistic simulations in nutrition intervention trials

Joint Test of Parametric and Nonparametric Effects in Partial Linear Models for Gene-environment Interaction

Identifying Gene-Environment Interactions with A Least Relative Error Approach