You have (unimportant lines omitted) > c2= survdiff(Surv(act.surv.time,censoring)~treatgrp ,data=b) > plot(c2) The problem is that you are using the wrong function. It is survfit that creates plottable survival curves, survdiff only does the log-rank test. Terry Therneau