继续推荐1篇2023年第6期《数量经济技术经济研究》上DID方法主题文章

文章：劳动保护与企业智能制造转型——来自机器人进口的证据

来源：陈勇兵，李辉，林雄立．劳动保护与企业智能制造转型——来自机器人进口 的证据[J/OL]．数量经济技术经济研究. https://doi.org/10.13653/j.cnki.jqte.20230417.003

以事前企业所在行业劳动密集度作为处理组与对照 组的划分依据，借助 2008 年《劳动合同法》实施这一自然实验，设定如下双重差分估计模 型：

其中𝑌𝑖𝑗𝑡表示企业𝑖在𝑡年的工业机器人应用情况，若企业当年机 器人存量大于零，我们定义该企业应用了机器人（Adoption），取值为 1，否则为 0。

---

DID模型分析

clear 
use maindata,clear
xtset firm year
gen labor_protect=labor_inst*post
global controls "tariff_ind4 lnexport hhi lnfirm_number lnfixed_asset soe_share credit_share lnoutput"


//表1-基准回归
xi:reghdfe Robot labor_protect , absorb(firm year) vce(cluster cic4)
est store s1
xi:reghdfe Robot labor_protect  , absorb(firm cic2#year) vce(cluster cic4)
est store s2
xi:reghdfe Robot labor_protect $contols , absorb(firm cic2#year) vce(cluster cic4)
est store s3
xi:reghdfe Robot1 labor_protect $contols , absorb(firm cic2#year) vce(cluster cic4)
est store s4
xi:reghdfe lnallquantity labor_protect $contols , absorb(firm cic2#year) vce(cluster cic4)
est store s5
xi:reghdfe lnallquantity1 labor_protect $contols , absorb(firm cic2#year) vce(cluster cic4)
est store s6

一共得到六个方程

• 在第（1）列中，控制了企业固定效应以及年份固定效应，核心 解释变量𝑙𝑎𝑏𝑜𝑟_𝑖𝑛𝑠𝑡𝑗07 × 𝑃𝑜𝑠𝑡08的系数显著为正，表明劳动保护加强显著促进了高劳动密 集度行业中企业应用机器人。
• 第（2）列中进一步控制了行业-年份固定效应
• 第（3）列控制了一 系列行业层面控制变量

---

平行趋势检验

采用双重差分模型估计处理组企业与对照组企业的机器人应用在《劳动合同法》 实施前后的平均差异，前提条件就是处理组企业与对照组企业具有良好的可比性，即在 《劳动合同法》实施之前，两组企业的机器人应用并没有显著差异。

//图3-平行趋势检验
xi:reghdfe Robot lp_pre5 lp_pre4 lp_pre3 lp_pre2 lp_post0 lp_post1 lp_post2 lp_post3 lp_post4 lp_post5 indtrend $contols , absorb(firm  year#cic2) vce(cluster cic4)
forvalue i=2/5{
generate b_`i' = _b[lp_pre`i']
generate se_b_`i' = _se[lp_pre`i']
generate b_`i'LB = b_`i' - t0*se_b_`i'
generate b_`i'UB = b_`i' + t0*se_b_`i'
}

forvalue i=0/5{
generate b`i' = _b[lp_post`i']
generate se_b`i' = _se[lp_post`i']
generate b`i'LB = b`i' - t0*se_b`i'
generate b`i'UB = b`i' + t0*se_b`i'
}

generate b = .
generate LB = .
generate UB = .

forvalue i=2/5{
replace b = b_`i'  if post_year == -`i'

replace LB = b_`i'LB if post_year == -`i'

replace UB = b_`i'UB if post_year == -`i'

}

forvalue i=0/5{
replace b = b`i'  if post_year == `i'
replace LB = b`i'LB  if post_year == `i'
replace UB = b`i'UB  if post_year == `i'
}
keep post_year b LB UB
duplicates drop
sort post_year
drop if post_year<-5
replace b=0 if post_year==-1
replace LB=0 if post_year==-1
replace UB=0 if post_year==-1
twoway (connected b post_year, sort lcolor(navy) mcolor(navy) msymbol(circle_hollow) cmissing(n)) (rcap LB UB post_year, lcolor(navy) lpattern(dash) msize(medium)),ytitle(企业机器人应用) ytitle(, size(small)) yline(0, lwidth(vthin) lpattern(dash) lcolor(teal)) ylabel(, labsize(small) angle(horizontal) nogrid) xtitle(相对于劳动合同法颁布年份) xtitle(, size(small)) xline(-1, lwidth(vthin) lpattern(dash) lcolor(teal)) xlabel(-5(1)5, labsize(small)) xmtick(-5(1)5, nolabels ticks) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) 

安慰剂检验

• 参照 Topalova（2010）的做 法构建安慰剂检验，通过保留《劳动合同法》实施之前的样本，并分别使用 2003 年、2004 年、2005 年作为虚假的《劳动合同法》实施时间。

//表2-安慰剂检验
use maindata.dta,clear
xi:reghdfe Robot placebo_lp1 $contols if year<2008, absorb(firm cic2#year) vce(cluster cic4)
est store s1
xi:reghdfe Robot placebo_lp2 $contols if year<2008, absorb(firm cic2#year) vce(cluster cic4)
est store s2
xi:reghdfe Robot placebo_lp3 $contols if year<2008, absorb(firm cic2#year) vce(cluster cic4)
est store s3 

推荐文章：