R/auxlindley.R

Defines functions auxlindley

auxlindley <-
function(beta,alpha,nAA,nBB,N) {
  lnum <- nAA*(alpha + beta) + nBB*(alpha - beta)
  lden <- N*log((1+exp(alpha)*cosh(beta)))
  y <- exp(lnum-lden)
  return(y)
}

Try the HardyWeinberg package in your browser

Any scripts or data that you put into this service are public.

HardyWeinberg documentation built on May 7, 2022, 5:05 p.m.