306878644 发表于 2016-3-11 16:42 
W*X的问题你解决了吗
T=30; % number of time periods
N=46; % number of regions
% row-normalize W
W=normw(W1); % function of LeSage
y=A(:,[3]); % column number in the data matrix that corresponds to the dependent variable
x=A(:,[4,6]); % column numbers in the data matrix that correspond to the independent variables
SDM模型也有程序代码的,链接https://bbs.pinggu.org/thread-2657434-1-1.html
里面解决了这个问题,WX不是W*X,是用以下这段代码处理的
for t=1:T
t1=(t-1)*N+1;t2=t*N;
wx(t1:t2,:)=W*x(t1:t2,:);
end
xconstant=ones(N*T,1);
[nobs K]=size(x);
这只是我个人理解,其实我也看不太懂程序。。。