 Research article
 Open access
 Published:
Use of partitioned GMM marginal regression model with timedependent covariates: analysis of Chinese Longitudinal Healthy Longevity Study
BMC Medical Research Methodology volume 20, Article number: 128 (2020)
Abstract
Background
Elderly population’s health is a major concern for most industrial nations. National health surveys provide a measure of the state of elderly health. One such survey is the Chinese Longitudinal Healthy Longevity Survey. It collects data on risk factors and outcomes on the elderly. We examine these longitudinal survey data to determine the changes in health and to identify risk factors as they impact health outcomes including the elderly’s ability to do a physical check.
Methods
We use a Partitioned GMM logistic regression model to identify risk factors. The model also accounts for the correlation between lagged timedependent covariates and the outcomes. It addresses present and past measures of timedependent covariates on simultaneous outcomes. The relation produces additional regression coefficients as byproduct of the Partitioned model, identifying the immediate, delayed effects (lag − 1), further delayed (lag2), etc. Therefore, the model presents the opportunity for decision makers to monitor the covariate over time. This technique is particularly useful in healthcare and health related research. We use the Chinese Longitudinal Health Longevity Survey data to identify those risk factors and to display the utility of the model.
Results
We found that one’s ability to make own decisions, frequently consuming vegetables, exercise frequently, one’s ability to transfer without assistance, having visual difficulties and being able to pick book from floor while standing had varying effects of significance on one’s health and ability to complete physical checks as they get older.
Conclusions
The partitioning of the covariates as immediate effect, delayed effect or further delayed effect are important measures in a declining population.
Background
Longitudinal studies in medicalrelated research are useful in identifying changes in outcomes as impacted by certain risk factors. While the repeated measurements on subjects generate correlated observations, they are of different types of correlation. There is correlation among the responses. There is correlation between the timedependent covariates and the response. These correlations have different impacts on the outcomes. Thus, any models fitted to these data need to address these correlations accordingly.
Modelling timedependent covariates when analyzing binary outcomes in longitudinal studies has drawn attention. There are methods due to Generalized Estimating Equations (GEE) and others based on Generalized Method of Moments (GMM) [1,2,3,4,5]. However, these methods do not separate out the impact of the timedependent covariates on the outcomes. In fact, they provide estimates that represent an average of the impacts. Obermeier et.al [6]. suggested that when modeling longitudinal data, one could not assume that the association between a timedependent covariate and the outcome was only direct and simultaneous. This is because the outcome might depend on past measurements of the covariate. Thus, an alternative approach is to separate the different impacts of the covariate. Heagerty [7] indicated that one way to properly model longitudinal outcomes with timedependent covariates is to include appropriate lagged values of such covariates. This approach requires additional regression coefficients for each segment of timedependent covariate. These additional coefficients allow parsing of the effect of the covariate on the response, rather than assuming that the association maintains the same strength and direction over time. It provides insight into the effects of timedependent covariates on present and future values of the outcomes.
Motivating example
Elderly population’s health is a major concern for most industrial nations. National health surveys provide a measure of the state of elderly health. One such survey uses the Chinese Longitudinal Healthy Longevity Study (CLHLS) [8]. It collects data on risk factors and outcomes on the elderly population. The CLHLS was designed to identify key factors contributing to healthy longevity among elderly adults in China. The survey was conducted over time but we concentrated on four waves 2005, 2008, 2011 and 2014. This survey is of particular interest in China, as their annual growth rate of the elderly population is approximately 4.4% and approximately 20% of the world’s oldest population live in China [8]. Gu, Zhang and Zeng [9] investigated the impact of adequate access to healthcare. Li, Zhang and Liang [10] used waves 1 & 2 to determine how living arrangements in 1998 impacted selfrated health in 2000. Zheng et.al [11]. studied the associations of environmental variables. Wu and Schimmele [12] tested how levels of psychological disposition in 1998 impacted selfrated health in 2000. Wang, Zheng, Kurosawa and Inaba [13] studied gender and age differences in health among elderly Chinese using data collected in 2002. However, in all of these studies only one or two waves of data were used and researchers were only able to determine crosssectional or lag1 effects of timedependent covariates on the outcomes.
In this paper, we made use of four waves to demonstrate the fit of Partitioned GMM for binary simultaneous outcomes, completion of a physical check and their health status. These responses were objectively measured by an interviewer. There are subjective measures but we concentrated on the objective measures. We focused our attention on the longitudinal aspect of the data and used all four waves. This increased number of waves used allows us to optimize the longitudinal nature of the data.
Data
The data consisted of elderly people 64 years and older living in 22 of 31 provinces in China. There were 8084 observations measured on 2021 individuals over the four waves. We fit models to interviewerrated health and completion of a physical check that included the timeindependent covariate gender. These models also included the timedependent covariates: able to make own decision, consumed vegetables frequently, exercised, transfer without assistance, visual difficulty and ability to pick up book from floor while standing. Descriptive statistics for the outcomes and timedependent covariates are given in Tables 1 and 2, respectively. Our initial observation suggested a steady decline in the percentage of interviewees considered healthy over time, Table 1.
Methods
We fit a partitioned GMM logistic regression model [14] to the Chinese Longitudinal Healthy Longevity Study data to determine the effects of timedependent covariates on the binary outcomes. The model measures the impact of time independent and timedependent covariates X on the outcome Y measured at four different time points. Thus, there are some relations between X and Y other than cross sectional that must be addressed, Fig. 1. Thus, the partitioned GMM logistic regression model [14] provides coefficient estimates for the effect of X on Y when both are measured at the same time, for when X is measured onetime period ahead of Y, for when X is measured twotime periods ahead to Y and for when X is measured threetime periods ahead to Y.
Partitioned GMM logistic regression models with time dependent covariates
Let y_{it} denote the binary observation for individual i (i = 1, …, N) at time t (t = 1, …, T). Let x_{it} = (x_{i1t}, …, x_{iJt}) be a vector of J timedependent covariates, where x_{ijt} is the j^{th} covariate observed at time t for individual i. Assume that observations y_{is} and y_{kt} are independent when i ≠ k but not necessarily when i = k and s ≠ t. The Partitioned GMM logistic regression model accounts for the relationships between the outcomes observed at time t, y_{i} = (y_{i1}, .., y_{iT}) and the j^{th} covariate observed at time s, x_{ijs} for s ≤ t. For each subject i and each timedependent covariate x_{ijt} measured at times t = 1, 2, …, T; the data matrix is reconfigured as a lower triangular matrix,
where the superscript denotes the difference, t − s in timeperiods between the response time t and the covariate time s. In this matrix, \( {\boldsymbol{x}}_{ij}^{\left[0\right]} \) contains values of the timedependent covariate observed at the same time as the outcome, \( {\boldsymbol{x}}_{ij}^{\left[1\right]} \) includes values of the timedependent covariate observed onetime period prior to outcomes, and so on such that \( {\boldsymbol{x}}_{ij}^{\left[T1\right]} \) consists of the values of the covariate measured T − 1 time periods prior to outcome. Thus, the model for the outcome at time t with one timeindependent covariate and one timedependent covariate is
while the model for all time periods in matrix form is
The coefficient \( {\beta}_j^{tt} \) denotes the effect of the covariate x_{ijt} on the response Y_{t} when both are observed in the same timeperiod, while the vector of coefficient β_{F} denotes the effect of the timeindependent covariate x_{F} on the response Y_{t}. When s < t, we denote the lagged effect of the covariate x_{js} on the response Y_{t} by the coefficients \( {\beta}_j^{\left[1\right]},{\beta}_j^{\left[2\right]},\dots, {\beta}_j^{\left[T1\right]} \). In general, each of the J timedependent covariates yield a maximum of T partitions of β_{j}. Thus, for a model with J covariates, the data matrix X has a maximum dimension of NT by (J × T) + 1, and β is a vector of maximum length (J × T) + 1.
This method of estimating regression coefficients relies on valid moment conditions resulting from the covariate values at different times on the outcome at other times. The moment conditions are valid at crosssectional measurements where covariates are measured at the same time as the outcome [2]. However, valid moment conditions between lagged covariates and the outcomes need to be tested. We do so through a test of bivariate correlation developed by Lalonde, Wilson and Yin [3]. Once the valid moments are identified, the regression parameters are estimated using a GMM approach [14]. We do not rehash the derivations here. We encourage the readers, who want to see that development to go to Lalonde, Wilson, and Yin [3], and Irimata, Broatch, and Wilson [14]. We fit these models through SAS 9.4 software using the %partitionedGMM macro (https://github.com/kirimata/PartitionedGMM) [15]. It includes the test for valid moment conditions [3].
In our analysis of data in CLHLS, we fit two partitioned GMM logistic regression models to model interviewerrated health and interviewees’ ability to complete a physical check separately.
Results
Health
Immediate impacts were identified for vegetables (OR = 1.70 with 95% CI: 1.30, 2.23), exercise (OR = 2.03 with 95% CI: 1.52, 2.71), transfer without assistance (OR = 3.65 with 95% CI: 2.39, 5.59), having visual difficulties (OR = 0.64 with 95% CI: 0.49, 0.84) and pick book from floor while standing (OR = 4.11 with 95% CI: 3.11, 5.43), Table 3. For a one timeperiod lag (i.e. delayed effect), exercise (OR = 1.39 with 95% CI: 1.03, 1.89) and transfer without assistance (OR = 1.76 with 95% CI: 1.05, 2.95) significantly impacts the outcome. Across a twotime period lag (further delayed effect), transfer without assistance (OR = 0.44 with 95% CI: 0.24, 0.81) had a significant impact on this outcome, Table 3. There were no significant effects across a threetime period lag (furthermost delayed effect).
Complete physical check
Immediate impacts were obtained for making own decisions (OR = 1.61 with 95% CI: 1.05, 2.48), transfer without assistance (OR = 13.83 with 95% CI: 8.23, 23.27), visual difficulties (OR = 0.39 with 95% CI: 0.25, 0.61) and pick up book from floor while standing (OR = 5.88 with 95% CI: 4.04, 8.54), Table 3, Fig. 3.
Further impacts were seen at lag2 for transfer without assistance (OR = 4.30 with 95% CI: 1.78, 10.43). An additional delayed impact at lag3 was seen for eating vegetables frequently (OR = 2.12 95% CI: 1.04, 4.33), Table 3, Fig. 3.
Discussion
The uniqueness of the partitioned GMM logistic regression models allows the immediate effect as well as future effects of timedependent covariates on outcomes to be measured. Unlike the previous studies, researchers analyzed the CLHLS data but were only able to estimate crosssectional or lag1 effects of timedependent covariates. However, we were able to determine both crosssectional and lag1 associations as well as lag2 and lag3 relationships between the timedependent covariates and our two binary outcomes, Table 4.
Figure 2 presents the relationships between the timedependent covariates and interviewerrated health, over time. We found that gender and the ability to make one’s own decision did not impact the probability of good health. Frequent consumption of vegetables increased good health immediately, but did not have any significant lagged effects. Exercising significantly increased the likelihood of being in good health immediately and in the next time period. The ability to transfer without assistance has a positive impact on good health immediately and in the next time period. Having visual challenges has an immediate negative impact on having good health. The ability to pick book from floor while standing has an immediate positive impact on good health.
Gender did not significantly impact the likelihood of completing a physical check. The ability to make one’s own decisions has an immediate positive impact on completing a physical check. Consumption of vegetables frequently in the first wave significantly increased the likelihood of completing a physical check in the last wave. Exercising did not impact the completion of a physical check at any point in time. Ability to transfer without assistance significantly increases the likelihood of completing a physical check immediately and across a two timeperiod lag. Having visual challenges negatively impacted completing a physical check immediately. Being able to pick up a book from floor while standing increases the probability of completing a physical check. Figure 3 presents the changing relationships between the timedependent covariates and the ability to complete a physical check.
Conclusions
Though we fitted the Partitioned GMM model to two binary outcomes, this model readily accommodates continuous outcomes. The partitioning of the data matrix with the use of additional coefficients provides an opportunity to measure the covariate on the responses at different periods.
Availability of data and materials
The dataset analyzed during the current study is available at the InterUniversity Consortium for Political and Social Research repository, https://www.icpsr.umich.edu/icpsrweb/ICPSR/studies/36692.
Abbreviations
 CLHLS:

Chinese Longitudinal Healthy Longevity Study
 GMM:

Generalized Method of Moments
 WPA:

World Population Aging
 OR:

Odds Ratio
 CI:

Confidence Interval
References
Liang KY, Zeger SL. Longitudinal data analysis using generalized linear models. Biometrika. 1986;73(1):13–22.
Lai TL, Small D. Marginal regression analysis of longitudinal data with timedependent covariates: a generalised method of moments approach. J R Stat Soc Ser B. 2007;69(1):79–99.
Lalonde TL, Wilson JR, Yin J. GMM logistic regression models for longitudinal data with timedependent covariates. Stat Med. 2014;33(27):4756–69.
Guerra M, Shults J, Amsterdam J, TenHav T. The analysis of binary longitudinal data with timedependent covariates. Stat Med. 2012;31(10):931–48.
Zhou Y, Lefante J, Rice J, Chen S. Using modified approaches on marginal regression analysis of longitudinal data with time dependent covariates. Stat Med. 2014;33(19):3354–64.
Obermeier V, Scheipl F, Heumann C, Wassermann J, Küchenhoff H. Flexible distributed lags for modelling earthquake data. J R Stat Soc Ser C Appl Stat. 2015;64(2):395–412.
Heagerty PJ. Marginalized transition models and likelihood inference for longitudinal cateogorical data. Biometrics. 2002;58(2):342–51.
Zeng Y, Vaupel J, Xiao Z, Liu Y, Zhang Z. Chinese Longitudinal Healthy Longevity Survey (CLHLS), 19982012. Ann Arbor, MI;.
Gu D, Zhang Z, Zheng Y. Access to healthcare services makes a difference in healthy longevity among older Chinese adults. Soc Sci Med. 2009;68(2):210–9.
Li L, Zhang J, Liang J. Health among oldest old in China: which living arrangements make a difference. Soc Sci Med. 2009;68(2):220–7.
Zheng Y, Gu D, Purser J, Hoenig H, Christakis N. Associations of environmental factors with elderly health and mortality in China. Am J Public Health. 2010;100(2):298–305.
Wu Z, Schimmele C. Psychological disposition and selfreported health among the oldestold in China. Ageing Soc. 2006;26(1):135–51.
Wang D, Zheng J, Kurosawa M, Inaba Y. Relationship between age and gender differentials in health among older people in China. Ageing Soc. 2009;29(7):1141–54.
Irimata KM, Broatch J, Wilson JR. Partitioned GMM logistic regression models for longitudinal data. Stat Med. 2019;38(12):2171–83.
Irimata KM, Wilson JR. Using SAS to estimate lagged coefficients with the %partitionedGMM macro. In: SAS Global Forum 2018 Conference Proceedings. Denver, CO; 2018.
Acknowledgements
Not Applicable.
Declarations
Not Applicable.
Funding
Not Applicable.
Author information
Authors and Affiliations
Contributions
EVA conducted literature review and final data analysis, created tables and figures, wrote methods, results, discussion and conclusion sections. DX identified outcomes of interest and conducted initial analysis and data cleaning. JRW revised and edited all drafts of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not Applicable.
Consent for publication
Not Applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
VazquezArreola, E., Xue, D. & Wilson, J.R. Use of partitioned GMM marginal regression model with timedependent covariates: analysis of Chinese Longitudinal Healthy Longevity Study. BMC Med Res Methodol 20, 128 (2020). https://doi.org/10.1186/s12874020010030
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12874020010030