cancerit / ascatNgs

Somatic copy number analysis using WGS paired end wholegenome sequencing
http://cancerit.github.io/ascatNgs/
GNU Affero General Public License v3.0
68 stars 17 forks source link

Force output to generate ASPCF data #115

Open tomvdbosch opened 2 years ago

tomvdbosch commented 2 years ago

Hi,

I am using ASCATngs on 40 WGS samples, and on 39 it works perfectly. The 40th I had to use the -f option to force the output, as no viable solution was found. I think this is due to the nature of the sample, which seems to be a mixture of a highly differentiated and less differentiated part of the tumour and thus is extremely heterogeneous which messes up the fitting to integer copy numbers.

However, I still get the ASPCF output plot which shows the logR and BAF values and still looks decent as shown below. I would like to have these LogR and BAF values and the segmentation for further analysis, but the .copynumber.segBAF.bw, .copynumber.logR.bw etc. which usually contain these values are empty. Is there a way to find the data in this plot?

CF48530 ASPCF