Commit c8bc12ab authored by Johannes Pfeifer 's avatar Johannes Pfeifer
Browse files

forcst.m: fix case with measurement error

parent cf37b515
...@@ -85,7 +85,7 @@ for i = 1:horizon ...@@ -85,7 +85,7 @@ for i = 1:horizon
end end
if nargout==3 if nargout==3
var_yf_ME=var_yf; var_yf_ME=var_yf;
[loc_H, loc_varlist] = ismember(options_.varobs', options_.varlist); [loc_H, loc_varlist] = ismember(options_.varobs, var_list);
loc_varlist(loc_varlist==0) = []; loc_varlist(loc_varlist==0) = [];
if ~isempty(loc_varlist) if ~isempty(loc_varlist)
var_yf_ME(:,loc_varlist) = var_yf(:,loc_varlist)+repmat(diag(M_.H(loc_H,loc_H))', horizon, 1); var_yf_ME(:,loc_varlist) = var_yf(:,loc_varlist)+repmat(diag(M_.H(loc_H,loc_H))', horizon, 1);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment