calculating martingale residual on new data using "predict.coxph"

Hi list,

I was trying to use "predict.coxph" to calculate martingale residuals on a test 
data, however, as pointed out before

http://tolstoy.newcastle.edu.au/R/e4/help/08/06/13508.html

predict(mycox1, newdata, type="expected") is not implemented yet.  Dieter 
suggested to use 'cph' and 'predict.Design', but from my reading so far, I'm not 
sure they can do that.

Do you other ways to calculate martingale residuals on a new data?

Thank you very much!

...Tao