BoevaLab / FREEC

Control-FREEC: Copy number and genotype annotation in whole genome and whole exome sequencing data
152 stars 49 forks source link

input for ABSOLUTE for CCF calculation #118

Open amootta opened 1 year ago

amootta commented 1 year ago

Hi there, I've run FREEC on all of my tumour-normal matched samples to call somatic CNAs. I'm wondering if there is a way to make the input for ABSOLUTE, which requires the columns "Chromosome", "Start", "End", "Num_Probes" and "Segment_Mean". I would like to try and use it to calculate the CCF of my variants. Please let me know if possible

Thanks very much

valeu commented 1 year ago

Hi, no, currently FREEC does not output such a file. However, it can created based on the _ratio.txt output. Could you please share the first lines of the format you need?

amootta commented 1 year ago

Hi, thanks for your response. It would look something like this (this is just an example i pulled out of someone else's query elsewhere): image

valeu commented 1 year ago

You can try this R script: FREEC_ratio2Absolute.zip You can run it as: cat FREEC_ratio2Absolute.R | R --slave --args <*_ratio.txt>

valeu commented 1 year ago

But keep in mind that you will get -inf if median ratio is 0. - Here FREEC outputs median and not mean. And the last column will be log2Median

valeu commented 1 year ago

I have also added it to https://github.com/BoevaLab/FREEC/tree/master/scripts