PLotting function now in R script. Better to use grid.arrange rather than messing about with facetting. Facetting made me output 2 extra columns from Marginals class but these may not be necessary in future
Added code to run Lydia's example