Displaying 1 result from an estimated 1 matches for "pptest2".
Did you mean:
kptest2
2013 Jan 24
1
predicted HR in coxph with psline
...nce value
for the predictions?
I did some test code below:
###
par(mfrow = c(1,3))
test = coxph(Surv(time, status) ~ age, cancer)
test2 = coxph(Surv(time, status) ~ pspline(age, df=2), cancer)
test3 = coxph(Surv(time, status) ~ pspline(age, df=5), cancer)
pptest3=predict(test3, type="lp")
pptest2= predict(test2, type="lp")
pptest = predict(test, type="lp")
plot(cancer$age, exp(pptest))
abline(v=mean(cancer$age))
abline(h=1)
plot(cancer$age, exp(pptest2))
abline(v=mean(cancer$age))
abline(h=1)
plot(cancer$age, exp(pptest3))
abline(v=mean(cancer$age))
abline(h=1)
In the f...