の例 計測プロセスフォームでのCoxモデルの適合

医学研究者は、急性白血病の治療として骨髄移植からの回復の成功率を決定したいと考えています。回復は リスクカテゴリ 、移植時の患者、その患者 病気の段階、および血小板数が正常なレベルに戻ったかどうかなどの要因に依存する。リスクカテゴリ 研究 病気の段階 全体を通して変化しないため、固定予測変数です。しかし、患者の血小板数は、回復プロセス中にカウントが変化する可能性があるため、時間依存の予測値です。

医学研究者は、移植を受けた後に137人の患者を研究し、彼らが無病である日数を記録します。血小板数が正常に戻る前に死亡した場合、または血小板数が正常に戻った後に白血病が戻った場合、患者は病気を起こすものではありません。の値は はい 、無病患者を示し、検閲観察である。打ち切られた観測は、観測時間の終わりまでにイベントが発生しない場合です。

データは棚卸プロセスフォームに入っているため、複数の行が各患者を表します。各行は、すべての変数の値が一定である時間間隔を表します。時間依存の予測変数は、行間で変化します。間隔は開始時刻の直後に開始され、終了時刻が含まれます。

たとえば、次の表には、1 の患者のデータが含まれています 識別 。観測値との観測値 リスクカテゴリ病気の段階 は、各行で同じです。通常の血小板数は研究中に変化する可能性があるため、各患者は、この予測変数が変化するたびに新しいデータ行を必要とします。最初の行は、患者が移植後13日後の最初の間隔で正常な血小板数を有していなかったことを示す。2行目は、患者が13日目から2,081日目の研究終了までの正常な血小板数を持っていたことを示しています。

識別 リスクカテゴリ 開始時刻 終了時間 病気のない 通常の血小板 病気の段階
1 1 0 13 はい いいえ ノーマル
1 1 13 2081 はい はい ノーマル

これらのデータは、クラインとモーシュバーガー(2003)1.

  1. サンプルデータ、骨髄.MTWを開きます。
  2. 統計 > 信頼性/生存時間 > Cox 回帰 > カウントプロセスフォームでのコックスモデルの適合 を選択します。
  3. 開始時間にを入力します。終了時間にを入力します。
  4. 打ち切り列 (オプション)にを入力します。打ち切りの値はいを選択します。
  5. カテゴリ予測変数 に、リスクカテゴリ通常の血小板、および 病気の段階 と入力します。
  6. オプション を選択します。ケース ID (対象残差の場合)に, 識別 を入力します。
  7. 各ダイアログボックスの OK をクリックします。

結果を解釈する

まず、研究者は適合度検定を使用して、モデルの全体的な適合度を評価します。3つのテストのp値はすべて0.05以下なので、研究者はモデルがデータによく適合すると結論付けます。その後、研究者はANOVA表を使用して、個々の用語の効果を評価します。回帰表では、デザインと重量のp値が0.05のα水準で有意です。したがって、医学研究者は、移植時の患者のリスクカテゴリー、疾患段階、および血小板数が正常レベルであるかどうかは、患者が骨髄移植から回復するかどうかに統計的に有意な影響を及ぼすと結論付けている。

研究者は、カテゴリ予測変数の相対リスク表を使用して、予測変数の異なるレベル間のリスクを評価します。例えば、正常血小板を有する患者の死亡または白血病の再発のリスクは、正常な血小板のない患者よりも0.37倍低い。さらに、信頼区間は、正常な血小板を有する患者の死亡または再発の真のリスクは、95%の信頼レベルで、正常血小板のない患者のリスクよりもわずか0.19倍または0.7倍も低くなる可能性があることを示している。信頼区間には1が含まれていないため、正常な血小板の有無にかかわらず死亡または再発のリスクの差は統計的に有意である。

方法

Coxモデルタイプ計測プロセスフォーム
カテゴリ予測変数のコーディング(1, 0)
同順位の調整エフロン(Efron)
未使用の行1

打ち切り情報

打ち切られていない単位打ち切られた単位合計打ち切られたパーセント
8317325667.58%
打ち切り値: 病気のない = はい

回帰式

リスクスコア=0.0 リスクカテゴリ_1 - 0.793 リスクカテゴリ_2 - 0.033 リスクカテゴリ_3 + 0.0 通常の血小板_いいえ - 1.004 通常の血小板_はい
+ 0.0 病気の段階_ハイリスク - 0.696 病気の段階_ノーマル

係数

係数係数の標準誤差Z-値p値
リスクカテゴリ       
  2-0.7930.321-2.470.014
  3-0.0330.325-0.100.919
通常の血小板       
  はい-1.0040.332-3.020.003
病気の段階       
  ノーマル-0.6960.275-2.530.011

カテゴリ予測変数の相対的リスク

水準A水準B相対的リスク95%信頼区間
リスクカテゴリ     
  210.4524(0.2409, 0.8495)
  310.9673(0.5116, 1.8290)
  322.1383(1.2487, 3.6616)
通常の血小板     
  はいいいえ0.3666(0.1912, 0.7029)
病気の段階     
  ノーマルハイリスク0.4986(0.2909, 0.8547)
水準Bと比べた水準Aのリスク

モデル要約

モデル対数尤度R二乗AICAICc(修正済み
赤池情報量基準)
BIC(ベイズ
情報量基準)
項なし-373.30746.59746.59746.59
項付き-358.6011.47%725.20725.71734.88

適合度検定

検定自由度カイ二乗p値
尤度比429.390.000
ワルド(Wald)432.470.000
スコア435.220.000

分散分析



ワルド検定
要因自由度カイ二乗p値
リスクカテゴリ29.770.008
通常の血小板19.130.003
病気の段階16.410.011
1 Klein, J.P. & Moeschberger, M.L.コペランからの公開データセットに基づいて適応されました。(2003). Semiparametric proportional hazards regression with fixed covariates. Survival Analysis: Techniques for Censored and Truncated Data (2nd ed., pp. 243-293). Springer