diff --git a/matlab/cli/prior.m b/matlab/cli/prior.m
index 703297d73dda4ca91c5aa60c9d961441dd231a50..21a6a8a734308d2c89e56520aaa18ac881045dcc 100644
--- a/matlab/cli/prior.m
+++ b/matlab/cli/prior.m
@@ -161,6 +161,7 @@ if ismember('moments(distribution)', varargin) % Prior simulations (BK).
                 dr = tmp.pdraws{j,3};
                 oo__ = oo_;
                 oo__.dr = dr;
+                Model=set_parameters_locally(Model,tmp.pdraws{j,1});% Needed to update the covariance matrix of the state innovations.
                 oo__ = disp_th_moments(oo__.dr, [], Model, options_, oo__);
                 FirstOrderMoments(:,iter) = oo__.mean;
                 SecondOrderMoments(:,:,iter) = oo__.var;