ppoig <- function(q, frac, rate, shape, lower.tail=TRUE, log.p=FALSE) {
if(length(frac) > 1 | length(rate) > 1 | length(shape) > 1) stop("Vectorization not implemented for the parameters")
nonwhole <- !is.wholenumber(q)
y <- suppressWarnings(cumsumW(dpoig, q, list(frac=frac, rate=rate, shape=shape), lower.tail, log.p, pad=FALSE))
if (any(nonwhole)) {
warning("non integer values in q")
if(log.p) y[nonwhole] <- -Inf
else y[nonwhole] <- 0
}
if(any(is.nan(y))) warning("NaNs produced")
return(y)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.