R/cdflap.R

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

"cdflap" <-
function(x,para) {
    if(! are.parlap.valid(para)) return()
    XI <- para$para[1]
    A  <- para$para[2]

    f <- vector(mode="numeric", length=length(x))
    Y <- (x - XI)/A
    f[x <= XI] <-     0.5*exp( Y[x <= XI])
    f[x >  XI] <- 1 - 0.5*exp(-Y[x >  XI])
    names(f) <- NULL
    return(f)
}