باستخدام سيم () مع لمر()

Question

لقد قمت بتشغيل اثنين من الانحدارات اللوجستية متعددة المستويات باستخدام نفس المتنبئين ، ولكن على ردين مختلفين:

fruitMLM <- lmer(InsuffFruit ~ Income + HDI + Income:HDI + (1 + Income | Country),family=binomial(link="logit"))  
fuelMLM <- lmer(Pollution ~ Income + HDI + Income:HDI + (1 + Income | Country),family=binomial(link="logit"))

Income منفصلة مع القيم c(-2,-1,0,1,2), HDI مستمر بين 0 و 1, Country هو قاطع ، والردود على حد سواء 1/0.

لرسم العصابات الثقة تشغيل محاكاة باستخدام سيم () وظيفة من حزمة أرم:

sim(fruitMLM,100)  
sim(fuelMLM,100)

أول واحد يحسب غرامة.الثانية ترجع الخطأ التالي:

Error in mvnorm(n.sims, bhat[j,], V.beta) :  
  'Sigma' is not positive definite

أنا فعلا أفعل هذا مع 8 ردود مختلفة.ستة منهم عملت بشكل جيد واثنين منهم عاد هذا الخطأ.

لا أحد يعرف كيفية تصحيح هذا?

Answer 1

يمكنك معرفة ما إذا كانت رسالة الخطأ صحيحة, أي, تحقق "باليد" ما إذا كانت مصفوفة التغاير من نموذج الفشل هو في الواقع إيجابية شبه محددة?لست متأكدا ، لكنني أتخيل أن هناك بعض الفحوصات التي يمكنك إجراؤها-ربما يوجد شيء هنا: http://en.wikipedia.org/wiki/Positive-definite_matrix