diff --git a/matlab/particle/auxiliary_particle_filter.m b/matlab/particle/auxiliary_particle_filter.m
index 47ff28135013259773231617626a1b1c08a6f6d8..eda6ff52709473295467c4d27e11cbb5c1da3037 100644
--- a/matlab/particle/auxiliary_particle_filter.m
+++ b/matlab/particle/auxiliary_particle_filter.m
@@ -104,7 +104,7 @@ for t=1:sample_size
     sum_tau_tilde = sum(tau_tilde) ;
     lik(t) = log(sum_tau_tilde) ;
     tau_tilde = tau_tilde/sum_tau_tilde;
-    indx_resmpl = resample(tau_tilde);
+    indx_resmpl = resample(tau_tilde,DynareOptions.particle.resampling.method1,DynareOptions.particle.resampling.method2);
     yhat = yhat(:,indx_resmpl);
     wtilde = wtilde(indx_resmpl);
     epsilon = Q_lower_triangular_cholesky*randn(number_of_structural_innovations,number_of_particles);