--- begin inclusion -- Hi everybody. I'm trying to fit a weibull survival model with a spline basis for the predictor, using the survival library. I've noticed that it doesn't seem to be possible to use the aic method to choose the degrees of freedom for the spline basis in a parametric regression (although it's fine with the cox model, or if the degrees of freedom are specified directly by the user), and I was wondering if there is some reason for this? ... --- end inclusion ---
A simple reason as it turns out: there was an incorrect variable name in the survpenal.fit function. A "variable not found" error is exactly what one would expect. I've now repaired the function for version 2.36.10 (which I hope to post shortly). The big surprise is that it has been wrong for about a decade, and no one had ever tried a situation that exercised that particular line of code -- or at least they never told me about the error. Thank you for the very clear description of the problem, which made it easy for me to track it down. 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.