學(xué)而時(shí)習(xí)之，不亦說乎,？之前在B站上看到了凱叔的視頻,，演示了如何利用stata計(jì)算公司注冊地周圍一定范圍里的的開發(fā)區(qū)數(shù)量。講解非常細(xì)致,，但是視頻較為模糊,。由此，自己聯(lián)想到了財(cái)會(huì)中的非正式制度領(lǐng)域——企業(yè)所面臨的地區(qū)文化指標(biāo),，是否可以以視頻代碼作為基礎(chǔ),，利用stata進(jìn)行計(jì)算呢？自己進(jìn)行了實(shí)驗(yàn),，并取得了結(jié)果,。
既有文獻(xiàn)對(duì)于文化的衡量一般采用計(jì)算以公司注冊地為圓心，以特定距離為半徑所涵蓋的道觀,、孔廟等的數(shù)量的方式。如,，杜興強(qiáng)等(2016)發(fā)表在《會(huì)計(jì)研究》的文章中,，“參考 Du(2013) 和陳冬華等 (2013)，采用地理近鄰性度量宗教的方法構(gòu)建宗教變量,。首先,，獲取公司注冊地的經(jīng)緯度; 其次，獲取佛教寺廟所在地的經(jīng)緯度,；第三,， 根據(jù)公司和寺廟的經(jīng)緯度，計(jì)算公司注冊地與寺廟的距離,。最終,，統(tǒng)計(jì)公司方圓100、200,、300公里內(nèi)全國重點(diǎn)寺廟的個(gè)數(shù),，并據(jù)此定義本文的解釋變量BUD100、BUD200 和 BUD300,?！?/span>徐細(xì)雄和李萬利（2019）發(fā)表在《金融研究》的文章中“若企業(yè)注冊地一定半徑范圍內(nèi)分布的儒家書院數(shù)量越多,， 則表明該企業(yè)受到的儒家文化影響越強(qiáng)。我們分別使用企業(yè)注冊地200(Confu_200) 和300公里(Confu_300) 半徑范圍內(nèi)儒家書院數(shù)量作為儒家文化強(qiáng)度的代理變量,，并加 1 取自然對(duì)數(shù)進(jìn)行標(biāo)準(zhǔn)化處理”,。……從中不難看出,，文化指標(biāo)計(jì)算的重點(diǎn)是如何計(jì)算公司注冊地與特定宗教建筑的距離,。以下對(duì)數(shù)據(jù)來源以及計(jì)算代碼進(jìn)行說明：有關(guān)于上市公司注冊地經(jīng)緯度的數(shù)據(jù)來源于CNRDS數(shù)據(jù)庫-公司特色庫-公司年度信息-CAID。本文暫時(shí)采用孔廟數(shù)量來代指儒家文化,?？讖R的經(jīng)緯度數(shù)據(jù)來源于CNRDS數(shù)據(jù)庫-經(jīng)濟(jì)特色庫-儒家文化。
首先用Stata對(duì)兩份數(shù)據(jù)進(jìn)行預(yù)處理：
cd C:\Users\Dell\Desktop\經(jīng)緯度 //指定存儲(chǔ)路徑
import excel 孔廟遺存.xlsx,first clear cellrange(A2) //導(dǎo)入下載的孔廟EXCEL數(shù)據(jù)destring _all,replace //將表中的經(jīng)緯度數(shù)據(jù)由字符型轉(zhuǎn)換為數(shù)值型ren 經(jīng)度 kongmiaolng //重命名經(jīng)度名稱,，以免同上市公司注冊的經(jīng)度變量名混淆ren 緯度 kongmiaolat //重命名緯度名稱,，以免同上市公司注冊的緯度變量名混淆save 孔廟遺存,replace //將EXCEL數(shù)據(jù)導(dǎo)出為DTA格式
import excel 上市公司年度信息.xlsx,first clear cellrange(A2) //導(dǎo)入下載的上市公司基本信息EXCEL數(shù)據(jù)keep 股票代碼 會(huì)計(jì)年度 經(jīng)度 緯度 所在省份 所在地級(jí)市 //保留主要關(guān)注的變量destring _all,replace //將表中的經(jīng)緯度數(shù)據(jù)由字符型轉(zhuǎn)換為數(shù)值型ren 經(jīng)度 firmlng //重命名經(jīng)度名稱ren 緯度 firmlat //重命名緯度名稱save 上市公司注冊地址,replace //將EXCEL數(shù)據(jù)導(dǎo)出為DTA格式
隨后，利用geodist命令,、cross命令以及循環(huán)語句計(jì)算出注冊地對(duì)于所有孔廟的距離,。cross命令主要是將每條上市公司-年度注冊地地址與所有的孔廟數(shù)據(jù)進(jìn)行合并。自己的數(shù)據(jù)中,，有關(guān)上市公司-年度注冊地經(jīng)緯度觀測值共有47251條,，孔廟數(shù)據(jù)共有200條，合并后共計(jì)47,251×200=9,450,200條觀測值,。合并后的樣本大致長這樣： 
采用上市公司-年度的形式而不只用上市公司形式的原因在于上市公司可能在樣本期內(nèi)發(fā)生注冊地變更,。與之相關(guān)的論文為全怡等（2019）發(fā)表在《會(huì)計(jì)研究》的《注冊地變更與企業(yè)費(fèi)用粘性——基于政策性優(yōu)惠和負(fù)擔(dān)的視角》，該文證實(shí)了上市公司存在注冊地變更行為并且影響企業(yè)成本管理,。 代碼的邏輯在于,，將合并后的數(shù)據(jù)中，抽取出每條觀測值進(jìn)行計(jì)算保存,，之后將計(jì)算結(jié)果合并,。 ssc install geodist,replace //計(jì)算之前首先下載/更新geodist命令
use 上市公司注冊地址,clearglobal I=_N //計(jì)算最大樣本數(shù)，并將其存放在I的全局元中g num=_n //計(jì)算每條樣本所對(duì)應(yīng)的行forvalues i=1(1) $I { preserve qui keep if num==`i' qui cross using 孔廟遺存.dta //將每條上市公司-年度注冊地地址都與所有的孔廟數(shù)據(jù)進(jìn)行合并 geodist firmlat firmlng kongmiaolat kongmiaolng,gen(dist) //采用geodist命令計(jì)算,，并將距離存放在dist變量中 sort dist qui save dist`i',replace //將每一條觀測值計(jì)算出的距離存放在dist中 restore}
clearforval i=1(1) $I{ append using dist`i'} //縱向合并所有的計(jì)算結(jié)果save 上市公司——孔廟距離,replace
clearforval i=1(1) $I{ erase dist`i'.dta} //清理所有的臨時(shí)文件
由于樣本量較大,，因此計(jì)算可能費(fèi)一些時(shí)間。合并后,，樣本長這樣：
最后一列的dist即為最終結(jié)果,，在geodist中，如果沒有特別設(shè)定,，默認(rèn)距離是以km作為單位的,。因此，現(xiàn)在就可以根據(jù)參考文獻(xiàn)的設(shè)定,，劃定半徑范圍以及計(jì)算孔廟個(gè)數(shù)了,。以方圓100km為例子：
use 上市公司——孔廟距離,clear
preservekeep if dist<=100bys 股票代碼 會(huì)計(jì)年度:gen N=_Ng Culture=ln(1+N)keep 股票代碼 會(huì)計(jì)年度 Cultureduplicates dropsave dist100,replace //首先先把方圓100km內(nèi)存在孔廟的樣本計(jì)算出culturerestore
keep 股票代碼 會(huì)計(jì)年度duplicates dropmerge 1:1 股票代碼 會(huì)計(jì)年度 using dist100replace Culture=0 if _m==1 //定義未匹配到的觀測值為Culture=0,，及方圓100km內(nèi)不存在孔廟drop if _m==2drop _msave Culture,replace
最終得到全部47,251份樣本的儒家文化測度值，樣本大致長這樣：
在公司注冊地經(jīng)緯度已知的情況下,，可以變更孔廟,，換成道觀、佛寺,、清真寺等同樣適合本文代碼,。
以上代碼較為通俗，但是對(duì)于更大樣本的數(shù)據(jù),，通常stata中的mata矩陣是一個(gè)更加快捷的方法,。
以上是建立在特定地址的經(jīng)緯度已知的情況，如果經(jīng)緯度未知,。則需參照Stata：批量獲取經(jīng)緯度數(shù)據(jù)-空間計(jì)量| 連享會(huì)主頁 (lianxh.cn)首先獲取經(jīng)緯度數(shù)據(jù),。 以上代碼建立在B站凱叔的視頻中，自己根據(jù)需要進(jìn)行改動(dòng),。非營利目的,，且在文后注明來源，故不涉及侵權(quán)行為,。
[1]徐細(xì)雄,李萬利.儒家傳統(tǒng)與企業(yè)創(chuàng)新:文化的力量[J].金融研究,2019(09):112-130.[2]全怡,嚴(yán)麗娜,劉磊.注冊地變更與企業(yè)費(fèi)用粘性——基于政策性優(yōu)惠和負(fù)擔(dān)的視角[J].會(huì)計(jì)研究,2019(08):47-54.[3]杜興強(qiáng),蹇薇,曾泉,?，摤?宗教影響、控股股東與過度投資:基于中國佛教的經(jīng)驗(yàn)證據(jù)[J].會(huì)計(jì)研究,2016(08):50-57+97.公司金融研究_STATA大批量計(jì)算兩組地址之間地理距離的幾種思路_嗶哩嗶哩_bilibili參考網(wǎng)站：Stata：我和她離多遠(yuǎn),？基于百度地圖API的地理距離計(jì)算| 連享會(huì)主頁 (lianxh.cn)
|
