inst/non_lin_models/Mod_negexpo.R

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

#NEGATIVE EXPONENTIAL (Holdridge et al. 1971)
model <- list(
  name=c("Negative exponential"),
  formula=expression(S == d*(1 - exp(-z*A) )),
  exp=expression(d*(1-exp(-z*A))),
  shape="convex",
  asymp=function(pars)pars["d"],
  #limits for parameters
  parLim = c("Rplus","Rplus"),
  custStart=function(data)c(max(data$S),.01),
  #initials values function
  init=function(data){
    d=max(data$S) 
    Z=(-log( (-data$S/(max(data$S)+1))+1))/data$A
    z = mean(Z)
    c(d,z)}
)