# Only set if the Rmd file is not in the parent directory (ie. 'projectname/')
knitr::opts_knit$set(root.dir = '../')

knitr::opts_chunk$set(collapse = TRUE, echo = FALSE, message = FALSE, warning = FALSE)
library(tidyverse)
devtools::load_all()
load_data(update = TRUE)
set_options()

ds <- project_data

Purpose: assess if there are any differences in the MDRD and CKD-Epi eGFR calculation equations. Calculations were done using the nephro package. Both formulas use serum creatinine as a surrogate measure of GFR.

Correlation between different eGFR calculation methods

There appears to be a strong correlation between eGFR calculated using the CKD-Epi formula (eGFR Epi) and eGFR calculated using MDRD (eGFR MDRD) as assessed using Spearman's rank correlation (r=0.95, p<0.001). However, a scatter plot of the association found three subjects with extremely high eGFR MDRD (eGFR > 300 ml/min). These subjects were removed from future analysis.

ds %>% 
  cor.test(formula = ~ eGFR + eGFR_mdrd, data = ., method = "spearman")
project_data %>% 
  scatter_plot("eGFR_mdrd", "eGFR",
               "eGFR MDRD", "eGFR CKDEpi")

Serum and Urine Creatinine

Serum and urine creatinine are measured using the same assay. Since serum creatinine is used to estimate eGFR, we need to assess if there are any correlations between serum and urine creatinine. There was no significant associations between serum and urine creatinine at baseline (r=0.20, p=0.003), 3-year follow-up (r=), or 6-year follow-up () assessed using Pearson's correlation.

ds %>% 
  dplyr::filter(VN == 3) %>% 
  cor.test(formula = ~ Creatinine + UrineCreatinine, data = ., method = "pearson")
ds %>% 
  scatter_plot("Creatinine", "UrineCreatinine",
               "Serum creatinine", "Urine Creatinine",
               facet = TRUE)


windyzn/urinaryDBP documentation built on May 4, 2019, 6:32 a.m.