How to extract hazard functions

I need to predict hazard functions with confidence intervals and I have tried many different ways without success. I have created a graph but I need the functions for each spline in numbers to be able to redo the graph in another program (SAS). Do someone have any idea?
The command I have used for the graph is the following:

rcsgen maxerystorage, df(5) gen(maxrcs) center(14)
stcox maxrcs* rbcspl* i.year sex agespl* i.indication previoustransfusion abo i.wdy crbcspl1 sexagespl1 Country#year platelets plasma, strata(hospital)
partpred hr_maxerystorage, for(maxrcs*) eform ci(hr_maxerystorage_lci3 hr_maxerystorage_uci3)

twoway (rarea hr_maxerystorage_lci hr_maxerystorage_uci maxerystorage, sort pstyle(ci)) ///
(line hr_maxerystorage maxerystorage, sort) ///
, yscale(log) ylabel(0.5 (0.1) 2)
, legend(off) xtitle(maxerystorage) ytitle(HR),///