Add automatic coverage determination

[johnlees]

We have a mixture of Poissons model from stan which should work. First just going to try optimisation without gradient