Guessing What Data Was Observed Based On The Model? Suppose I fit a Survival Cox-PH Regression Model in R and get the following results:
Call:
coxph(formula = Surv(time, status) ~ age + sex + ph.ecog, data = lung)

             coef exp(coef)  se(coef)      z        p
age      0.011067  1.011128  0.009267  1.194 0.232416
sex     -0.552612  0.575445  0.167739 -3.294 0.000986
ph.ecog  0.463728  1.589991  0.113577  4.083 4.45e-05

Likelihood ratio test=30.5  on 3 df, p=1.083e-06
n= 227, number of events= 164 
   (1 observation deleted due to missingness)

Based on these results, I can infer information such as:

*

*The number of observations

*The number of events

*The number of variables

*The estimate for the effect of each variable

My Question: Given this information, is it possible to simulate the covariate and response information for n = 227 such observations - such that if a similar Cox-PH model was fit to these newly simulated 227 observations, the resulting regression coefficients would approximately be equal to the original regression coefficients? Can I try to guess (and recreate)  observations might have been observed based on the regression model coefficients?
For example, I know that if I were to "fix" the covariate information for a group of n = 227 "arbitrary created" patients, I could then simulate their survival times (e.g. https://cran.r-project.org/web/packages/simsurv/index.html) - however, if I were to then fit a Cox-PH model to these observations, the model coefficients would not necessarily be close to the original model coefficients.
In general, is this possible to do? Only given the above model summary, could I try and somehow generate the original dataset that this model was trained on?
Thanks!
Note: I realize there are probably an infante number of n = 227 samples that can be randomly simulated such that a Cox-PH Model produces the same regression coefficient estimates as above.


In general, is this possible to do? Only given the above model summary, could I try and somehow generate the original dataset that this model was trained on?

You'll need additional information about the data to get the distribution of the covariates correct. Even then, you won't match the data exactly, but you'll get closer.
There are answer on Stack Exchange that explain how to simulate survival times from a Cox-PH model. See here, here, and here. Each of these answers use the method in Blender et al (2005) to simulate survival times: invert the cumulative hazard function and plug in uniform random variables.
Your simulation will include the correct effect of each covariate, but you won't be able to know the distribution of the covariates. You know the log hazard of sex is $-0.55$, but you don't know if this was estimated from 10, 50, 100, 200, etc men. You won't know the distribution of age either. Was the study restricted to 50-70 year olds? What were the ages of each person. You can usually find a demographics table that lists the distributions of important variables. If you include this in your simulations, you'll be able to more accurately reproduce survival times.
This additional information will only get you so far. You'll never be able to reproduce the original data exactly. Like you said, there are many datasets that will produce identical or nearly identical parameters.
