Get data into R
## ~time_90days + status_90days + age_admission + sex + ethnicity +
## liver_transplant_listed + final_type_of_aki + site + etiology_cirrhosis +
## albumin_given_admission + loop_diuretic + aldosterone_antagonist +
## lactulose + rifaximin + prophylactic_antibiotic + nsaids +
## beta_blockers + diabetes + cad + htn + ckd + ascites_admission +
## encephalopathy_admission + gi_bleed_admission + peritonitis_admission +
## hcc_admission + tips_admission + lvp_admission + alcoholic_hepatitis_admission +
## MELD_Na_baseline + creatinine_admission + log_na_admit +
## log_k_admit + log_cl_admit + log_co2_admit + log_bun_admit +
## log_alt_admit + log_alkphos_admit + log_tb_admit + log_alb_admit +
## log_inr_admit + log_wbc_admit + log_hct_admit + log_plt_admit +
## sbp_admission + dbp_admission# Backward selection based on the AIC
sfgAIC <- pec::selectFGR(ff, cause = 1, data =ckd_data, rule = "AIC", direction = "backward")
# Backward selection based on the BIC
sfgBIC <- pec::selectFGR(ff, cause = 1, data = ckd_data,rule = "BIC", direction = "backward")sfgAIC$fit ##
## Right-censored response of a competing.risks model
##
## No.Observations: 1949
##
## Pattern:
##
## Cause event right.censored
## 1 721 0
## 2 118 0
## unknown 0 1110
##
##
## Fine-Gray model: analysis of cause 1
##
## Competing Risks Regression
##
## Call:
## riskRegression::FGR(formula = Hist(time_90days, status_90days) ~
## age_admission + ethnicity2 + ethnicity3 + liver_transplant_listed1 +
## final_type_of_aki2 + final_type_of_aki3 + final_type_of_aki4 +
## final_type_of_aki5 + siteindiana + sitejacksonville +
## sitekentukey + sitemgh + siterochester + siteyale + etiology_cirrhosis2 +
## etiology_cirrhosis4 + etiology_cirrhosis6 + etiology_cirrhosis7 +
## albumin_given_admission1 + lactulose1 + rifaximin1 +
## beta_blockers1 + encephalopathy_admission1 + gi_bleed_admission1 +
## hcc_admission1 + alcoholic_hepatitis_admission1 + MELD_Na_baseline +
## log_k_admit + log_co2_admit + log_alt_admit + log_tb_admit +
## log_alb_admit + log_wbc_admit + log_plt_admit + sbp_admission,
## data = data, cause = cause)
##
## coef exp(coef) se(coef) z p-value
## age_admission 0.01851 1.019 0.00393 4.71 2.5e-06
## ethnicity2 0.35553 1.427 0.15193 2.34 1.9e-02
## ethnicity3 1.03336 2.810 0.23492 4.40 1.1e-05
## liver_transplant_listed1 -1.47320 0.229 0.15767 -9.34 0.0e+00
## final_type_of_aki2 0.90704 2.477 0.13337 6.80 1.0e-11
## final_type_of_aki3 0.93821 2.555 0.09594 9.78 0.0e+00
## final_type_of_aki4 0.36340 1.438 0.19583 1.86 6.4e-02
## final_type_of_aki5 0.85162 2.343 0.14862 5.73 1.0e-08
## siteindiana 0.24578 1.279 0.12825 1.92 5.5e-02
## sitejacksonville -0.66851 0.512 0.26501 -2.52 1.2e-02
## sitekentukey 0.34704 1.415 0.15488 2.24 2.5e-02
## sitemgh -0.18409 0.832 0.12382 -1.49 1.4e-01
## siterochester -0.25046 0.778 0.14313 -1.75 8.0e-02
## siteyale -0.23905 0.787 0.14374 -1.66 9.6e-02
## etiology_cirrhosis2 0.27871 1.321 0.12737 2.19 2.9e-02
## etiology_cirrhosis4 0.73499 2.085 0.28800 2.55 1.1e-02
## etiology_cirrhosis6 0.29045 1.337 0.14662 1.98 4.8e-02
## etiology_cirrhosis7 0.37367 1.453 0.13344 2.80 5.1e-03
## albumin_given_admission1 0.26088 1.298 0.11299 2.31 2.1e-02
## lactulose1 0.17588 1.192 0.09342 1.88 6.0e-02
## rifaximin1 -0.19423 0.823 0.09729 -2.00 4.6e-02
## beta_blockers1 -0.21601 0.806 0.08359 -2.58 9.8e-03
## encephalopathy_admission1 0.42819 1.534 0.09789 4.37 1.2e-05
## gi_bleed_admission1 0.12215 1.130 0.08228 1.48 1.4e-01
## hcc_admission1 0.52506 1.691 0.10961 4.79 1.7e-06
## alcoholic_hepatitis_admission1 -0.26099 0.770 0.12122 -2.15 3.1e-02
## MELD_Na_baseline 0.01929 1.019 0.00800 2.41 1.6e-02
## log_k_admit -0.28939 0.749 0.18832 -1.54 1.2e-01
## log_co2_admit -0.50422 0.604 0.15345 -3.29 1.0e-03
## log_alt_admit 0.12361 1.132 0.04836 2.56 1.1e-02
## log_tb_admit 0.19930 1.221 0.05647 3.53 4.2e-04
## log_alb_admit -0.70586 0.494 0.16268 -4.34 1.4e-05
## log_wbc_admit 0.20625 1.229 0.07740 2.66 7.7e-03
## log_plt_admit -0.18960 0.827 0.07207 -2.63 8.5e-03
## sbp_admission -0.00683 0.993 0.00240 -2.84 4.4e-03
##
## exp(coef) exp(-coef) 2.5% 97.5%
## age_admission 1.019 0.982 1.011 1.027
## ethnicity2 1.427 0.701 1.059 1.922
## ethnicity3 2.810 0.356 1.773 4.454
## liver_transplant_listed1 0.229 4.363 0.168 0.312
## final_type_of_aki2 2.477 0.404 1.907 3.217
## final_type_of_aki3 2.555 0.391 2.117 3.084
## final_type_of_aki4 1.438 0.695 0.980 2.111
## final_type_of_aki5 2.343 0.427 1.751 3.136
## siteindiana 1.279 0.782 0.994 1.644
## sitejacksonville 0.512 1.951 0.305 0.861
## sitekentukey 1.415 0.707 1.044 1.917
## sitemgh 0.832 1.202 0.653 1.060
## siterochester 0.778 1.285 0.588 1.031
## siteyale 0.787 1.270 0.594 1.044
## etiology_cirrhosis2 1.321 0.757 1.029 1.696
## etiology_cirrhosis4 2.085 0.480 1.186 3.667
## etiology_cirrhosis6 1.337 0.748 1.003 1.782
## etiology_cirrhosis7 1.453 0.688 1.119 1.887
## albumin_given_admission1 1.298 0.770 1.040 1.620
## lactulose1 1.192 0.839 0.993 1.432
## rifaximin1 0.823 1.214 0.681 0.996
## beta_blockers1 0.806 1.241 0.684 0.949
## encephalopathy_admission1 1.534 0.652 1.267 1.859
## gi_bleed_admission1 1.130 0.885 0.962 1.328
## hcc_admission1 1.691 0.592 1.364 2.096
## alcoholic_hepatitis_admission1 0.770 1.298 0.607 0.977
## MELD_Na_baseline 1.019 0.981 1.004 1.036
## log_k_admit 0.749 1.336 0.518 1.083
## log_co2_admit 0.604 1.656 0.447 0.816
## log_alt_admit 1.132 0.884 1.029 1.244
## log_tb_admit 1.221 0.819 1.093 1.363
## log_alb_admit 0.494 2.026 0.359 0.679
## log_wbc_admit 1.229 0.814 1.056 1.430
## log_plt_admit 0.827 1.209 0.718 0.953
## sbp_admission 0.993 1.007 0.989 0.998
##
## Num. cases = 1949
## Pseudo Log-likelihood = -4875
## Pseudo likelihood ratio test = 764 on 35 df,
##
## Convergence: TRUEsfgBIC$fit ##
## Right-censored response of a competing.risks model
##
## No.Observations: 1949
##
## Pattern:
##
## Cause event right.censored
## 1 721 0
## 2 118 0
## unknown 0 1110
##
##
## Fine-Gray model: analysis of cause 1
##
## Competing Risks Regression
##
## Call:
## riskRegression::FGR(formula = Hist(time_90days, status_90days) ~
## age_admission + ethnicity2 + ethnicity3 + liver_transplant_listed1 +
## final_type_of_aki2 + final_type_of_aki3 + final_type_of_aki4 +
## final_type_of_aki5 + siteindiana + sitejacksonville +
## sitekentukey + sitemgh + siterochester + siteyale + etiology_cirrhosis2 +
## etiology_cirrhosis4 + etiology_cirrhosis6 + etiology_cirrhosis7 +
## albumin_given_admission1 + lactulose1 + rifaximin1 +
## beta_blockers1 + encephalopathy_admission1 + gi_bleed_admission1 +
## hcc_admission1 + alcoholic_hepatitis_admission1 + MELD_Na_baseline +
## log_k_admit + log_co2_admit + log_alt_admit + log_tb_admit +
## log_alb_admit + log_wbc_admit + log_plt_admit + sbp_admission,
## data = data, cause = cause)
##
## coef exp(coef) se(coef) z p-value
## age_admission 0.01851 1.019 0.00393 4.71 2.5e-06
## ethnicity2 0.35553 1.427 0.15193 2.34 1.9e-02
## ethnicity3 1.03336 2.810 0.23492 4.40 1.1e-05
## liver_transplant_listed1 -1.47320 0.229 0.15767 -9.34 0.0e+00
## final_type_of_aki2 0.90704 2.477 0.13337 6.80 1.0e-11
## final_type_of_aki3 0.93821 2.555 0.09594 9.78 0.0e+00
## final_type_of_aki4 0.36340 1.438 0.19583 1.86 6.4e-02
## final_type_of_aki5 0.85162 2.343 0.14862 5.73 1.0e-08
## siteindiana 0.24578 1.279 0.12825 1.92 5.5e-02
## sitejacksonville -0.66851 0.512 0.26501 -2.52 1.2e-02
## sitekentukey 0.34704 1.415 0.15488 2.24 2.5e-02
## sitemgh -0.18409 0.832 0.12382 -1.49 1.4e-01
## siterochester -0.25046 0.778 0.14313 -1.75 8.0e-02
## siteyale -0.23905 0.787 0.14374 -1.66 9.6e-02
## etiology_cirrhosis2 0.27871 1.321 0.12737 2.19 2.9e-02
## etiology_cirrhosis4 0.73499 2.085 0.28800 2.55 1.1e-02
## etiology_cirrhosis6 0.29045 1.337 0.14662 1.98 4.8e-02
## etiology_cirrhosis7 0.37367 1.453 0.13344 2.80 5.1e-03
## albumin_given_admission1 0.26088 1.298 0.11299 2.31 2.1e-02
## lactulose1 0.17588 1.192 0.09342 1.88 6.0e-02
## rifaximin1 -0.19423 0.823 0.09729 -2.00 4.6e-02
## beta_blockers1 -0.21601 0.806 0.08359 -2.58 9.8e-03
## encephalopathy_admission1 0.42819 1.534 0.09789 4.37 1.2e-05
## gi_bleed_admission1 0.12215 1.130 0.08228 1.48 1.4e-01
## hcc_admission1 0.52506 1.691 0.10961 4.79 1.7e-06
## alcoholic_hepatitis_admission1 -0.26099 0.770 0.12122 -2.15 3.1e-02
## MELD_Na_baseline 0.01929 1.019 0.00800 2.41 1.6e-02
## log_k_admit -0.28939 0.749 0.18832 -1.54 1.2e-01
## log_co2_admit -0.50422 0.604 0.15345 -3.29 1.0e-03
## log_alt_admit 0.12361 1.132 0.04836 2.56 1.1e-02
## log_tb_admit 0.19930 1.221 0.05647 3.53 4.2e-04
## log_alb_admit -0.70586 0.494 0.16268 -4.34 1.4e-05
## log_wbc_admit 0.20625 1.229 0.07740 2.66 7.7e-03
## log_plt_admit -0.18960 0.827 0.07207 -2.63 8.5e-03
## sbp_admission -0.00683 0.993 0.00240 -2.84 4.4e-03
##
## exp(coef) exp(-coef) 2.5% 97.5%
## age_admission 1.019 0.982 1.011 1.027
## ethnicity2 1.427 0.701 1.059 1.922
## ethnicity3 2.810 0.356 1.773 4.454
## liver_transplant_listed1 0.229 4.363 0.168 0.312
## final_type_of_aki2 2.477 0.404 1.907 3.217
## final_type_of_aki3 2.555 0.391 2.117 3.084
## final_type_of_aki4 1.438 0.695 0.980 2.111
## final_type_of_aki5 2.343 0.427 1.751 3.136
## siteindiana 1.279 0.782 0.994 1.644
## sitejacksonville 0.512 1.951 0.305 0.861
## sitekentukey 1.415 0.707 1.044 1.917
## sitemgh 0.832 1.202 0.653 1.060
## siterochester 0.778 1.285 0.588 1.031
## siteyale 0.787 1.270 0.594 1.044
## etiology_cirrhosis2 1.321 0.757 1.029 1.696
## etiology_cirrhosis4 2.085 0.480 1.186 3.667
## etiology_cirrhosis6 1.337 0.748 1.003 1.782
## etiology_cirrhosis7 1.453 0.688 1.119 1.887
## albumin_given_admission1 1.298 0.770 1.040 1.620
## lactulose1 1.192 0.839 0.993 1.432
## rifaximin1 0.823 1.214 0.681 0.996
## beta_blockers1 0.806 1.241 0.684 0.949
## encephalopathy_admission1 1.534 0.652 1.267 1.859
## gi_bleed_admission1 1.130 0.885 0.962 1.328
## hcc_admission1 1.691 0.592 1.364 2.096
## alcoholic_hepatitis_admission1 0.770 1.298 0.607 0.977
## MELD_Na_baseline 1.019 0.981 1.004 1.036
## log_k_admit 0.749 1.336 0.518 1.083
## log_co2_admit 0.604 1.656 0.447 0.816
## log_alt_admit 1.132 0.884 1.029 1.244
## log_tb_admit 1.221 0.819 1.093 1.363
## log_alb_admit 0.494 2.026 0.359 0.679
## log_wbc_admit 1.229 0.814 1.056 1.430
## log_plt_admit 0.827 1.209 0.718 0.953
## sbp_admission 0.993 1.007 0.989 0.998
##
## Num. cases = 1949
## Pseudo Log-likelihood = -4875
## Pseudo likelihood ratio test = 764 on 35 df,
##
## Convergence: TRUE