is_stable: Evaluate evolutionary stability

In rscherrer/speciomx: Adaptive Dynamics Approximation of the Speciome Model

Evaluate evolutionary stability

Description

Assesses whether an equilibrium trait value is evolutionarily stable, or invasion-proof.

Usage

is_stable(xeq, pars, init, value = FALSE)

Arguments

xeq	Some equilibrium trait value
pars	An unevaluated parameter-list (e.g. as returned by get_default_pars)
init	A vector of two starting values for solving of the demographic equilibrium
value	Whether to return the actual value of the curvature of the fitness function (defaults to FALSE, i.e. only returns whether this value is negative)

Details

This functions computes an expression for the second derivative of the invasion fitness function with respect to the mutant trait value, evaluated at the resident trait value at equilibrium xeq. If this value is negative then the equilibrium is a fitness maximum, i.e. it is evolutionarily stable. It is unstable otherwise.

Value

A boolean

See Also

is_convergent

Examples

pars <- get_default_pars()
is_stable(0, pars, init = rep(1000, 2))

rscherrer/speciomx documentation built on March 28, 2023, 8:49 p.m.