Computes the excess kurtosis \(\gamma_2 = E[(X - \mu)^4] / \sigma^4 - 3\) for \(X \sim \mathrm{MHN}(\alpha, \beta, \gamma)\).
Details
Integrates the fourth central moment against the unnormalised kernel, in the variable centred on its peak. Building it from raw moments by the Lemma 2b recurrence is exact algebra but cancels once the standard deviation is small next to the mean, which is what a large tilt produces; that expansion is kept only for the parameter values where the integration window cannot be established.
References
Sun, J., Kong, M., & Pal, S. (2023). The Modified-Half-Normal distribution: Properties and an efficient sampling scheme. Communications in Statistics - Theory and Methods, 52(5), 1591–1613. (Lemma 2b)