Analysis of multivariate longitudinal substance use outcomes using multivariate mixed cumulative logit model
Ontology highlight
ABSTRACT: Abstract
Background
Longitudinal assessments of usage are often conducted for multiple substances (e.g., cigarettes, alcohol and marijuana) and research interests are often focused on the inter-substance association. We propose a multivariate longitudinal modeling approach for jointly analyzing the ordinal multivariate substance use data. Methods
We describe how the binary random slope logistic regression model can be extended to the multi-category ordinal outcomes. We also describe how the proportional odds assumption can be relaxed by allowing differential covariate effects on different cumulative logits for multiple outcomes. Data are analyzed from a P01 study that evaluates the usage levels of cigarettes, alcohol and marijuana repeatedly across 8 measurement waves during 7 consecutive years. Results
1263 subjects participated in the study with informed consent, among whom 56.6% are females. Males and females show significant differences in terms of the time trend for substance use. Specifically, males showed steeper trends on cigarette and marijuana use over time compared to females, while less so for alcohol. For all three substances, age effects appear to be different for different cumulative logits, indicating the violation of proportional odds assumption. Conclusions
The multivariate mixed cumulative logit model offers the most flexibility and allows one to examine the inter-substance association when proportional odds assumption is violated. Supplementary Information
The online version contains supplementary material available at 10.1186/s12874-021-01444-1.
SUBMITTER: Lin X
PROVIDER: S-EPMC8571881 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BACKGROUND:Renal transplant patients are mandated to have continuous assessment of their kidney function over time to monitor disease progression determined by changes in blood urea nitrogen (BUN), serum creatinine (Cr), and estimated glomerular filtration rate (eGFR). Multivariate analysis of these outcomes that aims at identifying the differential factors that affect disease progression is of great clinical significance. Thus our study aims at demonstrating the application of different joint modeling approaches with random coefficients on a cohort of renal transplant patients and presenting a comparison of their performance through a pseudo-simulation study. The objective of this comparison is to identify the model with best performance and to determine whether accuracy compensates for complexity in the different multivariate joint models. METHODS AND RESULTS:We propose a novel application of multivariate Generalized Linear Mixed Models (mGLMM) to analyze multiple longitudinal kidney function outcomes collected over 3 years on a cohort of 110 renal transplantation patients. The correlated outcomes BUN, Cr, and eGFR and the effect of various covariates such patient's gender, age and race on these markers was determined holistically using different mGLMMs. The performance of the various mGLMMs that encompass shared random intercept (SHRI), shared random intercept and slope (SHRIS), separate random intercept (SPRI) and separate random intercept and slope (SPRIS) was assessed to identify the one that has the best fit and most accurate estimates. A bootstrap pseudo-simulation study was conducted to gauge the tradeoff between the complexity and accuracy of the models. Accuracy was determined using two measures; the mean of the differences between the estimates of the bootstrapped datasets and the true beta obtained from the application of each model on the renal dataset, and the mean of the square of these differences. The results showed that SPRI provided most accurate estimates and did not exhibit any computational or convergence problem. CONCLUSION:Higher accuracy was demonstrated when the level of complexity increased from shared random coefficient models to the separate random coefficient alternatives with SPRI showing to have the best fit and most accurate estimates.
| S-EPMC4467678 | biostudies-other
Project description:Multilevel item response theory (MLIRT) models have been widely used to analyze the multivariate longitudinal data of mixed types (e.g., categorical and continuous) in clinical studies. The MLIRT models often have unidimensional assumption, that is, the multiple outcomes are clinical manifestations of a univariate latent variable. However, the unidimensional assumption may be unrealistic because some diseases may be heterogeneous and characterized by multiple impaired domains with variable clinical symptoms and disease progressions. We relax this assumption and propose a multidimensional latent trait linear mixed model (MLTLMM) to allow multiple latent variables and within-item multidimensionality (one outcome can be a manifestation of more than one latent variable). We conduct extensive simulation studies to assess the unidimensional MLIRT model and the proposed MLTLMM model. The simulation studies suggest that the MLTLMM model outperforms unidimensional model when the multivariate longitudinal outcomes are manifested by multiple latent variables. The proposed model is applied to two motivating studies of amyotrophic lateral sclerosis: a clinical trial of ceftriaxone and the Pooled Resource Open-Access ALS Clinical Trials database. Copyright © 2017 John Wiley & Sons, Ltd.
| S-EPMC5540878 | biostudies-literature
Project description:We consider analysis of Genetic Analysis Workshop 18 data, which involves multiple longitudinal traits and dense genome-wide single-nucleotide polymorphism (SNP) markers. We use a multivariate linear mixed model to account for the covariance of random effects and multivariate residuals. We divide the SNPs into groups according to the genes they belong to and score them using weighted sum statistics. We propose a penalized approach for genetic variant selection at the gene level. The overall modeling and penalized selection method is referred to as the penalized multivariate linear mixed model. Cross-validation is used for tuning parameter selection. A resampling approach is adopted to evaluate the relative stability of the identified genes. Application to the Genetic Analysis Workshop 18 data shows that the proposed approach can effectively select markers associated with phenotypes at gene level.
| S-EPMC4143695 | biostudies-literature
Project description:Many previous studies have identified associations between gene expression, measured using high-throughput genomic platforms, and quantitative or dichotomous traits. However, we note that health outcome and disease status measurements frequently appear on an ordinal scale, that is, the outcome is categorical but has inherent ordering. Identification of important genes may be useful for developing novel diagnostic and prognostic tools to predict or classify stage of disease. Gene expression data are usually high-dimensional, meaning that the number of genes is much larger than the sample size or number of patients. Herein we describe some existing frequentist methods for modeling an ordinal response in a high-dimensional predictor space. Following Tibshirani (1996), who described the LASSO estimate as the Bayesian posterior mode when the regression coefficients have independent Laplace priors, we propose a new approach for high-dimensional data with an ordinal response that is rooted in the Bayesian paradigm. We show that our proposed Bayesian approach outperforms existing frequentist methods through simulation studies. We then compare the performance of frequentist and Bayesian approaches using a study evaluating progression to hepatocellular carcinoma in hepatitis C infected patients.
| S-EPMC9153983 | biostudies-literature
Project description:The Therapeutic Community (TC) model is considered an effective treatment for substance dependence, particularly for individuals with complex presentations. While a popular approach for this cohort across a number of countries, few studies have focussed on biopsychosocial and longer-term outcomes for this treatment modality. This study reports on substance use, dependence, and biopsychosocial outcomes up to 9 months post-exit from two TC sites. A longitudinal cohort study (n = 166) with two follow-up time points. Measures included substance use, dependence, subjective well-being, social functioning, and mental and physical health. Generalized Linear Models were employed to assess change over time. At 9 months, 68% of participants reported complete 90-day drug abstinence. Alcohol frequency and quantity were reduced by over 50% at 9 months, with 32% of the sample recording 90-day abstinence at 9 months. Both alcohol and drug dependence scores were reduced by over 60%, and small to medium effect sizes were found for a range of psychosocial outcomes at 9 months follow-up, including a doubling of wellbeing scores, and a halving of psychiatric severity scores. Residents who remained in the TC for at least 9 months reported substantially better outcomes. With notably high study follow-up rates (over 90% at 9 months post-exit), these data demonstrate the value of the TC model in achieving substantial and sustained improvements in substance use and psychosocial outcomes for a cohort with severe substance dependence and complex presentations. Implications for optimal length of stay are discussed.
| S-EPMC7020066 | biostudies-literature
Project description:It is important to characterize the temporal trajectories of disease-related biomarkers in order to monitor progression and identify potential points of intervention. These are especially important for neurodegenerative diseases, as therapeutic intervention is most likely to be effective in the preclinical disease stages prior to significant neuronal damage. Neuroimaging allows for the measurement of structural, functional, and metabolic integrity of the brain at the level of voxels, whose volumes are on the order of mm(3). These voxelwise measurements provide a rich collection of disease indicators. Longitudinal neuroimaging studies enable the analysis of changes in these voxelwise measures. However, commonly used longitudinal analysis approaches, such as linear mixed effects models, do not account for the fact that individuals enter a study at various disease stages and progress at different rates, and generally consider each voxelwise measure independently. We propose a multivariate nonlinear mixed effects model for estimating the trajectories of voxelwise neuroimaging biomarkers from longitudinal data that accounts for such differences across individuals. The method involves the prediction of a progression score for each visit based on a collective analysis of voxelwise biomarker data within an expectation-maximization framework that efficiently handles large amounts of measurements and variable number of visits per individual, and accounts for spatial correlations among voxels. This score allows individuals with similar progressions to be aligned and analyzed together, which enables the construction of a trajectory of brain changes as a function of an underlying progression or disease stage. We apply our method to studying cortical β-amyloid deposition, a hallmark of preclinical Alzheimer's disease, as measured using positron emission tomography. Results on 104 individuals with a total of 300 visits suggest that precuneus is the earliest cortical region to accumulate amyloid, closely followed by the cingulate and frontal cortices, then by the lateral parietal cortex. The extracted progression scores reveal a pattern similar to mean cortical distribution volume ratio (DVR), an index of global brain amyloid levels. The proposed method can be applied to other types of longitudinal imaging data, including metabolism, blood flow, tau, and structural imaging-derived measures, to extract individualized summary scores indicating disease progression and to provide voxelwise trajectories that can be compared between brain regions.
| S-EPMC4912927 | biostudies-literature
Project description:BackgroundSubstance use and overdose deaths make up a substantial portion of injury-related deaths in the United States, with the state of Ohio leading the nation in rates of diagnosed substance use disorder (SUD). Ohio's growing epidemic has indicated a need to improve SUD care in a primary care setting through the engagement of multidisciplinary providers and the use of a comprehensive approach to care.ObjectiveThe purpose of this study was to assess the ability of the Weitzman Extension for Community Healthcare Outcomes (ECHO): Comprehensive Substance Use Disorder Care program to both address and meet 7 series learning objectives and address substances by analyzing (1) the frequency of exposure to the learning objective topics and substance types during case discussions and (2) participants' change in knowledge, self-efficacy, attitudes, and skills related to the treatment of SUDs pre- to postseries. The 7 series learning objective themes included harm reduction, team-based care, behavioral techniques, medication-assisted treatment, trauma-informed care, co-occurring conditions, and social determinants of health.MethodsWe used a mixed methods approach using a conceptual content analysis based on series learning objectives and substances and a 2-tailed paired-samples t test of participants' self-reported learner outcomes. The content analysis gauged the frequency and dose of learning objective themes and illicit and nonillicit substances mentioned in participant case presentations and discussions, and the paired-samples t test compared participants' knowledge, self-efficacy, attitudes, and skills associated with learning objectives and medication management of substances from pre- to postseries.ResultsThe results of the content analysis indicated that 3 learning objective themes-team-based care, harm reduction, and social determinants of health-resulted in the highest frequencies and dose, appearing in 100% (n=22) of case presentations and discussions. Alcohol had the highest frequency and dose among the illicit and nonillicit substances, appearing in 81% (n=18) of case presentations and discussions. The results of the paired-samples t test indicated statistically significant increases in knowledge domain statements related to polysubstance use (P=.02), understanding the approach other disciplines use in SUD care (P=.02), and medication management strategies for nicotine (P=.03) and opioid use disorder (P=.003). Statistically significant increases were observed for 2 self-efficacy domain statements regarding medication management for nicotine (P=.002) and alcohol use disorder (P=.02). Further, 1 statistically significant increase in the skill domain was observed regarding using the stages of change theory in interventions (P=.03).ConclusionsThese findings indicate that the ECHO program's content aligned with its stated learning objectives; met its learning objectives for the 3 themes where significant improvements were measured; and met its intent to address multiple substances in case presentations and discussions. These results demonstrate that Project ECHO is a potential tool to educate multidisciplinary providers in a comprehensive approach to SUD care.
| S-EPMC11019412 | biostudies-literature
Project description:Despite extensive research efforts, reconstruction of gene regulatory networks (GRNs) from transcriptomics data remains a pressing challenge in systems biology. While non-linear approaches for reconstruction of GRNs show improved performance over simpler alternatives, we do not yet have understanding if joint modelling of multiple target genes may improve performance, even under linearity assumptions. To address this problem, we propose two novel approaches that cast the GRN reconstruction problem as a blend between regularized multivariate regression and graphical models that combine the L2,1-norm with classical regularization techniques. We used data and networks from the DREAM5 challenge to show that the proposed models provide consistently good performance in comparison to contenders whose performance varies with data sets from simulation and experiments from model unicellular organisms Escherichia coli and Saccharomyces cerevisiae. Since the models' formulation facilitates the prediction of master regulators, we also used the resulting findings to identify master regulators over all data sets as well as their plasticity across different environments. Our results demonstrate that the identified master regulators are in line with experimental evidence from the model bacterium E. coli. Together, our study demonstrates that simultaneous modelling of several target genes results in improved inference of GRNs and can be used as an alternative in different applications.
| S-EPMC10414675 | biostudies-literature
Project description:In the statistical literature, the methods to understand the relationship of explanatory variables on each individual outcome variable are well developed and widely applied. However, in most health-related studies given the technological advancement and sophisticated methods of obtaining and storing data, a need to perform joint analysis of multivariate outcomes while explaining the impact of predictors simultaneously and accounting for all the correlations is in high demand. In this manuscript, we propose a generalized approach within a Bayesian framework that models the changes in the variation in terms of explanatory variables and captures the correlations between the multivariate continuous outcomes by the inclusion of random effects at both the location and scale levels. We describe the use of a spherical transformation for the correlations between the random location and scale effects in order to apply separation strategy for prior elicitation while ensuring positive semi-definiteness of the covariance matrix. We present the details of our approach using an example from an ecological momentary assessment study on adolescents.
| S-EPMC4768818 | biostudies-literature
Project description:BackgroundJoint modelling of longitudinal and time-to-event outcomes has received considerable attention over recent years. Commensurate with this has been a rise in statistical software options for fitting these models. However, these tools have generally been limited to a single longitudinal outcome. Here, we describe the classical joint model to the case of multiple longitudinal outcomes, propose a practical algorithm for fitting the models, and demonstrate how to fit the models using a new package for the statistical software platform R, joineRML.ResultsA multivariate linear mixed sub-model is specified for the longitudinal outcomes, and a Cox proportional hazards regression model with time-varying covariates is specified for the event time sub-model. The association between models is captured through a zero-mean multivariate latent Gaussian process. The models are fitted using a Monte Carlo Expectation-Maximisation algorithm, and inferences are based on approximate standard errors from the empirical profile information matrix, which are contrasted to an alternative bootstrap estimation approach. We illustrate the model and software on a real data example for patients with primary biliary cirrhosis with three repeatedly measured biomarkers.ConclusionsAn open-source software package capable of fitting multivariate joint models is available. The underlying algorithm and source code makes use of several methods to increase computational speed.
| S-EPMC6047371 | biostudies-literature