Computes the skewness \(\gamma_1 = E[(X - \mu)^3] / \sigma^3\) for \(X \sim \mathrm{MHN}(\alpha, \beta, \gamma)\).
Details
Uses the moment recurrence (Sun et al., 2023, Lemma 2b) to compute raw moments up to third order, then converts to central moments.
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), 1507–1536. (Lemma 2b)