[size=13.63636302947998px]在做空间误差模型的检验时,上述两个出错无法输出检验结果,代码如下:
[size=13.63636302947998px]load wmat.dat;
[size=13.63636302947998px]load anselin.dat;
[size=13.63636302947998px]y = anselin(:,1); nobs = length(y);
[size=13.63636302947998px]x = [ones(nobs,1) anselin(:,2:3)];
[size=13.63636302947998px]W = sparse(wmat( :,1),wmat( :,2),wmat( :,3));
[size=13.63636302947998px]vnames = strvcat('crime','const','income','house value');
[size=13.63636302947998px]res1 = moran(y,x,W);
[size=13.63636302947998px]prt(res1);
[size=13.63636302947998px]res2 = lmerror(y,x,W);
[size=13.63636302947998px]prt(res2);
[size=13.63636302947998px]z = speye(49) - 0.1*sparse(W);
[size=13.63636302947998px]p = colmmd(z);
[size=13.63636302947998px]res3 = lratios(y,x,W);
[size=13.63636302947998px]prt(res3);
[size=13.63636302947998px]res4 = walds(y,x,W);
[size=13.63636302947998px]prt(res4);
[size=13.63636302947998px]res5 = lmsar(y,x,W,W);
[size=13.63636302947998px]prt(res5);
[size=13.63636302947998px]res = sem(y,x,W);
[size=13.63636302947998px]prt(res,vnames);
[size=13.63636302947998px]显示:Undefined function 'colmmd' for input arguments of type 'double'.,换成colamd后也还是出错