R/Lemma5.R

In iva: Instrumental Variable Analysis in Case-Control Association Studies

Lemma5 <- function(par, dat){

  map <- dat$map
  n0 <- dat$n0
  n1 <- dat$n1

  Jn <- Lemma4(par, dat)
  Jna <- Jn[, map$a, drop = FALSE]
  In <- Jn - (n0 + n1)/n0/n1 * Jna %*% t(Jna)

  In

}