R/emma_delta_ML_LL_wo_Z.R

In Eagle: Multiple Locus Association Mapping on a Genome-Wide Scale

emma.delta.ML.LL.wo.Z <- function (logdelta, lambda, etas, xi)
{
    n <- length(xi)
    delta <- exp(logdelta)
    return(0.5 * (n * (log(n/(2 * pi)) - 1 - log(sum((etas *
        etas)/(lambda + delta)))) - sum(log(xi + delta))))
}