coxme比例ハザード仮定

Rのcoxme関数{coxme}を使用して混合効果Coxモデルを実行しています。比例ハザードの仮定を確認したいと思います。coxme比例ハザード仮定

PHの仮定は、cox.phphモデルのcox.zph関数{survival}で確認できます。

しかし、私はcoxmeモデルの同等物を見つけることができません。

2015年にも同様の質問が掲載されました。hereですが、答えはありませんでした。

私の質問は： 1）混合効果coxモデルcoxmeでPHの仮定をテストする方法1）？ 2）coxmeモデルのcox.zphに相当するものがない場合は、科学論文の出版物が混合効果coxmeモデルを実行するが、coxmeモデルと同一のcox.phモデルでPHの仮定をテストすることは有効ですか？効果？

ご回答いただきありがとうございます。よろしくお願いします。

出典

2016-12-30 user3637399

1）私はこの質問がCross Validatedに適していると思います。 – adibender

2）最も簡単な方法は、データをカウントプロセスフォーマットに変換することです（例えば 'survival :: survSplit'を参照し、PHの仮定をテストしたい変数と時間変数自身の変換をあなたは 'coxph'呼び出しに虚偽の言葉を含めることができます。例えば' x * log（t） '）を使用してモデルに入力してください。https://www.jstatsoft.org/article/view/v061c01/v61c01.pdf – adibender

3）虚弱用語。 '？survival :: frailty'の例を見てください。私はあなたが 'coxph'関数から返されたオブジェクトに対してcox.zph関数を弱い言葉で使うことができるはずだと思います。 – adibender

の機能でfrailtyオプションを使用できます。たとえば、ランダムエフェクト変数がB、固定効果変数がAであるとします。そして、あなたが今

myfit <- coxph(Surv(Time, Censor) ~ A + frailty(B) , data = mydata)

以下のようにあなたのモデルを当てはめるには、比例ハザード仮定をテストするためにcox.zph(myfit)を使用することができます。

出典

2017-07-20 20:50:05 user3373954

答えて

関連する問題