R/e.step.2normal.R

In idr: Irreproducible Discovery Rate

Defines functions e.step.2normal

Documented in e.step.2normal

e.step.2normal <-
function(z.1, z.2, mu, sigma, rho, p){

  e.z <- p/((1-p)*exp(d.binormal(z.1, z.2, 0, 1, 0)-d.binormal(z.1, z.2, mu, sigma, rho))+ p)
  
  invisible(e.z)
}

Try the idr package in your browser