Latent variable modeling improves AKI risk factor identification and AKI prediction compared to traditional methods
© The Author(s). 2017
Received: 13 October 2016
Accepted: 30 January 2017
Published: 8 February 2017
Acute kidney injury (AKI) is diagnosed based on postoperative serum creatinine change, but AKI models have not consistently performed well, in part due to the omission of clinically important but practically unmeasurable variables that affect creatinine. We hypothesized that a latent variable mixture model of postoperative serum creatinine change would partially account for these unmeasured factors and therefore increase power to identify risk factors of AKI and improve predictive accuracy.
We constructed a two-component latent variable mixture model and a linear model using data from a prospective, 653-subject randomized clinical trial of AKI following cardiac surgery (NCT00791648) and included established AKI risk factors and covariates known to affect serum creatinine. We compared model fit, discrimination, power to detect AKI risk factors, and ability to predict AKI between the latent variable mixture model and the linear model.
The latent variable mixture model demonstrated superior fit (likelihood ratio of 6.68 × 1071) and enhanced discrimination (permutation test of Spearman’s correlation coefficients, p < 0.001) compared to the linear model. The latent variable mixture model was 94% (−13 to 1132%) more powerful (median [range]) at identifying risk factors than the linear model, and demonstrated increased ability to predict change in serum creatinine (relative mean square error reduction of 6.8%).
A latent variable mixture model better fit a clinical cohort of cardiac surgery patients than a linear model, thus providing better assessment of the associations between risk factors of AKI and serum creatinine change and more accurate prediction of AKI. Incorporation of latent variable mixture modeling into AKI research will allow clinicians and investigators to account for clinically meaningful patient heterogeneity resulting from unmeasured variables, and therefore provide improved ability to examine risk factors, measure mechanisms and mediators of kidney injury, and more accurately predict AKI in clinical cohorts.
KeywordsAcute kidney injury Creatinine Latent variable Mixture model Prediction Risk factor
The diagnosis of acute kidney injury (AKI) relies primarily on changes in serum creatinine concentrations (ΔSCr) [1–3]. Changes in serum creatinine, however, can be insensitive and nonspecific for renal injury [4, 5] due to unmeasured confounders such as changes in creatinine production, myocyte injury, intravenous fluid administration, and renal functional reserve [6, 7]. Many of these unmeasured confounders, also known as latent variables , represent an important source of patient heterogeneity with respect to how and if AKI manifests, but are clinically impractical or impossible to measure. Failure to account for factors like these decreases power to identify risk factors for AKI and hinders accurate prediction of postoperative AKI.
Latent variable mixture modeling improves the ability to assess the associations between independent variables and an outcome by accounting for the effect of a latent variable. The model uses measured covariates to empirically stratify a cohort into subpopulations of patients, represented by a latent variable, which are more homogenous than the total cohort. These subpopulations are represented by component models which can be combined to form a comprehensive model to represent the entire cohort . We hypothesized that a latent variable mixture model would increase power to identify significant risk factors for AKI and improve accuracy in predicting a patient’s postoperative ΔSCr compared to a traditional linear model. To test this hypothesis, we built a traditional linear model and a two-component latent variable mixture model to predict ΔSCr in a well-phenotyped clinical trial of AKI following cardiac surgery and compared the models’ goodness-of-fit, power to identify established AKI risk factors, discrimination, and prediction of 48-h postoperative ΔSCr.
After Vanderbilt University Medical Center Institutional Review Board approval, we collected data from a 653-subject prospective clinical trial of perioperative statin use to prevent AKI following cardiac surgery (NCT00791648). The study was conducted according to the Declaration of Helsinki. Patients were eligible to participate in the trial if they were scheduled for elective coronary artery bypass grafting, valve surgery, or ascending aortic surgery requiring thoracotomy or sternotomy. Patients receiving preoperative renal replacement therapy, with liver dysfunction, acute coronary syndrome, pregnancy, current CYP3A4 inhibitor use, and a history of kidney transplant or statin intolerance were ineligible to participate. Six hundred fifty-three patients provided written informed consent. Thirty-eight patients were excluded for failing inclusion criteria or withdrew for personal reasons prior to study initiation, and one patient that completed the study received hemodialysis on postoperative day one and was excluded from 48-h ΔSCr model development since this patient’s ΔSCr no longer reflected renal injury or function. Thus 614 patients were included. No significant association between perioperative statin use and postoperative AKI was demonstrated in the clinical trial .
We chose maximum ΔSCr from baseline to postoperative day 2 to model AKI because serum creatinine is the most common and best characterized marker of renal injury, a 48-h interval is consistent with current consensus guidelines for AKI diagnosis, and a continuous scale rather than a binomial threshold for AKI preserves the measurement of AKI severity and provides the best opportunity to ascertain differences between linear and latent variable mixture modeling techniques. Baseline serum creatinine concentration was defined as the most recent preoperative creatinine measurement and was measured in inpatients on the morning of surgery and within a week prior to surgery in outpatients. Postoperative serum creatinine concentrations were measured at 2:00 am daily throughout hospitalization.
We selected model covariates a priori based on established predictors of post-cardiac surgery AKI and factors known to affect serum creatinine production or dilution [6, 9, 11–14]. Including well-established risk factors for AKI facilitates comparison of each model’s ability to identify significant AKI risk factors for the prediction of ΔSCr. Selected covariates were identical for both the linear model and the latent variable mixture model and included age, body mass index (BMI), baseline glomerular filtration rate estimated using the CKD-EPI formula (eGFR) , baseline serum creatinine, age●baseline serum creatinine interaction term, baseline hematocrit, presence of diabetes, presence of hypertension, duration of surgery, baseline pulse pressure, volume of hydroxyethyl starch administered during surgery, volume of urine output during surgery, duration of cardiopulmonary bypass, duration of aortic cross clamp, maximum intraoperative arterial lactate concentration, and average intraoperative mean blood pressure adjusted for baseline mean blood pressure. Dataset completion was excellent (100% of all serum creatinine data were complete; >99% of all covariate data were complete).
A linear model and a two-component latent variable mixture model were each fit to the maximum ΔSCr from baseline over the first 48 postoperative hours.
The latent variable mixture model is composed of two traditional linear models, known as component models. Each component model represents a subpopulation of patients formed by the latent variable. During fitting, the mixture model agnostically identifies two distinctive subpopulations based on covariate patterns with respect to observed 48-h postoperative ΔSCr. Given that there is uncertainty regarding individual patient subpopulation membership (i.e., subpopulation membership is determined by each patient’s unknown latent variable status, 0 or 1), a probability of being in each subpopulation is initially randomly assigned to each patient and then refined during the iterative model fitting process until convergence criteria are met. Therefore, at the conclusion of model fitting, a patient whose covariate pattern is very consistent with subpopulation 1, for example, may be assigned a 90% probability of subpopulation 1 membership and a 10% probability of subpopulation 2 membership. In this way, each patient’s data may contribute to both component models, improving overall model fit. Each component model represents a data-identified patient subpopulation. If two distinct subpopulations are not identified during the fitting process, the first component model would become identical to the traditional linear model and the coefficients for all the covariates of the second component model would be assigned a value of zero. After completion of model fitting, we developed a support vector machine algorithm to predict patient subpopulation allocation probabilities based on covariate patterns but independent of observed ΔSCr. This enables prediction of ΔSCr using the latent variable mixture model and allows us to compare ΔSCr prediction between latent variable mixture and linear models.
Patient characteristics were summarized with the 50th (10th, 90th) percentiles for continuous variables and percentages for categorical variables. To evaluate the latent variable mixture model relative to the linear model, we compared model: 1) goodness-of-fit, 2) average power to identify established risk factors for AKI, 3) discrimination (ability to rank subjects in order of predicted ΔSCr), and 4) accuracy to predict maximum 48-h postoperative ΔSCr.
Goodness-of-fit was assessed with calibration plots and r2 calculations of predicted ΔSCr versus observed ΔSCr for the latent variable mixture and linear models. To account for the increased flexibility of the latent variable mixture model with respect to differential model fit, Bayesian Information Criteria (BIC) values were calculated for each model and compared using a relative likelihood calculation.
The ability to identify AKI risk factors was assessed using a post hoc calculation of each model’s power to identify established risk factors as significant. For this calculation, 5000 new datasets were generated from our original dataset using standard parametric bootstrapping techniques, and both the latent variable mixture model and the linear model were refit in each new dataset. This produced a set of new risk factor coefficients and associated p-values for each model. Using these sets of new model coefficients, individual risk factor identification power comparisons between the two models were performed, taking our original fitted model coefficients as the power calculations’ alternative hypotheses. A sign test was used to determine the significance of the power comparison between the two models. Additionally, quantile-quantile (Q-Q) plots were used to assess the normalcy of each model’s errors in order to compare each model’s covariate coefficient accuracy.
Model discrimination was evaluated with a permutation test of each model’s Spearman’s correlation coefficients between predicted and observed ΔSCr.
To evaluate prediction of maximum 48-h postoperative ΔSCr, we compared the average of the square of the difference between the predicted and true ΔSCr (i.e., mean squared error relative difference [(predicted ΔSCr – true ΔSCr)2]) .
Models were bootstrapped with 200 replicates to assess for over-fitting and provide internal validation. Statistical analyses were performed in R (version 3.2.0, R Foundation, http://www.r-project.org) and included pROC and flexmix packages.
Subject characteristics and AKI
All subjects (n = 615)
67 (50, 81)
Body mass index, kg/m2
27.7 (22.5, 36.9)
Congestive heart failure
Left ventricular ejection fraction, %
60 (35, 60)
Prior cardiac surgery
Chronic obstructive pulmonary disease
Peripheral vascular disease
Preoperative medication use
Baseline laboratory data
1.01 (0.74, 1.60)
eGFR, ml/min/1.73 m2
72.8 (38.5 96.7)
34 (25, 43)
Perioperative atorvastatin treatment assignment
Cardiopulmonary bypass use
Cardiopulmonary bypass time, min
110.0 (0, 211.6)
Aortic cross clamp use
Aortic cross clamp time, min
0 (0, 139.6)
Intravenous crystalloid, mL
1600 (1000, 3000)
Intravenous hydroxyethyl starch, mL
0 (0, 0)a
Urine output, mL
430 (175, 946)
Arterial lactate, maximum intraoperative, mmol/L
1.7 (0.9, 3.8)
Length of surgery, hours
5.1 (3.6, 7.8)
The median serum creatinine and eGFR at baseline were 1.01 mg/dl (10th, 90th percentile: 0.74, 1.60) and 73 ml/min/1.73 m2 (38, 97). The median maximum ΔSCr within 48 h of surgery was 0.07 (−0.13, 0.52) for the total cohort, 0.50 mg/dl (0.27, 1.04) in patients that developed KDIGO AKI, and 0.03 (−0.15, 0.20) in patients that did not develop AKI.
Latent variable mixture model subpopulation assignments
AKI risk factor identification and estimation accuracy
Associations between established AKI risk factor covariates and maximum 48-h serum creatinine change from baseline using a linear model and each subpopulation of a two-component latent variable mixture model
Latent variable mixture model
Age (per 10 years)
0.037 (−0.044, 0.118)
0.040 (−0.021, 0.101)
0.042 (0.016, 0.068)**
BMI (per 5 kg/m2)
0.040 (0.013, 0.068)**
0.107 (0.081, 0.132)***
0.012 (0.005, 0.019)***
History of hypertension
0.005 (−0.046, 0.056)
0.080 (0.002, 0.158)*
−0.017 (−0.031, −0.003)*
History of diabetes
−0.024 (−0.082, 0.035)
−0.123 (−0.177, −0.068)***
−0.007 (−0.021, 0.008)
Baseline pulse pressure (per 10 mmHg)
0.003 (−0.009, 0.015)
−0.019 (−0.035, −0.003)*
0.004 (7.2e-5, 0.008)*
Baseline SCr (per mg/dL)
0.203 (−0.309, 0.715)
0.054 (−0.217, 0.326)
0.158 (−0.023, 0.339)
Baseline SCr:age interaction
−0.001 (−0.008, 0.007)
0.002 (−0.003, 0.007)
−0.001 (−0.003, 0.002)
Baseline eGFR (per 30 mL/min/1.73 m2)
0.081 (−0.012, 0.174)
0.045 (−0.066, 0.156)
0.099 (0.051, 0.147)***
Baseline hematocrit (per %)
−0.010 (−0.016, −0.005)***
−0.034 (−0.040, −0.028)***
−0.003 (−0.005, −0.001)***
Cardiopulmonary bypass time (per hour)
0.006 (−0.018, 0.030)
−0.072 (−0.108, −0.036)***
0.012 (2.0e-4, 0.024) **
Aortic cross clamp time (per hour)
0.036 (0.001, 0.072)*
0.156 (0.120, 0.192)***
−0.006 (−0.018, 0.006)
Intraoperative hydroxyethyl starch volume (per L)
0.200 (0.000, 0.400)
0.300 (0.100, 0.500)**
0.000 (−0.056, 0.094)
Intraoperative urine output (per L)
−0.100 (−0.200, −0.048)**
−0.300 (−0.400, −0.200)***
−0.100 (−0.094, −0.016)***
Mean intraoperative MAP adjusted for baseline MAP (per 10 mmHg)
0.023 (0.003, 0.042)*
0.064 (0.042, 0.086)***
0.005 (0.001, 0.009)*
Maximum intraoperative lactate (per mmol/L)
0.004 (−0.021, 0.028)
−0.009 (−0.033, 0.016)
0.013 (0.007, 0.019)***
Length of surgery (per hour)
0.034 (0.009, 0.059)**
0.113 (0.086, 0.140)***
0.027 (0.021, 0.034)***
Model discrimination and prediction of 48-h postoperative ΔSCr
The latent variable mixture model demonstrated superior discrimination for predicted ΔSCr compared to the linear model (Permutation test of Spearman’s correlation coefficients, p < 0.001). The relative mean squared error reduction for the latent variable mixture model comparative to the linear model was 6.8%, meaning that the latent variable mixture model predicted 48-h postoperative ΔSCr 6.8% more accurately.
In this study of perioperative AKI, a latent variable mixture model had markedly more power to identify established risk factors for AKI and improved ability to predict a patient’s postoperative ΔSCr than a traditional linear model. These benefits were likely due to superior goodness-of-fit, improved accuracy of covariate coefficient estimation, and enhanced discrimination of predicted postoperative ΔSCr. Latent variable mixture modeling may offer substantial benefits to the study of AKI, and future studies that seek to isolate risk factors for AKI, measure mechanisms of AKI, test therapies for AKI, or seek to predict AKI in clinical cohorts should consider using this methodology.
The improvement of AKI modeling with the latent variable mixture modeling technique indicates that substantial heterogeneity exists within the perioperative AKI population that is not accounted for by observed covariates, and that reliance on traditional linear modeling techniques which inherently assume observed covariates are the only relevant covariates obscures this heterogeneity. This unaccounted for patient heterogeneity within the AKI population may explain why numerous AKI prevention and intervention trials have failed to demonstrate efficacy despite promising preclinical trials.
While new to studies of AKI, latent variable mixture modeling is an established statistical methodology to account for patient heterogeneity in other clinical domains. It has long been used in psychology and genetics research [17, 18], and more recently in oncology. For example, the use of latent variable mixture modeling to model small cell lung cancer growth dynamics from serum biomarker data has improved the prediction of treatment outcomes and decreased reliance on sequential imaging . In acute lung injury, a latent variable mixture modeling technique recently identified patient phenotypes associated with differential treatment effects of high versus low positive end expiratory pressure where traditional modeling had failed . Identification of latent variable subpopulations in patients at risk for AKI may also lead to the identification of subpopulation-specific treatment benefits, enhanced risk stratification, and improved prediction of long-term outcomes.
In the current study, the latent variable mixture model displayed greater power to identify established risk factors for AKI. This improvement results in increased power to identify and characterize novel candidate risk factors, including baseline characteristics, intraoperative exposures, perioperative biomarkers, and patient management techniques that could be modified to reduce AKI. Candidate factors are easily evaluated using the latent variable mixture model by adding the candidate factor to both component models before mixture model fitting. The p-value associated with the candidate factor’s coefficient for each component model determines the significance of the candidate factor in each patient subpopulation. Using latent variable mixture modeling to assess candidate factors will increase discernment of their association with AKI and benefit the search for other non-latent, modifiable AKI risk factors, particularly in modestly sized patient cohorts where power may be low.
Development of the latent variable mixture model does not itself identify the latent variable or binomial pattern of variables, but can suggest potential candidates including renal functional reserve, genetic polymorphisms, clusters of disease exposure, fluid management strategies, or surgical treatments. For example, renal functional reserve is a potentially source of heterogeneity in susceptibility that leads to variation in the manifestation of AKI across patients [7, 21–26]. In our study, older age and higher comorbidity burden (e.g., diabetes, hypertension) is potentially consistent with a population with less renal reserve compared to subpopulation 2 in whom traditional risk modeling performed less well [7, 27, 28]. In the former, a potential lack of renal reserve might explain the larger model coefficients associated with established AKI risk factors such as history of hypertension and diabetes, BMI, baseline hematocrit, aortic cross clamp duration, and length of surgery. In contrast, the potential presence of renal functional reserve might contribute to smaller, and frequently statistically insignificant, model coefficients for established AKI risk factors in subpopulation 2. The latter subpopulation might represent patients in whom sensitive AKI biomarkers may better predict the potential long-term impact of AKI than currently emphasized risk factors. Irrespective of the identity of the latent variable, our results indicate that latent variable mixture modeling can identify subpopulations of patients that may be used to enrich outcomes in clinical trials, target monitoring and interventions, and shed novel insight into the pathophysiology of AKI.
Strengths of this study include the use of high-quality unbiased data collected as part of a prospective clinical trial with little to no missing data. We also retained serum creatinine as a continuous variable to enhance AKI discrimination and prediction [29, 30]. At the same time we acknowledge potential limitations. We did not evaluate latent variable mixture models with more than two subpopulations or perform latent class analysis to empirically determine the number of subpopulations to model. Given the goal of comparing latent variable mixture modeling to traditional linear modeling techniques, we selected the simplest latent variable mixture model for this initial assessment. We observed dramatic results, but increased latent variable flexibility could further improve AKI modeling. A second limitation was the small number of patients that developed moderate or severe AKI (100 or 200% ΔSCr – KDIGO stage II or III), which limited our power to compare latent variable mixture modeling to linear modeling techniques with high precision in patients with moderate or severe AKI. A majority of patients that develop postoperative AKI, however, develop mild AKI, and this outcome remains associated with major short and long-term morbidity [31–33]. Finally, it should be acknowledged that the predictive accuracy of the latent variable mixture model could be reduced in datasets missing model covariate information. However, due to increased model flexibility, the latent variable mixture model will always display equal to or better predictive accuracy compared to a linear model with the same covariates.
A latent variable mixture model increased power to identify established AKI risk factors, more accurately ranked the severity of patients’ 48-h ΔSCr, and more accurately predicted 48-h postoperative ΔSCr compared to a linear model. Latent variable mixture modeling may improve clinicians’ ability to identify novel risk factors and advance the understanding of AKI pathophysiology. Employment of this technique could also advance preoperative AKI risk stratification and provide opportunities to further phenotype and target higher risk patient subpopulations with specific monitoring, preventative strategies, and treatments. Latent variable mixture modeling may provide a powerful technique to advance the study of AKI.
Acute kidney injury
Bayesian information criteria
Body mass index
Estimated glomerular filtration rate
Kidney disease improving global outcomes
Change in serum creatinine concentration
We acknowledge Patricia Hendricks, RN, for nursing support.
We acknowledge funding support from the United States National Institutes of Health (LES GM108554, EDS DK92192-07, FTB GM102676 and GM112871, UL1 TR000445), the Veterans Health Affairs (EDS HSRD & Merit Award IIR 13-073-3), the Vanderbilt Center for Kidney Disease, and the Vanderbilt University Medical Center Department of Anesthesiology. None of these funding bodies had any role in the design of this study, or in the collection, analysis, or interpretation of the data, or in the writing of this manuscript.
Availability of data and materials
The data analyzed during the current study are available from the corresponding author upon reasonable request.
LES, DKS, and FTB contributed to study conception and design. All authors analyzed and interpreted the data and were major contributors to the writing of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics and approval and consent to participate
The Vanderbilt University Medical Center Institutional Review Board approved the study. All subjects provided written consent for the use of their data for research purposes including studies separate from the clinical trial. The study was conducted according to the Declaration of Helsinki.
Open AccessThis 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.
- Bellomo R, Ronco C, Kellum JA, Mehta RL, Palevsky P. Acute renal failure - definition, outcome measures, animal models, fluid therapy and information technology needs: the Second International Consensus Conference of the Acute Dialysis Quality Initiative (ADQI) Group. Crit Care. 2004;8:R204–12.View ArticlePubMedPubMed CentralGoogle Scholar
- KDIGO AKI guideline work group. KDIGO clinical practice guidelines for acute kidney injury. Kidney Int Suppl. 2012;2:1–141.View ArticleGoogle Scholar
- Mehta RL, Kellum JA, Shah SV, Molitoris BA, Ronco C, Warnock DG, Levin A. Acute Kidney Injury Network: report of an initiative to improve outcomes in acute kidney injury. Crit Care. 2007;11:R31.View ArticlePubMedPubMed CentralGoogle Scholar
- Waikar SS, Betensky RA, Emerson SC, Bonventre JV. Imperfect gold standards for kidney injury biomarker evaluation. J Am Soc Nephrol. 2012;23:13–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Waikar SS, Betensky RA, Bonventre JV. Creatinine as the gold standard for kidney injury biomarker studies? Nephrol Dial Transplant Off Publ Eur Dial Transpl Assoc - Eur Ren Assoc. 2009;24:3263–5.Google Scholar
- Bellomo R, Kellum JA, Ronco C. Defining acute renal failure: physiological principles. Intensive Care Med. 2004;30:33–7.View ArticlePubMedGoogle Scholar
- Sharma A, Mucino MJ, Ronco C. Renal functional reserve and renal recovery after acute kidney injury. Nephron Clin Pract. 2014;127:94–100.View ArticlePubMedGoogle Scholar
- Schmiege SJ, Meek P, Bryan AD, Petersen H. Latent variable mixture modeling: a flexible statistical approach for identifying and classifying heterogeneity. Nurs Res. 2012;61:204–12.View ArticlePubMedGoogle Scholar
- Berg KS, Stenseth R, Wahba A, Pleym H, Videm V. How can we best predict acute kidney injury following cardiac surgery?: A prospective observational study. Eur J Anaesthesiol. 2013;30:704–12.View ArticlePubMedGoogle Scholar
- Billings FT, Hendricks PA, Schildcrout JS, Shi Y, Petracek MR, Byrne JG, Brown NJ. High-dose perioperative atorvastatin and acute kidney injury following cardiac surgery: a randomized clinical trial. JAMA. 2016;315:877–88.View ArticlePubMedPubMed CentralGoogle Scholar
- Harel Z, Chan CT. Predicting and preventing acute kidney injury after cardiac surgery. Curr Opin Nephrol Hypertens. 2008;17:624–8.View ArticlePubMedGoogle Scholar
- Huen SC, Parikh CR. Predicting acute kidney injury after cardiac surgery: a systematic review. Ann Thorac Surg. 2012;93:337–47.View ArticlePubMedPubMed CentralGoogle Scholar
- Parolari A, Pesce LL, Pacini D, Mazzanti V, Salis S, Sciacovelli C, Rossi F, Alamanni F, Monzino Research Group on Cardiac Surgery Outcomes. Risk factors for perioperative acute kidney injury after adult cardiac surgery: role of perioperative management. Ann Thorac Surg. 2012;93:584–91.View ArticlePubMedGoogle Scholar
- Kim WH, Lee SM, Choi JW, Kim EH, Lee JH, Jung JW, Ahn JH, Sung KI, Kim CS, Cho HS. Simplified clinical risk score to predict acute kidney injury after aortic surgery. J Cardiothorac Vasc Anesth. 2013;27:1158–66.View ArticlePubMedGoogle Scholar
- Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF, Feldman HI, Kusek JW, Eggers P, Van Lente F, Greene T, Coresh J, Chronic Kidney Disease Epidemiology Collaboration. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150:604–12.View ArticlePubMedPubMed CentralGoogle Scholar
- James G, Witten D, Hastie T, Tibshirani R. An introduction to statistical learning: with applications in R. 1st ed. New York: Springer; 2013. p. 29–33.View ArticleGoogle Scholar
- Bentley MJ, Lin H, Fernandez TV, Lee M, Yrigollen CM, Pakstis AJ, Katsovich L, Olds DL, Grigorenko EL, Leckman JF. Gene variants associated with antisocial behaviour: a latent variable approach. J Child Psychol Psychiatry. 2013;54:1074–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu MK, Gaysina D, Barnett JH, Scoriels L, van de Lagemaat LN, Wong A, Richards M, Croudace TJ, Jones PB, LHA genetics group. Psychometric precision in phenotype definition is a useful step in molecular genetic investigation of psychiatric disorders. Transl Psychiatry. 2015;5:e593.View ArticlePubMedPubMed CentralGoogle Scholar
- Buil-Bruna N, López-Picazo J-M, Moreno-Jiménez M, Martín-Algarra S, Ribba B, Trocóniz IF. A population pharmacodynamic model for lactate dehydrogenase and neuron specific enolase to predict tumor progression in small cell lung cancer patients. AAPS J. 2014;16:609–19.View ArticlePubMedPubMed CentralGoogle Scholar
- Calfee CS, Delucchi K, Parsons PE, Thompson BT, Ware LB, Matthay MA, NHLBI ARDS Network. Subphenotypes in acute respiratory distress syndrome: latent class analysis of data from two randomised controlled trials. Lancet Respir Med. 2014;2:611–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Bosch JP. Renal reserve: a functional view of glomerular filtration rate. Semin Nephrol. 1995;15:381–5.PubMedGoogle Scholar
- DeSanto NG, Anastasio P, Coppola S, Barba G, Jadanza A, Capasso G. Age-related changes in renal reserve and renal tubular function in healthy humans. Child Nephrol Urol. 1991;11:33–40.PubMedGoogle Scholar
- Epstein M. Aging and the kidney. J Am Soc Nephrol. 1996;7:1106–22.PubMedGoogle Scholar
- Rodríguez-Iturbe B, Herrera J, García R. Response to acute protein load in kidney donors and in apparently normal postacute glomerulonephritis patients: evidence for glomerular hyperfiltration. Lancet. 1985;2:461–4.View ArticlePubMedGoogle Scholar
- Thomas DM, Coles GA, Williams JD. What does the renal reserve mean? Kidney Int. 1994;45:411–6.View ArticlePubMedGoogle Scholar
- Friedman EA, Woredekal Y. Diabetic nephropathy. In: Nissenson AR, Berns JS, Lerma EV, editors. Current diagnosis & treatment nephrology & hypertension. 1st ed. USA: McGraw-Hill; 2009. p. 483–91.Google Scholar
- Haase M, Devarajan P, Haase-Fielitz A, Bellomo R, Cruz DN, Wagener G, Krawczeski CD, Koyner JL, Murray P, Zappitelli M, Goldstein SL, Makris K, Ronco C, Martensson J, Martling C, Venge P, Siew E, Ware LB, Ikizler TA, Mertens PR. The outcome of neutrophil gelatinase-associated lipocalin-positive subclinical acute kidney injury: a multicenter pooled analysis of prospective studies. J Am Coll Cardiol. 2011;57:1752–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Ronco C, Kellum JA, Haase M. Subclinical AKI is still AKI. Crit Care. 2012;16:313–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Dawson NV, Weiss R. Dichotomizing continuous variables in statistical analysis: a practice to avoid. Med Decis Mak. 2012;32:225–6.View ArticleGoogle Scholar
- Fitzsimons GJ. Death to dichotomizing: figure 1. J Consum Res. 2008;35:5–8.View ArticleGoogle Scholar
- Elmistekawy E, McDonald B, Hudson C, Ruel M, Mesana T, Chan V, Boodhwani M. Clinical impact of mild acute kidney injury after cardiac surgery. Ann Thorac Surg. 2014;98:815–22.View ArticlePubMedGoogle Scholar
- Gameiro J, Neves JB, Rodrigues N, Bekerman C, Melo MJ, Pereira M, Teixeira C, Mendes I, Jorge S, Rosa R, Lopes JA. Acute kidney injury, long-term renal function and mortality in patients undergoing major abdominal surgery: a cohort analysis. Clin Kidney J. 2016;9:192–200.View ArticlePubMedPubMed CentralGoogle Scholar
- Lassnigg A. Minimal changes of serum creatinine predict prognosis in patients after cardiothoracic surgery: a prospective cohort study. J Am Soc Nephrol. 2004;15:1597–605.View ArticlePubMedGoogle Scholar