Generates the transition matrix for a given probability and an interaction network using the niche model.

buildMarkov(nbsp, colo, metaweb, basexti, mn, mx, shape)

Arguments

nbsp

an integer giving the number of species.

colo

a vector giving the colonization rates of all species.

metaweb

a matrix describing the interaction network.

basexti

extinction rate without interaction.

mn, mx

a real number giving the minimum and maximum rate of extinction.

shape

a positive real number controlling the impact of biotic interactions: the larger the value the stronger the impact of interactions.

Value

A list including:

  • makov: the transition of the Markov chain to be solved.

  • speciesRichness: species richness for all states.

  • speciesPresence: matrix of presence for every species in all states.

References

Cazelles, K., Mouquet, N., Mouillot, D. & Gravel, D. On the integration of biotic interaction and environmental constraints at the biogeographical scale.