Apakah pemasangan model Cox dengan interaksi strata dan strata-kovariat berbeda dari pemasangan dua model Cox?

13

Dalam Regresi Modeling Strategies oleh Harrell (edisi kedua) ada bagian (S. 20.1.7) membahas model Cox termasuk interaksi antara kovariat yang efek utamanya pada kelangsungan hidup yang ingin kita perkirakan juga (umur dalam contoh di bawah) dan kovariat yang efek utamanya tidak ingin kami perkirakan (jenis kelamin dalam contoh di bawah).

Secara konkret: anggaplah bahwa dalam suatu populasi bahaya (tidak diketahui, benar) h(t) mengikuti model

h(t)={hf(t)exp(β1age),for female patienshm(t)exp((β1+β2)age),for male patiens
di manahf ,hm tidak diketahui, benar, tidak dapat diperkirakan fungsi bahaya awal danβ1 ,β2 tidak diketahui, parameter sebenarnya yang diperkirakan dari data.

(Contoh ini diambil hampir secara harfiah dari buku ini.)

Sekarang Harrell menyatakan bahwa situasi di atas dapat ditulis ulang sebagai model model Cox bertingkat 1 :

h(t)=hgender(t)exp(β1age+β2X)
mana 'istilah interaksi'X sama dengan nol untuk wanita dan usia untuk pria. Ini nyaman karena itu berarti kita dapat menggunakan teknik standar untuk memperkirakanβ1 danβ2 .

Sekarang untuk pertanyaannya. Misalkan dua peneliti A dan B diberi sampel yang sama dari pasien yang diambil dari populasi yang dijelaskan di atas. Peneliti A cocok model 1, memperoleh perkiraan ß 1 , β 2 untuk parameter benar ß 1 , β 2 bersama-sama dengan interval kepercayaan.β^1β^2β1,β2

Peneliti B mengambil pendekatan yang lebih naif pas dua (yaitu unstratisfied) Cox-model biasa: Model 2a:

h(t)=hf(t)exp(γ1age)
pada pasien wanita dalam sampel saja dan Model 2b:
h(t)=hm(t)exp(γ2age)
pada pasien pria dalam sampel saja. Dengan demikian memperoleh estimasi γ1^ , γ2^dari parameter benar β1,β1+β2 masing-masing, bersama dengan interval kepercayaan.

Pertanyaan:

  • Apakah perkiraan ini tentu sama (dalam arti bahwa β 1 = γ 1 , ß 2 = γ 2 - γ 1 )? (Ingat bahwa kedua peneliti melihat data yang sama.)β^1=γ^1β^2=γ^2γ^1
  • Apakah interval kepercayaan harus sama?
  • Apakah masuk akal untuk mengatakan bahwa peneliti A memiliki keunggulan psikologis dibandingkan peneliti B dalam hal β2=0 , karena peneliti A kemudian lebih cenderung mencurigai hal itu dan beralih untuk memperkirakan model yang lebih pelit h(t)=hgender(t)exp(β1age) ?
Vincent
sumber

Jawaban:

4

YM.=αM.+βM.SebuahgeYF=αF+βFSebuahgeY=λ+λFF+γSebuahge+γFFSebuahgeαM.=λ,βM.=γ,αF-αM.=λF,βF-βM.=γF. Dalam hal itu, saya setuju dengan Anda bahwa model unik akan memungkinkan untuk memiliki ide langsung tentang perbedaan gender (diberikan oleh parameter interaksi,λF, karena perbedaan kemiringan memiliki interpretasi yang lebih jelas, dan pertanyaan Anda merujuk pada hal itu). Namun, dengan model Cox segalanya berbeda. Pertama-tama, jika kita tidak memasukkan gender dalam regresi mungkin ada alasan, yaitu bahwa hal itu tidak memenuhi asumsi bahaya proporsional. Juga, jika kita membangun model unik dengan gender sebagai istilah interaksi, kita mengasumsikan fungsi dasar bahaya yang umum (kecuali saya salah memahami arti darihjenis kelamin(t)), sementara pendekatan dua model terpisah memungkinkan untuk dua fungsi bahaya dasar yang terpisah, dengan demikian model yang berbeda tersirat.

Lihat, misalnya, Bab "Analisis Kelangsungan Hidup" dari Kleinbaum dan Klein, 2012, Bagian dari seri Statistik untuk Biologi dan Kesehatan.

Federico Tedeschi
sumber