use 地理信息_2022, clear
bys S: fillmissing L* R*
bys CITYCODE: fillmissing L* R*
save 0, replace
keep CITY
duplicates drop
g i = C+"政府办公室"
loc b = "********************"
cngcode, baidukey(`b') cit(C) add(i) long(lon) lat(lat)
save 1, replace
merge 1:m CITY using 0, nogen
order Stkcd year
so Stk y
geodist Lat Lng lat lon, g(dist)
compress
g GovInvol = ln(1+dist)
la var G 企业政府介入
la var lon 政府办公室经度
la var lat 政府办公室纬度
keep Stkcd year L* l* G
save 企业政府介入, replace
d
ta y
kdensity G
Contains data from 企业政府介入.dta
Observations: 58,073
Variables: 7 11 Jan 2024 09:49
---------------------------------------------------------------------------------------
Variable Storage Display Value
name type format label Variable label
---------------------------------------------------------------------------------------
Stkcd long %10.0g 股票代码
year int %10.0g
lon double %10.0g 政府办公室经度
lat double %10.0g 政府办公室纬度
Lng double %10.0g 办公地经度
Lat double %10.0g 办公地纬度
GovInvol float %9.0g 企业政府介入
---------------------------------------------------------------------------------------
Sorted by: Stkcd year
. ta y
year | Freq. Percent Cum.
------------+-----------------------------------
2000 | 1,088 1.87 1.87
2001 | 1,166 2.01 3.88
2002 | 1,223 2.11 5.99
2003 | 1,293 2.23 8.21
2004 | 1,379 2.37 10.59
2005 | 1,374 2.37 12.95
2006 | 1,457 2.51 15.46
2007 | 1,571 2.71 18.17
2008 | 1,625 2.80 20.97
2009 | 1,773 3.05 24.02
2010 | 2,128 3.66 27.68
2011 | 2,363 4.07 31.75
2012 | 2,494 4.29 36.05
2013 | 2,536 4.37 40.41
2014 | 2,652 4.57 44.98
2015 | 2,842 4.89 49.88
2016 | 3,281 5.65 55.52
2017 | 3,666 6.31 61.84
2018 | 3,767 6.49 68.32
2019 | 3,977 6.85 75.17
2020 | 4,427 7.62 82.80
2021 | 4,847 8.35 91.14
2022 | 5,144 8.86 100.00
------------+-----------------------------------
Total | 58,073 100.00