set.seed(23)
a<-aregImpute(formula1, data=somedata, n.impute=5)
f<-fit.mult.impute(Surv(ptime2, pstatus2) ~ formula2, cph, a, data=somedata #formula2 has the same terms as formula1, save 2 exceptions 1) ptime2 and pstatus2 are not in formula2 2) some continuous terms have rcs(*,3) functions around them.
p1<-plot(f,conf.int=T,ref.zero=T,fun=exp,age=c(18:70),xlab=c("Recipient age"),ylab=c("RR of Mortality"))
write.table(p1$x.xbeta,file=afilename,sep=",",dimnames.write=F)

set.seed(23)
a<-aregImpute(formula1, data=somedata, n.impute=5)
f<-fit.mult.impute(Surv(ptime2, pstatus2) ~ formula2, cph, a, data=somedata #formula2 has the same terms as formula1, save 2 exceptions 1) ptime2 and pstatus2 are not in formula2 2) some continuous terms have rcs(*,3) functions around them.
p2<-Predict(f ,age  ,conf.int=T ,ref.zero=T ,fun=exp )