रैखिक मिश्रित प्रभाव मॉडल के लिए उत्तोलन और कुक की दूरी को कैसे निकालना / गणना करना है


10

क्या किसी को पता है कि (या पैकेज के merमाध्यम से प्राप्त lme4) एक वर्ग वस्तु के लिए उत्तोलन और कुक की दूरी की गणना कैसे करें (निकालें )? मैं एक अवशिष्ट विश्लेषण के लिए इनकी साजिश करना चाहता हूं।

जवाबों:


13

आपको आर पैकेज पर एक नजर डालनी चाहिए influence.ME। यह आपको मिश्रित प्रभाव मॉडल द्वारा उत्पन्न प्रभावशाली डेटा के उपायों की गणना करने की अनुमति देता है lme4

एक उदाहरण मॉडल:

library(lme4)
model <- lmer(mpg ~ disp + (1 | cyl), mtcars)

फ़ंक्शन influenceसभी और चरणों के लिए आधार है:

library(influence.ME)
infl <- influence(model, obs = TRUE)

कुक की दूरी की गणना करें:

cooks.distance(infl)

प्लॉट कुक की दूरी:

plot(infl, which = "cook")

यहां छवि विवरण दर्ज करें


धन्यवाद! यह निश्चित रूप से मदद करता है। कुक की दूरी बनाम उत्तोलन की साजिश के लिए उत्तोलन की गणना कैसे करें?
Roey Angel

@RoeyAngel मुझे लगता है कि यह influence.MEपैकेज के साथ संभव नहीं है । दुर्भाग्य से, मेरे पास इस कार्य का हल नहीं है।
स्वेन हेनस्टीन

क्या यह नहीं होना चाहिए infl <- influence(model, group = "cyl"), क्योंकि आपने यादृच्छिक प्रभाव निर्दिष्ट किया है (1|cyl)? मुझे नहीं पता, मैं यह बिल्कुल नहीं समझता, मैंने सिर्फ प्रभाव स्थापित किया है ... लेकिन मैं वास्तव में नहीं जानता कि कब उपयोग करना है obs = TRUEऔर कब उपयोग करना है group...
जिज्ञासु

मैं निम्नलिखित जोड़ना चाहूंगा: यदि आप कुक की डी दूरियां होने वाली पंक्ति संख्या प्राप्त करना चाहते हैं - प्लॉटिंग के बिना प्लॉट में होने वाली एक ही संख्या है, तो आप कुक के साथ कुक डी दूरी संख्याओं के बारे में निम्नलिखित r सूत्र का उपयोग कर सकते हैं ०.१ के मूल्य से दूरcooksD_data<-as.data.frame(cooks.distance(ft1)) cooksD_data_select<-cooksd[cooksD_data>0.1,drop=FALSE,] cooksD_oultiers<-as.numeric(rownames(cooksD_data_select))]
एलियास एस्टैटिस्टिक्स

क्या यह यहां सुझाए गए hatvalues()फ़ंक्शन से बेहतर है ?
जिज्ञासु
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.