From b219a99a16d24378fafe7b22988b185655d1c613 Mon Sep 17 00:00:00 2001
From: Houtan Bastani <houtan@dynare.org>
Date: Thu, 24 Jan 2019 11:20:54 +0100
Subject: [PATCH] olsgibbs: fix bug in returned fitted value

---
 matlab/olsgibbs.m | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/matlab/olsgibbs.m b/matlab/olsgibbs.m
index 65a78ac13..fb51a7ee2 100644
--- a/matlab/olsgibbs.m
+++ b/matlab/olsgibbs.m
@@ -117,7 +117,7 @@ else
 end
 
 %% Parse equation
-[Y, ~, X, fp, lp] = common_parsing(ds, get_ast({eqtag}), true);
+[Y, lhssub, X, fp, lp] = common_parsing(ds, get_ast({eqtag}), true);
 lhsname = Y{1}.name;
 Y = Y{1}.data;
 X = X{1};
@@ -192,7 +192,7 @@ if ~isempty(fitted_names_dict)
         yhatname = fitted_names_dict{idx, 2};
     end
 end
-oo_.olsgibbs.(eqtag).Yhat = dseries(X*oo_.olsgibbs.(eqtag).posterior.mean.beta, fp, yhatname);
+oo_.olsgibbs.(eqtag).Yhat = dseries(X*oo_.olsgibbs.(eqtag).posterior.mean.beta, fp, yhatname) + lhssub;
 
 % Apply correcting function for Yhat if it was passed
 if any(idx) ...