層および層と共変量の相互作用を持つCoxモデルの適合は、2つのCoxモデルの適合とは異なりますか?


13

回帰モデリング戦略ハレル(第2版)によってセクション(S. 20.1.7)主効果生存に対する我々は(以下の例では年齢)も推定する共変量の間の相互作用を含むコックスモデルを検討し、あります主効果を推定したくない共変量(下の例では性別)。

具体的には、母集団では(未知、真)ハザードh(t)がモデルに従うと仮定します

h(t)={hf(t)exp(β1age),for female patienshm(t)exp((β1+β2)age),for male patiens
hfhm未知です、真の、ベースラインハザード関数と推定されるべきではなくβ1β2 未知の、データから推定される真のパラメーターです。

(この例はほとんど文字通り本から取られています。)

ハレルは、上記の状況を成層コックスモデルモデル1として書き直すことができると述べています。

h(t)=hgender(t)exp(β1age+β2X)
'対話用語'X男性女性、年齢はゼロに等しいです。それは我々が推定するための標準的な技術を使用できることを意味しますので、これは便利ですβ1およびβ2

さて、質問です。2人の研究者AとBに、上記の母集団から抽出された患者の同じサンプルが与えられたとします。研究員Aフィットモデル1、取得推定値はβ 1β 2β^1β^2、真のパラメータがためにβ1,β2信頼区間と一緒に。

研究者Bは、二つの通常の(すなわちunstratisfied)コックス・モデルをフィッティングのより単純なアプローチをとる:モデルA:

h(t)=hf(t)exp(γ1age)
のみの試料との女性患者にモデル2B:
h(t)=hm(t)exp(γ2age)
のみ試料中の男性患者に。こうして推定取得γ1^γ2^真のパラメータをβ1,β1+β2と共に信頼区間で、それぞれ。

質問:

  • これらの推定値は(という意味では必ずしも同じβ 1 = γ 1β 2 = γ 2 - γ 1)?(両方の研究者が同じデータを見ていることを思い出してください。)β^1=γ^1β^2=γ^2γ^1
  • 信頼区間は必ずしも同じですか?
  • それは研究者Aは、ケースには研究者Bの上に心理的な利点を持っていると言っては何の意味も持たないことにβ2=0、研究者Aは、そのスイッチより倹約モデルの推定に疑い、より多くの可能性が高いためh(t)=hgender(t)exp(β1age)

回答:


4

各パラメーターを推定する必要があるモデル(通常の最小二乗など)では、2つの別々のモデルが相互作用項を持つ単一モデルと同じ推定値を持つ状況を作成できます。例えば、我々は持っている可能性があり:Y F = α F + β F * Gの電子により要約:Y = λ + λ F * F + γ * グラムのE + γYM=αM+βMageYF=αF+βFage λ F。これにより、切片と勾配の両方で性差を直接推定できます。実際には: α M = λ β M = γ α F - α M = λ Fβ F - β M = γ F。その場合、ユニークなモデルにより、性別の違いについてすぐに考えられるようになることに同意します(相互作用パラメーター、Y=λ+λFF+γage+γFFageαM=λ,βM=γ,αFαM=λF,βFβM=γFλF、勾配の違いはより明確な解釈を持ち、あなたの質問はそれを参照しているため)。ただし、Coxモデルでは状況が異なります。まず、回帰に性別を含めない場合、理由があります。つまり、性別が比例ハザードの仮定を満たさないということです。また、性別を相互作用用語とする一意のモデルを構築する場合、共通のベースラインハザード関数を想定しています(意味を誤解していない限り)。ハザード関数、したがって異なるモデルが暗示されています。 hgender(t)

たとえば、Kleinbaum and Klein、2012年、「Statistics for Biology and Health」シリーズの一部の「Survival Analysis」の章を参照してください。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.