hap.py
hap.py copied to clipboard
How vary AF and DP for PR ROC?
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