R/dbs.R

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

dbs <- function(x, N, S, log = FALSE){
	N[ !is.finite(N) | N <= 0] <- NaN
	S[ !is.finite(S) | S <= 0] <- NaN
  y <- (S-1)*(1-x/N)^(S-2)/N
	if (any(is.nan(y))) warning ("NaNs produced")
	y[ x < 0 | x > N] <- 0
	if(log) return (log(y))
	else return(y)
}