Stata：最新横截面+工具变量+面板空间计量do文档必知必会

数量经济学专辑汇总！计量百科·资源·干货：

Stata  |Python  |Matlab  |Eviews  |R 

Geoda  |ArcGis  |GeodaSpace  |SPSS 

一文读懂   |数据资源   |回归方法  |网络爬虫                                                               

门限回归   |工具变量   |内生性   |空间计量 

因果推断   |合成控制法   |倾向匹配得分   |断点回归   |双重差分   

面板数据  | 动态面板数据

    

 
          *========================================
          *              高级计量经济学
          *========================================
    
          *            数量经济学
   
          *           
          *       ===================================
          *            空间计量及Stata应用
          *       ===================================

*-------------------------------------------------------------------------------               
*        参考资料：
*        《初级计量经济学及Stata应用：Stata从入门到进阶》             
*        《高级计量经济学及Stata应用：Stata回归分析与应用》
*        《高级计量经济学及Eviews应用》
*        《空间计量入门》 
*        《零基础|轻松搞定空间计量：空间计量及GeoDa、Stata应用》
*        《空间计量第二部：空间计量及Matlab应用课程》
*        《空间计量第三部：空间计量及Stata应用课程》
*        《空间计量第四部：《空间计量及ArcGis应用课程》
*        《空间计量第五部：空间计量经济学》
*        《空间计量第六部：《空间计量及Python应用》
*        《空间计量第七部：《空间计量及R应用》
*        《空间计量第八部：《高级空间计量经济学》
*-------------------------------------------------------------------------------
      

*-------------------------------------------------------------------------------
*一、//////空间计量经济学
*-------------------------------------------------------------------------------
  .
  *0、空间计量      
    set more off
  
  *1、导入数据
    use homicide1990
    
  *2、空间设定
      spset
    
  *3、绘制可视化的等值区域图（ choropleth ）    
    grmap hrate
  
  *4、生成空间权重矩阵
    spmatrix create contiguity W
    /*创建空间权重矩阵*/
    spmatrix create idistance M /*空间反距离矩阵*/
    spmatrix dir
    spmatrix create contiguity W /*空间邻近矩阵*/
    spmatrix dir

    /*输出空间权重矩阵*/
    spmatrix export W using W.txt
    spmatrix export M using M.txt
    
  *5、回归分析regress
      regress hrate
    
  *6、检验空间相关性 
    estat moran, errorlag(W)
    estat moran, errorlag(M)
    
  *7、SAR模型/*空间滞后模型*/
    spregress hrate ln_population ln_pdensity gini, gs2sls dvarlag(W)
    
    spregress hrate ln_population ln_pdensity gini, gs2sls dvarlag(W) errorlag(W)
    estimates store gs2sls_model
    
  *7、SEM模型/*空间误差模型*/
    spregress hrate ln_population ln_pdensity gini, gs2sls errorlag(W)
    spregress hrate ln_population ln_pdensity gini, ml errorlag(W)

  *8、距离空间权重矩阵
  
    spmatrix create idistance M
    spmatrix dir
  
    spregress hrate ln_population ln_pdensity gini, gs2sls dvarlag(W) errorlag(W) ivarlag(W: ln_population ln_pdensity gini) dvarlag(M)  
  
    spregress hrate ln_population ln_pdensity gini, gs2sls dvarlag(W) errorlag(M) ivarlag(W: ln_population ln_pdensity gini)
    estimates store model_ex1_last

  *9、反应空间效应
    estat impact
  
    spregress hrate ln_population ln_pdensity gini, gs2sls heteroskedastic dvarlag(W) errorlag(M) ivarlag(W: ln_population ln_pdensity gini)
  
  
  *10、汇总以及对比
    quietly spregress hrate ln_population ln_pdensity gini, ml dvarlag(W) errorlag(W)
    estimates store model_ml  
    quietly spregress hrate ln_population ln_pdensity gini, gs2sls dvarlag(W) errorlag(W)
    estimates store gs2sls_model
    estimates table gs2sls_model model_ml, b(%6.3f) se(%6.3f)

          *========================================
          *           高级计量经济学
          *========================================
    
      
          *        计量经济学服务中心
      

*-------------------------------------------------------------------------------               
*        参考资料：
*        《初级计量经济学及Stata应用：Stata从入门到进阶》             
*        《高级计量经济学及Stata应用：Stata回归分析与应用》
*        《高级计量经济学及Eviews应用》
*        《空间计量入门》 
*        《零基础|轻松搞定空间计量：空间计量及GeoDa、Stata应用》
*        《空间计量第二部：空间计量及Matlab应用课程》
*        《空间计量第三部：空间计量及Stata应用课程》
*        《空间计量第四部：《空间计量及ArcGis应用课程》
*        《空间计量第五部：空间计量经济学》
*        《空间计量第六部：《空间计量及Python应用》
*        《空间计量第七部：《空间计量及R应用》
*        《空间计量第八部：《高级空间计量经济学》
*-------------------------------------------------------------------------------



*-------------------------------------------------------------------------------
*一、横截面空间计量经济学
*-------------------------------------------------------------------------------


*1、数据查看
use http://www.stata-press.com/data/r15/texas_ue.dta
edit
describe 

**2、描述统计分析
summarize unemployment
grmap unemployment  

*3、OLS回归
regress unemployment college income  

*4、创建空间权重矩阵
spmatrix create contiguity W  //生成相邻权重矩阵W
spmatrix create idistance  M  //生成反距离空间权重矩阵M

*5、空间诊断
regress unemployment college income 
estat moran, errorlag(W) 

regress unemployment college income 
estat moran, errorlag(M) 
spmatrix export W using W.xls 


*6、横截面空间计量分析
*6.1 空间滞后效应（SAR模型）
regress unemployment college income
spregress unemployment college income, gs2sls  dvarlag(W)
spregress unemployment college income, ml    dvarlag(W) 


*6.2 空间误差效应（SEM模型）
regress unemployment college income
spregress unemployment college income, gs2sls  errorlag(W)
spregress unemployment college income, ml   error(W) 

*6.3空间混合回归模型
spregress unemployment college income, gs2sls dvarlag(W) errorlag(W) ivarlag(W:college income ) 
estat impact

*6.4空间杜宾回归模型（SDM）
spregress unemployment college income, gs2sls dvarlag(W)  ivarlag(W:college income ) 

*6.5广义嵌套空间计量模型（GNS）
spregress unemployment college income, gs2sls dvarlag(W) errorlag(W) ivarlag(W:college income ) 


*7、工具变量空间计量分析
copy http://www.stata-press.com/data/r15/dui_southern.dta .
copy http://www.stata-press.com/data/r15/dui_southern_shp.dta .
use dui_southern
spset

spmatrix create contiguity W
spivregress dui nondui vehicles i.dry (police = elect), dvarlag(W) errorlag(W)
spivregress dui nondui vehicles i.dry (police = elect), dvarlag(W) errorlag(W) ivarlag(W: i.dry)


*8、面板空间计量学分析

use homicide_1960_1990
xtset _ID year
spset
*生成空间权重矩阵
spmatrix create contiguity W if year == 1990

*8.1 随机效应
spxtregress hrate ln_population ln_pdensity gini i.year, re dvarlag(W)
spmatrix create idistance M if year == 1990

spxtregress hrate ln_population ln_pdensity gini i.year, re sarpanel dvarlag(M) errorlag(M)

*8.2 固定效应
spxtregress hrate ln_population ln_pdensity gini i.year, fe dvarlag(M) errorlag(M)