Rの時変共変量を使用したワイブル生存モデル

私はワイブルアプローチを使用して生存モデルを実行しようとしていますが、しわがあるのは、時変共変量があるためです。私はRでサバイバルパッケージを使用しています。私の呼び出しは次のとおりです。

output <- survreg(Surv(start, stop, fail) ~ gdppc + [...] + cluster(name), data = mydata, dist="weibull")

次のエラーが発生します。

Error in survreg(Surv(start, stop, fail) ~ gdppclag + : 
  Invalid survival type

coxph手順が正常に動作しますが、私はワイブルを使用します。

私の最初の質問は次のとおりです。ワイブルアプローチは、時変共変量を説明できますか？私はいくつかのテキストを見回しましたが、Cox PHアプローチは時変共変量に拡張できることがわかりました。ワイブルアプローチがそれを実行できるかどうかはあまり明確ではありません。

次に、実際にワイブルが機能する場合、それを処理できるRのパッケージは何ですか？

flexurvパッケージには、これを行うことができます：https://cran.r-project.org/web/packages/flexsurv/index.html

のflexsurvreg代わりに関数を呼び出すだけですsurvreg。

あなたSurvivalは次のようにパッケージでそれを行うことができます：

Srv <- Surv(start, stop, fail, type="interval" )

そして、あなたはSrvあなたのモデルで次のように使うことができます：

output <- survreg(Surv(start, stop, fail) ~ gdppc + [...] + cluster(name), data = mydata, dist="weibull")

これが役立つことを願っています:)

これはaftreg、ehaパッケージのコマンドでも実行できます。パッケージのドキュメントで説明されているように、survivalパッケージのコマンドとflexsurvパッケージのsurvregコマンドで通常使用されるのと同じパラメータ設定を取得するには、オプションをに設定する必要があります。したがって、コードの適応は次のようになります。flexsurvregparam"LifeExp"

output <- aftreg(Surv(start, stop, fail) ~ gdppc + [...] + cluster(name), data = mydata, dist="weibull", param="lifeExp")

このeha::aftregコマンドの利点の1つは、stargazerと互換性があるため、その出力をLaTeX形式に簡単にエクスポートできることです。