it’s not exactly log p(theta, x).

log [p(f(phi), x) |J_f| ] = log p(theta, x) + log det(J_f).

The sampling in Stan is done on the unconstrained parameters, phi,
which map through a transformation f to theta. It’s actually really neat.

returned when fit()

deprecated in 2.19 and included in trace