莫兰指数stata命令_一文读懂16个sp系列stata空间计量命令及应用
原标题:一文读懂16个sp系列stata空间计量命令及应用
在Stata 15中,推出了最新的空间计量官方命令,均以sp开头,表示 spatial data),可以处理横截面与面板形式的空间数据。本文主要为大家介绍空间计量命令。
一、命令汇总
1、Preparing data/数据准备
zipfile:Compress and uncompress files in zip archive format,表示压缩和解压缩zip存档格式的文档
spshape2dta:Translate shapefile to Stata format表示可以将 shapefile 文件转换为 Stata 格式
spset:Declare data to be Sp spatial data声明空间数据
spbalance:Make panel data strongly balanced使面板数据具有很强的平衡性
spcompress:Compress Stata-format shapefile压缩数据成stata格式的shapefile
2、Looking at data/观测数据
grmap:Graph choropleth maps绘制等值区域图
spdistance:Calculator for distance between places计算区域之间距离
3、Setting the spatial weighting matrix/设置空间权重矩阵
spmatrix:Create, manipulate, and import/export weighting matrices表示创建操作管理、导入和导出空间权重矩阵
spgenerate:Generate spatial lag variables生成空间滞后变量
4、Fitting models/模型估计
spregress:Fit cross-sectional SAR models估计截面SAR模型
spivregress:Fit cross-sectional SAR models with endogenous covariates 估计工具变量SAR模型
spxtregress:Fit panel-data SAR models估计面板数据SAR模型
5、Postestimation/估计检验命令
estat moran:Moran's test after regress莫兰指数检验
spregress postestimation:Postestimation tools for spregress截面SAR检验
spivregress postestimation:Postestimation tools for spivregress工具变量SAR检验
spxtregress postestimation:Postestimation tools for spxtregress面板SAR检验
2、操作应用介绍
首先导入并且下载homicide1990.dta和homicide1990_shp.dta ,其中homicide1990.dta表示的是具有相关变量的数据,而homicide1990_shp.dta是可以构建相关权重矩阵的数据,数据如下:
1、首先导入homicide1990数据,命令为:
use homicide1990
Convert the coordinates to latitude and longitude and measure the distance in miles
spset, modify coordsys(latlong, miles)
spmatrix命令简介
spmatrix可以用来定义、创建、导入、导出空间权重矩阵,基本语法格式为:
spmatrix create contiguity spmatname [if] [in] [, contoptions stdoptions] spmatrix create idistance spmatname [if] [in] [, idistoption stdoptions]
其中spmatrix create contiguity W,此命令将根据相邻关系(contiguity)来创建空间权重矩阵 ,而idistance 表示 inverse distance,将生成基于距离倒数的空间权重矩阵。
2、对于横截面数据,估计截面SAR模型命令为spregress,语法格式为:
Generalized spatial two-stage least squares表示进行 GS2SLS 估计
spregress depvar [indepvars] [if] [in], gs2sls [gs2sls_options]
Maximum likelihood表示进行 MLE 估计
spregress depvar [indepvars] [if] [in], ml [ml_options]
(1)在spregress depvar [indepvars] [if] [in], gs2sls [gs2sls_options]表达式中,选项包括:
dvarlag(spmatname):spatially lagged dependent variable; repeatable表示指定被解释变量(dependent variable)的空间权重矩阵
errorlag(spmatname):spatially lagged errors; repeatable,表示用于指定误差项(error term)的空间权重矩阵。
ivarlag(spmatname : varlist):spatially lagged independent variables; repeatable/表示用于指定自变量(independent variable)的空间权重矩阵。
noconstant:suppress constant term
heteroskedastic:treat errors as heteroskedastic
level(#):set confidence level; default is level(95)表示设置置信水平,默认为95%
另外选择项 robust 表示使用稳健标准误。
(2)在spregress depvar [indepvars] [if] [in], ml [ml_options]
表达式中,包括的选项基本一致,区别之处在于ml方法不同。
3、工具变量法SAR模型估计命令为:spivregress,语法格式为:
spivregress depvar [varlist_1] (varlist_2 = varlist_iv) [if] [in] [, options]
选择项包括ml gs2sls dvarlag(W) errorlag(W) ivarlag(W: x1 x2),含义与spregress一致。
4、面板SAR模型估计命令为:spxtregress,语法格式为:
Fixed-effects maximum likelihood
spxtregress depvar [indepvars] [if] [in], fe [fe_options]
Random-effects maximum likelihood
spxtregress depvar [indepvars] [if] [in], re [re_options]
选项包括如下:
dvarlag(spmatname):spatially lagged dependent variable
errorlag(spmatname):spatially lagged errors
ivarlag(spmatname : varlist):spatially lagged independent variables; repeatable
5、操作应用介绍
copy https://www.stata-press.com/data/r16/homicide_1960_1990.dta .
copy https://www.stata-press.com/data/r16/homicide_1960_1990_shp.dta
use homicide_1960_1990 //导入数据
xtset _ID year
spset
Create a contiguity weighting matrix with the default spectral normalization
spmatrix create contiguity W if year == 1990 //以1990年数据为基准构建空间权重矩阵
Fit a spatial autoregressive random-effects model
spxtregress hrate ln_population ln_pdensity gini i.year, re dvarlag(W)
//空间滞后模型回归
Create an inverse-distance weighting matrix with the default spectral normalization
spmatrix create idistance M if year == 1990
Same as above but use the alternative formulation of the estimator
spxtregress hrate ln_population ln_pdensity gini i.year, re sarpanel dvarlag(M) errorlag(M)
Fit a spatial autoregressive fixed-effects model
spxtregress hrate ln_population ln_pdensity gini i.year, fe dvarlag(M) errorlag(M)
责任编辑:
莫兰指数stata命令_一文读懂16个sp系列stata空间计量命令及应用相关推荐
- stata 求输出相关系数矩阵命令_一文读懂结果输出命令大全(上)
目录 描述统计量 help tabstat //Stata 官方命令 描述统计量组间均值差异检验 help ttest help ttable2 help estout 相关分析命令 help p ...
- stata 将数据集变量名称导出_一文读懂空间计量经济学及stata操作
在Stata 15中,推出了最新的空间计量官方命令,均以sp开头,表示 spatial data),可以处理横截面与面板形式的空间数据.本文主要为大家介绍空间计量命令之spregress的使用. 一. ...
- python空间计量经济学_一文读懂空间计量经济学及stata操作
spcompress:Compress Stata-format shapefile压缩数据成stata格式的shapefile 2.Looking at data/观测数据 2.Looking at ...
- psm倾向得分匹配法举例_一文读懂倾向得分匹配法(PSM)举例及stata实现(一)
原标题:一文读懂倾向得分匹配法(PSM)举例及stata实现(一) 一.倾向匹配得分应用之培训对工资的效应 政策背景:国家支持工作示范项目( National Supported Work,NSW ) ...
- python输入什么就输出什么_一文读懂Python的输入和输出
本文介绍了Python的输入和输出,既然是Python代码,那么就一定有输出量,那么,Python是如何输出的呢? 输出 用print()在括号中加上字符串,就可以向屏幕上输出指定的文字.比如输出'h ...
- java中date类型如何赋值_一文读懂java中的Reference和引用类型
简介 java中有值类型也有引用类型,引用类型一般是针对于java中对象来说的,今天介绍一下java中的引用类型.java为引用类型专门定义了一个类叫做Reference.Reference是跟jav ...
- gps导航原理与应用_一文读懂角速度传感器(陀螺仪)的应用场景
前文我们大致了解陀螺仪的来历,原理和种类,那么,它与我们的日常生活有怎样的关系呢? 陀螺仪器最早是用于航海导航,但随着科学技术的发展,它在航空和航天事业中也得到广泛的应用.陀螺仪器不仅可以作为指示仪表 ...
- hdfs读写流程_一文读懂HDFS分布式存储框架分析
一文读懂HDFS分布式存储框架分析 HDFS是一套基于区块链技术的个人的数据存储系统,利用无处不在的私人PC存储空间及便捷的网络为个人提供数据加密存储服务,将闲置的存储空间利用起来,服务于正处于爆发期 ...
- mysql 默认事务隔离级别_一文读懂MySQL的事务隔离级别及MVCC机制
回顾前文: <一文学会MySQL的explain工具> <一文读懂MySQL的索引结构及查询优化> (同时再次强调,这几篇关于MySQL的探究都是基于5.7版本,相关总结与结论 ...
- 超融合和服务器关系_一文读懂超融合服务器
原标题:一文读懂超融合服务器 1.什么叫超融合服务器 融合基础架构(Hyper-Converged Infrastructure)是一种集成了虚拟计算资源和存储设备的信息基础架构.在这样的架构环境中, ...
最新文章
- 技术图文:如何改进算法的运行效率?
- Eclipse中配置约束
- Hadoop基础-HDFS集群中大数据开发常用的命令总结
- 兰州大学萃英学院计算机,兰州大学萃英学院.PDF
- 查询速度提升200倍,ClickHouse到底有多快?
- Linux中防火墙命令笔记
- python扩展库不是用于科学计算的有_有哪些科学计算的开源库
- 基于matlab的数字下变频器的设计与仿真应用,基于MATLAB的数字下变频器的设计与仿真应用.pdf...
- C++模板-Traits
- MySQL之一篇读懂binlog redolog undolog
- Java中的String.hashCode()方法可能有问题?
- java怎么让遮罩层下面滚动_vue项目弹出层后禁止body底层的滚动事件
- Lottie-iOS的应用及部分源码分析
- 生成器 Generators
- python绘制条形图例题_[宜配屋]听图阁 - Python+matplotlib+numpy绘制精美的条形统计图...
- getch方法_c语言中中getch()函数为什么不起作用?
- Unity 内置渲染管线、SRP、URP、HDRP区别
- .net学科-杨中科-Unity3D视频教程
- bridging the gap between anchor-based and anchor-free detection via adaptive training sample select
- 论文图表录 出现几个 错误标签未定义的简单解决方法