diff --git a/matlab/backward/backward_model_irf.m b/matlab/backward/backward_model_irf.m
index 7760a2ffd7d42e85eb5cdfd90a32052b32e1d9ee..37952e6138d244d2f10ab6ced7070e6573f0092b 100644
--- a/matlab/backward/backward_model_irf.m
+++ b/matlab/backward/backward_model_irf.m
@@ -220,6 +220,7 @@ for i=1:length(listofshocks)
         name = listofshocks{i};
     end
     deviations.(name) = alldeviations{listofvariables{:}};
+    deviations.(name) = [deviations.(name) dseries(innovations, initialcondition.last+1, exonames)];
     if nargout>2
         irfs.(name) = allirfs{listofvariables{:}};
         irfs.(name) = [irfs.(name) dseries(innovations, initialcondition.last+1, exonames)];