# A test of nesting.  It makes sure the model.frame is built correctly
tfun <- function(fit, mydata) {
    survfit(fit, newdata=mydata)

myfit <- coxph(Surv(time, status) ~ age + factor(sex), lung)

temp1 <- tfun(myfit, lung[1:5,])
temp2 <- survfit(myfit, lung[1:5,])
indx <- match('call', names(temp1))  #the call components won't match

all.equal(unclass(temp1)[-indx], unclass(temp2)[-indx])

Try the survival package in your browser

Any scripts or data that you put into this service are public.

survival documentation built on March 31, 2023, 7:41 p.m.