I have been plotting my data with the rms package in R:

fit

Similar questions and discussions