We present two different methods to estimate parameters within a partial differential equation (PDE) model of cancer invasion. The model describes the spatio-temporal evolution of three variables -- tumour cell density, extracellular matrix density and matrix degrading enzyme concentration -- in a one-dimensional tissue domain. The first method is a likelihood-free approach associated with Approximate Bayesian Computation (ABC); the second is a two-stage gradient matching method based on smoothing the data with a Generalized Additive Model (GAM) and matching gradients from the GAM to those from the model. Both methods performed well on simulated data.  To increase realism, additionally we tested the gradient matching scheme with simulated measurement error and found that the ability to estimate some model parameters deteriorated rapidly as measurement error increased.

All related data presented in this repository regarding the manuscript Xiao et al. were simulated using R and MATLAB. The detailed methods of applying the two schemes can be found in the README files in the repository and the comments in each MATLAB or R file. 

The primary repository can be found directly on Github: https://github.com/ycx12341/Data-Code-Figures-ver-4. Updates to the datasets will first be made to this primary repository.