Re: [R] expected risk from coxph (survival)

From: Dieter Menne <>
Date: Sat, 07 Jun 2008 08:15:48 +0000 (UTC)

Reid Tingley <r_tingley <at>> writes:
> When I try to to obtain the expected risk for a new dataset using coxph in the
survival package I get an error.
> Using the example from ?coxph:

# Example rewritten by DM; please do not use HTML mail library(survival)

test1 <- list(time=  c(4, 3,1,1,2,2,3),
            x=     c(0, 2,1,1,1,0,0),
            sex=   c(0,0,0,0,1,1,1))

cox<-coxph( Surv(time, status) ~ x + strata(sex), test1) #stratified model
new<-list(time=  c(5, 1,1,2,2,4,3),
          x=     c(0, 2,1,1,1,0,0),
          sex=   c(0,0,0,0,1,1,1))


 # }

#    else if (type == "expected") {
#        if (missing(newdata)) 
#            pred <- y[, ncol(y)] - object$residuals
#        else stop("Method not yet finished")
Looks like this is "by design"; see the code above.

You might try to use cph and predict.Design from Frank Harrell's Design package instead.

Dieter mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Sat 07 Jun 2008 - 08:19:47 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 07 Jun 2008 - 21:31:00 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive