Description Usage Arguments Value Author(s) Examples
The estimation of the admixture rate is accomplished by fitting the input tumor somatic copy number (somatic ratio*2) of all segments with a Bayesian finite mixture model, with components centered at the discrete levels. Each segment was assigned with a discrete level based on corresponding posterior probability. Segments with ambiguous assignments will be classified as candidate subclonal events and excluded from admixture rate inference. The admixture rate will be estimated by an optimal solution contributed by explanation of tumor copy number with all remaining segments as integer level.
1 | admixtureRate(segmentwithratio, mcmc = 10000, burnin = 5000, p = 0.01, weight=FALSE)
|
segmentwithratio |
A GRanges object, segments with annotation of somatic ratio. |
mcmc |
number of MCMC iteration. |
burnin |
number of MCMC iteration for burnin. |
p |
posterior probability cutoff for ambiguous integer copy number assignments. |
weight |
whether penalize the segment based on its length. |
admix |
Admixture rate. |
mu |
Posterior mean for each discrete level. |
cluster |
A two-column matrix for input somatic ratio and corresponding integer somatic copy number level. |
Mengjie Chen
1 2 3 4 5 6 7 8 | data(segwithratio)
seg <- GRanges(seqnames=segwithratio$chromosome,
ranges=IRanges(start=segwithratio$start,
end=segwithratio$end),
medLAF=segwithratio$medLAF,
medgLAF=segwithratio$medgermlineLAF,
ratio=segwithratio$ratio)
bb <- admixtureRate(seg)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.