R/loglik.2binormal.R

In idr: Irreproducible Discovery Rate

Defines functions loglik.2binormal

Documented in loglik.2binormal

loglik.2binormal <-
function(z.1, z.2, mu, sigma, rho, p){

  l.m <- sum(d.binormal(z.1, z.2, 0, 1, 0)+log(p*exp(d.binormal(z.1, z.2, mu, sigma, rho)-d.binormal(z.1, z.2, 0, 1, 0))+(1-p)))
  
  return(l.m) 
}

Try the idr package in your browser