renCWR: Perform CWR simulation, Calculate Reny entropy (alfa=1,...

In DeadParrot69/CWERNI: Community Wide Evolutionary Rescue Neutral Implementation

Description

Perform CWR simulation, Calculate Reny entropy (alfa=1, alfa=2), re-simulate neutral communities using SADISA parameters estimated on CWR simulation, simulation neutral communities, calculate their Reny entropy (alfa=1, alfa=2) to generate a distribution.

Usage

renCWR(tmax, b1, d1, k1, interval, b2, d2, m12, m21, orgsimmaxcount,
  totalresamp)

Arguments

tmax	Arbitrary units of time the simulation should be run
b1	Birthrate of the mutants
d1	Deathrate of the mutants
k1	Carrying capacity
interval	After how many evenst should the population state be saved to the output matrix
b2	Birthrate of the residents
d2	Death rate of the residents
m12	Mutation rate of mutants to residents
m21	Mutation rate of residents to mutants
orgsimmaxcount	how many CWR simulations should be performed
totalresamp	how many neutral communities should be simulated for each single CWR simulation

Value

This function saves two .pdf plots for each CWR simulation performed, showing the Reny entropy of each CWR simulation and the distribution of Reny entorpy generated from the neutral re-simulations

Author(s)

Timo van Eldijk

Examples

renCWR(10, 0.6, 0.1, 16000,200, 0.05,0.1, 0, 0.0005, 1,1)

DeadParrot69/CWERNI documentation built on April 21, 2020, 3:10 a.m.