2016-07-02 25 views
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)) 

任意の考え?

+0

再現可能な例を追加することで、すべての作業があなたの質問を形にすることに感謝します。残念ながら、このエラーの解決方法はわかりませんが、問題の情報に基づいて再現できるようになりました。コミュニティの誰かが助けてくれることを願っています。 – josliber

+2

これはR-forgeで修正されました。バージョン1.15.8へのアップデート –

+0

@KamilBartoń、回答としてコメントを投稿できますか? (ちょうど古い 'lme4'の質問を整理する...) –

答えて

3

これはバージョン> = 1.15.8のfixed(すぐにCRAN、現在R-Forge上にある)のバグでした。

関連する問題