আমি একজন সহকর্মীর কাজের প্রতিলিপি তৈরির চেষ্টা করছি এবং স্টাটা থেকে আর-তে বিশ্লেষণটি সরিয়ে নিয়ে যাচ্ছি 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"))
এটি স্টাটায় বিশ্লেষণের ফলাফলগুলির প্রতিলিপি দেয় না যদিও সম্ভবত এটি ওএলএসে নেতিবাচক দ্বিপদী নয় কাজ করার জন্য ডিজাইন করা হয়েছে। সুতরাং অনুসন্ধান চলছে। আমি যেখানে ভুল করছি সে সম্পর্কে যে কোনও পয়েন্টার প্রশংসা করবে