R/Neg.LogLikelihood.VC.C1.R

In RVPedigree: Methods for Family-Based Rare-Variant Genetic Association Tests

Defines functions Neg.LogLikelihood.VC.C1

Neg.LogLikelihood.VC.C1 <- function(h, eps, U, S) # f 4.2 articl
{
    B <- h * (S - 1) + 1
    return(sum(log(B)) + t(eps) %*% U %*% diag(1/B) %*% t(U) %*% eps)
}

Try the RVPedigree package in your browser