"স্টাটা" বা "আর" -তে রিগ্রেশন বিচ্ছিন্নতার নকশায় গ্রাফগুলি

Question

লি এবং লেমিউক্স (পৃষ্ঠা 31, 2009) গবেষককে রিগ্রেশন বিচ্ছিন্নতা নকশা বিশ্লেষণ (আরডিডি) করার সময় গ্রাফগুলি উপস্থাপনের পরামর্শ দেন। তারা নিম্নলিখিত পদ্ধতিটি পরামর্শ দেয়:

"... কিছু ব্যান্ডউইথ জন্য , এবং বিন কিছু সংখ্যার জন্য কে 0 এবং বাম এবং কর্তক মান ঠিক আছে, যথাক্রমে, যা এই ধারণার বিন (গঠন করা হয় , ], জন্য + + , যেখানে "$h$$K_0$$K_1$$b_k$$b_{k+1}$$k = 1, . . . ,K = K_0$$K_1$$b_k = c−(K_0−k+1) \cdot h.$

c=cutoff point or threshold value of assignment variable
h=bandwidth or window width.

... তারপরে গড় ফলাফলগুলি কাটফয়েন্ট পয়েন্টের ঠিক বাম এবং ডানদিকে তুলনা করুন ... "

.. সব ক্ষেত্রে, আমরা কাটোফ পয়েন্টের প্রতিটি দিকে পৃথকভাবে অনুমান করা কোয়ার্টিক রিগ্রেশন মডেল থেকে values ​​tted মানগুলিও দেখাই ... (একই কাগজের পৃষ্ঠা 34)

আমার প্রশ্নটি হ'ল আমরা কীভাবে ধারালো আরডিডির জন্য অ্যাসাইনমেন্ট ভেরিয়েবলের সাথে ফলাফলের ভেরিয়েবলের (আত্মবিশ্বাসের ব্যবধানের সাথে) ফলাফলের গ্রাফ প্লট করার জন্য Stataবা সেই প্রক্রিয়াটিকে প্রোগ্রাম করব? এর Rএকটি নমুনা উদাহরণ এখানে এবং এখানেStata উল্লেখ করা হয়েছে (rd_obs দিয়ে rd প্রতিস্থাপন) এবং একটি নমুনা উদাহরণে Rহয় এখানে । যাইহোক, আমি মনে করি এই দুটিই পদক্ষেপ 1 কার্যকর করেনি implement দ্রষ্টব্য, উভয় প্লটের মধ্যে লাগানো লাইনগুলির সাথে কাঁচা ডেটা রয়েছে।

আত্মবিশ্বাসের পরিবর্তনশীল ছাড়াই নমুনা গ্রাফ [লি এবং লেমিউক্স, ২০০৯] আগাম আপনাকে ধন্যবাদ।

Comments

আপনার পতাকার প্রতিক্রিয়া হিসাবে, আপনার প্রশ্নকে পুনরজ্জীবিত করার একটি ভাল উপায় হ'ল এটি সম্পাদনা করা এবং একটি অনুগ্রহ অফার করা: এটি আপনার প্রশ্নকে ঘায়েল করবে এবং এতে আরও বেশি লোককে আগ্রহী করবে। আপনি যদি মনে করেন এই প্রশ্নটি স্ট্যাক ওভারফ্লোতে আরও ভাল পরিবেশিত হতে পারে তবে আমাদের জানান এবং আমরা আপনার জন্য এটি স্থানান্তর করতে পারি।

— chl

আমি এটি স্ট্যাক ওভারফ্লোতে স্থানান্তরিত হতে চাই।

— মেট্রিকস

দুর্ভাগ্যক্রমে, এই প্রশ্নটি স্ট্যাক ওভারফ্লোতে স্থানান্তরিত হওয়ার পক্ষে খুব পুরানো। আমি বিশ্বাস করি এটি ক্রস ভ্যালিডেটের সাথে সম্পর্কিত তবে আপনি যদি স্ট্যাক ওভারফ্লো সম্পর্কে জিজ্ঞাসা করতে চান ( প্রোগ্রামিং দিকটির উপর জোর দেওয়া এবং একটি নূন্যতম পুনরুত্পাদনযোগ্য উদাহরণ প্রদান করতে ) তবে আমাকে জানান এবং আমি এটি এখানে বন্ধ করব।

— chl

আপনি সেমিগ্রাম ব্যবহার করা উচিত । এটি আপনার প্রয়োজনীয় সমস্ত কিছু করে।

— ইয়ান সং

Answer 1

$K_i$

use votex // the election-spending data that comes with rd

tw 
(scatter lne d, mcolor(gs10) msize(tiny)) 
(lpolyci lne d if d<0, bw(0.05) deg(2) n(100) fcolor(none)) 
(lpolyci lne d if d>=0, bw(0.05) deg(2) n(100) fcolor(none)), xline(0)  legend(off)

