R/cmlmomco.R

In lmomco: L-Moments, Censored L-Moments, Trimmed L-Moments, L-Comoments, and Many Distributions

"cmlmomco" <- function(f,para) {
    if(! are.par.valid(para)) return()
    Qlo <- par2qua(0,para,paracheck=FALSE)
    Mu <- vector(mode="numeric", length=length(f))
    for(i in 1:length(f)) {
       if(f[i] == 1) { Mu[i] <- par2qua(1,para,paracheck=FALSE); next }
       tmp <- NULL
       "afunc" <- function(p) {
          return(par2qua(p,para,paracheck=FALSE))
       }
       try(tmp <- integrate(afunc, f[i], 1))
       Mu[i] <- ifelse(is.null(tmp), NA, tmp$value/(1-f[i]))
    }
    return(Mu)
}