randomNeedles: Simulate fine needle aspiration

In SITH: A Spatial Model of Intra-Tumor Heterogeneity

Description

Simulate a sampling procedure which takes a fine needle through the simulated tumor and reports the mutation allele frequency of the sampled cells.

Usage

randomNeedles(tumor, nsamples, threshold = 0.05, coverage = 0)

Arguments

tumor	A list which is the output of simulateTumor().
nsamples	The number of samples to take.
threshold	Only mutations with an allele frequency greater than the threshold will be included in the sample.
coverage	If nonzero then deep sequencing with specified coverage is performed.

Details

This sampling procedure is inspired by Chkhaidze et. al. (2019) and simulates fine needle aspiration. A random one-dimensional cross-section of the tumor is chosen, and the cells within this cross section are sampled, reporting mutation allele frequency.

Author(s)

Phillip B. Nicol

References

K. Chkhaidze, T. Heide, B. Werner, M. Williams, W. Huang, G. Caravagna, T. Graham, and A. Sottoriva. Spatially con- strained tumour growth affects the patterns of clonal selection and neutral drift in cancer genomic data. PLOS Computational Biology, 2019. https://doi.org/10.1371/journal.pcbi.1007243.

Examples

out <- simulateTumor(max_pop = 1000)
df <- randomNeedles(tumor = out, nsamples = 5)

SITH documentation built on Jan. 5, 2021, 5:08 p.m.