hap.py icon indicating copy to clipboard operation
hap.py copied to clipboard

How vary AF and DP for PR ROC?

Open KewinOgink opened this issue 4 years ago • 0 comments

Hi, I don't understand how the ROC Precision recall curve can be made with other metrics.

As I understand the --roc switch to feature on must be in the query vcf. This can be in the INFO field, but in my case this is empty. I would however like to vary the AF and FORMAT fields, but --roc FORMAT.AF only gives AF values for SNPs of 0.5714, 0.5049 and 0 in the roc.csv happy output file.

How can I vary the FORMAT AF (i.e. between 0.2 and 0.9 with steps of 0.1) for PR ROC plots and is there an example query VCF?

Best,

Kewin

KewinOgink avatar Jun 15 '21 08:06 KewinOgink