maxATAC
maxATAC copied to clipboard
Update threshold.py
Issue:
-
maxatac predictexpects the following column names in a cutoff file:- 'Monotonic_Avg_Precision', 'Monotonic_Avg_Recall', 'Standard_Thresh', 'Monotonic_Avg_log2FC', 'Avg_F1'.
-
However,
maxatac thresholdoutputs a file with the following column names:- 'Monotonic_Precision', 'Monotonic_Recall', 'Threshold', 'Monotonic_log2FC', 'Monotonic_F1'.
-
When these columns are not available,
maxatac predictwill not output a BED file of TFBS predictions, even when a cutoff file and type are specified by the user.
Fix:
- I changed the column names in line 143 to correspond to the expected column names for
maxatac predict.
I tested this change to ensure it works with v1.0.6 of maxATAC. A BED file of predictions is now outputted when a cutoff file and type are specified by the user.