Title: | Predict Regional Community Composition |
---|---|
Description: | Predict regional community composition at a fine spatial resolution using only sparse biological and environmental data. The package is based on the DynamicFOAM algorithm described in Mokany et al. (2011) <doi:10.1111/j.1461-0248.2011.01675.x>. |
Authors: | Craig Simpkins [aut, cre], Sebastian Hanss [aut], Maximilian Hesselbarth [aut], Matthias Spangenberg [aut], Jan Salecker [aut] |
Maintainer: | Craig Simpkins <[email protected]> |
License: | GPL-3 |
Version: | 1.0.2 |
Built: | 2024-10-27 04:19:41 UTC |
Source: | https://github.com/cran/spectre |
Matrix of predicted alpha diversity in each cell.
alpha_list
alpha_list
vector.
Calculate commonness error
calc_commonness_error(x, objective_matrix)
calc_commonness_error(x, objective_matrix)
x |
Results object from run_optimization_min_conf. |
objective_matrix |
Matrix from (modeled) alpha-diversity and Bray-Curtis dissimilarity |
Calculate mean absolute commonness error (MAE_c) and relative commonness error in percentage (RCE).
vector
Total (estimated) species in the system.
estimated_gamma
estimated_gamma
numeric
Creates a pairwise site by site commonness matrix from estimates of species richness and Bray-Curtis dissimilarity.
generate_commonness_matrix_from_gdm(gdm_predictions, alpha_list)
generate_commonness_matrix_from_gdm(gdm_predictions, alpha_list)
gdm_predictions |
a square pairwise |
alpha_list |
a |
generate_commonness_matrix_from_gdm
uses a vector of
estimated species richness per site and a pairwise matrix of site by site
Bray-Curtis dissimilarity (we recommend using the gdm-package
(Fitzpatrick et al. 2020) to generate this matrix) to produce a matrix of
the estimated species in common between site pairs (referred to as a
commonness matrix). The commonness between sites is calculated using
Where is the dissimilarity between sites,
is
the species in common between sites, and S is the number of species in
each site. For more details see Mokany et al 2011.
A pairwise site by site matrix
of the number of species in
common between each site pair, with dimensions equal to that of the
provided dissimilarity matrix.
Mokany, K., Harwood, T.D., Overton, J.M., Barker, G.M., &
Ferrier, S. (2011). Combining and
diversity models to fill
gaps in our knowledge of biodiversity. Ecology Letters, 14(10), 1043-1051.
gdm
packageList with example data created using the gdm
package
minimal_example_data
minimal_example_data
list
Plot commonness between observed and optimized data
plot_commonness(x, target)
plot_commonness(x, target)
x |
Results object of run_optimization_min_conf() |
target |
Pairwise matrix of species in common. |
Plot a heatmap of commonness between observed data and optimized data. This visual style allows for easier spatial understanding of commonness differences to be ascertained.
ggplot
Plot the absolute error
plot_error(x)
plot_error(x)
x |
Results object from run_optimization_min_conf |
Plot error over time
ggplot
xxx
Generate an optimized estimate of community composition (species presences and absences) for every site in the study area.
run_optimization_min_conf( alpha_list, total_gamma, target, max_iterations, partial_solution = NULL, fixed_species = NULL, seed = NA, verbose = TRUE, interruptible = TRUE )
run_optimization_min_conf( alpha_list, total_gamma, target, max_iterations, partial_solution = NULL, fixed_species = NULL, seed = NA, verbose = TRUE, interruptible = TRUE )
alpha_list |
|
total_gamma |
Total number of species present throughout the entire landscape. |
target |
Pairwise matrix of species in common between each site by site pair. Only the upper triangle of the matrix is actually needed. |
max_iterations |
The maximum number of iterations that the optimization algorithm may run through before stopping. |
partial_solution |
An initial |
fixed_species |
Fixed partial solution with species that are considered as given. Those species are not going to be changed during optimization. |
seed |
Seed for random number generator. Seed must be a positive integer value.
|
verbose |
If |
interruptible |
Allow a run to be interrupted before completion. |
run_optimization_min_conf
is the core function of the
spectre
package. The underlying algorithm of this function is
adapted from Mokany et al. (2011). A pairwise commonness matrix (having the
same structure as the target
matrix) is calculated from the
partial_solution
matrix and the value difference with the
target
determined. If a difference is present and depending on the
set stopping criteria the algorithm continues. A random site in the
presence/absence matrix is selected, and a random presence record at this
site replaced with an absence. Every absence in the selected site is then
individually flipped to a presence and the value difference with the
objective recorded. The presence record which resulted in the lowest value
difference (minimum conflict) is retained. This cycle continues, with a
random site selected every iteration, until the pairwise commonness and
objective matrices match or the algorithm runs beyond the
max_iterations
.
A species presence/absence matrix
of the study landscape.
Mokany, K., Harwood, T.D., Overton, J.M., Barker, G.M., &
Ferrier, S. (2011). Combining and
diversity models to fill
gaps in our knowledge of biodiversity. Ecology Letters, 14(10), 1043-1051.
The goal of spectre
is to provide an open source tool capable of predicting regional community composition at fine spatial resolutions using only sparse biological and environmental data.
Maintainer: Craig Simpkins [email protected]
Authors:
Sebastian Hanss
Maximilian Hesselbarth
Matthias Spangenberg
Jan Salecker
Pairwise matrix of species in common.
target_matrix
target_matrix
matrix