0

clogit完全なモデル(パッケージから)でドレッジを実行しようとしていますが、survival試行するたびに R がクラッシュし続けます。「R で致命的なエラーが発生しました。セッションが終了しました」というエラー メッセージが表示されます。

library(survival)
FullModel <- clogit(OBSERVED ~ Canopy + distgr_sca + Near_Hwy3 + strata(PID1), Compiled)
library(MuMIn)
dredge(FullModel)

これは私のコードです。これがクラッシュする理由、またはモデルで完全なドレッジを実行できる他の方法についてのアイデアはありますか?

ありがとうございました。

4

1 に答える 1

1

clogitモデル式で階層strata()が指定されていない場合、またはが唯一のモデル項である場合、R がクラッシュするようです。dredgeすべてのモデルを保持strata(PID1)し、サブセットを制限して最小 3 つの変数を持つように指示する必要があります。

dredge(fm0, eval = T, fixed = "strata(id)", m.min = 3)

再現可能な例:

library(survival)
library(MuMIn)

## from example(clogit)
resp <- levels(logan$occupation)
n <- nrow(logan)
indx <- rep(1:n, length(resp))
logan2 <- data.frame(logan[indx,], id = indx, tocc = factor(rep(resp, each=n)))
logan2$case <- (logan2$occupation == logan2$tocc)
fm <- clogit(case ~ tocc + tocc:education + strata(id), logan2, na.action = "na.fail")
##

dredge(fm, fixed = "strata(id)", m.min = 3)
于 2014-04-04T10:00:57.410 に答える