> lines(survexp( ~ trt, ratetable=pfit.wtd, data=pbc), col='purple') Error in eval(expr, envir, enclos) : object 'albumin' not found
The code does work, however, using the latest version of survival. (Use of case weights in survexp languished on my to-do list for a long long time.) Terry Therneau ______________________________________________ R-help@r-project.org mailing list 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.