seqem_tetraploid: SeqEM method to estimate parameters for tetraploid...

Description Usage Arguments Value Examples

View source: R/StatComp21017R.R

Description

EM algorithm is used to estimate nucleotide-read error and MAF

Usage

1
seqem_tetraploid(reads, iters, epsilon)

Arguments

reads

a data set for read sequences

iters

maximum iterations

epsilon

convergence criterion

Value

estimation of alpha and p

Examples

1
2
3
4
5
## Not run: 
sample <- sample_tetraploid(10,10,0.01,0.1)
seqem_tetraploid(sample, 100, 1e-8)

## End(Not run)

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