Update modelling.R