में मजबूत क्लस्टर्ड मानक त्रुटियों के साथ लॉजिस्टिक रिग्रेशन एक नौसिखिया प्रश्न: क्या कोई जानता है कि आर में क्लस्टर मानक त्रुटियों के साथ एक लॉजिस्टिक रिग्रेशन कैसे चलाया जाए? स्टाटा में यह सिर्फ logit Y X1 X2 X3, vce(cluster Z)
है, लेकिन दुर्भाग्य से मुझे पता नहीं चला है कि आर में एक ही विश्लेषण कैसे करें। अग्रिम धन्यवाद!आर
आर
उत्तर
आप rms
(प्रतिगमन मॉडलिंग रणनीतियों) पैकेज को देखना चाहते हैं। तो, lrm
रसद प्रतिगमन मॉडल है, और अगर fit
अपने उत्पादन का नाम है, तो आप कुछ इस तरह होगा:
fit=lrm(disease ~ age + study + rcs(bmi,3), x=T, y=T, data=dataf)
fit
robcov(fit, cluster=dataf$id)
bootcov(fit,cluster=dataf$id)
आपके पास मॉडल बयान में x=T
, y=T
निर्दिष्ट करने के लिए। rcs
3 समुद्री मील के साथ प्रतिबंधित क्यूबिक splines इंगित करता है।
बहुत बहुत धन्यवाद! यह चमत्कार काम किया है! मैं आरएमएस के मैनुअल को और अधिक बारीकी से पढ़ूंगा और देख सकता हूं कि देश और साल-दर-साल गुणांक को क्लस्टर करने का कोई तरीका है या नहीं। एक बार फिर आपका धन्यवाद! – danilofreire
यह उत्तर पहले से ही बहुत अच्छा है लेकिन अगर यह पूरी तरह से प्रतिकृति योग्य था तो इसे बेहतर किया जा सकता है। मुझे नहीं पता कि चर कहां से आते हैं, आउटपुट क्या है, और क्यों 'आरसीएस (बीएमआई, 3)' आवश्यक है। – MERose
मैं पिछले दो दिनों से इस समस्या के खिलाफ अपने सिर पर टक्कर लगी हूं; मैं जादुई रूप से पाया कि एक नया पैकेज प्रतीत होता है जो महान चीजों के लिए नियत लगता है - उदाहरण के लिए, मैं अपने विश्लेषण में कुछ क्लस्टर-मजबूत टोबीट मॉडल भी चला रहा हूं, और इस पैकेज में कार्यक्षमता भी बनाई गई है। सिंटैक्स का उल्लेख नहीं है कि मैंने देखा है कि अन्य सभी समाधानों की तुलना में अधिक क्लीनर है (हम स्वच्छ-स्टाटा के स्तर के बारे में बात कर रहे हैं)।
तो अपने खिलौना उदाहरण के लिए, मैं चलाने चाहते हैं:
library(Zelig)
logit<-zelig(Y~X1+X2+X3,data=data,model="logit",robust=T,cluster="Z")
एट देखा!
वाह, ऐसा लगता है कि मेरे आर कोड कभी ऐसा नहीं लगता है कि "बस काम करें" दिखाई देता है। क्या यह नई कार्यक्षमता है? यदि नहीं, तो ज़ेलिग आर में इसे हल करने का वैचारिक तरीका क्यों नहीं रहा है? – Philip
पता नहीं, लेकिन मुझे आशा है कि यह ऐसा हो जाएगा। [परियोजना] (http://zeligproject.org) निश्चित रूप से महत्वाकांक्षी लगता है! [Google समूह] (https://groups.google.com/forum/m/#!forum/zelig-statistic-software) हालांकि इतना सक्रिय प्रतीत नहीं होता है, इसलिए सुनिश्चित करें कि कितनी तेज़ी से प्रगति हुई है। – MichaelChirico
दुर्भाग्यवश, मुझे लगता है कि आदेश 'ज़ेलिग' (सीआरएएन पर) के नवीनतम संस्करण में काम नहीं करता है। मैंने बस 'क्लस्टर' तर्क के साथ और बिना कुछ मॉडल चलाए हैं और मानक त्रुटियां बिल्कुल वही हैं। मेरा मानना है कि यह संस्करण 4.0 के बाद से ऐसा ही रहा है, पिछली बार मैंने पैकेज का उपयोग किया था। – danilofreire
'सैंडविच' पैकेज में 'vcovHC()' फ़ंक्शन भी उपयोगी हो सकता है (सुनिश्चित नहीं है कि यह लॉजिस्टिक रिग्रेशन अनुमानों पर लागू होता है) –
यदि आप स्टाटा से माइग्रेट कर रहे हैं तो आपको '' plm'' नामक पैकेज मिल सकता है। उपयोगी। इसके अलावा, –
अनुमान के बाद भिन्नता मानक covariance मैट्रिक्स को जोड़कर पैनल को सही मानक त्रुटियों को कार्यान्वित करने के लिए '' pcse'' नामक पैकेज है, जो आपके उत्तरों, बेन और हबर्ट के लिए बहुत बहुत धन्यवाद। मैं आपके द्वारा सुझाए गए पैकेजों का भी परीक्षण करूंगा और देख सकता हूं कि वे लॉजिस्टिक अनुमानों के साथ काम करते हैं या नहीं। एक बार फिर धन्यवाद! – danilofreire