 Research note
 Open Access
 Published:
Bivariate longitudinal data analysis: a case of hypertensive patients at Felege Hiwot Referral Hospital, Bahir Dar, Ethiopia
BMC Research Notesvolume 10, Article number: 722 (2017)
Abstract
Objective
Longitudinal data are often collected to study the evolution of biomedical markers. The study of the joint evolution of response variables concerning hypertension over time was the aim of this paper. A hospital based retrospective data were collected from September 2014 to August 2015 to identify factors that affect hypertensive. The joint mixed effect model with unstructured covariance was fitted. A total of 172 patients screened for antihypertensive drugs treated were longitudinally considered from Felege Hiwot referral.
Results
The joint mixed effect model with unstructured covariance (AIC: 12,236.9 with \( \chi_{12}^{2} \) = 1007.8, P < 10^{−4}) was significantly best fit to the data. The correlation between the evolutions of DBP and SBP was 0.429 and the evolution of the association between responses overtime was found 0.257. Among all covariates included in jointmixedeffectmodels, sex, residence, related disease and time were statistically significant on evolution of systolic and diastolic blood pressure. The joint modeling of longitudinal bivariate responses is necessary to explore the association between paired response variables like systolic and diastolic blood pressure. Fitting joint model with modern computing method is recommended to address questions for association of the evolutions with better accuracy.
Introduction
Hypertension is a chronic condition resulting from high blood pressure in the arteries during circulation. Clinically, a person is said to be hypertensive if the individual’s systolic blood pressure (SBP) is greater than 140 mm Hg and or diastolic blood pressure (DBP) is greater than 90 mmHg and those who were already under medication [1, 2]. Monitoring hypertension is a necessity for patients to assess the progression of blood pressure. Typically two markers (systolic and diastolic blood pressure, SBP and DBP) are measured repeatedly over time after treatment to ensure that no signs of blood pressure problems. These markers are needed to ensure an accurate evaluation of the blood pressure since they are correlated and could be influenced by the socioeconomic status of the patient [3]. Given the interdependence of these outcomes in determining blood pressure, it is important to evaluate the factors that affect the rate of change in these outcomes in a joint manner [4].
When several markers are measured repeatedly, joint modeling was used in the context of jointly studying time to clinical event and repeated measures on surrogate outcomes [5,6,7,8]. Others included joint modeling of the multilevel item response theory (MLIRT) and Cox’s proportional hazard model for time to dependent terminal event with shared random effects to link the two models [9, 10]. As Jaffa et al. [5], Thiebaut et al. [11] and Fiewus and Verebeke [12,13,14] stated a joint modeling of biomarker data is necessary to quantify, the relationship between evolutions of different responses and the evolution of the relationship between different response variables over time. Thus, the aim of this study was using the joint mixed effect model with unstructured covariance structure to answer the evolution of systolic and diastolic blood pressure over time.
Main text
Methods
Hospital based retrospective studies were conducted among hypertensive patients attending antihypertensive clinic between September 2014 to August 2015 at Felege Hiwot referral hospital, Bahir Dar, Ethiopia. The area is selected because it is central referral hospital that provides organized hypertension follows up care. The sample size was calculated using 95% confidence level and the overall sample size was 172 which after adding 5% for nonresponse. A systematic random sampling method was adopted for selecting a representative sample from the list of the medical charts that contain the list of hypertensive patients’ name and identification number. Patients were selected randomly using their unique identification number. The study considered all hypertensive patients, whose age was above 15 years regardless of their treatment category during the study period in the referral hospital.
In this paper, a longitudinal data on two markers of hypertension (SBP and DBP), various socioeconomic status covariates and clinical factors were considered. SBP and DBP measurements were collected at baseline (SBP0 and DBP0), and every 3 months (SBP3 to SBP12 and DBP3 to DBP12) thereafter. In order to ensure sufficient available information, only a subsample of patients having at least four repeated with equally spaced measurements were included in the analyses. Thus, a followup data were included hypertensive patients had a maximum of five repeated measures on blood pressure. Data were entered and edited in SPSS v. 21 and analyzed by SAS v. 9.4 and the 95% confidence interval was used to determine significance tests. The missing data for both markers were replaced by last observation carried forward (LOCF) as [15].
Bivariate random mixed effect model
Multivariate longitudinal data arise when a set of different responses on the same unit are measured repeatedly over time. A joint modeling of such kind of data is necessary to quantify, firstly, the relationship between evolutions of different responses and, secondly, the evolution of the relationship between different response variables over time. Thus, a pair wise fitting approach has been used in this study as proposed in the literature [12, 16]. The joint modeling approach investigated in this study was the bivariate longitudinal mixed effect models that included both fixed and random effects. Thus, in this study DBP_{ij1} and SBP_{ij2} be the bivariate outcomes for the ith subject measured at jth times for outcomes 1 and 2.
We define a general bivariate linear mixed effect model including a random component, a variance–covariance process, and an independent error. Let \( Y^{k}_{i} = \left[ {\begin{array}{*{20}c} {DBP_{i} } \\ {SBP_{i} } \\ \end{array} } \right] \), the response vector for the ith subject with \( DBP_{i} \) and \( SBP_{i} \) having n_{i} sampled measurements of the marker k (k = 1, 2).
To take into account correlation between both markers the following bivariate linear mixed effect model was used.
The mixedeffect models assume that the regression coefficients are a random sample from some population of the possible coefficient and allow one to model variations between study units [17]. The random Coefficient mixed model gives that two random slopes (one for DBP and one for SBP) be fitted for each individual and the variances of measurement errors are different for different markers. The covariance matrix for the random slopes is \( G_{{}} = \left[ {\begin{array}{*{20}c} {\sigma_{1}^{2} } & {\sigma_{12}^{2} } \\ {\sigma_{21}^{2} } & {\sigma_{2}^{2} } \\ \end{array} } \right] \), and the two markers are independent of \( G_{{}} = \left[ {\begin{array}{*{20}c} {\sigma_{1}^{2} } & 0\\ 0& {\sigma_{2}^{2} } \\ \end{array} } \right] \).
In this paper, a bivariate linear mixed model including random effects and independent measurement error for both SBP and DBP was presented as outcome variables in each occasion. The two end points were longitudinally measured as a vector of responses, Y_{i}, at each occasion with this model: \( Y_{i} = X_{i} \beta + Z_{i} \gamma_{i} + W_{i} + \varepsilon_{i} \) with \( \left\{ {\begin{array}{*{20}c} {\varepsilon_{i}^{{}} \sim N(0,\varSigma_{i}^{{}} )} \\ {W_{i} \sim N(0,R_{i} )} \\ {\gamma_{i}^{{}} \sim N(0,G)} \\ \end{array} } \right. \), where \( X_{i} = \left[ {\begin{array}{*{20}c} {X_{i}^{1} } & 0 \\ 0 & {X_{i}^{2} } \\ \end{array} } \right] \), \( \beta = \left[ {\begin{array}{*{20}c} {\beta_{1} } \\ {\beta_{2} } \\ \end{array} } \right] \), \( Z_{i} = \left[ {\begin{array}{*{20}c} {Z_{i}^{1} } & 0 \\ 0 & {Z_{i}^{2} } \\ \end{array} } \right] \), \( \gamma_{i} = \left[ {\begin{array}{*{20}c} {\gamma_{i}^{1} } \\ {\gamma_{i}^{2} } \\ \end{array} } \right] \) and \( W_{i} = \left[ {\begin{array}{*{20}c} {W_{i}^{1} } \\ {W_{i}^{2} } \\ \end{array} } \right] \) is a \( 2n_{i} \)vector of realization of a variance–covariance process \( w_{i} (t) = \left[ {\begin{array}{*{20}c} {w_{i}^{1} (t)} \\ {w_{i}^{2} (t)} \\ \end{array} } \right] \) and \( \varepsilon_{i} = \left[ {\begin{array}{*{20}c} {\varepsilon_{i}^{1} } \\ {\varepsilon_{i}^{2} } \\ \end{array} } \right] \) represents independent measurement errors. Though many socioeconomic covariates and clinical factors were considered in the analysis, only covariates significantly associated with systolic and diastolic blood pressure were reported.
Results
The sample was composed of 172 patients with antihypertensive drugs treated, of which 38 (22.09%) and 21 (12.51%) patients had diabetic and stroke and 76 (44.19%) patients had no other related disease. The results showed that, 93 (54.09%) women, with a mean age of 51.87 (SD = 14.33) years. Whereas mean age of male patients was 56.82 (SD = 14.98) years and the mean systolic and diastolic blood pressure of patients declined from baseline to the next 3 months period, to the next 6 months period and soon up to the last 12 months period follow up time. The baseline SBP mean of patients was 153.58 (SD = 31.60) mmHg and declined to 129.48 (SD = 20.07) mmHg over time and similar history was found in DBO patients. These summaries suggest that the treatment has a significant effect on systolic and diastolic blood pressure over the followup time (Table 1).
Figure 1 shows mean estimates with 95% CI, for the SBP and DBP separately by different factors. Clearly, observe that decreasing trend in both SBP and DBP over time. The mean estimate of blood pressure for a female was more than male at baseline and significantly different (P < 0.0001) between the outcomes. Whereas, the mean estimate of diastolic and systolic blood pressure for rural and urban area seems to no difference The bivariate random mixedeffects model was fitted to show the association of correlate measurements to the 3monthinterval for blood pressure to identify the risk factors.
The bivariate random mixedeffects model was significantly best with AIC value of 12,236.9 and \( \chi_{12}^{2} \) = 1007.8, P < 10^{−4}. Output obtained by the bivariate random mixed effects model with UN provides estimations for the covariance matrix and fixed effects. The variance–covariance matrix for the random effects and the error components for each response were obtained as:
The correlation between the evolutions for the two random slopes is given by:
whereas, the marginal correlation between the two responses as a function of time is given by:
The result revealed that the correlation between the diastolic and systolic blood pressure for random intercepts at baseline was 0.257 and the remaining marginal correlation for 3 to 12 months were 0.368, 0.337, 0.327 and 0.320, respectively over the 3 month period follow up time.
The primary focus of this study was interested in modeling longitudinal measures of SBP and DBP patients as a function of baseline age, sex, residence and related disease. The bivariate random mixedeffects model revealed that baseline (P = 0.0001), 3 months period (P = 0.001) and 6 months period (P = 0.011) systolic and diastolic blood pressure were significantly different compared with 12 months period follow up time. The systolic and diastolic blood pressure of rural patients (P = 0.0308) was significantly different compared with urban patients. The systolic and diastolic blood pressure patients with diabetic showed a significant difference at baseline (P = 0.0077) and 3 months period (P = 0.0015) compared to patients with stroke (Table 2).
Discussion
The joint mixed effect model with unstructured covariance (AIC 12,236.9 with \( \chi_{12}^{2} \) = 1007.8, P < 10^{−4}) was significantly best fit to the covariates. This results in line with [18, 19] as they used a joint model for a longitudinal data and obtained the best final model with UN compared to others.
The correlation between the evolutions for the DBP and SBP was 0.429. This result was in line with a joint model fitted for female SBP and DBP that considered the before and after drug administration data [20]. Similarly, the evolution of the association between the responses over time was found 0.257. The remaining marginal correlation between DBP and SBP was 0.368, 0.337, 0.327 and 0.320, respectively over the followup time. As the marginal correlation converges to the random slopes, the relationship was underlined by the correlation between DBP and SBP at each followup time [12].
The mean difference in systolic and diastolic blood pressure for male and female had only significant (P < 0.034) between the baseline and after a 12month followuptime. This result in line with [21,22,23] as they stated hypertension is the leading cause of death in women than that of men and they concluded that antihypertensive treatment decreases the death of women. The current study revealed that there is an association between blood pressure at baseline and during followup among patients with stroke and diabetes. This result was similar to the study in the Louisiana State University Hospitalbased longitudinal data [24,25,26], and blood pressure at baseline and during followup and the risk of allcause mortality among patients with diabetes [27].
Conclusion and recommendation
The joint mixed effect model with unstructured covariance was preferred among others to fit the data. It can be generalized that, the two outcomes have a strong positive correlation and the joint mixed effect model was preferred. Thus, joint modeling of longitudinal bivariate responses is necessary to explore the association between paired response variables. The baseline mean of the two symptoms was out of the normal range for hypertensive patients but it declines through clinical treatment followup time of 3 month period intervals. As the joint model is able to address questions for the association of the evolutions with better accuracy, fitting joint mixed effect model is recommended.
Limitations
The joint modeling problem is failing to convergence because of a large number of parameter estimates. Moreover, some one also might want to look at modeling more than two response variables over time. Thus one might be able to be implemented using modern computing method for future work to model more than two response variables over time. Here in the joint mixed effect model, authors did not see the interaction effect of the predictors over time. Therefore, researchers should consider the contribution of the interaction effect of the predictors on to the joint mixed effect model.
Abbreviations
 AIC:

Akaike Information Criterion
 CI:

confidence interval
 DBP:

diastolic blood pressure
 LOCF:

last observation carried forward
 MLIRT:

multilevel item response theory
 SAS:

statistical analysis system
 SBP:

systolic blood pressure
 SD:

standard deviation
 SPSS:

Statistical package for Social Science
 UN:

unstructured
References
 1.
World Health Organization. Global status report on noncommunicable diseases 2010, WHO library cataloguinginpublication data. Geneva: World Health Organization; 2011. ISBN 9789241564229.
 2.
Choi YH, Chowdhury R, Swaminathan B. Prediction of hypertension based on the genetic analysis of longitudinal phenotypes: a comparison of different modeling approaches for the binary trait of hypertension. BMC Proc. 2014;8(Suppl 1):s78.
 3.
James GD, Sealey JE, Alderman M, Ljungman S, Mueller FB, Pecker MS, et al. A longitudinal study of urinary creatinine and creatinine clearance in normal subjects race, sex, and age differences. Am J Hypertens. 1988;1(2):124–31.
 4.
Beunckens C, Molenberghs G, Thijs H, Verbeke G. Incomplete hierarchical data. Stat Methods Med Res. 2007;16(5):457–92.
 5.
Jaffa MA, Gebregziabher M, Jaffa AA. Analysis of multivariate longitudinal kidney function outcomes using generalized linear mixed models. J Transl Med. 2015;13(1):192.
 6.
Xu J, Zeger SL. The evaluation of multiple surrogate endpoints. Biometrics. 2001;57(1):81–7.
 7.
Brown ER, Ibrahim JG, DeGruttola V. A flexible Bspline model for multiple longitudinal biomarkers and survival. Biometrics. 2005;61(1):64–73.
 8.
Brown ER, Ibrahim JG. Bayesian approaches to joint curerate and longitudinal models with applications to cancer vaccine trials. Biometrics. 2003;59(3):686–93.
 9.
He B, Luo S. Joint modeling of multivariate longitudinal measurements and survival data with applications to Parkinson’s disease. Stat Methods Med Res. 2016;25(4):1346–58.
 10.
Twisk JW. Applied longitudinal data analysis for epidemiology: a practical guide. Cambridge: Cambridge University Press; 2013.
 11.
Thiébaut R, JacqminGadda H, Chêne G, Leport C, Commenges D. Bivariate linear mixed models using SAS proc MIXED. Comput Methods Programs Biomed. 2002;69(3):249–56.
 12.
Fieuws S, Verbeke G. Joint modeling of multivariate longitudinal profiles: pitfalls of the randomeffects approach. Stat Med. 2004;23(20):3093–104.
 13.
Chakraborty H, Helms RW, Sen PK, Cohen MS. Estimating correlation by using a general linear mixed model: evaluation of the relationship between the concentration of HIV1 RNA in blood and semen. Stat Med. 2003;22(9):1457–64.
 14.
Bahçecitapar M, Karadag Ö, Aktas S. Specification of variancecovariance structure in bivariate mixed model for unequally timespaced longitudinal data. Int J Stat Med Res. 2015;4(4):370.
 15.
Verbeke G. Linear mixed models for longitudinal data. Linear mixed models in practice. Newyork: Springer; 1997. p. 63–153.
 16.
Fieuws S, Verbeke G. Pairwise fitting of mixed models for the joint modeling of multivariate longitudinal profiles. Biometrics. 2006;62(2):424–31.
 17.
Littell RC, Milliken GA, Stroup WW, Wolfinger RD, Schabenberger O. SAS for mixed models. Cary: SAS Institute; 2007.
 18.
Fissuh YH, Muletav G. A joint model for a longitudinal pulse rate and respiratory rate of congestive heart failure patients: at Ayder Referral Hospital of Mekelle University, Tigray, Ethiopia. J Biom Biostat. 2015;6(5):1.
 19.
Facts HA. State heart disease and stroke prevention programs addresses heart attack prevention. 2006.
 20.
Lambert P, Vandenhende F. A copulabased model for multivariate nonnormal longitudinal data: analysis of a dose titration safety study on a new antidepressant. Stat Med. 2002;21(21):3197–217.
 21.
Ong KL, Tso AW, Lam KS, Cheung BM. Gender difference in blood pressure control and cardiovascular risk factors in Americans with diagnosed hypertension. Hypertension. 2008;51(4):1142–8.
 22.
Gu Q, Burt VL, PauloseRam R, Dillon CF. Gender differences in hypertension treatment, drug utilization patterns, and blood pressure control among US adults with hypertension: data from the National Health and Nutrition Examination Survey 1999–2004. Am J Hypertens. 2008;21(7):789–98.
 23.
Keyhani S, Scobie JV, Hebert PL, McLaughlin MA. Gender disparities in blood pressure control and cardiovascular care in a national sample of ambulatory care visits. Hypertension. 2008;51(4):1149–55.
 24.
Zhao W, Katzmarzyk PT, Horswell R, Wang Y, Johnson J, Cefalu WT, et al. Blood pressure and stroke risk among diabetic patients. J Clin Endocrinol Metab. 2013;98(9):3653–62.
 25.
Li W, Wang Y, Chen L, Horswell R, Xiao K, Besse J, et al. Increasing prevalence of diabetes in middle or lowincome residents in Louisiana from 2000 to 2009. Diabetes Res Clin Pract. 2011;94(2):262–8.
 26.
Zhang Y, Li W, Wang Y, Chen L, Horswell R, Xiao K, et al. Increasing prevalence of hypertension in lowincome residents within Louisiana state university health care services division hospital system. Eur J intern Med. 2012;23(8):e179–84.
 27.
Li W, Katzmarzyk PT, Horswell R, Wang Y, Johnson J, Hu G. Blood pressure and allcause mortality among patients with type 2 diabetes. Int J Cardiol. 2016;1:206 (Epub 121).
Authors’ contributions
DL initiated the research and involved in the write up of the proposal, the data entry and data analysis, and write up of the manuscript. DT and HM initiated the research, wrote the research proposal, conducted the research, did data entry and analysis and wrote the manuscript. All authors read and approved the final manuscript.
Acknowledgements
The authors wish to thank the study participants and data collectors for making this study possible.
Competing interests
The authors declared that they have no competing interests.
Availability of data and materials
If needed the raw data in SPSS and SAS code for this article is available.
Consent for publication
Not applicable.
Ethics approval and consent to participate
This study was reviewed and approved by the Bahir Dar University of Science College Research Ethics Committee and Felege Hiwot Referal Hospital Patient’s Ethics Committee. As the study was based on retrospective data, informed consent was not needed.
Funding
This research was funded by Science College Research and Community Service, Bahir Dar University.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
About this article
Received
Accepted
Published
DOI
Keywords
 Joint mixed effect models
 Systolic blood pressure
 Diastolic blood pressure hypertension