Description Usage Arguments Details References Examples
Computes maximum loglikelihood and corresponding parameters for the requested model using the independent-species approach. For optimization it uses various auxiliary functions in the DDD package.
1 2 3 4 |
abund |
abundance vector or a list of abundance vectors. When a list is provided and mult = 'mg' (the default), it is assumed that the different vectors apply to different guilds. When mult = 'ms' then the different vectors apply to multiple samples. from the same metacommunity. In this case the vectors should have equal lengths and may contain zeros because there may be species that occur in multiple samples and species that do not occur in some of the samples. |
initpars |
a vector of initial values of the parameters to be optimized and fixed. See |
idpars |
a vector stating whether the parameters in |
labelpars |
a vector, a list of vectors or a list of lists of vectors indicating the labels
integers (starting at 1) of the parameters to be optimized and fixed. These integers correspond to
the position in |
model |
the chosen combination of metacommunity model and local community model as a vector, e.g. c('pm','dl') for a model with point mutation in the metacommunity and dispersal limitation. The choices for the metacommunity model are: 'pm' (point mutation), 'rf' (random fission), 'pr' (protracted speciation), 'dd' (density-dependence). The choices for the local community model are: 'dl' (dispersal limitation), 'dd' (density-dependence). |
mult |
When set to 'single' (the default), the loglikelihood for a single sample and single guild is computed. When set to 'mg', the loglikelihood for multiple guilds is computed. When set to 'ms' the loglikelihood for multiple samples from the same metacommunity is computed. |
tol |
a vector containing three numbers for the relative tolerance in the parameters, the relative tolerance in the function, and the absolute tolerance in the parameters. |
maxiter |
sets the maximum number of iterations |
optimmethod |
sets the optimization method to be used, either subplex (default) or an alternative implementation of simplex. |
num_cycles |
the number of cycles of opimization. If set at Inf, it will do as many cycles as needed to meet the tolerance set for the target function. |
Not all combinations of metacommunity model and local community model have been implemented yet. because this requires checking for numerical stability of the integration. The currently available model combinations are, for a single sample, c('pm','dl'), c('pm','rf'), c('dd','dl'), c('pr','dl'), c('pm','dd'), and for multiple samples, c('pm','dl').
Haegeman, B. & R.S. Etienne (2017). A general sampling formula for community structure data. Methods in Ecology & Evolution 8: 1506-1519. doi: 10.1111/2041-210X.12807
1 2 3 4 5 6 7 8 9 10 11 | utils::data(datasets);
utils::data(fitresults);
result <- SADISA_ML(
abund = datasets$dset1.abunvec[[1]],
initpars = fitresults$fit1a.parsopt[[1]],
idpars = c(1,1),
labelpars = c(1,2),
model = c('pm','dl'),
tol = c(1E-1, 1E-1, 1E-1)
);
# Note that tolerances should be set much lower than 1E-1 to get the best results.
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.