假設(shè)我的數(shù)據(jù)集是auto.dta
因變量為price
自變量為rep78
控制變量為headroom weight length
回歸模型為OLS神凑，即主回歸代碼為reg rep78 headroom weight length
下面是1000次placebo test的代碼
forvalue i=1/1000{sysuse auto,clear//調(diào)入數(shù)據(jù)*-思路：打亂rep78,即將rep78的全部取值拿出暫存紊浩，然后隨機(jī)賦給每一個樣本*-打亂rep78,即將rep78的全部取值拿出暫存? ? g obs_id=_n//初始樣本序號gen random_digit=runiform()//生成隨機(jī)數(shù)sort random_digit//按新生成的隨機(jī)數(shù)排序g random_id=_n//產(chǎn)生隨機(jī)序號preserve? ? ? ? keep random_id rep78//保留虛擬的rep78rename rep78 random_rep78? ? ? ? rename random_id id//重命名為id句携，以備與其他變量合并（merge）label var id 原數(shù)據(jù)與虛擬處理變量的唯一匹配碼? ? ? ? save random_rep78,replace? ? restore? ? ? ? drop random_digit random_id rep78//刪除原來的rep78rename obs_id id//重命名為id，以備與random_rep78合并（merge）label var id 原數(shù)據(jù)與虛擬處理變量的唯一匹配碼? ? ? ? save rawdata,replace*-合并融柬，回歸死嗦，提取系數(shù)userawdata,clear? ? ? ? merge1:1id using random_rep78,nogen? ? ? ? reg? price random_rep78 headroom weight length? ? ? ? g _b_random_rep78=_b[random_rep78]//提取x的回歸系數(shù)g _se_random_rep78=_se[random_rep78]//提取x的標(biāo)準(zhǔn)誤keep _b_random_rep78 _se_random_rep78? ? ? ? duplicates drop _b_random_rep78,force? ? ? ? save placebo`i',replace//把第i次placebo檢驗的系數(shù)和標(biāo)準(zhǔn)誤存起來}*-縱向合并1000次的系數(shù)和標(biāo)準(zhǔn)誤useplacebo1,clearforvalue i=2/1000{append using placebo`i'//縱向合并1000次回歸的系數(shù)及標(biāo)準(zhǔn)誤}gen tvalue=_b_random_rep78/_se_random_rep78kdensity tvalue,xtitle("t值")ytitle("分布")saving(placebo_test)*-刪除臨時文件forvalue i=1/1000{erase? placebo`i'.dta}*-mark：后續(xù)有時間再將它封裝成程序【完】
文獻(xiàn)范例
周茂、陸 毅粒氧、杜 艷越除、 姚 星，2018：《開發(fā)區(qū)設(shè)立于地區(qū)制造業(yè)升級》，《中國工業(yè)經(jīng)濟(jì)》第3期摘盆。
作者：虛童
鏈接：http://www.reibang.com/p/bad7471ab73b
來源：簡書
著作權(quán)歸作者所有翼雀。商業(yè)轉(zhuǎn)載請聯(lián)系作者獲得授權(quán)，非商業(yè)轉(zhuǎn)載請注明出處孩擂。
