R/dpareto.R

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

dpareto <- function(x, shape, scale = min(x), log = FALSE){
	shape[ !is.finite(shape) | shape <= 0] <- NaN
	scale[ !is.finite(scale) | scale <= 0] <- NaN
    lny <- log(shape) + shape*log(scale) - (shape+1)*log(x)
	if (any(is.nan(lny))) warning ("NaNs produced")
	lny [ x < scale ] <- -Inf
  if (log) return(lny)
  else return(exp(lny))
}