lambda: Expected migration rate per deme

In robjohnnoble/demonanalysis: Analysis and plotting of demon data

Description

Expected migration rate per deme

Usage

lambda(i, m, r1, r2, migration_type = 0, d = NA, K = NA)

Arguments

i	number of cells in the deme
m	migration rate per cell, relative to birth rate
r1	reference birth rate
r2	birth rate of migrating cells, relative to r1
migration_type	0 or 1 (see details)
d	migration distance relative to 1/sqrt(K) (see details)
K	deme carrying capacity (only required if d is not NA)

Details

If migration_type = 0 (meaning that migration is attempted only after a birth event) then the migration rate will be multiplied by a factor of r2 * r1. Otherwise no adjustment will be made.

If d and K are set then the calculation assumes that min(d / sqrt(K), 1) is the probability that a migration attempt will land outside the deme.

Importantly, this function does not account for the chance that the migrating cell will land in a deme that already occupied.

Value

The migration rate.

Examples

# if r1 * r2 = 1 then results are identical for migration_type = 0 and migration_type = 1: 
lambda(1, 0.1, 0.9, 1/0.9, 0)
lambda(1, 0.1, 0.9, 1/0.9, 1)

robjohnnoble/demonanalysis documentation built on June 30, 2020, 12:47 a.m.