空间杜宾模型检验

zxm045

16189

收藏 2017-03-24

求空间杜宾模型lr检验和wald检验的stata代码急求

扫码加我拉你入群

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

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

全部回复

飞飞法师

2017-3-24 12:24:47

zxm045 发表于 2017-3-24 10:28
求空间杜宾模型lr检验和wald检验的stata代码急求

matlab可以

扫码加我拉你入群

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

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

zxm045

2017-3-24 18:54:05

飞飞法师发表于 2017-3-24 12:24
matlab可以

只会stata

扫码加我拉你入群

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

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

PX0706

2017-4-2 16:23:46

很多命令都有检验

扫码加我拉你入群

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

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

tiamoyy

2018-9-5 10:06:44

同求！！！

扫码加我拉你入群

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

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

徐丽0

2018-9-8 19:51:56

zxm045 发表于 2017-3-24 18:54
只会stata

可以给一下stata的代码吗、

扫码加我拉你入群

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

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

点击查看更多内容…

啊啊啊啊哦哦哦哦

2018-11-25 14:19:05

飞飞法师发表于 2017-3-24 12:24
matlab可以

您好，能分享一下MATLAB代码吗？

扫码加我拉你入群

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

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

KrHt

2019-2-26 00:50:48

stata我做出来了，还需要吗

扫码加我拉你入群

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

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

杨大炮

2019-2-26 02:09:08

KrHt 发表于 2019-2-26 00:50
stata我做出来了，还需要吗

可以发一下代码吗？

扫码加我拉你入群

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

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

KrHt

2019-2-26 09:22:54

杨大炮发表于 2019-2-26 02:09
可以发一下代码吗？

见我发的帖子吧，里边有

扫码加我拉你入群

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

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

minyunsun

2019-3-7 11:22:37

KrHt 发表于 2019-2-26 00:50
stata我做出来了，还需要吗

大神，你有stata的代码吗？

扫码加我拉你入群

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

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

KrHt

2019-3-10 10:49:25

minyunsun 发表于 2019-3-7 11:22
大神，你有stata的代码吗？

见我的帖子

扫码加我拉你入群

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

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

beiluo08

2019-4-29 17:23:05

xsmle y x1 x2 x3, wmat(W) model(sdm) fe type(both) robust hausman nolog
***Wald Test for SAR
test [Wx]x1=[Wx]x2=[Wx]x3=0
***Wald Test for SEM
testnl ([Wx]x1=-[Spatial]rho*[Main]x1) ([Wx]x2=-[Spatial]rho*[Main]x2) ([Wx]x3=-[Spatial]rho*[Main]x3)

***LR Test for SAR or SEM
qui xsmle y x1 x2 x3, wmat(W) model(sdm) fe type(both) nolog
est sto sdm
qui xsmle y x1 x2 x3, wmat(W) model(sar) fe type(both) nolog
est sto sar
qui xsmle y x1 x2 x3, emat(W) model(sem) fe type(both) nolog
est sto sem
lrtest sdm sar
lrtest sdm sem

扫码加我拉你入群

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

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

jiguishen1

2019-5-4 20:21:45

beiluo08 发表于 2019-4-29 17:23
xsmle y x1 x2 x3, wmat(W) model(sdm) fe type(both) robust hausman nolog
***Wald Test for SAR
test ...

请问为什么都选择both呢？时间固定个体固定还是双固定直接estat ic看aic值和bic值可以吗，aic和bic即使是负数也是越小越好吗

扫码加我拉你入群

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

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

beiluo08

2019-5-5 10:15:58

jiguishen1 发表于 2019-5-4 20:21
请问为什么都选择both呢？时间固定个体固定还是双固定直接estat ic看aic值和bic值可以吗，aic和bic即使是 ...

可以选择ind或者time。both、ind、time三者的比较可以结合AIC、BC进行

扫码加我拉你入群

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

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

jiguishen1

2019-5-5 11:19:09

beiluo08 发表于 2019-5-5 10:15
可以选择ind或者time。both、ind、time三者的比较可以结合AIC、BC进行

谢谢，还想继续请教您几个问题。
1.我通过 lr检验这三个效应，进行个体和双固定的检验会失败，提示 lrtest sdm_sp sdm_both
df(unrestricted) = df(restricted) = 4，但是检验个体和时间，时间和双固定都没问题，而且为什么两两检验lrtet sdm_sp sdp_time和lrtest sdm_time sdm_sp结果是一样的呢？
2.AIC和BIC值即使是负数，也是越小越好吗，我做个体效应estat ic,AIC是-260，BIC是-240
3.如果没有合理解释只要不是比率变量都得取对数吗，可不可以部分取部分不取
问题有点多不好意思，希望您有时间可以回答一下，如果您需要论坛币也可以发帖发个文件定个二三百，我去下载，还不太会用这论坛

扫码加我拉你入群

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

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

杨大炮

2019-5-19 13:03:02

