R/lrtout.R

lrtout <- function(betaout, nsbp, alpha) {
  
  over <- betaout[, 1] / betaout[, 2]
  oversq <- over ^ 2
  cc <- min(sum(oversq * (over > 0)), sum(oversq * (over < 0)))
  P_value <- alphaf(nsbp, cc)
  c_value <- cvaluef(nsbp, alpha)
  Powerv <- powerval(betaout, nsbp, test = "LRT", alpha)
  
  list(index = c_value, pvalue = P_value, power = Powerv)
  
}

Try the QualInt package in your browser

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

QualInt documentation built on May 1, 2019, 11 p.m.