Skip to content

Kaplan Meier Post Hoc?

1 message · Terry Therneau

#
No, you didn't miss anything in the survival package.  I've never found post-hoc tests 
interesting so have little motivation to add such (and a very long "to do" list of things 
I would like to add).

If you simply must have them, why not do all pairwise tests?

chisq <- matrix(0., 4,4)
for (i in 1:4) {
     for (j in (1:4)[-i]) {
          temp <- survdiff(Surv(time, status) ~ group, data=mydata,
                              subset=(group %in% (unique(group))[c(i,j)]))
          chisq[i,j] <- temp$chisq
          }
     }

Terry Therneau
On 10/25/2012 05:00 AM, r-help-request at r-project.org wrote: