CKD-EPI计算eGFR，用R语言如何实现？

hushidi

CKD-EPI计算eGFR，用R语言如何实现？麻烦提供一下代码，谢谢

hellowolrd

伸手要代码。。

不难啊！

Vinnish

hellowolrd 6

tctcab

随便一搜就有结果嘛……

https://www.msdmanuals.com/medical-calculators/gfr_ckd_epi-zh.htm?ivk_sa=1024320u

calc_egfr = function(性别, 年龄, 种族, 血清肌酐) {
  kappa = 0.7
  alpha = -0.329
  sex = 1.018
  race = 1
  
  if (性别 == "男") {
    sex = 1
    alpha = -0.411
    kappa = 0.9
  }
  if (种族 == "黑人") {
    race = 1.159
  }
  egfr = 141 * min(血清肌酐 / kappa, 1) ^ alpha * max (血清肌酐 / kappa, 1) ^ (-1.209) * 0.993 ^
    (年龄) * sex * race
  
  return(egfr)
}

calc_egfr(性别 = "女",
          年龄 = 30,
          种族 = "亚洲",
          血清肌酐 = 100)
#> [1] 0.2885147

参考
Levey AS, Stevens LA, Schmid CH, et al. for the CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration). A New Equation to Estimate Glomerular Filtration Rate. Ann Intern Med. 2009 May 5; 150:604-612. PubMed ID: 19414839 PubMed Logo

dryezl

tctcab
之前就听说谢益辉老师说你很有耐心，现在是真的感觉到了

稍微改一下，用上向量，支持计算多个egfr

calc_egfr = function(性别, 年龄, 种族, 血清肌酐) {
  kappa = ifelse(性别 == "男", 0.9, 0.7) 
  alpha = ifelse(性别 == "男", -0.411, -0.329) 
  sex =  ifelse(性别 == "男", 1, 1.018)
  race =  ifelse(种族 == "黑人", 1.159, 1)
  egfr = 141 * pmin(血清肌酐 / kappa, 1) ^ alpha * pmax (血清肌酐 / kappa, 1) ^ (-1.209) * 0.993 ^
    (年龄) * sex * race
  
  return(egfr)
}

calc_egfr(性别 = "女",
          年龄 = 30,
          种族 = "亚洲",
          血清肌酐 = 100)