Displaying 1 result from an estimated 1 matches for "qpred".
Did you mean:
pred
2008 Mar 03
1
Problem plotting curve on survival curve
...al curves using the predict function. Here is a
simple example:
> library(survival)
> tfit <- survreg(Surv(time, status) ~ factor(ph.ecog), data=lung)
> table(lung$ph.ecog)
0 1 2 3 <NA>
63 113 50 1 1
> tdata <- data.frame(ph.ecog=factor(0:3))
> qpred <- predict(tfit, newdata= tdata, type='quantile', p=1:99/100)
> matplot(t(qpred), 99:1/100, type='l')
The result of predict is a matrix with one row per group and one column per
quantile. The final plot uses "99:1" so as to show 1-F(t) = S(t) rather than F.
Don&...