r - coxph オブジェクトの虚弱推定

Question

obj=coxph(... + frailty(id) ) を使用すると、オブジェクトは各個人の (log) frailty 推定値も返します。これは obj$frail で抽出できます。

これらの見積もりがどのように得られているか知っている人はいますか? それらは経験的ベイズ推定値ですか?

ありがとう！

テオドール

Answer 1

のデフォルトの分布frailtyは、?frailtyページで「ガンマ」であることがわかります。frailty 関数 (非表示ではない) を見ると、ディストリビューションの名前が単に "frailty" に貼り付けられていることがわかります。get()適切な関数を取得するために使用します。frailty.gammaあなたの質問への答えを見つけるために（隠されていない）を見てください。ヘルプ ページをもう一度見直すと、ページの一番上にコードが表示されているので、コードを見なくてもすべてを理解できたはずであることがわかります。しかし、R を使用して知識を得るには多くの方法があります (それらは「経験的ベイズ」推定ではなく ML です)。

ヘルプ ページは、著者 (Therneau) が、コードを読んでも明らかでない詳細については、Therneau と Grambsch に相談することを期待していることを示唆しています。R でサバイバル モデルを扱う真面目な仕事をしているなら、持っておくと非常に便利な本です。「サバイバル」パッケージの基盤を理解する上で非常に明確で役立ちます。