Hi, I am fitting an AFT model assuming a Weibull distribution and I would like to check the residuals compared to the Kaplan Meier residuals, but when I try to create the Kaplan Meier residuals I get an error: Time and status are different lengths.
I am using the following script: # Fitting the AFT model fit.weib <- survreg(Surv(TimeDeath, event) ~ age + sex + mutation + ventilation + BM1 + BM2, data = DF, dist = "weibull") fits <- fit.weib$linear.predictors resids <- (fit.weib$y[, 1] - fits)/fit.weib$scale resKM <- survfit(Surv(resids, event) ~ 1, data = DF) I get the error from the last line of the script. I tried some things that didn't work and I was wondering if anyone could help me. If you need more information please let me know. Thanks in advance, M [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.