#' @keywords internal zgamma <- function (x, y, gamma_power) if (x^2 < y) out <- abs(x)^gamma_power else if (gamma_power == 1) out <- 1.094 * sqrt(y)if (gamma_power == 2) out <- 1.207 * yif (gamma_power == 4/3) out <- 1.129 * y^(2/3)return(out)
Useful links