survreg
Object
survreg
function is deprecated; please use
instead.
Predicted values for a
survreg
object
predict.survreg(object, newdata, type="response" se.fit=F, terms=labels.lm(object), p=c(0.1, 0.9), ripley=F)
survreg
function.
"response"
),
the linear predictor (
"linear"
, with
"lp"
as an allowed abbreviation),
a predicted quantile on the original scale of the data (
"quantile"
),
a quantile on the linear predictor scale (
"uquantile"
),
or the matrix of terms for the linear predictor ("
terms"
).
At this time
"link"
and linear predictor (
"lp"
) are identical.
TRUE
, include the standard errors of the prediction in the result.
terms
is a matrix with
one column for every term (excluding the intercept) in the model.
Escobar, L. A. and Meeker, W. Q. (1992). Assessing influence in regression analysis with censored data. Biometrics 48, 507-528.
# Draw figure 1 from Escobar and Meeker fit <- survreg(Surv(time,status) ~ age + age^2, data=stanford2, dist='lognormal') plot(stanford2$age, stanford2$time, xlab='Age', ylab='Days', xlim=c(0,65), ylim=c(.01, 10^6), log='y') pred <- predict(fit, newdata=list(age=1:65), type='quantile', p=c(.1, .5, .9)) matlines(1:65, pred, lty=c(2,1,2), col=1)