باستخدام سيم () مع لمر()
سؤال
لقد قمت بتشغيل اثنين من الانحدارات اللوجستية متعددة المستويات باستخدام نفس المتنبئين ، ولكن على ردين مختلفين:
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 ردود مختلفة.ستة منهم عملت بشكل جيد واثنين منهم عاد هذا الخطأ.
لا أحد يعرف كيفية تصحيح هذا?
المحلول
يمكنك معرفة ما إذا كانت رسالة الخطأ صحيحة, أي, تحقق "باليد" ما إذا كانت مصفوفة التغاير من نموذج الفشل هو في الواقع إيجابية شبه محددة?لست متأكدا ، لكنني أتخيل أن هناك بعض الفحوصات التي يمكنك إجراؤها-ربما يوجد شيء هنا: http://en.wikipedia.org/wiki/Positive-definite_matrix
لا تنتمي إلى StackOverflow