R/mra.default.R

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

mra.default <- function(oformula, odata, eformula, edata){

  dat <- init(oformula, odata,
              eformula, edata,
              scale = TRUE)

  tsr <- TSR(dat)
  wald <- Wald.test(dat)

  i <- dat$map$bet
  lm <- LM.test(dat, wald$par[i], sqrt(wald$vcov[i, i]))

  ct <- c.test(dat, wald$par, wald$vcov)

  #pt1 <- pleio.test1(dat, wald$par)
  pt <- pleio.test(dat, wald$par, wald$vcov)
  
  return(pt)
  
  fit <- reorganize(dat, wald, lm, ct, pt, tsr)
  fit$call <- match.call()
  class(fit) <- 'mra'
  
  fit

}