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)

  fit <- reorganize(dat, wald, lm, ct, tsr)
  fit$call <- match.call()
  class(fit) <- 'mra'



Try the iva package in your browser

Any scripts or data that you put into this service are public.

iva documentation built on May 2, 2019, 3:25 a.m.