draw.clone: Graphic representation of clone growing upon a finite time.

In rcqls/flan: Fluctuation Analysis of Mutant Cell Counts

Description

Simulates a clone up to a time t and represents the clone as a binary tree.

Usage

  draw.clone(t,mutprob=1e-2,fitness=1,death=0,
	     dist=list("lnorm",meanlog=-0.3795851,sdlog=0.3016223))

Arguments

t	time of end of experiment .
mutprob	mutation probability: numeric between 0 and 1. By default 1e-2.
fitness	fitness parameter: numeric positive.
death	death probability: numeric between 0 and 0.5.
dist	lifetime distribution for mutant cells. See Details.

Details

This function does not provide a representation of a realistic realization of a mutation model (mutation probability too high, time of end of experiment to small).

The argument dist is a list beginning with the distribution name followed by its parameters, and must be one the 4 following distributions: "dirac"(location), "exp"(rate), "lnorm"(meanlog, sdlog), "gamma"(shape, scale).

See Also

rflan

Examples

# Luria-Delbrück model, mutation probability 1e-2, fitness 1
draw.clone(t=9,dist=list(name="exp",rate=1))

# Luria-Delbrück model, mutation probability 0.1, fitness 0.6
draw.clone(t=9,mutprob=0.1,fitness=0.6,dist=list(name="exp",rate=1))

# Haldane model, mutation probability 1e-2, fitness 1
draw.clone(t=7,dist=list(name="dirac",location=1))

# Lognormal lifetime distribution 
draw.clone(t=7,fitness=0.5,death=0.1)

# Luria-Delbrück model with positive cell death probability
draw.clone(t=7,dist=list(name="exp",rate=1),death=0.2)

rcqls/flan documentation built on May 27, 2019, 3:05 a.m.