R/dominanceDecay.R

In MarioJose/nicheApport: Tokeshi's Niche Apportionment Species Abundance Distributions Models

# Dominance Decay
dominanceDecay <- function (N, S, count = FALSE){
  k <- runif(n = S, min = 0, max = 1)
  if(count){
    r <- floor(N * k[1])
    r[2] <- N - r[1]
    for(i in 3:S){
      r <- sort(r, decreasing = FALSE)
      r[i] <- floor(r[i - 1] * k[i])
      r[i - 1] <- abs(r[i - 1] - r[i])
    }
  }
  else{
    r <- N * k[1]
    r[2] <- N - r[1]
    for(i in 3:S){
      r <- sort(r, decreasing = FALSE)
      r[i] <- r[i - 1] * k[i]
      r[i - 1] <- abs(r[i - 1] - r[i])
    }
  }
  return(sort(r, decreasing = TRUE))
}