动态空间杜宾模型 Matlab程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

动态面板空间杜宾模型.rar (337.47 KB)

2017-9-10 09:36:19 上传

该程序根据J.Paul Elhorst...书改写，比静态空间杜宾模型的好的地方在于除了可以分解直接间接效应外，还可以对长短期效应进行分解。
部分代码：
    parms = chol(varcov)'*randn(size(btemp)) + btemp;
    deltasim = parms(npar-1,1); % coef WY(t)
    betasim = parms(3:npar-2,1);
    tausim = parms(1,1); % Coef Y(t-1)
    etasim = parms(2,1); % Coef WY(t-1)
    simresults(:,sim)=[tausim;etasim;betasim;deltasim];
    SS=(eye(N)-deltasim*W)\eye(N);
    SC=SS*((tausim-1)*eye(N)+(deltasim+etasim)*W);
    simdirc(1,sim)=sum(diag(SC))/N; % average direct effect
    simindc(1,sim)=sum(sum(SC,2)-diag(SC))/N; % average indirect effect
    simtotc(1,sim)=simdirc(1,sim)+simindc(1,sim);

扫码加我 拉你入群

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

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

分享0 收藏33 回帖

关键词：matlab 空间计量

longterm.M文件为程序文件，不会的可以用我里面的数据操作一下，使用的时候将数据保存路径复制上，T为（总时期-1）。
A=xlsread('D:\....\data.xlsx');%导入数据
W1=xlsread('D:\....\matrix.xlsx');%导入矩阵，根据自己文件存储位置修改%
T=13; % number of time periods
N=30; % number of regions

longterm.M文件为程序文件，不会的可以用我里面的数据操作一下，使用的时候将数据保存路径复制上，T为（总时期-1）。
A=xlsread('D:\....\data.xlsx');%导入数据
W1=xlsread('D:\....\matrix.xlsx');%导入矩阵，根据自己文件存储位置修改%
T=13; % number of time periods
N=30; % number of regions

楼主我买了，可以具体指导一下吗

楼主我在用，但是总是出现错误，还请您帮助，谢谢！

楼主，我刚购买了，但无法下载。能把程序发给我吗？37473119@qq.com

楼主你好，我在进行SDEM检验时出现这个错误：
Error using  *
Inner matrix dimensions must agree.

Error in direct_indirect_effects_estimates (line 86)
    parms = chol(results.cov)'*randn(size(parm)) + parm;

Error in model1 (line 34)
direct_indirect_effects_estimates(results,W,spat_model);
请问是什么原因造成的。

楼主，你这么做即违法，也不道德！

dkq1027 发表于 2018-6-9 10:52
楼主，你这么做即违法，也不道德！

我这不是照抄原作的，有自己的修改和创新，何为不道德，还谈什么违法。