sample_tetraploid: Generate nucleotide read counts for tetraploids in SeqEM.

In liuhyustc/StatComp21017: Project for 'statistical computing' course

Description

Generate read depths and variant counts with a symmetric error assumption.

Usage

sample_tetraploid(n, lambda, alpha, p)

Arguments

n	sample size
lambda	assume read depth distributed as Poisson with mean lambda
alpha	nucleotide-read error
p	minor allele frequency

Value

a list for reads of size n

Examples

## Not run: 
sample_tetraploid(10,10,0.01,0.1)

## End(Not run)

liuhyustc/StatComp21017 documentation built on Dec. 23, 2021, 11:16 p.m.