diff --git a/GPflow/model.py b/GPflow/model.py
index <HASH>..<HASH> 100755
--- a/GPflow/model.py
+++ b/GPflow/model.py
@@ -167,14 +167,15 @@ class Model(Parameterized):
     def compute_log_likelihood(self):
         return self.build_likelihood()
 
-    def sample(self, num_samples, Lmax=20, epsilon=0.01, verbose=False, return_logprobs=False, RNG=np.random.RandomState(0)):
+    def sample(self, num_samples, Lmax=20, epsilon=0.01, thin=1, burn=0,
+               verbose=False, return_logprobs=False, RNG=np.random.RandomState(0)):
         """
         Use Hamiltonian Monte Carlo to draw samples from the model posterior.
         """
         if self._needs_recompile:
             self._compile()
         return hmc.sample_HMC(self._objective, num_samples,
-                              Lmax, epsilon,
+                              Lmax=Lmax, epsilon=epsilon, thin=thin, burn=burn,
                               x0=self.get_free_state(), verbose=verbose,
                               return_logprobs=return_logprobs, RNG=RNG)
adding extra args to model.sample to pass through to hmc