2016-07-29 32 views
0

でTRUEはエラーではありません私はglmerをモデルに合うようにしようとしていると私はエラーを取得しておいてください。た(p < - NcoI部位(X))==をNcoI(Y)はglmer

Error : (p <- ncol(X)) == ncol(Y) is not TRUE

にです8つの機能と数千点の大きなモデルです。このバグを調査すると、それはあまりにも多くのNAによって引き起こされることがわかった。私は自分の機能を見て、そのうちの1人だけがNAを持っていました。モデルからその機能を除外しましたが、まだエラーがあります。データが大きすぎて転記できません。

http://r.789695.n4.nabble.com/Error-from-lme4-quot-Error-p-lt-ncol-X-ncol-Y-is-not-TRUE-quot-td4712617.html

モデルは次のとおりです。

covariates=c("Sex", "PC1", "PC2", "PC3", "PC4", "label", "Alive") 

fix.eff=paste("outcome","~",xi) 

if (!is.null(covariates)) {for (covi in covariates) fix.eff=paste(fix.eff,"+",covi) } 

fix.eff=formula(paste(fix.eff, "+(1|",ind.family,")")) 

fit <- try(glmer(fix.eff,family=binomial(link='logit'), data=x)) 
+0

再現可能な例がなければトラブルシューティングが難しくなります。あなたが今のように式を構築するのではなく、単純に 'glmer'でモデルを"手書き "で書くと同じエラーが出ますか? – aosmith

答えて

1

問題は次の行にあった:それはされている必要があります

fix.eff=formula(paste(fix.eff, "+(1|",ind.family,")")) 

fix.eff=formula(paste(fix.eff, "+(1| ind.family)")) 

ind.familyはデータセットの列x.元の方法で記述されたもので、列の代わりに現在の値として読み込まれました。式には(1|ind.family)の代わりに(1|1)という語句が含まれていました。それは固定されているので、それは動作します!

他の誰かがこのエラーを受け取った場合は、あなたの数式とデータをよく見て、それがあなたの考えであることを確認することをお勧めします。

関連する問題