run_neutral_trajectory: Runs a simple single compartment neutral simulation with a...

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

Runs a simple single compartment neutral simulation with a specified trajectory

Description

Runs a simple single compartment neutral simulation with a specified trajectory

Usage

run_neutral_trajectory(simpop, initial_division_rate, trajectory)

Arguments

simpop	Rate of symmetric cell division during development
initial_division_rate	Rate of symmetric cell division once population equilibrium is reached.
trajectory	data.frame - with fields ts(timestamp in days),target_pop_size,division_rate
nyears	Total number of years to run the simulation

Value

simpop object.

Examples

trajectory=data.frame(ts=365*(1:80),target_pop_size=5e4+100*(1:80),division_rate=1/(2*190))
trajectory$target_pop_size[5:10]=2*trajectory$target_pop_size[5:10]
trajectory$target_pop_size[11:15]=0.2*trajectory$target_pop_size[11:15]
sp=run_neutral_trajectory(NULL,0.5,trajectory)
plot(sp)

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