beiluo08 发表于 2019-4-29 17:23
xsmle y x1 x2 x3, wmat(W) model(sdm) fe type(both) robust hausman nolog
***Wald Test for SAR
test ...

请问一下 WX 代表什么意思，还有就是显示我的变量发现不了什么意思

扫码加我拉你入群

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

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

1241590250

2019-6-21 09:35:46

杨大炮发表于 2019-5-19 13:03
请问一下 WX 代表什么意思，还有就是显示我的变量发现不了什么意思

同问

扫码加我拉你入群

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

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

niceSimon7

2019-9-27 12:13:08

jiguishen1 发表于 2019-5-5 11:19
谢谢，还想继续请教您几个问题。
1.我通过 lr检验这三个效应，进行个体和双固定的检验会失败，提示 lrt ...

你好，想问一下你第一个问题解决了吗

扫码加我拉你入群

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

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

Sensual

2020-2-24 10:49:05

KrHt 发表于 2019-2-26 00:50
stata我做出来了，还需要吗

需要想要一下stat做LR检验的代码感谢感谢

扫码加我拉你入群

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

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

Sensual

2020-2-24 10:58:54

Sensual 发表于 2020-2-24 10:49
需要想要一下stat做LR检验的代码感谢感谢

不用啦看到了你发的帖子啦感谢！

扫码加我拉你入群

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

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

2020000001008

2021-12-7 08:47:24

***LR 检验       //用来检验SDM模型能否退化为SEM、SAR模型

xsmle y x1 x2 x3,fe model(sdm)wmat(w)type(both)nolog effects
est store sdm
xsmle y x1 x2 x3,fe model(sar)wmat(w)type(both)nolog effects
est store sar
xsmle y x1 x2 x3,fe model(sem)emat(w)type(both)nolog effects
est store sem
lrtest sdm sar //H0：空间杜宾模型可以简化为空间滞后模型
lrtest sdm sem  //H0：空间杜宾模型可以简化为空间误差模型

***WALD检验 //也是用来检验模型的适配性
clear all
cd"C:\Program Files (x86)\Stata14\ado"//定义路径
use shuju                            //调用数据
spatwmat using D.dta,name(w) standardize //调用数据矩阵D，并将其还行标准化
encode province,gen(pro)
xtset pro year
gen lny=ln(y)
gen lnx1=ln(x1)
gen lnx2=ln(x2)
gen lnx3=ln(x3)
xsmle lny lnx1 lnx2 lnx3,fe model(sdm) wmat(w)type(both)nolog noeffects

test [Wx]x1 = [Wx]x2 = [Wx]x3=0
testnl ([Wx]x1 = -[Spatial]rho*[Main]x1)([Wx]x2 = -[Spatial]rho*[Main]x2)([Wx]x3= -[Spatial]rho*[Main]x3)

扫码加我拉你入群

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

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

2020000001008

2021-12-7 09:48:00

Wald 检验
clear all
use data
spatwmat using W.dta,name(W) standardize
xtset id year
xsmle y x* a*, fe model(sdm) wmat(W) type(both) nolog noeffects
·Wald Test for SAR
test [Wx]x1 = [Wx]x2 = [Wx]x3 = [Wx]a1= [Wx]a2= [Wx]a3= [Wx]a4=0
·Wald Test for SEM
testnl ([Wx]x1 = -[Spatial]rho*[Main]x1)([Wx]x2 = -[Spatial]rho*[Main]x2)([Wx]x3= -
[Spatial]rho*[Main]x3) ([Wx]a1 = -[Spatial]rho*[Main]a1)([Wx]a2= -[Spatial]rho*[Main]a2) ([Wx]a3 =
-[Spatial]rho*[Main]a3)([Wx]a4= -[Spatial]rho*[Main]a4)

LR 检验
xsmle y x* a*, fe model(sdm) wmat(W) type(both) nolog noeffects
est store sdm
xsmle y x* a*, fe model(sar) wmat(W) type(both) nolog noeffects
est store sar
xsmle y x* a*, fe model(sem) emat(W) type(both) nolog noeffects
est store sem
lrtest sdm sar //H0：空间杜宾模型可以简化为空间滞后模型（SAR）
lrtest sdm sem //H0：空间杜宾模型可以简化为空间误差模型（SEM）

扫码加我拉你入群

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

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

东走西顾的兔

2022-1-2 16:38:09

jiguishen1 发表于 2019-5-5 11:19
谢谢，还想继续请教您几个问题。
1.我通过 lr检验这三个效应，进行个体和双固定的检验会失败，提示 lrt ...

你好，请问您的第一个问题解决了吗，代码如何实现呢

扫码加我拉你入群

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

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

wxzh1997

2022-4-22 10:12:13

beiluo08 发表于 2019-4-29 17:23
xsmle y x1 x2 x3, wmat(W) model(sdm) fe type(both) robust hausman nolog
***Wald Test for SAR
test ...

请问里面的Wx是什么呀？W是矩阵对吗？X是什么？

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群