全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
74298 192
2013-10-01
elhorst.rar
大小:(6.29 KB)

只需: 2 个论坛币  马上下载

本附件包括:

  • demopanelscompare.m
  • direct_indirect_effects_estimates.m
  • f2_sarpanel.m
  • f_sarpanel.m
  • sar_panel_FE.m

最近写论文用到空间杜宾模型SDM,由于jplv7包中的原有的SDM模型只能分析截面数据,而做面板数据需要jplv7的panel里面给出的SAR模型再自己改编,难度很高。附件是已经直接可以做面板的SDM代码,当然,也可以做SAR,SER等空间计量模型,非jplv7的原始版本(Elhorst在此基础上的修改),套用这个代码,做空间计量模型(SDM,SAR,SEM),除了可以得到模型系数估计值,还可以计算自变量对因变量的直接效应和间接效应。找的很辛苦,分享给大家。另外,Matlab的包具体如何使用,很难手把手教,最主要的还是靠自己去一句一句读代码,通过对DEMO的运行(这个在Elhorst文件中不含,建议下载Elhorst model文件,不仅有函数,还有实例和原数据,点击DEMO可以看到运行效果),查看它是如何对原始数据运用读取的,然后结合FUNCTION的编写解释,自行摸索,这是我的一点体会。
附件列表

elhorst model.rar

大小:98.3 KB

只需: 2 个论坛币  马上下载

这个是附带实例版本的,比上面那个多一些文件

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-10-1 21:47:31
有什么问题欢迎在此留言,大家一起交流讨论
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-10-2 11:44:13
已经下载第一个文件的朋友,把邮箱发给我,我把补充版本发给你们
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-10-2 11:46:49
附上程序操作的部分结果,看看你是否需要用到
Results obtained in Matlab version 7.10.0 after running demopanelscompare directly after Matlab was opened

Results reported in Matlab paper are indicated

>> demopanelscompare

Pooled model with spatially lagged dependent variable, no fixed effects
Dependent Variable =        logcit   
R-squared          =    0.3504   
corr-squared       =    0.2995   
sigma^2            =    0.0327
Nobs,Nvar,#FE      =   1380,     4,     3  
log-likelihood     =        394.57544
# of iterations    =      1   
min and max rho    =   -1.3924,   1.0000
total time in secs =    1.4840
time for optimiz   =    0.2660
time for lndet     =    0.0470
time for eigs      =    0.7340
time for t-stats   =    0.0930
No lndet approximation used
***************************************************************
Variable        Coefficient  Asymptot t-stat    z-probability
intercept          2.649149        16.690302         0.000000
logp              -0.759025       -19.195318         0.000000
logy               0.253923        10.127432         0.000000
W*dep.var.         0.189563         6.173944         0.000000

    direct    t-stat   indirect    t-stat   total    t-stat

ans =

   -0.7642  -19.9941   -0.1707   -5.7618   -0.9349  -22.1216
    0.2562   10.6002    0.0572    5.1699    0.3134   10.7289

Direct         Coefficient           t-stat           t-prob         lower 05         upper 95
logp             -0.766077       -20.452653         0.000000        -0.841696        -0.695877
logy              0.256334        10.659193         0.000000         0.209184         0.304201

Indirect       Coefficient           t-stat           t-prob         lower 05         upper 95
logp             -0.170935        -5.853636         0.000000        -0.228872        -0.114329
logy              0.057215         5.255482         0.000004         0.036976         0.079186

Total          Coefficient           t-stat           t-prob         lower 05         upper 95
logp             -0.937012       -22.000409         0.000000        -1.030269        -0.856131
logy              0.313549        10.786241         0.000000         0.256581         0.370397


Pooled model with spatially lagged dependent variable, no fixed effects
Dependent Variable =        logcit   
R-squared          =    0.4346   
corr-squared       =    0.3600   
sigma^2            =    0.0285
Nobs,Nvar,#FE      =   1380,     6,     5  
log-likelihood     =         475.5419
# of iterations    =      1   
min and max rho    =   -1.3924,   1.0000
total time in secs =    0.1250
time for optimiz   =    0.0160
time for lndet     =    0.0160
time for eigs      =    0.0460
No lndet approximation used
***************************************************************
Variable        Coefficient  Asymptot t-stat    z-probability
intercept          2.631341        15.820860         0.000000
logp              -1.250674       -21.795938         0.000000
logy               0.554194        14.957778         0.000000
W*logp             0.780485        11.148678         0.000000
W*logy            -0.444388       -10.114568         0.000000
W*dep.var.         0.336576        11.090923         0.000000

    direct    t-stat   indirect    t-stat   total    t-stat

ans =

   -1.2161  -23.3866    0.5075    7.2653   -0.7086  -13.3295
    0.5303   15.4833   -0.3664   -7.4704    0.1639    4.0665

Direct         Coefficient           t-stat           t-prob         lower 05         upper 95
logp             -1.214549       -22.699019         0.000000        -1.317824        -1.110875
logy              0.528213        15.097459         0.000000         0.458332         0.592698

Indirect       Coefficient           t-stat           t-prob         lower 05         upper 95
logp              0.505157         7.364146         0.000000         0.367755         0.635877
logy             -0.362380        -7.485479         0.000000        -0.458527        -0.266987

Total          Coefficient           t-stat           t-prob         lower 05         upper 95
logp             -0.709393       -13.544731         0.000000        -0.809527        -0.613405
logy              0.165833         4.055646         0.000191         0.080061         0.239749


Wald_spatial_lag =

   14.8284


prob_spatial_lag =

  6.0264e-004


LR_spatial_lag =

   15.7526


prob_spatial_lag =

  3.7964e-004


Wald_spatial_error =

    8.9801


prob_spatial_error =

    0.0112


LR_spatial_error =

    8.2292


prob_spatial_error =

    0.0163

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-10-18 14:17:39
楼主你的意思是第二个文件是第一个文件的完整版?只要下第二个文件就好么?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-10-18 14:18:13
hefengmocha 发表于 2013-10-2 11:44
已经下载第一个文件的朋友,把邮箱发给我,我把补充版本发给你们
楼主你的意思是第二个文件是第一个文件的完整版?只要下第二个文件就好么?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群