twoway$x$$s$$se$$ul$$ll$

lpoly lne d if d<0, bw(0.05) deg(2) n(100) gen(x0 s0) ci se(se0)
lpoly lne d if d>=0, bw(0.05) deg(2) n(100) gen(x1 s1) ci se(se1)

/* Get the 95% CIs */
forvalues v=0/1 {
    gen ul`v' = s`v' + 1.95*se`v' 
    gen ll`v' = s`v' - 1.95*se`v' 
};

tw 
(line ul0 ll0 s0 x0, lcolor(blue blue blue) lpattern(dash dash solid)) 
(line ul1 ll1 s1 x1, lcolor(red red red) lpattern(dash dash solid)), legend(off)  

আপনি দেখতে পাচ্ছেন যে প্রথম চক্রান্তের রেখাগুলি দ্বিতীয়টির মতোই।

Comments

@ ডিমিট্রি: সমাধানের জন্য +1 করুন। যাইহোক, আমি কাঁচা মানগুলি দেখানোর জন্য ছড়িয়ে পড়া প্লটের পরিবর্তে প্রতিটি বিনের (দয়া করে উপরের স্টাটা উদাহরণটি চালনা করুন) এর গড় মূল্য পেতে চাই। সিআই দুর্দান্ত।

— মেট্রিকস

আপনার অর্থ কী তা আমি নিশ্চিত নই am আমি কোডে যুক্ত করে দেখিয়েছি যে কীভাবে আপনি হাতের মুঠোয় প্রতিটি বিনের স্মুথযুক্ত অর্থ পাবেন। যদি আপনি যা খুঁজছেন তা যদি না হয় তবে দয়া করে আপনার মনে কী রয়েছে তা আরও বিশদে ব্যাখ্যা করুন। আমি যতদূর বলতে পারি, এই গ্রাফগুলি সাধারণত কাঁচা ডেটা এবং স্মুথড মানে দেখায়।

— দিমিত্রি ভি। মাস্টারভ

লি এবং লেমিউক্সের উদ্ধৃতি (পৃষ্ঠা 31, ২০০৯): "ডেটা গ্রাফিক করার একটি স্ট্যান্ডার্ড পদ্ধতি হ'ল অ্যাসাইনমেন্ট ভেরিয়েবলকে (ডি এখানে) একটি বিংশে বিভক্ত করা, কাট অফের প্রতিটি পাশের দুটি পৃথক বিন রয়েছে কিনা তা নিশ্চিত করে পয়েন্ট (চিকিত্সা এবং চিকিত্সাবিহীন পর্যবেক্ষণগুলি একই বিনাতে একত্রে মিশ্রণ এড়ানোর জন্য) তারপরে, ফলাফলের পরিবর্তনশীলের গড় মান প্রতিটি বিনের জন্য গণনা করা যায় এবং বিনের মাঝের পয়েন্টগুলির সাথে আঁকানো যায় "। সুতরাং, যদি 50 টি বিন থাকে, তবে আমাদের বাম এবং ডানদিকে কেবল 25 ডেটা পয়েন্ট থাকবে এবং সমস্ত কাঁচা তথ্য নেই (যেমন, রেফারেন্সের গ্রাফ 6 (খ): প্রশ্নে আপডেট হয়েছে)

— মেট্রিক্স

এখন এটা পরিষ্কার! আমি কার্নেলের সাথে একমত তবে আপনি কি নিশ্চিত যে এটি এখন 0 ডিগ্রি নয়? এটি সমান ওজনযুক্ত মানে মসৃণতার সাথে মিল রাখে।

— দিমিত্রি ভি। মাস্টারভ

আমি বিশ্বাস করি যে নিয়মিত কার্নেল এবং ডিগ্রি 0 বহুবর্ষের সাথে এলপোলির সাথে মিলে যায়

— দিমিত্রি ভি। মাস্টারভ

Answer 2

এখানে একটি ক্যানড অ্যালগরিদম। ক্যালোনিকো, কাতানোও এবং তিতুনিক সম্প্রতি শক্তিশালী ব্যান্ডউইথ নির্বাচনের জন্য একটি পদ্ধতি প্রস্তাব করেছিলেন। তারা তাদের তাত্ত্বিক কাজ স্টাটা এবং আর উভয়ের জন্যই বাস্তবায়িত করেছিল এবং এটি প্লট কমান্ডের সাথে আসে। এখানে আর এর একটি উদাহরণ রয়েছে:

# install.packages("rdrobust")
library(rdrobust)
set.seed(26950) # from random.org
x<-runif(1000,-1,1)
y<-5+3*x+2*(x>=0)+rnorm(1000)
rdplot(y,x)

এটি আপনাকে এই গ্রাফটি দেবে:

Comments

হাই, সিআই কিভাবে যুক্ত করবেন?

— ক্র্যান্টজ