[R] anova() method for coxph objects with robust standard errors

2009-11-23 Thread Mura Tamakou
Hi all, what is the reason why the anova() method for coxph objects does not work when robust standard errors have been requested, e.g., fit - coxph(Surv(futime, fustat) ~ resid.ds *rx + ecog.ps, data = ovarian, robust = T) anova(fit) any pointers will be much appreciated. Sincerely, Mura

[R] coxph() and survfit()

2009-10-23 Thread Mura Tamakou
Dear All, I have a question regarding the output of survfit() when I supply a Cox model. Lets say for example: library(survival) fit - coxph(Surv(time, status == 2) ~ factor(spiders), data = pbc) fit # HR for spiders is significant newdata - data.frame(spiders = factor(0:1)) sf -