simulateBiasedDriftData: Simulate data from biased drift model for testing inference.

In MorrisseyLab/CryptDriftR: Analyse crypt lineage tracing data

Description

simulateBiasedDriftData will use the analytical biased stem cell drift model solutions to similate count data.

Usage

simulateBiasedDriftData(
  lambda,
  Ns,
  tau,
  Pr,
  time_points,
  size_measure,
  num_crypts
)

Arguments

lambda	Stem cell replacement rate.
Ns	Number of functional stem cells per crypt.
tau	Time delay until drift starts.
Pr	Bias in stem cell replacement Pr = 0.5 is neutral.
time_points	vector of time point values.
size_measure	Fractions in which the crypts are measured in (e.g. 8ths).
num_crypts	Number of crypts measured per time point.

Value

Will retrun simulated crypt clone count data.

Examples

time_points = c(4, 7, 10, 14, 21)
WT_params = list(lambda = 0.1, N  = 5, tau = 1)
x = simulateBiasedDriftData(WT_params$lambda, WT_params$N, WT_params$tau, Pr = 0.7, time_points, 8, 300)

MorrisseyLab/CryptDriftR documentation built on Aug. 7, 2021, 8:22 p.m.