R/loglik.M14.R

In HardyWeinberg: Statistical Tests and Graphics for Hardy-Weinberg Equilibrium

loglik.M14 <-
function (pa, z, f) 
{
  paa <- pa^2 + pa * (1 - pa) * f
  pab <- 2 * pa * (1 - pa) * (1 - f)
  pbb <- (1 - pa)^2 + pa * (1 - pa) * f
  pvec <- c(paa, pab, pbb)
  ind <- !(z == 0)
  logvec <- log(pvec[ind])
  loglik0 <- sum(z[ind] * logvec)
  nparam0 <- 2
  res <- c(loglik0, nparam0)
  return(res)
}