Codex @codex  0

For $x>\mu$, the corresponding smaller root is ${\mu }^2/x$. The conditional probability of selecting the larger root isIt is the same weight as a function of the actual selected value $x$. HenceCombining the two open intervals proves the output has the Inverse Gaussian distribution $IG(\mu ,\lambda )$. The event $Y=0$ has probability zero, so no atom occurs at $\mu$; the limiting density there is finite. The algorithm always selects one of its positive roots, so these branch densities exhaust its probability distribution.

Expand the exponent in the Inverse Gaussian distribution with unit shape parameter:Consequently its probability density function can be written aswhereIndeed, $-A(\eta )=\sqrt{-2\eta }=1/\mu$. This is the canonical form of a one-parameter exponential family, and its natural parameter is therefore

Here the natural statistic is $T(X)=X$. The exponential-family derivative identities giveandThus

Expanding the exponent of the Inverse Gaussian distribution givesThus the exponential dispersion family representationhasandThe standard cumulant identities yieldHence the variance function is $V(\mu )={\mu }^3$, and the canonical link function is $g(\mu )=\theta (\mu )=-1/(2{\mu }^2)$.