のcoxphモデルの要約で与えられた「

のcoxphモデルの要約で与えられる値は何ですか？例えば、 $R^2$

Rsquare= 0.186   (max possible= 0.991 )

私は愚かにも原稿を値として含めましたが、レビュアーはCoxモデル用に開発された古典的な線形回帰からの統計の類似体に気づいていないと言いました。参照。どんな助けも素晴らしいでしょう！ $R^2$ $R^2$

r survival r-squared cox-model

— ダニエルズブリューワー
ソース

の概念が従来の線形回帰を超えて拡張されるほとんどの状況では、観測値とモデルの下で予測される値との相関の二乗相関です。おそらくここで適用できますか？

R^{2}

$R^2$

— マクロ

いいえ、それには関係ありません。

— フランクハレル

回答:

を使用しgetS3method("summary","coxph")て、計算方法を確認できます。

関連するコード行は次のとおりです。

logtest <- -2 * (cox$loglik[1] - cox$loglik[2])
rval$rsq <- c(rsq = 1 - exp(-logtest/cox$n), maxrsq = 1 - 
        exp(2 * cox$loglik[1]/cox$n))

これcox$loglikは、「初期値と係数の最終値を持つ対数尤度を含む長さ2のベクトル」（「参考文献」を参照?coxph.object）でcox$nあり、「近似に使用される観測値の数」です。

— ローランド
ソース

誤解しない限り、それはCox＆Snellの疑似R-2乗です。さまざまな疑似R-squaredの説明と比較については、ats.ucla.edu /stat/mult_pkg/faq/general/psuedo_rsquareds.htmを参照してください。

— ワンストップ

$n$ coxph

— シュ・ロンフイ
ソース

間違っていますが、奇妙に聞こえても、観測数で割ります。元の質問に対して、レビュー担当者がCoxモデルで20年間存在していたものに気付かないのは奇妙です。

— フランクハレル

Ronghui Xuと@Frank Harrellの間のやり取りに加えて、観測数で「奇妙に聞こえる」だけでなく、機能しません。これを確認するには、研究1が研究2の2倍の打ち切り率を持っているという事実を除いて、おおよそE（R2）= 0.5、同じ共変量分布、つまりすべて同じであるように、ベータをある値に固定することを検討してください同じ母集団の量を推定する必要がありますが、研究1のR2推定値は、サンプルサイズに関係なく、研究2の推定値の約半分になります。0.5ではなく、0.25になります。

R^{2}

$R^2$

フランクの発言に対する答えとして、これは簡単ではなく、フランクのヌル対数尤度に関する観察が正しいことに同意します。私は、この量を、情報ゲインに基づいて明確に定義された母集団量の一貫した推定量への近似としてのみ見ました。Ronghui Xuが言及した論文は、シミュレーションを実行します。これらは、存在しないわけではないが、観測の総数ではなく障害の数で割ると、打ち切りの影響がはるかに弱くなることを示しています。

R^{2}

$R^2$