Learn R Programming

STEPCAM (version 1.2)

STEPCAM_ABC: ABC-SMC inference of the STEPCAM model

Description

ABC-SMC inference of the relative contribution of stochasticy (or random removal), filtering and limiting similarity (e.g. of shared resources (competition) or predators) through the use of stepwise commmunity assembly and fitting Functional Diversity patterns of simulated communities with those of observed communities. It takes quite a lot of computing time to run this function, so it is recommended to first run it for a very simple (species poor, artificial) dataset (see example).

Usage

STEPCAM_ABC(data_abundances, data_species, numParticles, n_traits, plot_number, stopRate, stop_at_iteration, continue_from_file)

Arguments

data_abundances
data_abundances is a matrix containing the abundances of all species over all plots. rows every row corresponds to a different plot or community columns every column corresponds to a different species. For an example of a suitable data_abundances matrix, run the generate.Artificial.Data function
data_species
data_species is a matrix containing the species' specific traits. rows every row corresponds to a different species columns the first column contains the name of the species, every consecutive column corresponds to a different trait It is important to realize that ideally, the number of traits is lower than the number of species in the most species poor plot of the 'data_abundances' matrix, as FRic calculation requires at least as many species as their are traits, and will autmatically omit some traits if this assumption is not met. For an example of a suitable data_species matrix, run the generate.Artificial.Data function
numParticles
The number of particles to be used in the ABC-SMC analysis.
n_traits
The number of traits used in the analysis.
plot_number
Row number of the plot for which we want to infer the contributions of stochasticity, filtering and competition.
stopRate
The SMC procedure stops as soon as the acceptance rate of newly proposed parameter combinations drops below this level. Beware that low values for this parameter might lead into long computation times.
stop_at_iteration
If the user instead wants to perform a fixed number of iterations, this parameter sets the maximum number of iterations (the SMC stops when either the maximum number of iterations is reached, or the acceptance rate drops below the stopRate, whatever happens first)
continue_from_file
Should the ABC-SMC algorithm continue from previous output of the STEPCAM_ABC function?

Value

References

Fons van der Plas, Thijs Janzen, Alejandro Ordonez, Wimke Fokkema, Josephine Reinders, Rampal S. Etienne, and Han Olff 2015. A new modeling approach estimates the relative importance of different community assembly processes. Ecology 96:1502-1515.

Examples

Run this code
## Not run:  
# Artificial.Data <- generate.Artificial.Data(n_species = 40, n_traits = 3,
#                         n_communities = 5, occurence_distribution = 0.2,
#                         average_richness = 0.5, sd_richness = 0.2,
#                         mechanism_random=FALSE)
# O <- STEPCAM_ABC(Artificial.Data$abundances, Artificial.Data$traits,
#                  numParticles = 10, n_traits = 3, plot_number = 1, stopRate = 0.8)
# ## End(Not run)

Run the code above in your browser using DataLab