Consider a more scalable inference algorithm

[jnothman]

There are problems with both memory and time consumption of the implementation. We could consider implementing the variational inference described at https://arxiv.org/pdf/1711.05150.pdf. I've previously found an R implementation of this, but have now lost it.