नकारात्मक द्विपद प्रतिगमन का उपयोग करते समय क्लस्टर विकल्प के बराबर


10

मैं एक सहयोगी के काम को दोहराने की कोशिश कर रहा हूं और स्टैटा से आर तक विश्लेषण को आगे बढ़ा रहा हूं। वह मॉडल जो मानक त्रुटियों को क्लस्टर करने के लिए nbreg फ़ंक्शन के भीतर "क्लस्टर" विकल्प को लागू करता है।

Http://repec.org/usug2007/crse.pdf इस विकल्प के क्या और क्यों के एक पूर्ण विवरण के लिए देखें

मेरा सवाल यह है कि आर के भीतर नकारात्मक द्विपद प्रतिगमन के लिए इसी विकल्प को कैसे लागू किया जाए?

हमारे पेपर में प्राथमिक मॉडल निम्नानुसार निर्दिष्ट है

 xi: nbreg cntpd09 logpop08 pcbnkthft07 pccrunion07 urbanpop pov00 pov002 edu4yr ///
 black04 hispanic04 respop i.pdpolicy i.maxloan rollover i.region if isser4 != 1,   
 cluster(state)

और मैंने इसे बदल दिया है

pday<-glm.nb(cntpd09~logpop08+pcbnkthft07+pccrunion07+urbanpop+pov00+pov002+edu4yr+
black04+hispanic04+respop+as.factor(pdpolicy)+as.factor(maxloan)+rollover+
as.factor(region),data=data[which(data$isser4 != 1),])

जो स्पष्ट रूप से संकुल त्रुटियों का अभाव है।

क्या एक सटीक प्रतिकृति करना संभव है? यदि हां, तो कैसे? यदि नहीं, तो कुछ उचित विकल्प क्या हैं?

धन्यवाद

[संपादित करें] जैसा कि टिप्पणियों में कहा गया है, मैं एक समाधान की उम्मीद कर रहा था जो मुझे बहुस्तरीय मॉडल के दायरे में नहीं ले गया था। जबकि मेरा प्रशिक्षण मुझे यह देखने की अनुमति देता है कि इन चीजों को संबंधित होना चाहिए, यह एक छलांग से अधिक है जितना मैं अपने आप को लेने में सहज हूं। जैसे मैंने खुदाई की और यह लिंक पाया: http://landroni.wordpress.com/2012/06/02/fama-macbeth-and-cluster-robust-by-firm-and-time-standard-errors-in- आर /

जो मैं चाहता हूँ करने के लिए कुछ सीधे कोड की ओर इशारा करता है:

library(lmtest)
pday<-glm.nb(cntpd09~logpop08+pcbnkthft07+pccrunion07+urbanpop+pov00+pov002+edu4yr+
 black04+hispanic04+respop+as.factor(pdpolicy)+as.factor(maxloan)+rollover+
 as.factor(region),data=data[which(data$isser4 != 1),])
summary(pday)

coeftest(pday, vcov=function(x) vcovHC(x, cluster="state", type="HC1"))

यह स्टैटा में विश्लेषण से परिणामों की नकल नहीं करता है, संभवतः क्योंकि यह ओएलएस पर काम करने के लिए डिज़ाइन किया गया है न कि नकारात्मक द्विपद। तो खोज चलती है। मैं गलत हो रहा हूँ, जहां पर किसी भी संकेत बहुत सराहना की जाएगी


3
आपको बेन बोलकर के नोट यहां उपयोगी हो सकते हैं।
शाम

और यह भी देखने इस पहले प्रश्न
fmark

FYI यहां स्टाटा की मजबूत क्लस्टरिंग मानक त्रुटियों की एक परिभाषा है। उन्हें लगता है कि लागू करने के लिए कठिन नहीं है। IMO आप बूटस्ट्रैप्ड या कटे-फटे मानक त्रुटियों के साथ वैसे भी बेहतर हो सकते हैं ( vce पर सहायता देखें )। मैं किसी भी आर संकुल का सुझाव नहीं दे सकता। एक प्रतिस्थापन खोजने पर शुभकामनाएँ!
एंडी डब्ल्यू

धन्यवाद @fmark - बहुत उपयोगी टिप्पणियाँ, मेरे "उत्तर" से बहुत बेहतर है और मैंने इसे तदनुसार अपडेट किया है।
पीटर एलिस

सभी को धन्यवाद। मुझे लगता है कि मेरे प्रश्न का संक्षिप्त उत्तर यह है कि कोई सीधा प्रतिस्थापन नहीं है (उदाहरण के लिए, एक पूर्व-निर्मित फ़ंक्शन जो बिल्कुल क्लस्टर को बदलता है)। स्पष्ट रूप से अधिक अनुभव वाला कोई व्यक्ति बेन बोल्कर के नोट्स के माध्यम से रास्ता देख सकता है, लेकिन यह मुझे नए क्षेत्र में ले जाता है जहां मुझे यकीन नहीं हो रहा था कि मुझे सूत्र कथन सही मिल रहे हैं। मुझे यकीन नहीं है कि एक उत्तर को स्वीकार किए बिना "धन्यवाद" कहने का उपयुक्त तरीका क्या है, लेकिन आपके पास मेरा धन्यवाद है, और कमियां मेरी अपनी हैं।
csfowler

जवाबों:


4

यह दस्तावेज़ दिखाता है कि कैसे एक glm प्रतिगमन के लिए क्लस्टर SEs प्राप्त करें:

http://dynaman.net/R/clrob.pdf


मुझे इसे स्टैटा के परिणामों के साथ एक तुलनात्मक परीक्षा देनी होगी, लेकिन यह बिल्कुल वैसा ही है जैसा मैं उम्मीद कर रहा था।
csfowler

1

यह पूरी तरह से संतोषजनक जवाब नहीं है ...

मैंने इसे स्वयं आज़माया नहीं है, लेकिन ऐसा लगता है कि glmmADMB पैकेज ऐसा कर सकता है जो आप चाहते हैं।

मैं इस सवाल पर @ fmark की टिप्पणी से बेशर्मी से चुटकी लूंगा और उससे सहमत हूं कि बेन बोल्कर के नोट्स उपयोगी हैं, जैसा कि यह पहले वाला सवाल है , जो बिल्कुल सटीक डुप्लिकेट नहीं है लेकिन बहुत समान मुद्दों को शामिल करता है।

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