上市公司环保费改税政策DID变量

///更多stata学习资料尽在丘知数据分析会员之家，网址：www.qqrdata.com
//导入原始数据
use 原始数据,clear
duplicates drop stkcd year, force
label data "更多stata学习资料尽在丘知数据分析会员之家：www.qqrdata.com"

///更多stata学习资料尽在丘知数据分析会员之家，网址：www.qqrdata.com
//生成虚拟变量，是否重污染行业企业（行业是否属于B06、B07、B08、B09、C17、C19、C22、C25、C26、C28、C29、C30、C31、C32、D44，是为1；否为0）
gen 是否重污染行业企业 = 0
foreach i in B06 B07 B08 B09 C17 C19 C22 C25 C26 C28 C29 C30 C31 C32 D44{
replace 是否重污染行业企业= 1 if regexm(行业代码, "`i'")
}
label var 是否重污染行业企业 "行业是否属于B06、B07、B08、B09、C17、C19、C22、C25、C26、C28、C29、C30、C31、C32、D44，是为1；否为0"
ren 是否重污染行业企业 pollute

///更多stata学习资料尽在丘知数据分析会员之家，网址：www.qqrdata.com 
//生成treat变量，企业所在地区为河北、江苏、山东、河南、湖南、四川、重庆、贵州、海南、广西、山西、北京的，为税费负担提标地区，treat变量赋值为1。而企业所在地区为湖北、浙江、辽宁、福建、吉林、安徽、江西、陕西、甘肃、新疆、西藏、宁夏、青海、内蒙古、黑龙江、云南、天津、上海和广东的，为税费负担平移地区，treat变量赋值为0。
g treat = 0
foreach i in 河北 江苏 山东 河南 湖南 四川 重庆 贵州 海南 广西 山西 北京 {
replace treat = 1 if regexm(所属省份, "`i'")
}
lab var treat "税费负担提标地区，treat变量赋值为1；税费负担平移地区，treat变量赋值为0"

//生成after变量，表示《环保税法》实施前后的虚拟变量，即2018年当年及之后取值为1，否则为0。
g after = year >= 2018
lab var after "2018年当年及之后取值为1，否则为0"

///更多stata学习资料尽在丘知数据分析会员之家，网址：www.qqrdata.com 
//生成DID变量，即交互项treat×after为核心解释变量。
g did = treat * after