নেতিবাচক দ্বিপদী রিগ্রেশন ব্যবহার করার সময় ক্লাস্টার বিকল্পের সমতুল্য


10

আমি একজন সহকর্মীর কাজের প্রতিলিপি তৈরির চেষ্টা করছি এবং স্টাটা থেকে আর-তে বিশ্লেষণটি সরিয়ে নিয়ে যাচ্ছি she

এই বিকল্পটি কী এবং কেন এর পুরোপুরি বিবরণের জন্য 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-erferences-in- R /

যা আমি চাই তা করার জন্য কিছুটা সহজ সরল কোডকে নির্দেশ করে:

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


এফওয়াইআই এখানে স্টাতার শক্তিশালী ক্লাস্টার স্ট্যান্ডার্ড ত্রুটির সংজ্ঞা রয়েছে। তারা বাস্তবায়ন করতে যে কঠিন মনে হয় না। আইএমও আপনি যেভাবেই হোক বুটস্ট্র্যাপড বা জ্যাকনিফাইড স্ট্যান্ডার্ড ত্রুটিগুলি ( ভেসে সহায়তা দেখুন ) দিয়ে ভাল হতে পারেন । যদিও আমি কোনও আর প্যাকেজ প্রস্তাব করতে পারি না। প্রতিস্থাপনের জন্য সৌভাগ্য!
অ্যান্ডি ডাব্লু

ধন্যবাদ @ ফারমার্ক - খুব দরকারী মন্তব্য, আমার "উত্তর" এর চেয়ে অনেক ভাল এবং আমি সে অনুযায়ী আপডেট করেছি।
পিটার এলিস

সবাইকে ধন্যবাদ. আমি মনে করি আমার প্রশ্নের সংক্ষিপ্ত উত্তরটি হ'ল কোনও সরল প্রতিস্থাপন নেই (উদাহরণস্বরূপ, একটি প্রাক-তৈরি ফাংশন যা ক্লাস্টারের বিকল্পটিকে হুবহু প্রতিস্থাপন করে)। স্পষ্টতই আরও অভিজ্ঞতার সাথে বেন বলকারের নোটগুলির পথ দেখতে পাবে তবে এটি আমাকে নতুন অঞ্চলে নিয়ে যায় যেখানে আমি নিশ্চিত হতে পারি না যে আমি সূত্রের বক্তব্যগুলি সঠিকভাবে পেয়েছি। কোনও উত্তর গ্রহণ না করেই "থ্যাঙ্কস" বলার উপযুক্ত উপায়টি কী তা আমি নিশ্চিত নই, তবে আপনার আমার ধন্যবাদ আছে এবং ত্রুটিগুলি আমার নিজস্ব।
csfowler

উত্তর:


4

এই দস্তাবেজটি দেখায় যে কীভাবে কোনও গ্ল্যাম রিগ্রেশন জন্য ক্লাস্টার এসই পাবেন:

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


স্ট্যাটার ফলাফলের সাথে আমাকে এটির পরীক্ষার তুলনা দিতে হবে, তবে এটি ঠিক যা প্রত্যাশা করেছিলাম তার মতোই।
csfowler

1

এটি সম্পূর্ণ সন্তোষজনক উত্তর নয় ...

আমি নিজে চেষ্টা করে দেখিনি তবে মনে হচ্ছে glmmADMB প্যাকেজটি আপনি যা চান তা করতে পারে।

আমি প্রশ্নটিতে @ ফারমার মন্তব্য থেকে নির্লজ্জভাবে চিমটি টানব এবং তার সাথে একমত হব যে বেন বলকারের নোটগুলি দরকারী, এটি পূর্ববর্তী প্রশ্নটি , যা একেবারে যথাযথ সদৃশ নয় তবে এটি খুব অনুরূপ ইস্যুগুলিকে অন্তর্ভুক্ত করে।

আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.