该程序根据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);