library(survival)
#
# 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])
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.