Tweak the likelihood functions to follow a functional code paradigm

[adadams]

I think to prevent headaches in the future, it may be good at some point to have the likelihood+prior+posterior functions be purely functional, i.e. to not use outside variables.