学习空间计量经济学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值。但是，第二个结果每运行一次，结果都和上一次不一致，而且和教材里不一样。请问大家遇到过这个问题吗？怎么解决的？谢谢。