上市公司儒家文化影响指标（孔庙数量法）

///更多数据尽在丘知数据分析会员之家： www.qqrdata.com
//导入原始数据1
use 原始数据1,clear
label variable 丘知数据分析会员之家 "www.qqrdata.com"
mata
mata clear
//转化为弧度
KLongit=st_data(.,("经度")):*(pi()/180)
KLatit=st_data(.,("纬度")):*(pi()/180)
end

///更多数据尽在丘知数据分析会员之家 ：www.qqrdata.com
//导入原始数据2
use 原始数据2,clear
destring 注册地经度 注册地纬度,replace
label variable 丘知数据分析会员之家 "www.qqrdata.com"
mata
//转化为弧度
FLongit=st_data(.,("注册地经度")):*(pi()/180)  
FLatit=st_data(.,("注册地纬度")):*(pi()/180)
//计算公司与孔庙之间距离
distance=acos(sin(FLatit)*sin(KLatit')+(cos(FLatit)*cos(KLatit'):*(cos(FLongit)*cos(KLongit')+sin(FLongit)*sin(KLongit')))):*6371    

//更多数据尽在丘知数据分析会员之家 ：www.qqrdata.com
//A、B、C为矩阵，如果distance对应位置的元素小于100km、200km、300km，则A、B、C对应位置的元素取1，否则为0
A=distance:<=100 
B=distance:<=200
C=distance:<=300

//计算矩阵A、B、C的行总和，即距离小于100km、200km、300km的孔庙个数