假如dynare模块程序如下：

var xxxxx;

varexo  xxxxxxx  ;

parameters xxxxxx;

model(linear);

    y = xxxxx

    ……

end;

steady;

check;

shocks;

    var …..

end;

estimated_params;                                               

  xxxx, normal_pdf, 0.03, 0.01;

…..

    stderr epsilon_r, inv_gamma_pdf, 0.01, inf;

    ….

end;

varobs y  pi;

stoch_simul(irf=40, periods=600) y  pi r;

estimation(datafile=data201208, prefilter=1, first_obs=1, mh_drop=0.45, mode_compute=4, mode_check, mh_replic=5000);

上述模型，经运行，可以出估计结果，如脉冲响应函数，后验分布等，但就是没有包括方差分解的结果。

请问高人：
如何在上面的程序中做出方差分解结果？比如：什么命令？放在什么地方？
亲，要具体点哦~我理解能力比较差........