R/ppower.R

In sads: Maximum Likelihood Models for Species Abundance Distributions

ppower <- function(q, s, lower.tail = TRUE, log.p = FALSE){
  if (length(s) > 1) stop("vectorization of parameters is not implemented")
	if( !is.finite(s) | s <= 1 ) return(rep(NaN, length(q)))
	y <- c()
	for (i in 1:length(q)){
		if (is.nan(q[i])) y[i] <- NaN
		else y[i] <- log(sum(1/(1:q[i])^s)) - log(zeta(s))
	}
	y <- exp(y)
	if (any(!is.wholenumber(q))) warning("non integer values in q")
	y[ ! is.wholenumber(q) | q < 1 ] <- 0
	if(!lower.tail) y <- 1 - y
	if(log.p) y <- log(y)
	if (any(is.nan(y))) warning ("NaNs produced")
	return(y)
}