The image below displays an approach of using an auxiliary variable to map the parametric curves of a standard normal pdf and cdf.
In Equation (1), z as r.v. is clearly one-dimensional. However, after defining z as an auxiliary variable, it is 2-dimensional
I'm a little bit confused on why they insert the auxiliary variable into the pdf function $f_Z$.
Why is the parametric equation defined as $f_Z(z, f_Z(z))$ and not $(z, f_Z(z))$?
The authors later use this approach to define parametric curves of the pdf after transforming quantiles of the standard normal distribution with the tukey g- and h- transformation, because the pdf does not exist in closed form.
Reference: Headrick, T. C., Kowalchuk, R. K., & Sheng, Y. (2008). Parametric probability densities and distribution functions for Tukey g-and-h transformations and their use for fitting data. Applied Mathematical Sciences, 2(9), 449-462.