R/pdfgum.R

"pdfgum" <-
function(x,para) {
   if(! are.pargum.valid(para)) return()
   U <- para$para[1]
   A <- para$para[2]

   Y <- -(x - U)/A
   f <- A^(-1)*exp(Y)*exp(-exp(Y))

   names(f) <- NULL
   f[! is.finite(f)] <- NA
   f[is.na(f)] <- 0 # decision Dec. 2015
   return(f)
}

Try the lmomco package in your browser

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

lmomco documentation built on May 29, 2024, 10:06 a.m.