空间Durbin模型(.lagsarlm代码):我无法使用回归值的z值和p值对影响进行汇总

时间:2019-08-08 16:50:30

标签: r spatial-regression

我没有太多时间来运行我的空间回归模型,因此,我将非常感谢您的帮助。

我不知道如何运行空间durbin模型(.lagsarlm)的碰撞函数,这样我才能看到回归变量总效应的p值。

我使用的是Queen连续性矩阵(.gal文件),我在下面针对“空间误差Durbin模型”或“空间滞后X模型”运行了完全相同的代码,并且工作正常。但是对于使用.lagsarlm函数的空间Durbin模型,它不起作用

# The model:
regSDM.5=lagsarlm(reg.eq5,data=spat.data,listw1,type="mixed")

# Impacts, not working:
impacts (regSDM.5, listw=listw1)
summary(impacts(regSDM.5,listw=listw.wts,R=500),zstats=TRUE)

我得到的只是一个错误:

> summary(impacts(regSDM.5,listw=listw1,R=500),zstats=TRUE)



Error in eigen(Sigma, symmetric = TRUE) : 
infinite or missing values in 'x'
In addition: Warning message:
Method impacts.sarlm moved to the spatialreg package

我需要获得这样的结果(此处显示了SLX示例),最重要的是该结果最底部的总体效果的p值。

> summary(impacts(regSLX.5,listw=listw1,R=500),zstats=TRUE)
Impact measures (SLX, estimable, n-k):
                   Direct      Indirect         Total
ROI_Act    -0.79420742945  1.0099291341  0.2157217047
Home_Aff   -0.00731795312 -0.0588033114 -0.0661212645
Emmissions -0.01044269992 -0.0280940042 -0.0385367041
Pop_Dens   -0.00002670612 -0.0001296893 -0.0001563955
lnEl_Cons  -0.67494998732 -0.3364665653 -1.0114165526
lnDisp_Y   -1.28264355334 -0.3288232887 -1.6114668421
Own_Occ    -0.00564662179 -0.0495448745 -0.0551914962
Detached    0.01517019074  0.0123018795  0.0274720702
========================================================
Standard errors:
                 Direct      Indirect         Total
ROI_Act    0.1669871799 0.18392867854 0.08612939462
Home_Aff   0.0206790208 0.04270741415 0.04326666551
Emmissions 0.0043738633 0.00817558337 0.00888052007
Pop_Dens   0.0000358991 0.00005264358 0.00004279778
lnEl_Cons  0.4147250666 0.66547133438 0.59146099150
lnDisp_Y   0.3609449462 0.56077965701 0.54018324532
Own_Occ    0.0064277564 0.01138136501 0.01089799081
Detached   0.0089160665 0.01223168037 0.00945775531
========================================================
Z-values:
               Direct   Indirect     Total
ROI_Act    -4.7560982  5.4908736  2.504623
Home_Aff   -0.3538830 -1.3768877 -1.528226
Emmissions -2.3875232 -3.4363302 -4.339465
Pop_Dens   -0.7439218 -2.4635356 -3.654289
lnEl_Cons  -1.6274637 -0.5056064 -1.710031
lnDisp_Y   -3.5535712 -0.5863681 -2.983186
Own_Occ    -0.8784748 -4.3531575 -5.064374
Detached    1.7014443  1.0057391  2.904714

p-values:
           Direct       Indirect       Total        
ROI_Act    0.0000019737 0.000000039995 0.0122582    
Home_Aff   0.72342660   0.16854698     0.1264563    
Emmissions 0.01696233   0.00058965     0.00001428302
Pop_Dens   0.45692378   0.01375743     0.0002579    
lnEl_Cons  0.10363865   0.61313299     0.0872602    
lnDisp_Y   0.00038004   0.55762819     0.0028527    
Own_Occ    0.37968613   0.000013419067 0.00000040975
Detached   0.08885959   0.31454108     0.0036759  

我发现运行它(“ tr”)需要某些电源曲线,但我不知道它是否适用于皇后型矩阵。

很抱歉,有很多细节,谢谢您提前获得帮助!

0 个答案:

没有答案