R/nloglik.null.R

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

nloglik.null <- function(par.null, dat, bet.null){

  ## I assume map$bet == 1
  if(dat$map$bet != 1){
    stop('debug nloglik.null')
  }
  par <- c(bet = bet.null, par.null)

  nloglik(par, dat)

}