学习空间计量经济学Elhost第二章的时候,尝试相应的MATLAB命令。有一个问题弄不明白:
当检验直接效应和简介效应的时候,用了以下代码:
nvar=2;
beta=result.beta(2:3);
rho=result.rho;
for p=1:nvar
C=zeros(N,N);
for i=1:N
for j=1:N
if (i==j) C(i,j)=beta(p);
% else C(i,j)=beta(nvar+p)*W(i,j);
end
end
end
S=(eye(N)-rho*W)\C;
EAVD(p,1)=sum(diag(S))/N; % average direct effect
EAVI(p,1)=sum(sum(S,2)-diag(S))/N; % average indirect effect
EAVC(p,1)=sum(sum(S,1)'-diag(S))/N; % average indirect effect
EAVCtot(p,1)=EAVD(p,1)+EAVI(p,1);
end
fprintf(1,' direct indirect total \n');
[EAVD EAVI EAVCtot]
result.N=result.nobs;
result.cflag=1;
result.xwith=x;
spat_model=0;
direct_indirect_effects_estimates(result,W,spat_model);
然后呢,会出两个结果:一个结果只有变量的直接效益、简介间接效应和总效应,结果和教材结果是一致的,但是没有t值;另一个结果包括了直接效应、简介效应和总效应,还包括各自的t值。但是,第二个结果每运行一次,结果都和上一次不一致,而且和教材里不一样。请问大家遇到过这个问题吗?怎么解决的?谢谢。


雷达卡



京公网安备 11010802022788号







