From 98591af8462ace94879fc9a5a6cd0392e6f3c099 Mon Sep 17 00:00:00 2001
From: Houtan Bastani <houtan@dynare.org>
Date: Tue, 12 Feb 2019 19:57:39 +0100
Subject: [PATCH] sur: update M_.Sigma_e using estimated parameters

---
 matlab/ols/sur.m | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/matlab/ols/sur.m b/matlab/ols/sur.m
index 6dab7fbf9..6589249e3 100644
--- a/matlab/ols/sur.m
+++ b/matlab/ols/sur.m
@@ -137,6 +137,10 @@ kLeye = kron(chol(inv(M_.Sigma_e)), eye(oo_.sur.(model_name).dof));
 [q, r] = qr(kLeye*X.data, 0);
 oo_.sur.(model_name).beta = r\(q'*kLeye*Y.data);
 
+resid = Y.data - X.data * oo_.sur.(model_name).beta;
+resid = reshape(resid, oo_.sur.(model_name).dof, neqs);
+M_.Sigma_e = resid'*resid/oo_.sur.(model_name).dof;
+
 M_.params(opidxs) = oo_.sur.(model_name).beta;
 
 % Yhat