cancerTiming | R Documentation |

Fixed E-step of EM algorithm (in internal function

`.estimateQ`

) to provide numerical stability for calculating probability of true allele frequency given the observed counts. Previously the function could give NaN probabilities due to the dividing of the sum of probabilities which might be small, and therefore cause the timing to crash.Added

`plotPi0`

to create plots of confidence intervals of estimates of pi0

Added functions

`eventTimingOverList`

and`getPi0Summary`

to aid in applying`eventTiming`

over many segmented regions at one time.In the EM algorithm, the common event histories of (0,2) and (1,2) are given an algebraic solution in the M-Step to avoid calling the optim function, which should speed calculations. Identification is by matching of the history matrix.

Added argument

`mutationId`

for individual mutations in`eventTiming`

, which allows for tracking of assignments of mutationsMade argument

`returnAssignments`

and`returnData`

for function`eventTiming`

a single option of`returnAssignments`

and the assignments are returned in one data.frame object that has the original data (`x`

and`m`

) in addition to the probabilities assigned to them.-
`GenomicRanges`

is no longer a package listed as 'Depends' but merely 'Suggests'. It is only used for a tiny helper function (`labelSeg`

) useful for labeling segments as to the arm of the chromosome, e.g. '3q' or '17p'. Requiring GenomicRanges was a large installation overhead for such a minor function since the package is part of bioconductor.