This vectorises the prior calculation done at the end of the run. I was finding this was running super-slowly for large numbers of samples and parameters as it was basically a nested for loop.
