Invert gamma = 1 / (alpha (1 - alpha)) for alpha

Uses the conjugate form to keep precision when gamma is large.

Usage

morie_alpha_from_gamma(gamma)

Arguments

gamma

Threshold, must be at least 4.

Value

A numeric value in (0, 1/2].

Examples

morie_alpha_from_gamma(4)     # gamma = 4 gives alpha = 1/2
#> [1] 0.5
morie_alpha_from_gamma(100)   # alpha (1 - alpha) = 1/100
#> [1] 0.01010205