powered by
This function computes the maximum likelihood estimates of the parameters of the Neutral model, using the Etienne Sampling Formula
maxLikelihood.ESF(init_vals, abund, verbose = FALSE)
the output is a list containing the following:
a vector containing the parameter values at the maximum likelihood c(theta, m)
the likelihood at the corresponding parameter values
gives a message on convergence of optimization; conv = 0 means convergence
A vector of initial starting values, of the format c(theta, m)
Vector containing a record of the number of individuals per species
TRUE/FALSE flag, indicates whether intermediate output is shown on screen
Thijs Janzen
Etienne, R.S. (2005). A new sampling formula for neutral biodiversity. Ecology Letters, 8(3), 253-260.
A <- c(1, 1, 1, 3, 5, 8) maxLikelihood.ESF( c(7, 0.1), abund = A)
Run the code above in your browser using DataLab