Computes the logit (log-odds) function \(\mathrm{logit}(p) = \log\left(p /
(1-p)\right)\), the canonical link function for binomial/logistic-family
models throughout this package. p is first clamped to
\([\code{zero\_one\_logit\_clamp}, 1 - \code{zero\_one\_logit\_clamp}]\)
before transforming, so exact 0 or 1 inputs (which would
otherwise map to \(-\infty\)/\(\infty\)) instead return a large but
finite value; this is what lets proportion/fractional responses with mass
exactly at the boundary be used as pseudo-continuous inputs to logit-scale
machinery elsewhere in the package (e.g. fast_ols_cpp-backed
logit-transform-then-OLS shortcuts) without producing non-finite values.
See also
inv_logit for the inverse transform.
