From 24d41f66c8fe1052fb57470c961a34f2623a8c9f Mon Sep 17 00:00:00 2001
From: Marco Ratto <marco.ratto@jrc.ec.europa.eu>
Date: Fri, 5 Mar 2010 11:26:29 +0100
Subject: [PATCH] fixed bug with tolerance for steady state kalman filter
---
matlab/kalman/likelihood/univariate_diffuse_kalman_filter.m | 2 +-
.../kalman/likelihood/univariate_diffuse_kalman_filter_corr.m | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/matlab/kalman/likelihood/univariate_diffuse_kalman_filter.m b/matlab/kalman/likelihood/univariate_diffuse_kalman_filter.m
index 925b2bbeb7..e6024a22f5 100644
--- a/matlab/kalman/likelihood/univariate_diffuse_kalman_filter.m
+++ b/matlab/kalman/likelihood/univariate_diffuse_kalman_filter.m
@@ -134,7 +134,7 @@ while t < smpl
Zi = Z(i,:);
prediction_error = Y(i,t) - Zi*a;
Fi = Zi*Pstar*Zi'+H(i);
- if Fi > crit
+ if Fi > kalman_tol
Ki = Pstar*Zi';
a = a + Ki*prediction_error/Fi;
Pstar = Pstar - Ki*Ki'/Fi;
diff --git a/matlab/kalman/likelihood/univariate_diffuse_kalman_filter_corr.m b/matlab/kalman/likelihood/univariate_diffuse_kalman_filter_corr.m
index 0370b89ca0..fc588898d5 100644
--- a/matlab/kalman/likelihood/univariate_diffuse_kalman_filter_corr.m
+++ b/matlab/kalman/likelihood/univariate_diffuse_kalman_filter_corr.m
@@ -183,7 +183,7 @@ while t < smpl
Zi = ZZ(i,:);
prediction_error = Y(i,t) - Zi*a;
Fi = Zi*Pstar*Zi'+H(i);
- if Fi > crit
+ if Fi > kalman_tol
Ki = Pstar*Zi';
a = a + Ki*prediction_error/Fi;
Pstar = Pstar - Ki*Ki'/Fi;
--
GitLab