This function generates initial values for NPGL distribution using the method of moments estimators given in equations (19) and (20) of Altun (2021).
estim_mu_sigma_NPGL_MM(y)y <- rNPGL(n=100, mu=0.1, sigma=2)
estim_mu_sigma_NPGL_MM(y=y)
#> mu_hat sigma_hat
#> -0.3163345 -3.7711112