R/dls.R

In piLaboratory/sads: Maximum Likelihood Models for Species Abundance Distributions

dls <- function(x, N, alpha, log = FALSE){
	N[ !is.finite(N) | N <= 0] <- NaN
	alpha[ !is.finite(alpha) | alpha <= 0] <- NaN
	X <- N/(N+alpha)
	gama <- function(y) {1/log(1/(1-y))}
	y <- gama(X)*(X^x)/x
	if (any(is.nan(y))) warning ("NaNs produced")
	if (any(!is.wholenumber(x))) warning("non integer values in x")
	y[ ! is.wholenumber(x) | x < 1] <- 0
	if(log) return (log(y))
	else return(y)
}