Description Usage Arguments Value
Computes MLE for single type species under a clade specific scenario
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23  DAISIE_ML1(
datalist,
initparsopt,
idparsopt,
parsfix,
idparsfix,
idparsnoshift = 6:10,
res = 100,
ddmodel = 0,
cond = 0,
eqmodel = 0,
x_E = 0.95,
x_I = 0.98,
tol = c(1e04, 1e05, 1e07),
maxiter = 1000 * round((1.25)^length(idparsopt)),
methode = "lsodes",
optimmethod = "subplex",
CS_version = 1,
verbose = 0,
tolint = c(1e16, 1e10),
island_ontogeny = NA,
jitter = 0
)

datalist 
Data object containing information on colonisation and
branching times. This object can be generated using the DAISIE_dataprep
function, which converts a userspecified data table into a data object,
but the object can of course also be entered directly.
It is an R list object with the following elements. 
initparsopt 
The initial values of the parameters that must be optimized, they are all positive. 
idparsopt 
The ids of the parameters that must be optimized. The ids
are defined as follows: 
parsfix 
The values of the parameters that should not be optimized. 
idparsfix 
The ids of the parameters that should not be optimized, e.g. c(1,3) if lambda^c and K should not be optimized. 
idparsnoshift 
For datatype = 'single' only: The ids of the parameters that should not be different between two groups of species; This can only apply to ids 6:10, e.g. idparsnoshift = c(6,7) means that lambda^c and mu have the same values for both groups. 
res 
Sets the maximum number of species for which a probability must be computed, must be larger than the size of the largest clade. 
ddmodel 
Sets the model of diversitydependence: 
cond 
cond = 0 : conditioning on island age 
eqmodel 
Sets the equilibrium constraint that can be used during the
likelihood optimization. Only available for datatype = 'single'. 
x_E 
Sets the fraction of the equlibrium endemic diversity above which the endemics are assumed to be in equilibrium; only active for eqmodel = 13 or 15. 
x_I 
Sets the fraction of the equlibrium nonendemic diversity above which the system is assumed to be in equilibrium; only active for eqmodel = 15. 
tol 
Sets the tolerances in the optimization. Consists of: 
maxiter 
Sets the maximum number of iterations in the optimization. 
methode 
Method of the ODEsolver. See package deSolve for details. Default is "lsodes". 
optimmethod 
Method used in likelihood optimization. Default is "subplex" (see subplex package). Alternative is 'simplex' which was the method in previous versions. 
CS_version 
a numeric or list. Default is 1 for the standard DAISIE model, for a relaxedrate model a list with the following elements:

verbose 
In simulation and dataprep functions a logical,

tolint 
Vector of two elements containing the absolute and relative tolerance of the integration. 
island_ontogeny 
In 
jitter 
Numeric for 
The output is a dataframe containing estimated parameters and maximum loglikelihood.
lambda_c 
gives the maximum likelihood estimate of lambda^c, the rate of cladogenesis 
mu 
gives the maximum likelihood estimate of mu, the extinction rate 
K 
gives the maximum likelihood estimate of K, the carryingcapacity 
gamma 
gives the maximum likelihood estimate of gamma, the immigration rate 
lambda_a 
gives the maximum likelihood estimate of lambda^a, the rate of anagenesis 
loglik 
gives the maximum loglikelihood 
df 
gives the number of estimated parameters, i.e. degrees of feedom 
conv 
gives a message on convergence of optimization; conv = 0 means convergence 
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.