This code has been adapted from the software by Frank Dondelinger uploaded for their publication in ICML 2014: Frank Dondelinger, Maurizio Filippone, Simon Rogers and Dirk Husmeier ODE parameter inference using adaptive gradient matching with Gaussian processes Proceedings of the 16th International Conference on Artificial Intelligence and Statistics, Journal of Machine Learning Research - Workshop and Conference Proceedings Volume 31, pp. 216-228, 2013 http://jmlr.org/proceedings/papers/v31/dondelinger13a.pdf Code adapted by Benn Macdonald, for the simulations carried out in the paper: Benn Macdonald, Mu Niu, Simon Rogers, Maurizio Filippone and Dirk Husmeier "Approximate parameter inference in systems biology using gradient matching: a comparative evaluation" Biomedical Engineering Online http://eprints.gla.ac.uk/113366/ In print The compressed file contains: LB10 Scheme.R is the script that performs parameter estimation on the Fitz-Hugh Nagumo model. Running this script locates and runs all necessary functions e.g. fitting the Gaussian process, performing the MCMC sampling etc. Run this file to see an example of the code working. It is set-up to use log base 10 increments of the gradient mismatch parameter. This can be changed in the file if so desired. Fitz-Hugh Nagumo Model.R specifies the ODEs for the Fitz-Hugh Nagumo model. FhN SD Noise 0.38 Average SNR 10 1.R is an example dataset for the Fitz-Hugh Nagumo model. Dataset Creation Fitz-Hugh Nagumo Model.R allows for simulating datasets for the Fitz-Hugh Nagumo model. V&G Model 1.R specifies the ODEs for the protein signalling transduction pathway. Dataset Creation V&G.R allows for simulating datasets for the protein signalling transduction pathway. The folder MCMC Source contains the files for the MCMC. Namely: fitGP.R fits a Gaussian process functionsMCMC.R specifies functions to allow the MCMC to run init_utilTemperingMismatch.R is a script that initialises things before the MCMC begins population_mcmcTemperingMismatch.R is a script that performs adaptive gradient matching with tempering of the gradient mismatch parameter. It is run from the script LB10 Scheme.R The files contain comments to ease understanding. Read through LB10 Scheme.R before running for the first time.