run_selection_sim: Runs a single compartment/single driver scenario

In NickWilliamsSanger/rsimpop: Simulates Cellular Evolution using a Continous Time Birth Death Model

Runs a single compartment/single driver scenario

Description

The driver is subject to stochastic extinction and so is repeatedly dropped in until it "takes". The prevailing state when the driver is introduced is saved and reinstated with each attempted introduction.

Usage

run_selection_sim(
  initial_division_rate = 0.1,
  final_division_rate = 1/365,
  target_pop_size = 1e+05,
  nyears_driver_acquisition = 15,
  nyears = 40,
  fitness = 0.2,
  minprop = 0.001,
  mindriver = 1,
  maxtry = 40
)

Arguments

initial_division_rate	Rate of symmetric cell division during development
final_division_rate	Rate of symmetric cell division once population equilibrium is reached.
target_pop_size	Size of target population
nyears_driver_acquisition	When driver is acquired.
nyears	Total number of years to run the simulation
fitness	relative fitness advantage. Cells carrying this divide at a rate=(1+fitness)*baserate
minprop	Minimum aberrant cell fraction at nyears for driver to be regarded as have "taken"
maxtry	Maximum number of attempts to introduce a driver

Value

simpop object.

Examples

selsim=run_selection_sim(0.05,1/365,target_pop_size = 5e4,nyears = 50,fitness=0.3)

NickWilliamsSanger/rsimpop documentation built on Sept. 6, 2024, 12:42 a.m.