Use unconstrained distributions internally
This version is almost entirely concerned with updating the underlying Stan code. The goal is to yield decreased computation time for examples with large numbers of samples or taxa.
Rather than using constrained distributions (e.g., lognormal), we now model the parameters on the original (e.g., log) scale and then back-transform.