求助一份matlab空间面板动态模型命令，本人在论坛上找过，但是发现运行结果和stata的结果完全不一样。而且stata还可以构建被解释变量在时间上、空间上、空间和时间上的滞后项可以建三个模型，matlab只见有一个模型的命令，好奇怪。

以下是命令，但总感觉有点不太对，有没有大神给一份新的呢。


T=9;
N=31;
W=normw(juli);
y=lunwen(:,1);
x=lunwen(:,[2,3,4,5,6,7,8]);
for t=1:T+1
    t1=(t-1)*N+1;t2=t*N;
    wx(t1:t2,:)=W*x(t1:t2,:);
    Wy(t1:t2,1)=W*y(t1:t2,1);
end
[nobs K]=size(x);
xconstant=ones(nobs,1);
%--------------------------------------------------------------------

info.lflag=0;
info.tl=1;
info.stl=1;
info.ted=1;
info.dyn=1;
info.model=3;
info.fe=0;
info.bc=0;
results=sar_panel_FE(y(N+1:end),[y(1:end-N) Wy(1:end-N) x(N+1:end,:) wx(N+1:end,:)],W,T,info);
vnames=char('cxcy','timelag cxcy','spacetimelag cxcy', 'cywh','lncfsp','hlw','lncsfj','lnjyzc','W*cywh','W*lncfsp','W*hlw','W*lncsfj','W*lnjyzc');
results1=sar_jihai_time(y(1:nobs),[x(N+1:nobs,:) wx(N+1:nobs,:)],W,info);
results.beta=results1.theta1(1:end-2);
results.rho=results1.theta1(end-1);
results.tstat=results1.tstat1(1:end-1);
prt_sp(results,vnames,1);