library(survival)
data(flchain)
# complete data analysis.
flchain = flchain[complete.cases(flchain), ]
# short-hand variable names.
r = cbind(time = flchain$futime + 1, status = flchain$death)
s = flchain[, c("age", "sex")]
p = flchain[, c("sample.yr", "kappa", "lambda", "flc.grp", "creatinine", "mgus",
"chapter")]
if (FALSE) {
m = fgrrm(response = r, sensitive = s, predictors = p, family = "cox",
unfairness = 0.05)
summary(m)
}
Run the code above in your browser using DataLab