Simple R code to generate plots in the paper "Systematic structural discrepancy assessment for computer models". Code is for demonstration purposes and has not been optimised in any way.

1. Create three directories: "Observed Data", "plots_health" and "runs" in the working directory. Alternatively comment out all the "save" commands.

2. Evaluate the file "Stochastic_LotVol_Real_World.R" first, which creates observed data from a more complex "real world" model which uses a stochastic version of the LV model.

3. Then evaluate the file "Deter_LotVol_Assess_Mod_Disc_HealthBio.R" which performs the analysis in the paper. It may take a few mins to run. 

The plots in the paper will be stored in the "plots_health" directory.