search for: modsf3

...) and x = 1 on (1,2) datnew <- data.frame(x=c(0,1,0,1), z=c(0,0,0,1), y=c(2,2,1,2), e=rep(0,4), id=c(1,2,3,3)) datnew # as expected modsf1 <- survfit(mod, newdata=datnew, id=id) modsf1 # not as expected modsf2 <- survfit(mod, newdata=datnew, id=id, censor=FALSE) modsf2 # for comparison modsf3 <- survfit(mod, newdata=datnew[1:2,]) modsf3 # appears to work when individual=FALSE (id not specified) modsf4 <- survfit(mod, newdata=datnew[1:2,], censor=FALSE) modsf4 # visually par(mfrow=c(2,2)) plot(modsf1, col=1:3, lty=1:3, conf.int=FALSE) plot(modsf2, col=1:3, lty=1:3, conf.int=FALSE...