Read Data

  • Import data in R and create an R dataframe named newdata

Read GRAIL Model model as R object from a file

  • Print ANOVA Table (if desired)

Logistic (Proportional Odds) Ordinal Regression Model

orm(formula = MGFR ~ Age + Sex + rcs(logAlb, 3) + rcs(logCr, 
    3) + rcs(logBUN, 3), data = grail_data, x = TRUE, y = TRUE, 
    scale = TRUE)
Model Likelihood
Ratio Test		 Discrimination
Indexes		 Rank Discrim.
Indexes
Obs 1235 LR χ2 1278.13 R2 0.645 ρ 0.771
Distinct Y 151 d.f. 8 R21235 0.645
Y0.5 96 Pr(>χ2) <0.0001 R28,1235 0.642
max |∂log L/∂β| 0.0004 Score χ2 1561.71 R21234.3 0.645
Pr(>χ2) <0.0001 R28,1234.3 0.643
g 2.696
gr 14.825
|Pr(Y ≥ median)-½| 0.287
β S.E. Wald Z Pr(>|Z|)
Age  -0.0467  0.0057 -8.13 <0.0001
Sex=Male   1.3326  0.1121 11.88 <0.0001
logAlb   1.9265  0.4870 3.96 <0.0001
logAlb'  -0.8810  0.5481 -1.61 0.1080
logCr  -4.7246  0.4113 -11.49 <0.0001
logCr'  -0.2968  0.6386 -0.46 0.6421
logBUN  -0.4928  0.2665 -1.85 0.0644
logBUN'  -1.2407  0.3703 -3.35 0.0008

Calculate the Score from the Model

  • Applied to saved model to the Harmony_data to calculate the score

KM Curves

KM Curves by GRAIL and CKD-Epi Event: Kidney Ouucome

KM Curves by GRAIL and CKD-Epi Event: Kidney/Death

Cumulative Incidence Rates of Kidney/Death Events by GRAIL and CKD-Epi

Characteristic 15 Days 30 Days 45 Days 60 Days 75 Days 90 Days
GRAIL_CKDEpi
    GRAIL<30 CKD-Epi 30+ 16% (12%, 20%) 27% (23%, 32%) 35% (29%, 39%) 39% (33%, 43%) 39% (34%, 44%) 42% (37%, 47%)
    GRAIL<30 CKD-Epi<30 20% (17%, 22%) 28% (25%, 31%) 34% (31%, 37%) 39% (35%, 42%) 42% (38%, 45%) 46% (42%, 49%)
    GRAIL30+ CKD-Epi 30+ 16% (13%, 18%) 24% (21%, 27%) 29% (26%, 32%) 32% (29%, 35%) 34% (30%, 37%) 37% (34%, 41%)

Cumulative Incidence Rates of Kidney/Death Events by GRAIL and CKD-Epi

Characteristic 15 Days 30 Days 45 Days 60 Days 75 Days 90 Days
GRAIL_CKDEpi
    GRAIL<30 CKD-Epi 30+ 16% (12%, 20%) 27% (23%, 32%) 35% (29%, 39%) 39% (33%, 43%) 39% (34%, 44%) 42% (37%, 47%)
    GRAIL<30 CKD-Epi<30 20% (17%, 22%) 28% (25%, 31%) 34% (31%, 37%) 39% (35%, 42%) 42% (38%, 45%) 46% (42%, 49%)
    GRAIL30+ CKD-Epi 30+ 16% (13%, 18%) 24% (21%, 27%) 29% (26%, 32%) 32% (29%, 35%) 34% (30%, 37%) 37% (34%, 41%)