I'm trying to figure out where the actual box-cox transformation comes from. I've looked at the original paper, and some of it's references, but for the most part, it seems that they just drop the transformation without saying why or how they arrived at this particular transformation. In particular, how do we arrive at the transformed variable $z_i,$ where $\tilde{y}^{\lambda-1}$ is the geometric mean?
$$ z_i = \frac{y_i^{\lambda}-1}{\lambda\tilde{y}^{\lambda-1}} \text{ when }\lambda \neq 0 \\ z_i = \frac{\log(y_i)}{\tilde{y}^{\lambda-1}} \text{ when } \lambda=0 $$
It seems that parts of this transformation come from the normal log-likelihood, but I'm not sure how to move from one step to the other.