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检验代码。我用jplv7包。在时间,空间,混合估计下结果可以跑出来,但是在无固定效应下显示*内部矩阵维度出错。请问有高手可以解决吗?