* 图2
use $source/数据.dta, clear
gen tt=-4 if b4==1|b5==1|b6==1|b7==1|b8==1
replace tt=-3 if b3==1
replace tt=-2 if b2==1
replace tt=-1 if b1==1
replace tt=0 if t0==1
replace tt=1 if a1==1
replace tt=2 if a2==1
replace tt=3 if a3==1
replace tt=4 if a4==1|a5==1
tab tt, g(tt)
forvalue i=1/9{
replace tt`i'=0 if tt==.
}
xtreg Patent tt1 tt2 tt3 tt4 tt6 tt7 tt8 tt9 Size Age Lev Capex Agdp Pop i.year, fe cluster(stkcd)
cap drop b
cap drop se b_up b_low id
gen b=.
gen se=.
forvalue i=1/4{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=1/4{
replace se=_se[tt`i'] if _n==`i'
}
replace b=0 if _n==5
replace se=0 if _n==5
forvalue i=6/9{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=6/9
{
replace se=_se[tt`i'] if _n==`i'
}
gen b_up=b+1.64*se
gen b_low=b-1.64*se
gen idd=_n
scatter b idd if idd<=9, connect(l) lcolor(black) lwidth(*0.6) lp(solid) msymbol(circle) msize(*0.4) mcolor(black) || (rcap b_low b_up idd if idd<=9, lcolor(black) lwidth(*0.6))
,yline(0, lcolor(black) lwidth(0.01)) xline(5, lwidth(0.01) lcolor(black) lp(shortdash)) ylabel(-0.16 "-0.16" 0 "0" 0.16 "0.16" 0.32 "0.32", labsize(small) tl(0.8) nogrid)
xlabel(1 "≤-4" 2 "-3" 3 "-2" 4 "-1" 5 "0" 6 "1" 7 "2" 8 "3" 9 "≥4", labsize(small) tl(0.8))
xtitle("政策实施相对时间""a", box margin(small) size(small) fcolor(white) lc(white))
ytitle("估计系数大小", box margin(small) size(small) fcolor(white) lc(white))
graphregion(color(white))
legend(off)
name(Patent, replace)
use $source/数据.dta, clear
gen tt=-4 if b4==1|b5==1|b6==1|b7==1|b8==1
replace tt=-3 if b3==1
replace tt=-2 if b2==1
replace tt=-1 if b1==1
replace tt=0 if t0==1
replace tt=1 if a1==1
replace tt=2 if a2==1
replace tt=3 if a3==1
replace tt=4 if a4==1|a5==1
tab tt, g(tt)
forvalue i=1/9{
replace tt`i'=0 if tt==.
}
xtreg Citation tt1 tt2 tt3 tt4 tt6 tt7 tt8 tt9 Size Age Lev Capex Agdp Pop i.year, fe cluster(stkcd)
cap drop b
cap drop se b_up b_low id
gen b=.
gen se=.
forvalue i=1/4{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=1/4{
replace se=_se[tt`i'] if _n==`i'
}
replace b=0 if _n==5
replace se=0 if _n==5
forvalue i=6/9{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=6/9{
replace se=_se[tt`i'] if _n==`i'
}
gen b_up=b+1.64*se
gen b_low=b-1.64*se
gen idd=_n
scatter b idd if idd<=9, connect(l) lcolor(black) lwidth(*0.6) lp(solid) msymbol(circle) msize(*0.4) mcolor(black) || (rcap b_low b_up idd if idd<=9, lcolor(black) lwidth(*0.6))
,yline(0, lcolor(black) lwidth(0.01)) xline(5, lwidth(0.01) lcolor(black) lp(shortdash)) ylabel(-0.15 "-0.15" 0 "0" 0.15 "0.15" 0.30 "0.30", labsize(small) tl(0.8) nogrid)
xlabel(1 "≤-4" 2 "-3" 3 "-2" 4 "-1" 5 "0" 6 "1" 7 "2" 8 "3" 9 "≥4", labsize(small) tl(0.8))
xtitle("政策实施相对时间""b", box margin(small) size(small) fcolor(white) lc(white))
ytitle("估计系数大小", box margin(small) size(small) fcolor(white) lc(white))
graphregion(color(white))
legend(off)
name(Citation, replace)
use $source/数据.dta, clear
gen tt=-4 if b4==1|b5==1|b6==1|b7==1|b8==1
replace tt=-3 if b3==1
replace tt=-2 if b2==1
replace tt=-1 if b1==1
replace tt=0 if t0==1
replace tt=1 if a1==1
replace tt=2 if a2==1
replace tt=3 if a3==1
replace tt=4 if a4==1|a5==1
tab tt, g(tt)
forvalue i=1/9{
replace tt`i'=0 if tt==.
}
xtreg Fcitation1 tt1 tt2 tt3 tt4 tt6 tt7 tt8 tt9 Size Age Lev Capex Agdp Pop i.year, fe cluster(stkcd)
cap drop b
cap drop se b_up b_low id
gen b=.
gen se=.
forvalue i=1/4{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=1/4{
replace se=_se[tt`i'] if _n==`i'
}
replace b=0 if _n==5
replace se=0 if _n==5
forvalue i=6/9{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=6/9{
replace se=_se[tt`i'] if _n==`i'
}
gen b_up=b+1.64*se
gen b_low=b-1.64*se
gen idd=_n
scatter b idd if idd<=9, connect(l) lcolor(black) lwidth(*0.6) lp(solid) msymbol(circle) msize(*0.4) mcolor(black) || (rcap b_low b_up idd if idd<=9, lcolor(black) lwidth(*0.6))
,yline(0, lcolor(black) lwidth(0.01)) xline(5, lwidth(0.01) lcolor(black) lp(shortdash)) ylabel(-0.23 "-0.23" 0 "0" 0.23 "0.23" 0.46 "0.46", labsize(small) tl(0.8
) nogrid)
xlabel(1 "≤-4" 2 "-3" 3 "-2" 4 "-1" 5 "0" 6 "1" 7 "2" 8 "3" 9 "≥4", labsize(small) tl(0.8))
xtitle("政策实施相对时间""c", box margin(small) size(small) fcolor(white) lc(white))
ytitle("估计系数大小", box margin(small) size(small) fcolor(white) lc(white))
graphregion(color(white))
legend(off)
name(Fcitation1, replace)
use $source/数据.dta, clear
gen tt=-4 if b4==1|b5==1|b6==1|b7==1|b8==1
replace tt=-3 if b3==1
replace tt=-2 if b2==1
replace tt=-1 if b1==1
replace tt=0 if t0==1
replace tt=1 if a1==1
replace tt=2 if a2==1
replace tt=3 if a3==1
replace tt=4 if a4==1|a5==1
tab tt, g(tt)
forvalue i=1/9{
replace tt`i'=0 if tt==.
}
xtreg Fcitation2 tt1 tt2 tt3 tt4 tt6 tt7 tt8 tt9 Size Age Lev Capex Agdp Pop i.year, fe cluster(stkcd)
cap drop b
cap drop se b_up b_low id
gen b=.
gen se=.
forvalue i=1/4{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=1/4{
replace se=_se[tt`i'] if _n==`i'
}
replace b=0 if _n==5
replace se=0 if _n==5
forvalue i=6/9{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=6/9{
replace se=_se[tt`i'] if _n==`i'
}
gen b_up=b+1.64*se
gen b_low=b-1.64*se
gen idd=_n
scatter b idd if idd<=9, connect(l) lcolor(black) lwidth(*0.6) lp(solid) msymbol(circle) msize(*0.4) mcolor(black) || (rcap b_low b_up idd if idd<=9, lcolor(black) lwidth(*0.6))
,yline(0, lcolor(black) lwidth(0.01)) xline(5, lwidth(0.01) lcolor(black) lp(shortdash)) ylabel(-0.20 "-0.20" 0 "0" 0.20 "0.20" 0.40 "0.40", labsize(small) tl(0.8) nogrid)
xlabel(1 "≤-4" 2 "-3" 3 "-2" 4 "-1" 5 "0" 6 "1" 7 "2" 8 "3" 9 ":≥4", labsize(small) tl(0.8))
xtitle("政策实施相对时间""d", box margin(small) size(small) fcolor(white) lc(white))
ytitle("估计系数大小", box margin(small) size(small) fcolor(white) lc(white))
graphregion(color(white))
legend(off)
name(Fcitation2, replace)
use $source/数据.dta, clear
gen tt=-4 if b4==1|b5==1|b6==1|b7==1|b8==1
replace tt=-3 if b3==1
replace tt=-2 if b2==1
replace tt=-1 if b1==1
replace tt=0 if t0==1
replace tt=1 if a1==1
replace tt=2 if a2==1
replace tt=3 if a3==1
replace tt=4 if a4==1|a5==1
tab tt, g(tt)
forvalue i=1/9{
replace tt`i'=0 if tt==.
}
xtreg RD tt1 tt2 tt3 tt4 tt6 tt7 tt8 tt9 Size Age Lev Capex Agdp Pop i.year, fe cluster(stkcd)
cap drop b
cap drop se b_up b_low id
gen b=.
gen se=.
forvalue i=1/4{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=1/4{
replace se=_se[tt`i'] if _n==`i'
}
replace b=0 if _n==5
replace se=0 if _n==5
forvalue i=6/9{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=6/9{
replace se=_se[tt`i'] if _n==`i'
}
gen b_up=b+1.64*se
gen b_low=b-1.64*se
gen idd=_n
scatter b idd if idd<=9, connect(l) lcolor(black) lwidth(*0.6) lp(solid) msymbol(circle) msize(*0.4) mcolor(black) || (rcap b_low b_up idd if idd<=9, lcolor(black) lwidth(*0.6))
,yline(0, lcolor(black) lwidth(0.01)) xline(5, lwidth(0.01) lcolor(black) lp(shortdash)) ylabel(-0.18 "-0.18" 0 "0" 0.18 "0.18" 0.36 "0.36", labsize(small) tl(0.8) nogrid)
xlabel(1 "≤-4" 2 "-3" 3 "-2" 4 "-1" 5 "0" 6 "1" 7 "2" 8 "3" 9 "≥4", labsize(small) tl(0.8))
xtitle("政策实施相对时间""e", box margin(small) size(small) fcolor(white) lc(white))
ytitle("估计系数大小", box margin(small) size(small) fcolor(white) lc(white))
graphregion(color(white))
legend(off)
name(RD, replace)
use $source/数据.dta, clear
gen tt=-4 if b4==1|b5==1|b6==1|b7==1|b8==1
replace tt=-3 if b3==1
replace tt=-2 if b2==1
replace tt=-1 if b1==1
replace tt=0 if t0==1
replace tt=1 if a1==1
replace tt=2 if a2==1
replace tt=3 if a3==1
replace tt=4 if a4==1|a5==1
tab tt, g(tt)
forvalue i=1/9{
replace tt`i'=0 if tt==.
}
xtreg Subsidy tt1 tt2 tt3 tt4 tt6 tt7 tt8 tt9 Size Age Lev Capex Agdp Pop i.year, fe cluster(stkcd)
cap drop b
cap drop se b_up b_low id
gen b=.
gen se=.
forvalue i=1/4{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=1/4{
replace se=_se[tt`i'] if _n==`i'
}
replace b=0 if _n==5
replace se=0 if _n==5
forvalue i=6/9{
replace b=_b[tt`i'] if _n==`i'
}
forvalue i=6/9{
replace se=_se[tt`i'] if _n==`i'
}
gen b_up=b+1.64*se
gen b_low=b-1.64*se
gen idd=_n
scatter b idd if idd<=9, connect(l) lcolor(black) lwidth(*0.6) lp(solid) msymbol(circle) msize(*0.4) mcolor(black) || (rcap b_low b_up idd if idd<=9, lcolor(black) lwidth(*0.6))
,yline(0, lcolor(black) lwidth(0.01)) xline(5, lwidth(0.01) lcolor(black) lp(shortdash)) ylabel(-0.01 "-0.01" 0 "0" 0.01 "0.01" 0.02 "0.02", labsize(small) tl(0.8) nogrid)
xlabel(1 "≤-4" 2 "-3" 3 "-2" 4 "-1" 5 "0" 6 "1" 7 "2" 8 "3" 9 "≥4", labsize(small) tl(0.8))
xtitle("政策实施相对时间""f", box margin(small) size(small) fcolor(white) lc(white))
ytitle("估计系数大小", box margin(small) size(small) fcolor(white) lc(white))
graphregion(color(white))
legend(off)
name(Subsidy, replace)
graph combine Patent Citation Fcitation1 Fcitation2 RD Subsidy, graphregion(color(white))