Description Usage Arguments Value Author(s) References See Also Examples

This function computes the maximum likelihood estimates of the parameters of
the DAISIE model with clade-specific diversity-dependence for data from
lineages colonizing an island. It also outputs the corresponding
loglikelihood that can be used in model comparisons.
The result of sort(c(idparsopt, idparsfix, idparsnoshift)) should be
identical to c(1:10). If not, an error is reported that the input is
incoherent. The same happens when the length of initparsopt is different
from the length of idparsopt, and the length of parsfix is different from
the length of idparsfix.

Including the 11th parameter (p_f) in either
idparsopt or idparsfix (and therefore initparsopt or parsfix) is optional.
If this parameter is not specified, then the information in the data is
used, otherwise the information in the data is overruled.

1 2 3 4 5 6 7 | ```
DAISIE_ML_CS(datalist, datatype = "single", initparsopt, idparsopt,
parsfix, idparsfix, idparsnoshift = 6:10, idparsmat = NULL,
res = 100, ddmodel = 0, cond = 0, island_ontogeny = NA,
eqmodel = 0, x_E = 0.95, x_I = 0.98, tol = c(1e-04, 1e-05,
1e-07), maxiter = 1000 * round((1.25)^length(idparsopt)),
methode = "lsodes", optimmethod = "subplex", CS_version = 1,
verbose = 0, tolint = c(1e-16, 1e-10))
``` |

`datalist` |
Data object containing information on colonisation and
branching times. This object can be generated using the DAISIE_dataprep
function, which converts a user-specified 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. |

`datatype` |
Sets the type of data: 'single' for a single island or archipelago treated as one, and 'multiple' for multiple archipelagoes potentially sharing the same parameters |

`initparsopt` |
The initial values of the parameters that must be optimized |

`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 |

`idparsmat` |
For datatype = 'multiple' only: Matrix containing the ids
of the parameters, linking them to initparsopt and parsfix. Per island
system we use the following order: |

`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 diversity-dependence: |

`cond` |
cond = 0 : conditioning on island age |

`island_ontogeny` |
type of island ontonogeny. If NA, then constant ontogeny is assumed |

`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 non-endemic 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 ODE-solver. 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` |
For internal testing purposes only. Default is 1, the original DAISIE code. |

`verbose` |
sets whether parameters and likelihood should be printed (1) or not (0) |

`tolint` |
Vector of two elements containing the absolute and relative tolerance of the integration |

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 carrying-capacity |

`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 |

`lambda_c2` |
gives the maximum likelihood estimate of lambda^c2, the rate of cladogenesis for the optional second group of species |

`mu2` |
gives the maximum likelihood estimate of mu2, the extinction rate for the optional second group of species |

`K2` |
gives the maximum likelihood estimate of K2, the carrying-capacity for the optional second group of species |

`gamma2` |
gives the maximum likelihood estimate of gamma2, the immigration rate for the optional second group of species |

`lambda_a2` |
gives the maximum likelihood estimate of lambda^a2, the rate of anagenesis for the optional second group of species |

`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 |

Rampal S. Etienne

Valente, L.M., A.B. Phillimore and R.S. Etienne (2015). Equilibrium and non-equilibrium dynamics simultaneously operate in the Galapagos islands. Ecology Letters 18: 844-852. <DOI:10.1111/ele.12461>.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 | ```
cat("
### When all species have the same rates, and we want to optimize all 5 parameters,
# we use:
utils::data(Galapagos_datalist)
DAISIE_ML(
datalist = Galapagos_datalist,
initparsopt = c(2.5,2.7,20,0.009,1.01),
ddmodel = 11,
idparsopt = 1:5,
parsfix = NULL,
idparsfix = NULL
)
### When all species have the same rates, and we want to optimize all parameters
# except K (which we set equal to Inf), we use:
utils::data(Galapagos_datalist)
DAISIE_ML(
datalist = Galapagos_datalist,
initparsopt = c(2.5,2.7,0.009,1.01),
idparsopt = c(1,2,4,5),
parsfix = Inf,
idparsfix = 3
)
### When all species have the same rates except that the finches have a different
# rate of cladogenesis, and we want to optimize all parameters except K (which we
# set equal to Inf), fixing the proportion of finch-type species at 0.163, we use:
utils::data(Galapagos_datalist_2types)
DAISIE_ML(
datalist = Galapagos_datalist_2types,
initparsopt = c(0.38,0.55,0.004,1.1,2.28),
idparsopt = c(1,2,4,5,6),
parsfix = c(Inf,Inf,0.163),
idparsfix = c(3,8,11),
idparsnoshift = c(7,9,10)
)
### When all species have the same rates except that the finches have a different
# rate of cladogenesis, extinction and a different K, and we want to optimize all
# parameters, fixing the proportion of finch-type species at 0.163, we use:
utils::data(Galapagos_datalist_2types)
DAISIE_ML(
datalist = Galapagos_datalist_2types,
ddmodel = 11,
initparsopt = c(0.19,0.09,0.002,0.87,20,8.9,15),
idparsopt = c(1,2,4,5,6,7,8),
parsfix = c(Inf,0.163),
idparsfix = c(3,11),
idparsnoshift = c(9,10)
)
### When all species have the same rates except that the finches have a different
# rate of extinction, and we want to optimize all parameters except K (which we
# set equal to Inf), and we also# want to estimate the fraction of finch species
# in the mainland pool. we use:
utils::data(Galapagos_datalist_2types)
DAISIE_ML(
datalist = Galapagos_datalist_2types,
initparsopt = c(2.48,2.7,0.009,1.01,2.25,0.163),
idparsopt = c(1,2,4,5,7,11),
parsfix = c(Inf,Inf),
idparsfix = c(3,8),
idparsnoshift = c(6,9,10)
)
### When we have two islands with the same rates except for immigration and anagenesis rate,
# and we want to optimize all parameters, we use:
utils::data(Galapagos_datalist)
DAISIE_ML(
datalist = list(Galapagos_datalist,Galapagos_datalist),
datatype = 'multiple',
initparsopt = c(2.5,2.7,20,0.009,1.01,0.009,1.01),
idparsmat = rbind(1:5,c(1:3,6,7)),
idparsopt = 1:7,
parsfix = NULL,
idparsfix = NULL
)
### When we consider the four Macaronesia archipelagoes and set all parameters the same
# except for rates of cladogenesis, extinction and immigration for Canary Islands,
# rate of cladogenesis is fixed to 0 for the other archipelagoes,
# diversity-dependence is assumed to be absent
# and we want to optimize all parameters, we use:
utils::data(Macaronesia_datalist)
DAISIE_ML(
datalist = Macaronesia_datalist,
datatype = 'multiple',
initparsopt = c(1.053151832,0.052148979,0.512939011,0.133766934,0.152763179),
idparsmat = rbind(1:5,c(6,2,3,7,5),1:5,1:5),
idparsopt = c(2,4,5,6,7),
parsfix = c(0,Inf),
idparsfix = c(1,3)
)
")
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.