question: public static double quantile ( double p , double k , double theta , double shift ) { return Math . exp ( GammaDistribution . quantile ( p , k , theta ) ) + shift ; }
target:  Compute probit (inverse cdf) for LogGamma distributions.