Cox Regression on the Plane

技術要約：平面上のCox回帰

問題提起
Cox比例ハザードモデルは単変量生存解析の標準であるが、二変量生存データ（相関のある故障時間）への拡張は依然として限定的である。既存のアプローチには、以下のような重大なトレードオフが存在する：

• 共通脆弱性（Shared frailty）およびコピュラモデル： これらは多くの場合、依存構造に関するパラメトリックな仮定（特定のコピュラ族や脆弱性分布など）に依存する。もし依存構造が誤って指定された場合、結合生存レベルにおける共変量の効果は、解釈不能または不正確になる可能性がある。
• 乗法的ハザード拡張： 単一故障率と二重故障率の3つの個別の比例ハザードモデルを指定しようとする過去の試みは、結果として得られるハザード関数が必須の構造的制約を満たさないため、有効な二変量生存分布を導き出せないことが多い。

著者らは、Coxモデルの半パラメトリックな性質を維持し、有効な二変量生存分布をもたらし、かつコピュラや脆弱性分布の指定を必要とせずに透明性の高い共変量の解釈を提供する手法の欠如を指摘している。

手法
本論文では、生存関数の一般化レーマン表現に基づく2つの半パラメトリックな拡張を提案している。

1. モデル：

   • 単純レーマンモデル（Simple Lehmann Model）： 二変量生存関数を $S(t_1, t_2 | Z) = [S_0(t_1, t_2)]^{\exp(b^T Z)}$ とする直接的な拡張である。これは、両方の周辺分布および結合依存性に対して共通の比例効果を課すことで、ベースラインの依存構造を維持する。
   • 一般化レーマンモデル（Generalized Lehmann Model）： 2つの周辺ハザード（$\alpha, \beta$）と依存成分（$\gamma$）に対して異なる回帰パラメータを許容する、より柔軟な定式化である。これは次のように定義される：
     $$S(t_1, t_2 | Z) = [S_{0T1}(t_1)]^{\exp(\alpha^T Z)} [S_{0T2}(t_2)]^{\exp(\beta^T Z)} \left[ \frac{S_0(t_1, t_2)}{S_{0T1}(t_1)S_{0T2}(t_2)} \right]^{\exp(\gamma^T Z)}$$
     このモデルは、依存成分をベースラインの依存比に対する乗法的因子として扱う。著者らは、対応するハザード関数を導出し、結果として得られる関数が有効で非負の生存分布を定義するための十分条件（例：$\alpha$ および $\beta$ に対する $\gamma$ の制約）を確立している。

2. 推定：

   • 著者らは、ベースラインハザードが相殺されない二変量の設定において部分尤度の複雑さを避けるため、疑似観測値（pseudo-observations）アプローチを利用している。
   • 単純モデル： 二変量の疑似観測値を応答変数とし、補完対数対数（cloglog）リンクを用いたGLMを用いたGEE（一般化推定方程式）を通じて推定される。
   • 一般化モデル： 2ステップの手順を用いて推定される：
     1. 周辺パラメータ（$\alpha, \beta$）を、単変量の疑似観測値と標準的なGEEを用いて推定する。
     2. 依存成分（結合生存と周辺の積の比）のための修正された応答量を構築し、第2のGEEステップを用いて $\gamma$ を推定する。
   • 著者らは、独立な検閲条件下におけるこれらの推定量の一致性と漸近正規性を確立している。

主な結果

• シミュレーション研究： 著者らは、正の依存関係（PQD）および負の依存関係（NQD）を代表する様々なコピュラ（Frank, Clayton, Gumbel-Barnett）から生成されたデータに対して、提案手法を評価した。
  • 推定値は、サンプルサイズ $n=800$ において、低いバイアスと公称レベルに近い被覆率を示した。
  • 提案モデルは、単純レーマンモデルが誤指定されている場合（すなわち $\alpha \neq \beta \neq \gamma$ の場合）でも、柔軟なコピュラベースの回帰フレームワーク（Marra and Radice, 2020）と比較して競争力のある予測精度（平均絶対誤差）を維持した。
  • 手法は、レーマンモデルもコピュラモデルも厳密には正しくない対数正規分布から生成されたデータに対しても、モデル誤指定に対して頑健であることを示した。
• データ適用（Global Retinoblastoma Outcome Study - GROS）：
  • 本手法を、両眼の摘出までの時間を分析するために、977例の両側網膜芽細胞腫症例に適用した。
  • 依存構造： 非パラメトリックな推定は、両眼の間に負の四分位依存性（NQD）があることを示した。しかし、標準的なコピュラベースのアプローチは、正の依存（PQD）を意味するClaytonコピュラを選択し、その結果、矛盾する条件付き生存確率を導いた。
  • 解釈可能性： レーマンモデルは、パラメトリックなコピュラを強制することなくNQD構造を正しく捉えた。結果は、片方の眼が早期に摘出された場合、もう片方の眼がより長く生存する確率は、最初の眼が生存した場合よりも高くなることを示しており、これは臨床的な期待と一致しているが、コピュラモデルでは見逃されていた。
  • 共変量の効果： モデルは、提示時の年齢、腫瘍のステージ、および所得水準が生存確率に有意な影響を与える一方で、性別や家族歴は影響を与えないことを特定した。

意義と主張
本論文は、透明性が高く、完全に半パラメトリックなフレームワークを提供し、両方の周辺生存確率およびその依存性に対する共変量の効果を評価できると主張している。その主な貢献は以下の通りである：

1. コピュラの誤指定の回避： 共変量の効果をパラメトリックなコピラ族を通じてではなく、二変量生存関数に直接モデリングすることで、依存構造の誤った仮定のリスクを回避している。これは条件付き確率を解釈する上で極めて重要である。
2. 妥当性と解釈可能性： 提案されたモデルは、特定の条件下で有効な二変量生存分布を定義することが証明されており、回帰係数は生存スケール上で直接的な解釈が可能である（例：$\gamma$ は共変量がベースラインの依存性をどのように増幅または減衰させるかを制御する）。
3. 新規の推定法： 著者らは、二変量生存回帰のための初の二段階疑似観測値ベースの手法を提示し、その理論的特性を確立した。
4. 臨床的関連性： 網膜芽細胞腫データへの適用は、標準的なモデル選択基準（BICなど）が、誤った依存構造（PQD 対 NCD）を持つコピュラを好む可能性がある一方で、提案されたフレームワークは経験的証拠およびドメイン知識と一致する、臨床的に解釈可能な結果をもたらすことを示している。

著者らは、疑似観測値のための生存推定量の再計算に伴う計算コストや、独立な検閲の仮定を含む限界についても言及しているが、本手法が、依存構造が未知または複雑な二変量生存解析において、堅牢な代替案を提供することを強調している。