求助：matlab的空间杜宾模型的wald test检验的代码和程序

求助：matlab的空间杜宾模型的wald test检验的代码和程序

网上可以下载得到吧

楼主解决没？？是哪一个？？我也没搞懂

wald test spatial lag
% Wald test for spatial Durbin model against spatial lag model
btemp=results.parm;
varcov=results.cov;
Rafg=zeros(K,2*K+2);
for k=1:K
    Rafg(k,K+k)=1; % R(1,3)=0 and R(2,4)=0;
end
Wald_spatial_lag=(Rafg*btemp)'*inv(Rafg*varcov*Rafg')*Rafg*btemp
prob_spatial_lag=1-chis_cdf (Wald_spatial_lag, K)

wald test spatial error
% Wald test spatial Durbin model against spatial error model
R=zeros(K,1);
for k=1:K
    R(k)=btemp(2*K+1)*btemp(k)+btemp(K+k); % k changed in 1, 7/12/2010
%   R(1)=btemp(5)*btemp(1)+btemp(3);
%   R(2)=btemp(5)*btemp(2)+btemp(4);
end
Rafg=zeros(K,2*K+2);
for k=1:K
    Rafg(k,k)    =btemp(2*K+1); % k changed in 1, 7/12/2010
    Rafg(k,K+k)  =1;
    Rafg(k,2*K+1)=btemp(k);
%   Rafg(1,1)=btemp(5);Rafg(1,3)=1;Rafg(1,5)=btemp(1);
%   Rafg(2,2)=btemp(5);Rafg(2,4)=1;Rafg(2,5)=btemp(2);
end   
Wald_spatial_error=R'*inv(Rafg*varcov*Rafg')*R
prob_spatial_error=1-chis_cdf (Wald_spatial_error,K)

压电振动 发表于 2018-5-15 16:13
wald test spatial lag
% Wald test for spatial Durbin model against spatial lag model
btemp=results ...

为什么我做的那个results.parm那个变量找不到啊，求解

先锋哈哈哈 发表于 2018-5-17 13:39
为什么我做的那个results.parm那个变量找不到啊，求解

设置路径时把elhorst文件放在jplv7文件夹之前

压电振动 发表于 2018-5-17 17:31
设置路径时把elhorst文件放在jplv7文件夹之前

怎么设置，找不到elhorst？

642253181an 发表于 2018-12-9 15:11
怎么设置，找不到elhorst？

软件设置路径，你要找到你的安装目录，你去经管之间首页搜索空间计量模型的操作步骤

https://bbs.pinggu.org/forum.php?mod=viewthread&tid=6886401&page=1#pid56308714

可以试试这个，我自己亲测有效，没有出现parm问题。有需要的可以点连接，里面有我qq联系方式，尽力解决大家的问题

% Wald test spatial Durbin model against spatial error model
R=zeros(K,1);
for k=1:K
    R(k)=btemp(2*K+1)*btemp(k)+btemp(K+k); % k changed in 1, 7/12/2010
%   R(1)=btemp(5)*btemp(1)+btemp(3);
%   R(2)=btemp(5)*btemp(2)+btemp(4);
end
Rafg=zeros(K,2*K+2);
for k=1:K
    Rafg(k,k)    =btemp(2*K+1); % k changed in 1, 7/12/2010
    Rafg(k,K+k)  =1;
    Rafg(k,2*K+1)=btemp(k);
%   Rafg(1,1)=btemp(5);Rafg(1,3)=1;Rafg(1,5)=btemp(1);
%   Rafg(2,2)=btemp(5);Rafg(2,4)=1;Rafg(2,5)=btemp(2);
end   
Wald_spatial_error=R'*inv(Rafg*varcov*Rafg')*R
prob_spatial_error=1-chis_cdf (Wald_spatial_error,K)