L'ajustement du modèle de Cox avec des strates et l'interaction strate-covariable diffère-t-il de l'ajustement de deux modèles de Cox?

13

Dans Regression Modeling Strategies de Harrell (deuxième édition), il y a une section (S. 20.1.7) traitant des modèles de Cox, y compris une interaction entre une covariable dont nous voulons également estimer l'effet principal sur la survie (âge dans l'exemple ci-dessous) et un covariable dont nous ne voulons pas estimer l'effet principal (sexe dans l'exemple ci-dessous).

Concrètement: supposons que dans une population le danger (inconnu, vrai) h(t) suit le modèle

h(t)={hf(t)exp(β1age),for female patienshm(t)exp((β1+β2)age),for male patiens
hf ,hm sont inconnus, vrai, à ne pas estimer les fonctions de risque de base etβ1 ,β2 sont des paramètres réels inconnus à estimer à partir des données.

(Cet exemple est tiré presque littéralement du livre.)

Harrell remarque maintenant que la situation ci-dessus peut être réécrite sous la forme du modèle de modèle de Cox stratifié 1 :

h(t)=hgender(t)exp(β1age+β2X)
où le «terme d'interaction»X est égal à zéro pour les femmes et à l'âge pour les hommes. C'est pratique car cela signifie que nous pouvons utiliser la technique standard pour estimerβ1 etβ2 .

Maintenant pour la question. Supposons que deux chercheurs A et B reçoivent le même échantillon de patients issus de la population décrite ci-dessus. Chercheur Un modèle unique 1, obtenant des estimations ß 1 , β 2 pour les vrais paramètres ß 1 , β 2 avec des intervalles de confiance.β^1β^2β1,β2

Le chercheur B adopte l'approche la plus naïve de l'ajustement de deux modèles de Cox ordinaires (c'est-à-dire non satisfaits): modèle 2a:

h(t)=hf(t)exp(γ1age)
sur les patientes de l'échantillon uniquement et modèle 2b:
h(t)=hm(t)exp(γ2age)
sur les patients masculins de l'échantillon uniquement. On obtient ainsi des estimations γ1^ , γ2^des vrais paramètres β1,β1+β2 respectivement, ainsi que des intervalles de confiance.

Question:

  • Sont ces estimations nécessairement les mêmes (dans le sens où β 1 = γ 1 , B 2 = γ 2 - γ 1 )? (Rappelons que les deux chercheurs regardent les mêmes données.)β^1=γ^1β^2=γ^2γ^1
  • Les intervalles de confiance sont-ils nécessairement les mêmes?
  • Est-il logique de dire que le chercheur A a un avantage psychologique sur le chercheur B dans le cas où β2=0 , car le chercheur A est alors plus susceptible de soupçonner cela et de passer à l'estimation du modèle le plus parcimonieux h(t)=hgender(t)exp(β1age) ?
Vincent
la source

Réponses:

4

YM=αM+βMageYF=αF+βFage summarized by: Y=λ+λFF+γage+γFFage, so that you could directly estimate the gender difference both in intercept and in slope. In fact: αM=λ,βM=γ,αFαM=λF,βFβM=γF. In that case, I agree with you that the unique model would allow to have an immediate idea on the gender difference (given by the interaction parameters, λF, since the slope difference has a clearer interpretation, and your question refers to that). However, with the Cox model things are different. First of all, if we don't include gender in the regression there may be a reason, i.e. that it doesn not fulfill the proportional hazard assumption. Also, if we build a unique model with gender as an interaction term, we are assuming a common baseline hazard function (unless I misunderstood the meaning of hgender(t)), while the two-separate-models approach allow for two separate baseline hazard functions, thus different models are implied.

See, for example, the Chapter "Survival Analysis" from Kleinbaum and Klein, 2012, Part of the series Statistics for Biology and Health.

Federico Tedeschi
la source