trials_tetraploid: Simulations to assess the accuracy of SeqEM when individuals...
In liuhyustc/StatComp21017: Project for 'statistical computing' course

Parameters are estimated using SeqEM

1	trials_tetraploid(m, n, lambda, alpha, p)

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

A matrix of estimates

## Not run: 
tr <- trials_tetraploid(1000,100,8,0.01,0.1)
apply(tr,1,mean)
apply(tr,1,sd)

## End(Not run)

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

liuhyustc/StatComp21017 index

Note that we can't provide technical support on individual packages. You should contact the package authors for that.

Description