求一份matlab空间面板动态模型命令

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助一份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);

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：MATLAB matla atlab 动态模型 空间面板

你把info.model=3;调成info.model=1;或info.model=2;试试