You report an error message: > plot(zph.revasFit[1]) Error in plot.window(xlim, ylim, log, asp, ...) : need finite 'ylim' values I have never seen this error before, and I cannot guess what causes it. You need to provide more information, and likely a small data set that produces the problem. Perhaps you have an x variable that is a constant? Terry Therneau