R/are.paremu.valid.R

"are.paremu.valid" <-
function(para, nowarn=FALSE) {
    if(! is.emu(para)) return(FALSE)
    if(any(is.na(para$para))) return(FALSE)

    E <- para$para[1]
    M <- para$para[2]
    op <- options()
    GO <- TRUE
    if(nowarn == TRUE) options(warn=-1)
    if(E < 0) {
      warning("Parameter ETA is not >= 0, invalid")
      GO <- FALSE
    }
    if(E > 1) {
      warning("Parameter ETA is not <= 1, invalid")
      GO <- FALSE
    }
    if(M <= 0) {
      warning("Parameter MU is not > 0, invalid")
      GO <- FALSE
    }
    options(op)
    if(GO) return(TRUE)
    return(FALSE)
}
wasquith/lmomco documentation built on Feb. 24, 2024, 1:31 a.m.