温馨提示: 定期 清理浏览器缓存,可以获得最佳浏览体验。
New!
lianxh
命令发布了:
随时搜索连享会推文、Stata 资源,安装命令如下:
. ssc install lianxh
使用详情参见帮助文件 (有惊喜):
. help lianxh
作者: 冷萱 (西南财经大学)
邮箱: lx@smail.swufe.edu.cn
目录
Stata 功能的强大不只是在于经济学中的因果推断,同时也能兼顾一些其他软件的功能,如用 ArcGIS 计算不同企业 (学校) 之间的距离。本文主要介绍如何通过 Stata 计算不同学校之间的距离。在此之前,也需要一些准备工作,如申请百度秘钥和数据整理。如果学会了计算不同学校的距离,那么将工业企业地理位置转换为经纬度,以及进行相关距离测算也就变得容易了。
在正式分析之前,我们需要准备如下材料:
cngcode
命令cross
命令geodist
命令
首先,我们需要在「百度地图开放平台」申请秘钥,具体操作详见「百度地图 api 密钥获取」。秘钥的主要作用是通过调取百度地图将地理信息转换为经纬度,并且个人用户每天有一定的额度限制,但也能满足一般科研需求。感兴趣的读者可自行申请。
本节数据整理主要分为两部分内容:
为了后续方便,本文分析对象仅为教育部部属高校。具体代码如下:
*数据地址:https://gitee.com/arlionn/data/blob/master/data01/W020200709292792106069.xls
clear all
import excel "W020200709292792106069.xls", sheet("sheet1") clear
nrow 3 //将第3行设置为变量名
drop if 学校名称 ==""
keep if 主管部门=="教育部" //仅保留部属高校
drop 备注
destring _all ,replace
format 学校标识码 %15.0f
compress
在原始数据整理好后,通过 cngcode
命令生成高校对应的经纬度。命令语法如下:
cngcode, baidukey(string) [options]
baidukey(string)
指定申请的百度秘钥;[option]
主要选项又包括:
address(varname)
指定包含地址字符串的变量;lat(varname) long(varname)
指定输出的经纬度变量名。具体而言,本文分析对应的命令如下:
. ssc install cngcode, replace
. cngcode, baidukey(DfsFvAXxxEVCElaB69iLXEen7OKOt9gB) ///
lat(startlat0) long(startlng0) address(学校名称 )
. save university.dta,replace
执行上述命令后会生成两个变量分别为 startlat0 和 startlng0。需要注意的是,经纬度是根据变量 学校名称 生成。
在原始数据清洗后,得到 76 所部属高校。本文的目的在于分析第 1 所高校与其他 75 所高校之间的距离,第 2 所高校与其他 75 所高校的距离。因此,需要将数据进行转换,使 76×1 的数据集变为 76×76 的数据集,使用到的命令为 cross
。
clear all
use "university.dta",clear
rename (学校名称 startlng0 startlat0) (name startlng0 startlat0)
save name, replace
use "university.dta",clear
rename (学校名称 startlng0 startlat0) (学校名称 startlng1 startlat1)
cross using "name.dta"
最终,我们将得到 76×76 个样本值,具体如下图所示:
在完成以上数据整理工作后,距离计算就显得轻而易举了。我们使用 geodist
命令即可实现,并且可以指定新生成的距离变量为 distance。
. ssc install geodist, replace
. geodist startlat1 startlng1 startlat0 startlng0, gen(distance)
. save distance, replace
clear all
import excel "W020200709292792106069.xls", sheet("sheet1") clear
nrow 3 //将第3行设置为变量名
drop if 学校名称 ==""
keep if 主管部门=="教育部" //仅保留部属高校
drop 备注
destring _all ,replace
format 学校标识码 %15.0f
compress
cngcode, baidukey(DfsFvAXxxEVCElaB69iLXEen7OKOt9gB) lat(startlat0) long(startlng0) address(学校名称 ) //生成经纬度
save university.dta,replace
clear all
use university.dta,clear
rename (学校名称 startlng0 startlat0) (name startlng0 startlat0)
save name, replace
use university.dta,clear
rename (学校名称 startlng0 startlat0) (学校名称 startlng1 startlat1)
cross using "name.dta"
geodist startlat1 startlng1 startlat0 startlng0, gen(distance) //计算距离
save distance, replace
Note: 产生如下推文列表的命令为:
lianxh 地理 API 地图, m
安装最新版lianxh
命令:
ssc install lianxh, replace
连享会-直播课 上线了!
http://lianxh.duanshu.com
免费公开课:
直击面板数据模型 - 连玉君,时长:1小时40分钟,课程主页 Stata 33 讲 - 连玉君, 每讲 15 分钟. Stata 小白的取经之路 - 龙志能,时长:2 小时,课程主页 部分直播课 课程资料下载 (PPT,dofiles等)
支持回看
专题 | 嘉宾 | 直播/回看视频 |
---|---|---|
⭐ 最新专题 | 因果推断, 空间计量,寒暑假班等 | |
⭕ 数据清洗系列 | 游万海 | 直播, 88 元,已上线 |
研究设计 | 连玉君 | 我的特斯拉-实证研究设计,-幻灯片- |
面板模型 | 连玉君 | 动态面板模型,-幻灯片- |
面板模型 | 连玉君 | 直击面板数据模型 [免费公开课,2小时] |
Note: 部分课程的资料,PPT 等可以前往 连享会-直播课 主页查看,下载。
关于我们
课程, 直播, 视频, 客服, 模型设定, 研究设计, stata, plus, 绘图, 编程, 面板, 论文重现, 可视化, RDD, DID, PSM, 合成控制法
等
连享会小程序:扫一扫,看推文,看视频……
扫码加入连享会微信群,提问交流更方便
✏ 连享会学习群-常见问题解答汇总:
✨ https://gitee.com/arlionn/WD
New!
lianxh
命令发布了:
随时搜索连享会推文、Stata 资源,安装命令如下:
. ssc install lianxh
使用详情参见帮助文件 (有惊喜):
. help lianxh