https://github.com/tansey/smoothfdr
Tip revision: c5b693d0a66e83c9387433b33c0eab481bd4a763 authored by Wesley Tansey on 08 May 2020, 15:42:20 UTC
Fixed bug in easy that created too large a support for the alternative distribution
Fixed bug in easy that created too large a support for the alternative distribution
Tip revision: c5b693d
oracle_run.py
import sys
import numpy as np
from smoothfdr.utils import calc_fdr
posteriors = np.loadtxt(sys.argv[1] + 'oracle_posteriors.csv', delimiter=',')
discoveries = calc_fdr(posteriors, 0.1)
np.savetxt(sys.argv[1] + 'oracle_discoveries.csv', discoveries, delimiter=',', fmt='%d')