powered by
Semiparametric Proportional odds model, that has the advantage that $$ logit(S(t|x)) = \log(\Lambda(t)) + x \beta $$ so covariate effects give OR of survival.
logitSurv(formula, data, offset = NULL, weights = NULL, ...)
formula with 'Surv' outcome (see coxph)
coxph
data frame
offsets for exp(x beta) terms
weights for score equations
Additional arguments to lower level funtions
This is equivalent to using a hazards model $$ Z \lambda(t) \exp(x \beta) $$ where Z is gamma distributed with mean and variance 1.
# NOT RUN { data(TRACE) dcut(TRACE) <- ~. out1 <- logitSurv(Surv(time,status==9)~vf+chf+strata(wmicat.4),data=TRACE) summary(out1) gof(out1) bplot(out1) # }
Run the code above in your browser using DataLab