最近学习了一下用matlab做时变系数var方法,因为时间比较紧没有详细读懂里面的每一个算法,可是发现同样的数据同样的代码,运行两遍,出来的结果都不相同,应该是MCMC上的问题。如何能获取可复制的结果?应该是要设定一个种子值,但我之前没接触过matlab这款软件,有比较了解的请麻烦帮忙看看应该怎么弄。
代码如下(已经调入数据矩阵my)
asvar = {'RPI'; 'R'; 'CC'}; % variable names
nlag = 2; % lags
setvar('data', my, asvar, nlag); % set data
setvar('ranseed', 1);
setvar('fastimp', 1); % fast computing of response
mcmc(10000); % MCMC
drawimp([1], 1);
drawimp([25 41 58], 0); % draw impulse response


雷达卡






京公网安备 11010802022788号







