GavinHaLab / ichorCNA

Estimating tumor fraction in cell-free DNA from ultra-low-pass whole genome sequencing.
GNU General Public License v3.0
12 stars 10 forks source link

In what situations should maxGenomeSubclonal and maxCnaSubclonal be reconsidered? #18

Open ilykos opened 11 months ago

ilykos commented 11 months ago

We are working with noisy (MAD < 0.15, avg. 0.10) low pass WGS data from FFPE blocks where the DNA was sampled from a histologically delineated tumor region.

IchorCNA frequently suggests results with whole genome amplification (dark red). We have done some work on understanding the algorithm for solution selection and have following questions

  1. Is it correct that the solution is ALWAYS the highest log likelihood after filtering by maxGenomeSubclonal and maxCnaSubclonal?
  2. Do you have any examples of situations where we should consider changing the default values (maxGenomeSubclonal=0.5, maxCnaSubclonal=0.7) considering the propreties of our data?