R/Get_Liu_Params_Mod.R

In KMgene: Gene-Based Association Analysis for Complex Traits

Get_Liu_Params_Mod <- function (c1)
{
    muQ <- c1[1]
    sigmaQ <- sqrt(2 * c1[2])
    s1 = c1[3]/c1[2]^(3/2)
    s2 = c1[4]/c1[2]^2
    beta1 <- sqrt(8) * s1
    beta2 <- 12 * s2
    type1 <- 0
    if (s1^2 > s2) {
        a = 1/(s1 - sqrt(s1^2 - s2))
        d = s1 * a^3 - a^2
        l = a^2 - 2 * d
    }
    else {
        type1 <- 1
        l = 1/s2
        a = sqrt(l)
        d = 0
    }
    muX <- l + d
    sigmaX <- sqrt(2) * a
    re <- list(l = l, d = d, muQ = muQ, muX = muX, sigmaQ = sigmaQ,
        sigmaX = sigmaX)
    return(re)
}