【stata代码】地方政府驱动企业参与乡村振兴的机制——乡村振兴注意力视角的分析
地方政府驱动企业参与乡村振兴的机制——乡村振兴注意力视角的分析
主要介绍一下相关stata代码
完整代码
cd "F:\BaiduSyncdisk\论文\中国农村经济"
set more off
clear all
use "publish20250126.dta", clear
tsset stkid year
global y "Participation"
global x "RuralEr RuralPr"
global med "ESGm Em Sm Gm lnArgfirm"
global inter "lnTotal lnPositive lnNegtitve"
global cv "lnSize Leverage lnAge Tangible Cashflow ROA SOE Top1 lnBroad MShare"cap drop dummy*
tab ind, gen (dummy_ind)
tab year, gen (dummy_year)
drop dummy_ind1 dummy_year1*附表1 相关系数检验
global x "RuralEr"
pwcorr_a $y RuralPr $cv, star1(0.01) star5(0.05) star10(0.1) bonferroni*表2 描述性统计
tabstat Participation_org RuralPr ESGm Em Sm Gm Argfirm_org Total_org Positive_org Negtitve_org Size_org Leverage Age_org Tangible Cashflow ROA SOE Top1 Broad_org MShare, s(N mean sd min p75 max) f(%6.3f) c(s)*表3 基准
reghdfe Participation RuralPr, absorb(ind year) vce(r)
est store m1
reghdfe Participation RuralPr lnSize Leverage lnAge Tangible Cashflow ROA, absorb(ind year) vce(r)
est store m2
reghdfe Participation RuralPr $cv, absorb(ind year) vce(r)
est store m3local model "m1 m2 m3"
local mtitle "Participation Participation Participation"
local order "RuralPr $cv"
esttab `model', mtitle(`mtitle') compress nogap star(* 0.1 ** 0.05 *** 0.01) ///
b(%6.3f) se(%6.3f) obslast scalars(r2_a F) order(`order')*表4 稳健性 x/word y:0-1
reghdfe Participation RuralEr $cv, absorb(ind year) vce(r)
est store m1
reghdfe Participation lnWordP $cv, absorb(ind year) vce(r)
est store m2
reghdfe Participation lnWordE $cv, absorb(ind year) vce(r)
est store m3
cap drop Participation_dum
gen Participation_dum=Participation>0
logit Participation_dum RuralPr $cv dummy*, vce(r)
est store m4
poisson Participation_dum RuralPr $cv dummy*, vce(r)
est store m5
reghdfe Participation RuralEr $cv, absorb(ind year provid) vce(r)
est store m6local model "m1 m2 m3 m4 m5 m6"
local mtitle "Participation Participation Participation Participation_dum Participation_dum Participation"
local order "RuralEr lnWordP lnWordE RuralPr $cv"
esttab `model', mtitle(`mtitle') compress nogap star(* 0.1 ** 0.05 *** 0.01) ///
b(%6.3f) se(%6.3f) obslast scalars(r2_a r2_p F) order(`order') ///
indicate("Industry FE=dummy_ind*" "Year FE=dummy_year*")*表5 替换研究样本
use "cpes.dta", clear
local cv "Connect Gender Age Edu Staff Firmage Family Party Roa"reghdfe Participation PRuralPr, absorb(Ind) vce(r)
est store m1
reghdfe Participation PRuralPr Connect Gender Age Edu, absorb(Ind) vce(r)
est store m2
reghdfe Participation PRuralPr `cv', absorb(Ind) vce(r)
est store m3local model "m1 m2 m3"
local mtitle "Participation Participation Participation"
local order "PRuralPr `cv'"
esttab `model', mtitle(`mtitle') compress nogap star(* 0.1 ** 0.05 *** 0.01) ///
b(%6.3f) se(%6.3f) obslast scalars(r2_a F) order(`order')use "publish20250126.dta", clear
*表6 内生性:地表起伏度
gen RdlsIV=Rdls*pop
reghdfe RuralPr RdlsIV $cv, absorb(ind year) vce(r)
est store m1
ivreghdfe Participation (RuralPr=RdlsIV) $cv, absorb(ind year) r
est store m1local model "m1 m2"
local mtitle "RuralPr Participation"
local order "RdlsIV RuralPr $cv"
esttab `model', mtitle(`mtitle') compress nogap star(* 0.1 ** 0.05 *** 0.01) ///
b(%6.3f) se(%6.3f) obslast scalars(r2_a F) order(`order')*表7 机制
reghdfe ESGm RuralPr $cv, absorb(ind year) vce(r)
est store m1
reghdfe Em RuralPr $cv, absorb(ind year) vce(r)
est store m2
reghdfe Sm RuralPr $cv, absorb(ind year) vce(r)
est store m3
reghdfe Gm RuralPr $cv, absorb(ind year) vce(r)
est store m4
reghdfe lnArgfirm RuralPr $cv, absorb(ind year) vce(r)
est store m5local model "m1 m2 m3 m4 m5"
local mtitle "ESGm Em Sm Gm lnArgfirm"
local order "RuralPr $cv"
esttab `model', mtitle(`mtitle') compress nogap star(* 0.1 ** 0.05 *** 0.01) ///
b(%6.3f) se(%6.3f) obslast scalars(r2_a F) order(`order')*表8 调节
egen stdRuralPr=std(RuralPr)
egen stdlnTotal=std(lnTotal)
egen stdlnPositive=std(lnPositive)
egen stdlnNegtitve=std(lnNegtitve)reghdfe Participation c.stdRuralPr##c.stdlnTotal $cv, absorb(ind year) vce(r)
est store m1
reghdfe Participation c.stdRuralPr##c.stdlnPositive $cv, absorb(ind year) vce(r)
est store m2
reghdfe Participation c.stdRuralPr##c.stdlnNegtitve $cv, absorb(ind year) vce(r)
est store m3local model "m1 m2 m3"
local mtitle "Participation Participation Participation Participation"
local order "stdRuralPr stdlnTotal stdlnPositive stdlnNegtitve c.stdRuralPr* $cv"
esttab `model', mtitle(`mtitle') compress nogap star(* 0.1 ** 0.05 *** 0.01) ///
b(%6.3f) se(%6.3f) obslast scalars(r2_a F) order(`order')*表9 异质性
**产权异质 是否国有企业
reghdfe Participation RuralPr $cv if SOE==1, absorb(ind year) vce(r)
est store m1
reghdfe Participation RuralPr $cv if SOE==0, absorb(ind year) vce(r)
est store m2
**行业异质 农业与非农业
reghdfe Participation RuralPr $cv if ind=="A", absorb(ind year) vce(r)
est store m3
reghdfe Participation RuralPr $cv if ind!="A", absorb(ind year) vce(r)
est store m4
**行业异质 是否强制披露
reghdfe Participation RuralPr $cv if Force==1, absorb(ind year) vce(r)
est store m5
reghdfe Participation RuralPr $cv if Force==0, absorb(ind year) vce(r)
est store m6local model "m1 m2 m3 m4 m5 m6"
local mtitle "SOE NSOE Arg NArg Force NForce"
local order "RuralPr $cv"
esttab `model', mtitle(`mtitle') compress nogap star(* 0.1 ** 0.05 *** 0.01) ///
b(%6.3f) se(%6.3f) obslast scalars(r2_a F) order(`order')**表10 地区异质
**东部10省(市)包括北京、天津、河北、上海、江苏、浙江、福建、山东、广东和海南;中部6省包括山西、安徽、江西、河南、湖北和湖南;西部12省(区、市)包括内蒙古、广西、重庆、四川、贵州、云南、西藏、陕西、甘肃、青海、宁夏和新疆;东北3省包括辽宁、吉林和黑龙江。
reghdfe Participation RuralPr $cv if Region==1, absorb(ind year) vce(r)
est store m1
reghdfe Participation RuralPr $cv if Region==2, absorb(ind year) vce(r)
est store m2
reghdfe Participation RuralPr $cv if Region==3, absorb(ind year) vce(r)
est store m3
reghdfe Participation RuralPr $cv if Region==4, absorb(ind year) vce(r)
est store m4local model "m1 m2 m3 m4"
local mtitle "East Central West Northeast"
local order "RuralPr $cv"
esttab `model', mtitle(`mtitle') compress nogap star(* 0.1 ** 0.05 *** 0.01) ///
b(%6.3f) se(%6.3f) obslast scalars(r2_a F) order(`order')