Weighting of mutations

[lima1]

Hi,

thanks for a great package!

I'm using it for tumor-only sequencing with our https://github.com/lima1/PureCN tool. PureCN essentially predicts somatic status by adjusting allele frequencies for tumor purity and allele-specific copy number. This works pretty well and I get a decent correlation of deconstructSigs scores with matched tumor/normal. But I was wondering if you see an easy way to incorporate the uncertainty of the somatic prediction. Essentially weighting mutations by somatic posterior probability would be perfect. (for example especially copy number losses can result in cases where germline and somatic have very similar expected allele frequencies).

I understand if you don't have the bandwidth for this. But if you see an easy solution, I'm happy to dig into the code and method.

Thanks in advance,
Markus