--- begin included message ---
I am experimenting with the function predict() in two versions of R and
the R extension package "survival".
library(survival)
set.seed(123)
testdat=data.frame(otime=rexp(10),event=rep(0:1,each=5),x=rnorm(10))
testfm=as.formula('Surv(otime,event)~x')
testfun=function(dat,fm)
{
predict(coxph(fm,data=dat),type='lp',newdata=dat)
}
-- end inclusion ----
The question was: this works under survival 2.35-8, but not survival
2.36-5
Answer: The predict and underlying model.frame functions for coxph were
brought into line with lm and other models. The major advantage is that
I now deal with factors and data dependent predictors like ns()
correctly.
You've shown a disadvantage of which I was not aware. Using your
example but replacing coxph() by lm() with otime ~x as the model I get a
similar failure. I'd like to ask a wider audience of R-devel since it
is bigger than coxph.
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.
