R/riglmomco.R

"riglmomco" <- function(f,para) {
    if(! are.par.valid(para)) return()
    Qu <- par2qua(f,para,paracheck=FALSE)
    Theta1 <- sapply(f, function(u) { return(rreslife.lmoms(u,para, nmom=1)$lambdas[1]) })
    Gu <- 1 - (Theta1/Qu)
    return(Gu)
}
wasquith/lmomco documentation built on April 20, 2024, 7:20 p.m.