run_transient_selection: Runs a single compartment/single driver scenario with...

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

Runs a single compartment/single driver scenario with transient selection

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_transient_selection(
  initial_division_rate,
  final_division_rate,
  target_pop_size = 1e+05,
  nyears_driver_acquisition = 15,
  nyears_transient_end = 30,
  nyears = 40,
  fitness = 0.2,
  minprop = 0.05
)

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_transient_end	When selective advantage is set to 0.
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

tselsim=run_transient_selection(0.05,1/365,target_pop_size = 5e4,nyears_driver_acquisition=15,
nyears_transient_end=30,nyears=50,fitness=0.5)

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