Add mh samples option to fsl_mrs.

CHANGELOG.rst

@@ -8,6 +8,7 @@ This document contains the FSL-MRS release history in reverse chronological orde
 - Dynamic fitting now handles multiple different basis sets.
 - Fix mapped parameter uncertainties in dynamic MRS results.
 - Dynamic fitting results can now be saved to and loaded from a directory.
+- Added MH sample option to fsl_mrs, matches interactive python interface.
 
 1.1.8 (Tuesday 5th October 2021)
 -------------------------------

fsl_mrs/scripts/fsl_mrs

@@ -85,6 +85,10 @@ def main():
                                    ' independently of other basis spectra.')
     fitting_args.add_argument('--disable_MH_priors', action="store_true",
                               help="Disable MH priors.")
+    fitting_args.add_argument('--mh_samples', type=int, default=500,
+                              help="Number of Metropolis Hastings samples,"
+                                   " every tenth sample is kept."
+                                   " Default = 500")
 
     # ADDITIONAL OPTIONAL ARGUMENTS
     optional.add_argument('--t1', type=str, default=None, metavar='IMAGE',
@@ -274,14 +278,16 @@ def main():
                    'baseline_order': args.baseline_order,
                    'metab_groups': metab_groups,
                    'model': 'lorentzian',
-                   'disable_mh_priors': args.disable_MH_priors}
+                   'disable_mh_priors': args.disable_MH_priors,
+                   'MHSamples': args.mh_samples}
     else:
         Fitargs = {'ppmlim': ppmlim,
                    'method': args.algo,
                    'baseline_order': args.baseline_order,
                    'metab_groups': metab_groups,
                    'model': 'voigt',
-                   'disable_mh_priors': args.disable_MH_priors}
+                   'disable_mh_priors': args.disable_MH_priors,
+                   'MHSamples': args.mh_samples}
 
     if args.verbose:
         print(mrs)