r - predict lme4 for type "response" (binary choices) -- Confidence intervals -

मैंने lme4 (प्रकार = "प्रतिक्रिया") के साथ एक मिश्रित प्रभाव लॉगिस्टग प्रतिगमन चलाया। अब मैंने भविष्यवाणी की सुविधा का इस्तेमाल किया और आत्मविश्वास अंतराल निर्धारित करना चाहता था।

मुझे भविष्यवाणियों के लिए यह कोड मिला (), यह काम करता है, लेकिन सीआई बाइनरी रिपॉन्स के लिए उपयुक्त नहीं हैं (मेरी भविष्यवाणियां 0 और 1 के बीच हैं, और अचानक सीआई 3 और 3 के बीच होती हैं। किसी को भी पता है कि यह कैसे समायोजित करें?

  लाइब्रेरी (lme4) लाइब्रेरी (ggplot2) # प्लॉटिंग एफएम 1 & lt; - ग्लमर (पसंद ~ रेटिंग + इंडिया + रेटिंग * इंडी + (1 | आईडी), डेटा = z , पारिवारिक = "द्विपद") newdat & lt; -data.frame (इंडी = फैक्टर (सी (1,1,1,1,1,1,2,2,2,2,2,2,2)), रेटिंग = फैक्टर ( 1: 6), आईडी = फैक्टर (प्रतिनिधि (सी (1:30), प्रत्येक = 12)), पसंद = 0) नईडेट $ समस्या & lt; -predict (fm1, newdata = newdat, re.form = null, type = " प्रतिक्रिया ") मिमी & lt; - model.matrix (terms (fm1), newdat) newdat $ विकल्प & lt; - भविष्यवाणी (fm1, newdat) ## या newdat $ विकल्प & lt; - mm% *% fixef (fm1) pvar1 & lt; - diag (मिमी% *% tcrossprod (vcov (fm1), मिमी)) tvar1 & lt; - pvar1 + VarCorr (एफएम 1) $ आईडी [1] ## अनुकूलित किया जाना चाहिए tvar1 & lt; - newdat & lt; - data.frame (newdat, plo = Newdat $ choice-2 * sqrt (pvar1), फाई = नईडेट $ पसंद + 2 * sqrt (pvar1), tlo = newdat $ पसंद- 2 * sqrt (tvar1), थी = newdat $ पसंद + 2 * sqrt (Tvar1)) #plot विश्वास g0 & lt; - ggplot (newdat, aes (x = रेटिंग, y = विकल्प, रंग = इंडी)) + geom_point () g0 + geom_errorbar (एईएस (ymin = plo, ymax = phi)) + Opts (शीर्षक = "फिक्स्ड-असर अनिश्चितता के आधार पर केवल सीआई") # प्लॉट भविष्यवाणी g0 + geom_errorbar (एईएस (ymin = tlo, ymax = thi)) + opts (शीर्षक = "एफई अनिश्चितता + आरई विचरण पर आधारित सीआई") < / Code> 

बहुत बहुत धन्यवाद!

आपको उलटा- लिंक फ़ंक्शन, plogis () इस मामले में:

  newdat & lt; - रूपांतरण (newdat, plo = plogis (plo), phi = plogis (phi), tlo = Plogis (tlo), thi = plogis (thi))  

अधिक सामान्यतः यदि आपने एक मॉडल gm1 को लगाया है तो व्युत्क्रम-लिंक फ़ंक्शन में < कोड> जीएम 1 resp $ family $ linkinv (यद्यपि इस तरह के किसी वस्तु के आंतरिक के साथ घूमना, भविष्य के संस्करणों के साथ संगत रहने की गारंटी नहीं है)।