O ajuste do modelo Cox com estratos e interação estrato-covariável difere do ajuste de dois modelos Cox?

13

Em Estratégias de modelagem de regressão, de Harrell (segunda edição), há uma seção (S. 20.1.7) discutindo modelos de Cox, incluindo uma interação entre uma covariável cujo principal efeito sobre a sobrevivência também queremos estimar (idade no exemplo abaixo) e uma covariável cujo principal efeito não queremos estimar (gênero no exemplo abaixo).

Concretamente: suponha que em uma população o risco (desconhecido, verdadeiro) h(t) siga o modelo

h(t)={hf(t)exp(β1age),for female patienshm(t)exp((β1+β2)age),for male patiens
, ondehf ,hm são desconhecidas, verdadeiro, não deve ser estimado as funções de risco de linha de base eβ1 ,β2 são desconhecidos, parâmetros verdadeiros a serem estimados a partir dos dados.

(Este exemplo é retirado quase literalmente do livro.)

Agora Harrell observa que a situação acima pode ser reescrita como o modelo estratificado modelo Cox 1 :

h(t)=hgender(t)exp(β1age+β2X)
que o "termo de interação"X é igual a zero para mulheres e idade para homens. Isso é conveniente porque significa que podemos usar a técnica padrão para estimarβ1 eβ2 .

Agora para a pergunta. Suponha que dois pesquisadores A e B recebam a mesma amostra de pacientes da população descrita acima. Pesquisador Ajusta um modelo 1, as estimativas obtenção p 1 , β 2 para os verdadeiros parâmetros p 1 ,β^1β^2β1,β2 , juntamente com intervalos de confiança.

Pesquisador B leva a abordagem mais ingénuo de encaixe de dois (isto é, unstratisfied) COX-modelos comuns: modelo 2a:

h(t)=hf(t)exp(γ1age)
nos pacientes do sexo feminino na amostra única e modelo 2b:
h(t)=hm(t)exp(γ2age)
em pacientes do sexo masculino em apenas a amostra. Assim, obtendo estimativas γ1^ , γ2^dos parâmetros verdadeiros β1,β1+β2 respectivamente, juntamente com intervalos de confiança.

Questão:

  • São estas estimativas necessariamente o mesmo (no sentido de que β 1 = γ 1 , p 2 = γ 2 - γ 1 )? (Lembre-se de que os dois pesquisadores analisam os mesmos dados.)β^1=γ^1β^2=γ^2γ^1
  • Os intervalos de confiança são necessariamente os mesmos?
  • Faz algum sentido dizer que o pesquisador A tem uma vantagem psicológica sobre o pesquisador B no caso de β2=0 , porque é mais provável que o pesquisador A suspeite disso e mude para estimar o modelo mais parcimonioso h(t)=hgender(t)exp(β1age) ?
Vincent
fonte

Respostas:

4

YM=αM+βMumage, YF=αF+βFumage resumido por: Y=λ+λFF+γumage+γFFumage, para que você possa estimar diretamente a diferença de gênero na interceptação e na inclinação. De fato:α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
fonte