simulateNeutralDriftData: Simulate data from neutral drift model for testing inference.

Description Usage Arguments Value Examples

View source: R/fit_NeutralDrift.R

Description

simulateNeutralDriftData will use the analytical neutral drift model solutions to similate count data.

Usage

1
2
3
4
5
6
7
8
simulateNeutralDriftData(
  lambda,
  Ns,
  tau,
  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.

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

1
2
time_points = c(4, 7, 10, 14, 21)
x = simulateNeutralDriftData(0.1, 5, 2, time_points, 8, 300)

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