मैं एक सहयोगी के काम को दोहराने की कोशिश कर रहा हूं और स्टैटा से आर तक विश्लेषण को आगे बढ़ा रहा हूं। वह मॉडल जो मानक त्रुटियों को क्लस्टर करने के लिए 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"))
यह स्टैटा में विश्लेषण से परिणामों की नकल नहीं करता है, संभवतः क्योंकि यह ओएलएस पर काम करने के लिए डिज़ाइन किया गया है न कि नकारात्मक द्विपद। तो खोज चलती है। मैं गलत हो रहा हूँ, जहां पर किसी भी संकेत बहुत सराहना की जाएगी