- Research note
- Open Access
Predictive performance of the estimating equations of renal function in Sri Lankan subjects
BMC Research Notes volume 12, Article number: 655 (2019)
This study validates two popular predictive equations of renal function firstly, Modifications of Diet in Renal Disease and secondly, Chronic Kidney Disease Epidemiology Collaboration equations for Sri Lankan cohort. We used data of the patients referred to Renal Research lab of University of Colombo for creatinine clearance measurement.
Predictive performances varied with the gender. Creatinine clearance and predicted renal functions were compared. Both fared unsatisfactorily with R2 ranging from 0.632 to 0.652, and overestimated renal function by 6–15%. The proportion chronic kidney disease staging 1 and 2 returned by Chronic Kidney Disease Epidemiology Collaboration equation showed significant difference, in females. Modifications of Diet in Renal Disease equation significantly under-estimated advanced chronic kidney disease in females. Chronic Kidney Disease Epidemiology Collaboration equation had better accuracy. The study sample had more females, Asian and lower body size and better renal functions than historic cohorts. Thai and Pakistani studies show both equations and their Asian adaptations fare poorly. Chronic kidney disease stages differ significantly with the equation used. Predictive equations have fared unsatisfactorily by overestimating renal functions. We recommend further studies using gold standards of measuring renal function.
The prevalence and incidence of renal diseases is continuously increasing in US and other countries [1,2,3,4], and similar trend is expected in Sri Lanka as well. To level off these incident rates, various initiatives, such as the Kidney Disease Outcomes Quality Initiative (K/DOQI), have provided physicians with guidelines to optimize the care of patients with chronic kidney disease (CKD). These guidelines emphasize the need to assess kidney function using predictive equations rather than serum creatinine alone . It is also highlighted the need of the equations to be validated in large samples of subjects, in particular that they should be tested in non-US populations and in individuals with mild decrease in kidney function or normal renal functions . Validation of the predictive formulas is also particularly important for patients aged 65 and older, who by far have the highest incident rates of end stage renal disease (ESRD) [4, 6, 7].
Estimated glomerular filtration rate (eGFR) is an important parameter in clinical practice as well. The renal dose adjustments of the drugs are based on CKD stages, which in turn dependent of eGFR. This makes predictive performance of the equations extremely important, to prevent over or under dosing of the drugs.
The formulae most widely used to estimate kidney function, as well as being recommended in adults by the K/DOQI guidelines,  are the Cockcroft–Gault (CG) formula  and the recently developed  and later simplified  Modification of Diet in Renal Disease (MDRD) formula. The CG formula is an estimate of creatinine clearance originally developed in a population of 236 Canadian patients, 209 of which were male. The MDRD formula have been developed as an estimation of 125I-Iothalamate renal clearance based GFR measurement in a population of 1628 patients with previously diagnosed CKD [9,10,11]. The Chronic Kidney Disease Epidemiology Collaboration (CKD EPI) formula was published in May 2009. It was developed in an effort to create formula more accurate than the MDRD formula, especially when actual GFR is greater than 60 ml/min/1.73 m2. Researchers pooled data from multiple studies to develop and validate this new equation. The CKDEPI equation performed better than the MDRD equation, especially at higher GFR, with less bias and greater accuracy. National Health and Nutrition Examination Survey (NHANES) data, the median GFR was 94.5 ml/min/1.73 m2 vs. 85.0 ml/min/1.73 m2, and the prevalence of chronic kidney disease was 11.5% versus 13.1% , using the two equations.
This is a study conducted retrospectively. The data source is from Renal Research Lab (RRL) of the Department of Clinical Medicine, Faculty of Medicine, University of Colombo, where patients’ height, weight, gender, age, serum creatinine and creatinine clearance are recorded. Ethical clearance was obtained through the Review Committee of National Hospital of Sri Lanka.
The creatinine measurements of the lab, is regularly validated against standard strength solutions, and have shown no deviations around the levels of clinical significance. Jaffe’s method uses sodium 2,4,6 trinitrophenolate and then spectrophotometry was the method of creatinine estimation. The inward health personnel supervised the urine collection at all times. In catheterized patients, urine was directly collected to the collection bottle. Unless contra-indicated, every subject was asked to keep an intake of 2.5–3.5 l during collection period, and to stick to their usual dietary habits. The urine bottles were transported to RRL immediately after the collection process, where the transportation time was less than half an hour.
We included all patients aged more than 18 years in the study. We excluded records of all patients with acute renal failure, rhabdomyolysis, aged < 18, BMI < 17 kg/m2, BMI > 40 kg/m2, pregnancy and history of muscular disorders from the study.
We adopted abbreviated MDRD equation from the original work  and CKD EPI from original work of authors . We corrected creatinine clearance for BSA, using the Dubois’ formula  and then, adjusted for over estimation by multiplying from 0.81  (please see Additional file 1: Table S1 for the equations). The bias of the measurement was defined as the mean of the difference between estimated and measured GFR, while precision is defined as the standard deviation of the bias. We constructed Combined Root Mean Square (CRMSE) using the formula of (bias2 + precision2)0.5. The recorded data was analyzed using SPSS 21.0 statistical program.
Study population consisted of 475 subjects (n = 475), of which 212 were males (44.6%).
Table 1 shows basic characteristics of the study group.
The eGFR using MDRD formula and CKD EPI formula overestimated the GFR 9.3 and 8.1 ml/min in average in males (Additional file 1: Figure S1), while it was 23.7 and 0.04 ml/min for the females respectively. While CCr and CKD EPI have not produced any outlier data, MDRD has produced number of outliers, in both male and female subgroups, at the high GFR end.
Precision and bias of the estimation
We compared the two methods of estimating GFR against the CCr using linear regression. Males and females have already shown wide differences, we carried out analysis for genders separately. The results of the regression analysis are displayed in Fig. 1.
The R2 values for four subgroups (gender vs. MDRD/CKD EPI) was between 0.626 and 0.652, showing moderately strong relationship.
The bias and the precision of the equations, as defined in methodology section, were used to calculate CRMSE, which assesses the accuracy of the equations. Additionally data was observed for the accuracy of prediction within 15%, 30% and 50% (p15, p30 and p50) of the actual value. Table 2 summarizes the above data on accuracy, and precision of measurement.
In the case of males, about 40% and 65% of all cases were predicted with 15% and 30% accuracy respectively by both equations. However the difference was marked in females. CKD EPI had 44.5% and 79.5% accuracies (15% and 30% respectively) in prediction; MDRD was much poorer with it with 25.5% and 47.1%.
The study sample was heterogeneous for the KDOQI CKD stage. Out of the 263 females in the study population, 62.4% (n = 164) had their eGFR > 90 ml/min. Among males this was 44.8% (n = 95). This is statistically significant difference (Χ2 = 14.63, p < 0.001). Advanced CKD (Stages 4 and 5) was seen in 8.0% among males and 10.2% among females. The CKD staging performed according to MDRD show striking similarity to that done with CCr in both genders. The trend is same with CKD EPI method other than in females. CKD EPI results in a significantly lower proportion of CKD1 in females (41.1% vs. 62.2%, p < 0.0001) and significantly higher proportion of CKD 2 (14.8% vs. 33.8%, p < 0.0001). Additionally MDRD equation tend to underestimate the advanced CKD in females (10.2% vs. 5.0%, p = 0.0246) (Additional file 1: Table S2).
The study sample of the study has accurate representation of a cross section of patients who are presenting to an adult nephrology facility. However, our study population is different from many previous studies due to plethora of reasons. Firstly, it consists of an exclusive Asian population. Caucasian and Black patients have been evaluated in setting up  the equations as well as its validation [15, 16], but there have been hardly any studies  that involved Asian patients, and even such studies had very limited number of subjects. A recent study in Thailand  showed a remarkable degree of variability in CKD prevalence and risk estimates when current equations developed in Caucasians and Asians were applied to the Thai population. MDRD and CKD-EPI overestimated the prevalence of CKD two to four fold in Thai patients, and adjusted CKD EPI for Japanese, sevenfold, illustrating the fact that predictions have been quite poor. A cross-sectional study done in Pakistan , has demonstrated similar results to our study, with P30 values of 68% and 76.1% for MDRD and CKD-EPI equations respectively, and a modified CKD-EPI equation for the Pakistani subjects had better results. Similar study done with participation of Chinese subjects  reaffirms the findings with CKD EPI performing better than MDRD over all GFR ranges and subgroups.
Secondly, our sample had a majority of women, all of them are non pregnant, which is unique among study samples, where the highest percentage has been 53.8% . Women have been traditionally under-represented despite being majority in many countries. Thirdly, the body habitus of the patients have been different to that of earlier studies. The average height (1.537 m), weight (57.7 kg) and BSAs (1.544 m2) all are significantly lower [9, 16, 17] than previous studies (p < 0.001).
Fourthly, our cohort is younger and has more subjects with better renal functions compared to the cohort that used to develop MDRD equation. The average CCr for males was 69.29 ml/min. and 76.84 ml/min for females. In contrast, cohort that was used to develop MDRD equation had mean GFR of 39.8 ± 21.2 ml/min, and a mean age of 50.6 ± 12.7 years. This is not surprising given the large number of potential live kidney donors that were tested in the RRL.
We used creatinine clearance as the measured GFR to compare with others. Creatinine clearance has been shown to overestimate GFR, has large inter- and intra-subject variability , but it was the method used as the gold standard in many of the original works. Inulin, chromium-51 EDTA or Iothalamine clearance would have been the better methods to use; sadly such resources are very costly for a nephrology unit which is specialized on transplant management. Creatinine clearance overestimates GFR about 19% , hence the actual GFRs should be lower than the former, and the predictive equations should return values lower than that of CCr, and we included this correction in calculating CCr for comparison with other methods that estimate GFR.
Both equations predicted an average value equal or greater than CCr, regardless of the gender. Furthermore, the CKD staging was significantly different when CKD EPI used to stage, despite it being more accurate than MDRD. The most significant difference was seen in CKD Stage 2, where the distinction between the healthy and the diseased is not very clear. The difference of the staging approaches an 20% in females in CKD Stage 2. A recent Thai study , shows similar results, where CKD stage changes significantly with the equation used to estimate GFR.
This is particularly important as many of the Sri Lankans are suffering from chronic interstitial nephritis of agricultural communities (CINAC) . The lower protein intakes, lesser body weights, it is expected that serum creatinine to be lower than that of an age and gender matched Caucasian population. Under such situations, overestimation is inevitable. Hence the actual incidence of CKD would be much higher than eGFR based CKD staging would suggest. This theory is vindicated by change of CKD stages in statistically significant proportions when different eGFR formulae are used. Furthermore, drug dosing in renal impairment should be performed with caution, as the KDOQI CKD staging significantly changes with the method used to estimate GFR.
The R2 values for all the subgroups are between 0.626 and 0.652 indicating that there is a large residual factor at play. The CCr accounts for an R2 of 0.869 , when compared to Iothalamine clearance, and MDRD accounts for a R2 of 0.879 to 0.903 in various studies. These shows MDRD and CCr are comparable methods in estimating renal functions. We used the following equation to derive R2 between CCr and MDRD mathematically.
The term R stands for correlation between CCr and MDRD, whereas R1 and R2 stands for that of Iothalamine vs. CCr and Iothalamine vs. MDRD. This gives an estimation of 0.528 < R2 < 0.997, which is within the observed value of 0.626 < R2 < 0.656. No such information is available CKD EPI equation.
The bias and the precision and CRMSE are another ways to look at the accuracy of the equations . CKDEPI had better accuracy than MDRD, especially in females. The p15, p30 and p50 values were considerably lesser than that of the comparative studies . The comparative studies using mono- and multiethnic cohorts shows CKD-EPI to be a better tool for screening of CKD as well as estimating GFR [19, 23, 24].
Conventional GFR estimating equations do not perform well in Sri Lankan population, and performance differs with gender. Most of them overestimate GFR. We recommend development of novel formulae to estimate GFR for the local population, using Inulin or Iothalamine clearance as the gold standard.
Most of the patient population is from the wet zone of the country, where the leading cause of CKD is diabetes and hypertension, and the representation of CKD of unknown origin is lesser. Secondly, we have used a method (creatinine clearance) to estimate GFR which is not regarded as the gold standard in current practice. Creatinine clearance shows intra- and inter-subject variability than inulin or iothalamine clearance.
Availability of data and materials
The datasets used and/or analysed during the current study is uploaded as a additional material (See Additional file 2).
body mass index
body surface area
chronic kidney disease
- CKD EPI:
chronic kidney disease epidemiology collaboration
combined root mean square error
estimated glomerular filtration rate
end stage renal disease
glomerular filtration rate
kidney disease outcomes quality initiative
modifications of diet into renal diseases
measured glomerular filtration rate
national health and nutrition examination survey
Geddes CC, Paul C, van Dijk W, McArthur S, Metcalfe W, Jager KJ, Zwinderman AH, Mooney M, Fox JG, Simpson K. The ERA-EDTA cohort study—comparison of methods to predict survival on renal replacement therapy. Nephrol Dial Transplant. 2005;21(4):945–56.
Schaubel DE, Morrison HI, Desmeules M, Parsons DA, Fenton SS. End-stage renal disease in Canada: prevalence projections to 2005. CMAJ. 1999;2005(160):1557–63.
Stengel B, Billon S, Van Dijk PC, Jager KJ, Dekker FW, Simpson K, Briggs JD. Trends in the incidence of renal replacement therapy for end-stage renal disease in Europe, 1990–1999. Nephrol Dial Transplant. 2003;18:1824–33.
USRDS. Incidence and prevalence of ESRD. In: US renal data system 2003 annual data report. Bethesda: National Institutes of Health, National Institutes of Diabetes, National Institutes of Digestive and Kidney Diseases; 2003. p. 47–60.
Foundation NK. K/DOQI Clinical practice guideline to define chronic kidney disease: evaluation, classification and stratification. Am J Kidney Dis. 2002;39(2):S1–266.
Jager KJ, van Dijk PC, Dekker FW, Stengel B, Simpson K, Briggs JD. The epidemic of aging in renal replacement therapy: an update on elderly patients and their outcomes. Clin Nephrol. 2003;30:352–60.
Jungers P, Chauveau P, Descamps-Latscha B, Labrunie M, Giraud E, Man NK, Grunfeld JP, Jacobs C. Age and gender related incidence of chronic renal failure in a French urban area: a prospective epidemiologic study. Nephrol Dial Transplant. 1996;11:1542–6.
Cockcroft DW, Gault MH. Prediction of creatinine clearance from serum creatinine. Nephron. 1976;16(1):31–41.
Levey AS, Bosch JP, Lewis JB, Greene T, Rogers N, Roth D. A more accurate method to estimate glomerular filtration rate from serum creatinine: a new prediction equation. Modification of Diet in Renal Disease Study Group. Ann Intern Med. 1999;130:461–70.
Levey AS, Greene T, Kusek JW, Beck GJ. A simplified equation to predict glomerular filtration rate from serum creatinine (Abstract). J Am Soc Nephrol. 2000;11:155A.
Klahr S, Levey AS, Beck GJ, Caggiula AW, Hunsicker L, Kusek JW, Striker G. The effects of dietary protein restriction and blood pressure control on the progression of chronic renal disease. N Engl J Med. 1994;330:877–84.
Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF, Feldman HI, et al. Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI), a new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–12.
Levey AS, Stevens LA, Schmid CH, Zhang YL. Castro AF. Ann Intern Med: A new equation to estimate glomerular filtration rate; 2009. p. 150.
Du Bois D, Du Bois EF. A formula to estimate the approximate surface area if height and weight be known. Arch Intern Med. 1916;17(6):863–71.
Fontseré N, Bonal J, Salinas I, Arellano MRD, Rios J, Ferran Torres F, et al. Is the new Mayo Clinic quadratic equation useful for the estimation of glomerular filtration rate in type 2 diabetic patients? Diabetes Care. 2008;31(12):2265–7.
Froissart M, Rossert J, Jacquot C, Paillard M, Houillier P. Predictive performance of the modification of diet in renal disease and Cockcroft–Gault equations for estimating renal function. J Am Soc Nephrol. 2005;16:763–73.
Al Wakeel JS, Hammad D, Al Suwaida A, Tarif N, Chaudhary A, Isnani A, Albedaiwi WA, Mitwalli AH, Ahmad SS. Validation of predictive equations for glomerular filtration rate in the Saudi population. Saudi J Kidney Dis Transpl. 2009;20(6):1030–7.
Kitiyakara C, Yamwong S, Vathesatogkit P, Chittamma A, Cheepudomwit S, Vanavanan S, Hengprasith B, Sritara P. The impact of different GFR estimating equations on the prevalence of CKD and risk groups in a Southeast Asian cohort using the new KDIGO guidelines. BMC Nephrol. 2012;6(13):1.
Jessani S, Levey AS, Bux R, Inker LA, Islam M, Chaturvedi N, Mariat C, Schmid CH, Jafar TH. Estimation of GFR in South Asians: a study from the general population in Pakistan. Am J Kidney Dis. 2014;63(1):49–58.
Li JT, Xun C, Cui CL, Wang HF, Wu YT, Yun AH, Jiang XF, Ma J. Relative performance of two equations for estimation of glomerular filtration rate in a Chinese population having chronic kidney disease. Chin Med J. 2012;125(4):599–603.
Pozuelos EG, Molina ML, Perera RJJ, Herrera DN, Barroso CL, Ramírez BF. Prevalence of hidden renal failure calculated through formulas on the degree of renal function in hypertense patients over 60 referred to out-patients for blood pressure monitoring. Aten Primaria. 2007;39(5):247–53.
Wanigasuriya KP, Peiris-John RJ, Wickremasinghe R. Chronic kidney disease of unknown aetiology in Sri Lanka: is cadmium a likely cause? BMC Nephrol. 2011;12:32.
Teo BW, Xu H, Koh YY, Li J, Subramanian S, Sinha AK, Shuter B, Toh QC, Sethi S, Lee EJ. Glomerular filtration rates in healthy Asians without kidney disease. Nephrology. 2014;19(2):72–9.
Teo BW, Xu H, Wang D, Li J, Sinha AK, Shuter B, Sethi S, Lee EJ. GFR estimating equations in a multiethnic Asian population. Am J Kidney Dis. 2011;58(1):56–63.
We like to acknowledge all the employees of the Renal Research lab of University of Colombo.
This is a self funded project.
Ethics approval and consent to participate
Ethics approval was taken from the Ethics review committee of the National Hospital of Sri Lanka.
Consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Weerakkody, R.M., Sheriff, M.H.R. Predictive performance of the estimating equations of renal function in Sri Lankan subjects. BMC Res Notes 12, 655 (2019). https://doi.org/10.1186/s13104-019-4692-3
- Glomerular filtration rate
- Sri Lanka
- Creatinine clearance