In Regression Modeling Strategies by Harrell (second edition) there is a section (S. 20.1.7) discussing Cox models including an interaction between a covariate whose main effect on survival we want to estimate as well (age in the example below) and a covariate whose main effect we do not want to estimate (gender in the example below).
Concretely: suppose that in a population the (unknown, true) hazard $h(t)$ follows the model
$$h(t) = \begin{cases} h_f(t) \exp(\beta_1 \textrm{age}), & \textrm{for female patiens} \\ h_m(t) \exp((\beta_1 + \beta_2) \textrm{age}), & \textrm{for male patiens} \end{cases} $$ where $h_f$, $h_m$ are unknown, true, not to be estimated baseline hazard functions and $\beta_1$, $\beta_2$ are unknown, true parameters to be estimated from the data.
(This example is taken almost literally from the book.)
Now Harrell remarks that the above situation can be rewritten as the stratified Cox model model 1:
$$h(t) = h_{\textrm{gender}}(t) \exp(\beta_1 \textrm{age} + \beta_2 X)$$ where the 'interaction term' $X$ equals zero for females and age for males. This is convenient because it means we can use the standard technique for estimating $\beta_1$ and $\beta_2$.
Now for the question. Suppose that two researchers A and B are given the same sample of patients drawn from the population described above. Researcher A fits model 1, obtaining estimates $\hat{\beta}_1$, $\hat{\beta}_2$ for the true parameters $\beta_1, \beta_2$ together with confidence intervals.
Researcher B takes the more naive approach of fitting two ordinary (i.e. unstratisfied) Cox-models: model 2a: $$h(t) = h_f(t)\exp(\gamma_1 \textrm{age})$$ on the female patients in the sample only and model 2b: $$h(t) = h_m(t)\exp(\gamma_2 \textrm{age})$$ on the male patients in the sample only. Thus obtaining estimates $\hat{\gamma_1}$, $\hat{\gamma_2}$ of the true parameters $\beta_1, \beta_1 + \beta_2$ respectively, together with confidence intervals.
Question:
- Are these estimates necessarily the same (in the sense that $\hat{\beta}_1 = \hat{\gamma}_1$, $\hat{\beta}_2 = \hat{\gamma}_2 - \hat{\gamma}_1$)? (Recall that both researchers look at the same data.)
- Are the confidence intervals necessarily the same?
- Does it make any sense to say that researcher A has a psychological advantage over researcher B in the case that $\beta_2 = 0$, because researcher A is then more likely to suspect that and switch to estimating the more parsimonious model $h(t) = h_{\textrm{gender}}(t)\exp(\beta_1 \textrm{age})$?