3

survival次のようにRパッケージを使用して、状態1から状態2への遷移のCox比例ハザードモデルを推定しようとしています:

Altman <- coxph(Surv(Tstart, Tstop, to == 2) ~ wWCTA + wRETA + wEBITTA + wMETL + 
            wSTA, data=Multistate, subset = from == 1, 
          frailty(x=id, distribution= "gamma"))

frailty引数なしでモデルを推定している場合、問題なく動作します。しかし、それを含めると、次のエラーが発生します。

Error in if (any(ord > 1)) stop("Penalty terms cannot be in an interaction") : 
missing value where TRUE/FALSE needed

誰が私が間違っているのか説明してもらえますか?

よろしく、

J

4

1 に答える 1

4

私は同じ問題を抱えていましたが、最終的に脆弱性という用語は相加的であることに気付きました! つまり、コンマの後ではなく、プラスの後に来る必要があります。

したがって、コードを次のように変更すると、問題が解決するはずです。

Altman <- coxph(Surv(Tstart, Tstop, to == 2) ~ wWCTA + wRETA + wEBITTA + wMETL + wSTA + frailty(x=id, distribution= "gamma"), data=Multistate, subset = from == 1)
于 2015-04-15T15:40:35.707 に答える