inst/non_lin_models/Mod_p1.R

In txm676/mmSAR2: Fit and Compare Species-Area Relationship Models Using Multimodel Inference

# PERSISTENCE FUNCTION 1 (TJORVE 2009)
model <- list(
  name = c("Persistence function 1"),
  formula = expression(S == c*A^z * exp(-d*A)),
  exp = expression(c*A^z * exp(-d*A)),
  mod = s~c*a^z * exp(-d*a),
  shape = "convex/sigmoid",
  asymp = function(pars)FALSE,
  custStart = function(data)c(5,.25,.15),
  #limits for parameters
  parLim = c("Rplus","Rplus","Rplus"),
  #initials values function
  init = function(data) {
    if (any(data$S == 0)) {
      log.data <- data.frame(A = log(data$A), S = log(data$S + .5))
    } else {
      log.data <- log(data)
    }
    res <- stats::lm(S~A,log.data)$coefficients
    res <- c(res[1],res[2],1)
    names(res) <- c("c","z","d")
    return(res)
    }
)