动态空间杜宾模型的程序实现

看elhorst资料已解决静态空间杜宾模型，程序也能跑出来，但动态空间杜宾模型程序运行出不来，出现??? Wrong # of arguments to prt_sp，不知道程序那个地方出现错误，请各位大神指教！！

连图你都不贴出来，你让人家怎么帮你解决

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);
>> ylevel=y;fprintf(1,'Y minus item 6 \n');
ylevel=100*ylevel;
Y minus item 6
>> fprintf(1,'W obv distance capitals \n');p=eig(W);eigmax=p(2);peig=1;
W obv distance capitals
>> vnames2=strvcat('FL(t)','FL(t-1)','W*FL(t-1)','tc','ftc','tftc','IR*tc','IR*ftc');
>> vnames=vnames2;fprintf(1,'x var 1-7 + 9-12 \n');
x var 1-7 + 9-12
>> for t=1:T+1
    t1=1+(t-1)*N;t2=t*N;
    Wylevel(t1:t2,1)=W*ylevel(t1:t2,1);
end
y=ylevel(N+1:end)-ylevel(1:end-N);
??? Index exceeds matrix dimensions.

>> for t=1:T+1
    t1=1+(t-1)*N;t2=t*N;
    Wylevel(t1:t2,1)=W*ylevel(t1:t2,1);
end
y=ylevel(N+1:end)-ylevel(1:end-N);
就是运行这段时出错的

这样是看不出来怎么出错的，你需要一步一步的步进程序来找出错误。

恩，已经解决了，谢谢！

林静树 发表于 2015-10-27 20:42
恩，已经解决了，谢谢！

请问动态SDM的程序在哪里可以下载？如果可以，能否发一份给我？谢谢
343103419@qq.com

林静树 发表于 2015-10-27 20:42
恩，已经解决了，谢谢！

请问问题是怎么解决的啊，谢谢！

林静树 发表于 2015-10-27 20:42
恩，已经解决了，谢谢！

你好！请问你的动态空间杜宾模型用的什么估计方法？谢谢！

楼主你怎么解决的，我也是??? Index exceeds matrix dimensions.