设为首页收藏本站

EPS数据狗论坛

搜索
快捷导航
查看: 650|回复: 3

stata进行空间计量分析

  [复制链接]

30

主题

272

金钱

445

积分

入门用户

发表于 2019-7-9 16:06:41 | 显示全部楼层 |阅读模式

第一步,打开是stata14,安装xsmle(本文使用的是面板数据);
  1. ssc install xsmle
复制代码

第二步,打开要分析的文件;
首先,单击file—import—选择导入的文件形式(本文导入的是.xls)
1.png
然后,点击Browse,找到所需要的文件,点击OK;
2.png
3.png
第三步,将变量取对数;
  1. gen lnGDP = log(GDP)
  2. gen lnpersonel = log(personel)
  3. gen lnincome = log(income)
  4. gen lnRdfare = log(Rdfare)
  5. gen lnexport = log(export)
  6. gen ln(location) = log(location)
  7. gen lnscienceExpense = log(scienceExpense)
复制代码

第四步,导入权重矩阵;
首先,将权重矩阵.xls转换为.dta格式,并保存(本文命名为weight.dta)
4.png
然后,在stata中打开(a1-a31是所有的变量的变量名,注意W0是新的名字,可以换成任何的名字)
  1. spmat dta W0  a1-a31
复制代码

接着,储存W0(引号里是存储位置,新的spmat的名字叫W3,可以换成任何名字)
  1. spmat save W0 using "F:\Desktop\空间计量的论文\W3.spmat"
复制代码

最后,打开刚生成的spmat文件(W33也是新的名字,打开的时候命名)
  1. spmat use W33 using "F:\Desktop\空间计量的论文\W3.spmat"
复制代码

第五步,使用xtset设置region和year(格式为xtset region year);
  1. xtset state A
复制代码

第六步,将权重矩阵W33标准化;
  1. spatwmat using weight.dta,n(W1) standardize
复制代码

第七步,使用聚类稳健的标准误估计随机效应的SDM模型;
  1. xsmle lnGDP lnpersonel lnincome lnexport lnRdfare lnlocation lnscienceExpense,wmat(W1) model(sdm)  robust nolog
  2. //个体效应
  3. xsmle lnGDP lnpersonel lnincome lnexport lnRdfare lnlocation lnscienceExpense,wmat(W1) model(sdm)  robust nolog type(ind)
  4. //时间效应
  5. xsmle lnGDP lnpersonel lnincome lnexport lnRdfare lnlocation lnscienceExpense,wmat(W1) model(sdm)  robust nolog type(time)
  6. //双效应
  7. xsmle lnGDP lnpersonel lnincome lnexport lnRdfare lnlocation lnscienceExpense,wmat(W1) model(sdm)  robust nolog type(both)
复制代码

如果使用SAR模型,则输入以下形式:
  1. xsmle lnGDP lnpersonel lnincome lnexport lnRdfare lnlocation lnscienceExpense,wmat(W1) model(sar)  robust nolog
复制代码

如果使用SEM模型,则输入以下形式:
  1. xsmle lnGDP lnpersonel lnincome lnexport lnRdfare lnlocation lnscienceExpense,emat(W1) model(sem) robust nolog
复制代码

第八步,进行固定效应的估计;
固定效应也分为时间固定、个体固定和双向固定,本文仅列举固定效应中的空间杜宾模型。若要加上时间固定、个体固定、双固定,在后面加上type(time),type(ind),type(both)即可。
  1. xsmle lnGDP lnpersonel lnincome lnexport lnRdfare lnlocation lnscienceExpense,wmat(W1)  model(sdm) robust nolog fe
复制代码

第九步,进行固定效应AIC和BIC检验
  1. estat ic
复制代码

最后,进行豪斯曼检验,判断使用随机效应模型还是固定效应模型;
  1. qui xsmle lnGDP lnpersonel lnincome lnexport lnRdfare lnlocation lnscienceExpense,wmat(W1)  model(sdm)   nolog
  2. est sto re
  3. qui xsmle lnGDP lnpersonel lnincome lnexport lnRdfare lnlocation lnscienceExpense,wmat(W1)  model(sdm)   nolog fe
  4. est sto fe
  5. // 豪斯曼检验
  6. hausman fe re
复制代码

如果豪斯曼统计量小于零,则接受随机效应的原假设

7

主题

2861

金钱

5308

积分

高级用户

发表于 2019-7-12 08:44:17 | 显示全部楼层
老大,若附上数据文件,参照参照一下,就更好了
回复 支持 反对

使用道具 举报

2

主题

311

金钱

460

积分

入门用户

发表于 2019-8-4 21:41:49 | 显示全部楼层
建议共享一下数据。谢谢!
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

关闭

站长推荐上一条 /1 下一条

客服中心
关闭
在线时间:
周一~周五
9:00-18:00
QQ群:
542280196
联系电话:
010-85786021-8014
在线咨询
客服中心

Powered by BFIT! X3.2© 2008-2018 BFIT Inc.

快速回复 返回顶部 返回列表
082 840 819 796