R/nderiv.null.R

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

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

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

  nderiv(par, dat)[-1]

}