diff --git a/matlab/evaluate_likelihood.m b/matlab/evaluate_likelihood.m
index cfb702701633569a17f2d9b5de27ea3fab41024b..7d5efbb1c6fa68539f64e029db8721d00ec3953f 100644
--- a/matlab/evaluate_likelihood.m
+++ b/matlab/evaluate_likelihood.m
@@ -102,7 +102,7 @@ if isempty(load_data)
     end
     data = transpose(rawdata);
     % Handle the missing observations.
-    [data_index,number_of_observations,no_more_missing_observations] = describe_missing_data(data,gend,n_varobs);
+    [data_index,number_of_observations,no_more_missing_observations] = describe_missing_data(data);
     missing_value = ~(number_of_observations == gend*n_varobs);
     % Determine if a constant is needed.
     if options_.steadystate_flag% if the *_steadystate.m file is provided.
diff --git a/matlab/imcforecast.m b/matlab/imcforecast.m
index 726b09815c43807b72ddc8bb028e4a1d327bead4..7fa6340c3bb49a13805be109ad6787a3efebee41 100644
--- a/matlab/imcforecast.m
+++ b/matlab/imcforecast.m
@@ -132,7 +132,7 @@ if estimated_model
     end
     data = transpose(rawdata);
     % Handle the missing observations.
-    [data_index,number_of_observations,no_more_missing_observations] = describe_missing_data(data,gend,n_varobs);
+    [data_index,number_of_observations,no_more_missing_observations] = describe_missing_data(data);
     missing_value = ~(number_of_observations == gend*n_varobs);
 
     [atT,innov,measurement_error,filtered_state_vector,ys,trend_coeff] = DsgeSmoother(xparam,gend,data,data_index,number_of_observations);