N=11; % number of regions
% row-normalize W
W=normw(w11); % function of LeSage
y=shuju1(:,[1]); % column number in the data matrix that corresponds to the dependent variable
x=shuju1(:,[2,3,4,5,6]); % column numbers in the data matrix that correspond to the independent variables
xconstant=ones(N*T,1);
[nobs K]=size(x);
% ----------------------------------------------------------------------------------------
% spatial fixed effects + spatially lagged dependent variable无固定效应
info.lflag=0; % required for exact results
info.model=0;
results1=sar_panel_FE(y,[xconstant x],W,T,info);
vnames1=strvcat('y','intercept','zyfg','lq','rjgdp','gyh','tr');
prt_spnew(results1,vnames1,1);
%空间固定效应
info.lflag=0; % required for exact results
info.model=1;
info.fe=0; % no print intercept and spatial fixed effects
results2=sar_panel_FE(y,x,W,T,info);
vnames2=strvcat('y','zyfg','lq','rjgdp','gyh','tr');
prt_spnew(results2,vnames2,1);
T=7; % number of time periods
N=11; % number of regions
% row-normalize W
W=normw(w11); % function of LeSage
y=shuju1(:,[1]); % column number in the data matrix that corresponds to the dependent variable
x=shuju1(:,[2,3,4,5,6]); % column numbers in the data matrix that correspond to the independent variables
xconstant=ones(N*T,1);
[nobs K]=size(x);
% ----------------------------------------------------------------------------------------
% spatial fixed effects + spatially lagged dependent variable无固定效应
info.lflag=0; % required for exact results
info.model=0;
results1=sar_panel_FE(y,[xconstant x],W,T,info);
vnames1=strvcat('y','intercept','zyfg','lq','rjgdp','gyh','tr');
prt_spnew(results1,vnames1,1);
%空间固定效应
info.lflag=0; % required for exact results
info.model=1;
info.fe=0; % no print intercept and spatial fixed effects
results2=sar_panel_FE(y,x,W,T,info);
vnames2=strvcat('y','zyfg','lq','rjgdp','gyh','tr');
prt_spnew(results2,vnames2,1);
哪位大侠帮我看看这段程序对吗?运行时候总是出错: