Has anyone tried plotting survival probabilities for a cox regression using Firth's penalized likelihood? I am using the coxphf() function in R. Any advice would be appreciated!
Thanks for this! Unfortunately, in the CRAN manual, they only provide an example of how to plot the penalized likelihood as compared to the wald.
I am interested in creating a graphic that illustrates survival probabilities over time (like you would for a regular cox ph regression). I would also love to figure out how to write code that predicts survival probabilities for new observations (i.e. at specified time points) based on the model. So, for example, I would like to pull the predicted survival for a specific covariate at 1 year, 2 years, etc. Does this make sense? Can you help? Any advice would be appreciated!