Multivariate analysis of longitudinal rates of change.
Ontology highlight
ABSTRACT:
SUBMITTER: Bryan M
PROVIDER: S-EPMC5097016 | biostudies-literature | 2016 Dec
REPOSITORIES: biostudies-literature
Publications
Multivariate analysis of longitudinal rates of change.
Bryan Matthew M Heagerty Patrick J PJ
Statistics in medicine 20160714 28
Longitudinal data allow direct comparison of the change in patient outcomes associated with treatment or exposure. Frequently, several longitudinal measures are collected that either reflect a common underlying health status, or characterize processes that are influenced in a similar way by covariates such as exposure or demographic characteristics. Statistical methods that can combine multivariate response variables into common measures of covariate effects have been proposed in the literature. ...[more]
Similar Datasets
Project description:PurposeTo investigate spatiotemporal correlations among ganglion cell complex (GCC) superpixel thickness measurements and explore underlying patterns of longitudinal change across the macular region.DesignLongitudinal cohort study.SubjectsOne hundred eleven eyes from 111 subjects from the Advanced Glaucoma Progression Study with ≥ 4 visits and ≥ 2 years of follow-up.MethodsWe further developed our proposed Bayesian hierarchical model for studying longitudinal GCC thickness changes across macular superpixels in a cohort of glaucoma patients. Global priors were introduced for macular superpixel parameters to combine data across superpixels and better estimate population slopes and intercepts.Main outcome measuresBayesian residual analysis to inspect cross-superpixel correlations for subject random effects and residuals. Principal component analysis (PCA) to explore underlying patterns of longitudinal macular change.ResultsAverage (standard deviation [SD]) follow-up and baseline 10-2 visual field mean deviation were 3.6 (0.4) years and -8.9 (5.9) dB, respectively. Superpixel-level random effects and residuals had the greatest correlations with nearest neighbors; correlations were higher in the superior than in the inferior region and strongest among random intercepts, followed by random slopes, residuals, and residual SDs. PCA of random intercepts showed a first large principal component (PC) across superpixels that approximated a global intercept, a second PC that contrasted the superior and inferior macula, and a third PC, contrasting inner and nasal superpixels with temporal and peripheral superpixels. PCs for slopes, residual SDs, and residuals were remarkably similar to those of random intercepts.ConclusionsIntroduction of cross-superpixel random intercepts and slopes is expected to improve estimation of population and subject parameters. Further model enhancement may be possible by including cross-superpixel random effects and correlations to address spatiotemporal relationships in longitudinal data sets.
| S-EPMC9559093 | biostudies-literature
Project description:In the context of multivariate multilevel data analysis, this paper focuses on the multivariate linear mixed-effects model, including all the correlations between the random effects when the dimensional residual terms are assumed uncorrelated. Using the EM algorithm, we suggest more general expressions of the model's parameters estimators. These estimators can be used in the framework of the multivariate longitudinal data analysis as well as in the more general context of the analysis of multivariate multilevel data. By using a likelihood ratio test, we test the significance of the correlations between the random effects of two dependent variables of the model, in order to investigate whether or not it is useful to model these dependent variables jointly. Simulation studies are done to assess both the parameter recovery performance of the EM estimators and the power of the test. Using two empirical data sets which are of longitudinal multivariate type and multivariate multilevel type, respectively, the usefulness of the test is illustrated.
| S-EPMC4990185 | biostudies-literature
Project description:BackgroundLongitudinal 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.MethodsWe 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.Results1263 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.ConclusionsThe multivariate mixed cumulative logit model offers the most flexibility and allows one to examine the inter-substance association when proportional odds assumption is violated.
| S-EPMC8571881 | biostudies-literature
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:Two-phase studies are crucial when outcome and covariate data are available in a first-phase sample (e.g., a cohort study), but costs associated with retrospective ascertainment of a novel exposure limit the size of the second-phase sample, in whom the exposure is collected. For longitudinal outcomes, one class of two-phase studies stratifies subjects based on an outcome vector summary (e.g., an average or a slope over time) and oversamples subjects in the extreme value strata while undersampling subjects in the medium-value stratum. Based on the choice of the summary, two-phase studies for longitudinal data can increase efficiency of time-varying and/or time-fixed exposure parameter estimates. In this manuscript, we extend efficient, two-phase study designs to multivariate longitudinal continuous outcomes, and we detail two analysis approaches. The first approach is a multiple imputation analysis that combines complete data from subjects selected for phase two with the incomplete data from those not selected. The second approach is a conditional maximum likelihood analysis that is intended for applications where only data from subjects selected for phase two are available. Importantly, we show that both approaches can be applied to secondary analyses of previously conducted two-phase studies. We examine finite sample operating characteristics of the two approaches and use the Lung Health Study (Connett et al. (1993), Controlled Clinical Trials, 14, 3S-19S) to examine genetic associations with lung function decline over time.
| S-EPMC9392467 | biostudies-literature
Project description:Parallel evolution of phenotypic traits is regarded as strong evidence for natural selection and has been studied extensively in a variety of taxa. However, we have limited knowledge of whether parallel evolution of host organisms is accompanied by parallel changes of their associated microbial communities (i.e., microbiotas), which are crucial for their hosts' ecology and evolution. Determining the extent of microbiota parallelism in nature can improve our ability to identify the factors that are associated with (putatively adaptive) shifts in microbial communities. While it has been emphasized that (non)parallel evolution is better considered as a quantitative continuum rather than a binary phenomenon, quantitative approaches have rarely been used to study microbiota parallelism. We advocate using multivariate vector analysis (i.e., phenotypic change vector analysis) to quantify direction and magnitude of microbiota changes and discuss the applicability of this approach for studying parallelism, and we compiled an R package for multivariate vector analysis of microbial communities ('multivarvector'). We exemplify its use by reanalyzing gut microbiota data from multiple fish species that exhibit parallel shifts in trophic ecology. We found that multivariate vector analysis results were largely consistent with other statistical methods, parallelism estimates were not affected by the taxonomic level at which the microbiota is studied, and parallelism might be stronger for gut microbiota function compared to taxonomic composition. This approach provides an analytical framework for quantitative comparisons across host lineages, thereby providing the potential to advance our capacity to predict microbiota changes. Hence, we emphasize that the development and application of quantitative measures, such as multivariate vector analysis, should be further explored in microbiota research in order to better understand the role of microbiota dynamics during their hosts' adaptive evolution, particularly in settings of parallel evolution.
| S-EPMC9797641 | biostudies-literature
Project description:To examine facility-level pressure ulcer (PrU) development rates and variations in these rates after a system-wide adoption of culture change in Veterans Health Administration (VHA) nursing homes.Four-year retrospective longitudinal design.VHA facilities (N=109) representing 132 nursing homes known as community living centers (CLCs).VHA nursing home residents.PrUs were identified using fiscal year (FY) 2008-11 Minimum Data Set (MDS) data. PrU development was defined as a Stage 2 or larger PrU on an MDS assessment with no PrU on the previous assessment. A risk adjustment model was developed using 105,274 MDS observations to predict the likelihood of PrUs (c-statistic=0.72). A Bayesian hierarchical model that adjusted for differences in the precision of PrU rates from different-size facilities was used to calculate smoothed risk-adjusted (SRA) rates for each facility. The statistical significance of the trend over the 4 years was determined by examining the 95% interval estimate for the slope.Over the 4-year period, the beginning of which coincided with the VHA's system-wide adoption of culture change as a performance measure, median SRA facility PrU development rates were fairly consistent at approximately 4%. The range in SRA rates declined over the years, from a 14.8-percentage point spread to 10.1-percentage point spread. Some facilities had significantly improving SRA rates (e.g., declined steadily from 5.5% to 3.9%) and some had significantly worsening SRA rates (e.g., increased steadily from 5.1% to 7.9%). Seven sites had significantly improving rates (P<.001) that were below the median across all 4 years.A system-wide culture change implementation did not affect CLC PrU development rates, but there was significant variation in facility rates, and some facilities had sustained high performance.
| S-EPMC4719574 | biostudies-literature
Project description:Longitudinal intervention studies with repeated measurements over time are an important type of experimental design in biomedical research. Due to the advent of “omics”-sciences (genomics, transcriptomics, proteomics, metabolomics), longitudinal studies generate increasingly multivariate outcome data. Analysis of such data must take both the longitudinal intervention structure and multivariate nature of the data into account. The ASCA+-framework combines general linear models with principal component analysis and can be used to separate and visualize the multivariate effect of different experimental factors. However, this methodology has not yet been developed for the more complex designs often found in longitudinal intervention studies, which may be unbalanced, involve randomized interventions, and have substantial missing data. Here we describe a new methodology, repeated measures ASCA+ (RM-ASCA+), and show how it can be used to model metabolic changes over time, and compare metabolic changes between groups, in both randomized and non-randomized intervention studies. Tools for both visualization and model validation are discussed. This approach can facilitate easier interpretation of data from longitudinal clinical trials with multivariate outcomes. Author summary Clinical trials are increasingly generating large amounts of complex biological data. Examples can include measuring metabolism or gene expression in tissue or blood sampled repeatedly over the course of a treatment. In such cases, one might wish to compare changes in not one, but hundreds, or thousands of variables simultaneously. In order to effectively analyze such data, both the study design and the multivariate nature of the data should be considered during data analysis. ANOVA simultaneous component analysis+ (ASCA+) is a statistical method which combines general linear models with principal component analysis, and provides a way to separate and visualize the effects of different factors on complex biological data. In this work, we describe how repeated measures linear mixed models, a class of models commonly used when analyzing changes over time and treatment effects in longitudinal studies, can be used together with ASCA+ for analyzing clinical trials in a novel method called repeated measures-ASCA+ (RM-ASCA+).
| S-EPMC8604364 | biostudies-literature
Project description:BackgroundImmuno-epidemiologists are often faced with multivariate outcomes, measured repeatedly over time. Such data are characterised by complex inter- and intra-outcome relationships which must be accounted for during analysis. Scientific questions of interest might include determining the effect of a treatment on the evolution of all outcomes together, or grouping outcomes that change in the same way. Modelling the different outcomes separately may not be appropriate because it ignores the underlying relationships between outcomes. In such situations, a joint modelling strategy is necessary. This paper describes a pairwise joint modelling approach and discusses its benefits over more simple statistical analysis approaches, with application to data from a study of the response to BCG vaccination in the first year of life, conducted in Entebbe, Uganda.MethodsThe study aimed to determine the effect of maternal latent Mycobacterium tuberculosis infection (LTBI) on infant immune response (TNF, IFN-γ, IL-13, IL-10, IL-5, IL-17A and IL-2 responses to PPD), following immunisation with BCG. A simple analysis ignoring the correlation structure of multivariate longitudinal data is first shown. Univariate linear mixed models are then used to describe longitudinal profiles of each outcome, and are then combined into a multivariate mixed model, specifying a joint distribution for the random effects to account for correlations between the multiple outcomes. A pairwise joint modelling approach, where all possible pairs of bivariate mixed models are fitted, is then used to obtain parameter estimates.ResultsUnivariate and pairwise longitudinal analysis approaches are consistent in finding that LTBI had no impact on the evolution of cytokine responses to PPD. Estimates from the pairwise joint modelling approach were more precise. Major advantages of the pairwise approach include the opportunity to test for the effect of LTBI on the joint evolution of all, or groups of, outcomes and the ability to estimate association structures of the outcomes.ConclusionsThe pairwise joint modelling approach reduces the complexity of analysis of high-dimensional multivariate repeated measures, allows for proper accounting for association structures and can improve our understanding and interpretation of longitudinal immuno-epidemiological data.
| S-EPMC8449434 | biostudies-literature
Project description:ObjectiveTo estimate the association between homicide and suicide rates in Brazilian municipalities over a period of 7 years.DesignWe conducted a longitudinal ecological study using annual mortality data from 5507 Brazilian municipalities between 2008 and 2014. Multivariable negative binomial regression models were used to examine the relationship between homicide and suicide rates. Robustness of results was explored using sensitivity analyses to examine the influence of data quality, population size, age and sex on the relationship between homicide and suicide rates.SettingA nationwide study of municipality-level data.ParticipantsMortality data and corresponding population estimates for municipal populations aged 10 years and older.Primary and secondary outcome measuresAge-standardised suicide rates per 100 000.ResultsMunicipal suicide rates were positively associated with municipal homicide rates; after adjusting for socioeconomic and demographic factors, a doubling of the homicide rate was associated with 22% increase in suicide rate (rate ratio=1.22, 95% CI: 1.13 to 1.33). A dose-response effect was observed with 4% increase in suicide rates at the third quintile, 9% at the fourth quintile and 12% at the highest quintile of homicide rates compared with the lowest quintile. The observed effect estimates were robust to sensitivity analyses.ConclusionsMunicipalities with higher homicide rates have higher suicide rates and the relationship between homicide and suicide rates in Brazil exists independently of many sociodemographic and socioeconomic factors. Our results are in line with the hypothesis that changes in homicide rates lead to changes in suicide rates, although a causal association cannot be established from this study. Suicide and homicide rates have increased in Brazil despite increased community mental health support and incarceration, respectively; therefore, new avenues for intervention are needed. The identification of a positive relationship between homicide and suicide rates suggests that population-based interventions to reduce homicide rates may also reduce suicide rates in Brazil.
| S-EPMC7643512 | biostudies-literature