3
"lme4" Rパッケージのglmer()を使ってGLMMを構築していますが、 "MuMIn"パッケージのr.squaredGLMM()を使ってR2値を見積もる際にエラーが発生することがあります。r.squaredGLMM()のエラー
私は合うようにしようとしていますモデルはこの1つのsimmilarです:次に
library(lme4)
lmA <- glmer(x~y+(1|w)+(1|w/k), data = data1, family = binomial(link="logit"))
、ESTIME R2に、私が使用:
library(MuMIn)
r.squaredGLMM(lmA)
そして、私はこの取得:
The result is correct only if all data used by the model has not changed since model was fitted. Error in .rsqGLMM(fam = family(x),
varFx = var(fxpred), varRe = varRe, : 'names' attribute [2] must be the same length as the vector [0]
を
このエラーが表示される理由をご存知ですか?たとえば、単一のランダム要素(この場合は(1|w)
)のみを使用すると、このエラーは表示されません。ここで
は私のデータセットです。
data1 <-
structure(list(w = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L,
1L, 2L, 1L), .Label = c("CA", "CB"), class = "factor"), k = structure(c(4L,
4L, 3L, 3L, 3L, 4L, 1L, 3L, 2L, 3L, 2L), .Label = c("CAF01-CAM01",
"CAM01", "CBF01-CBM01", "CBM01"), class = "factor"), x = c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L), y = c(-0.034973549,
0.671720643, 4.557044729, 5.347170897, 2.634240583, -0.555740207,
4.118277809, 2.599825716, 0.95853864, 4.327804344, 0.057331718
)), .Names = c("w", "k", "x", "y"), class = "data.frame", row.names = c(NA,
-11L))
任意の考え?
再現可能な例を追加することで、すべての作業があなたの質問を形にすることに感謝します。残念ながら、このエラーの解決方法はわかりませんが、問題の情報に基づいて再現できるようになりました。コミュニティの誰かが助けてくれることを願っています。 – josliber
これはR-forgeで修正されました。バージョン1.15.8へのアップデート –
@KamilBartoń、回答としてコメントを投稿できますか? (ちょうど古い 'lme4'の質問を整理する...) –