Algorithm to implement Fraction and Allele specific Copy number Estimate from Tumor/normal Sequencing.
The new version estimates the log-ratio level corresponding to the diploid state. It is embedded into the procSample call. In terms of using the package you can now do:
xx <- preProcSample(...)
# specify cval you like
oo <- procSample(xx, cval = 300)And go straight to
emcncf(oo)The output of procSample now has 4 elements:
jointseg– same as beforeout– as before with 3 additional columns: cf, tcn, lcndipLogR– the estimated location of diploid log-ratio valueflags– this gives an indication of whether the dipLogR is estimated well.
If flags is NULL then no obvious problem with the dipLogR estimate. It can have two other comments: "mafR not sufficiently small" and "could be polyclonal 1 copy loss". The first one means that there aren't segments with sufficiently balanced alleles and so the estimate may not be great. The second one means that it looks like genome doubling; the only lower level segments are 2+1 and 2+0 (from 2+2); a model without genome doubling but single copy loss (1+0 from 1+1) with two different cellular fraction could fit better.
Note: I am not claiming that I have covered all possible scenarios that lead to bad dipLogR estimate. If you come across anything that doesn't seem right but flags is NULL let me know.