auc_mtx[regulon_name].hist(bins=bins)
if threshold is None:
threshold = _derive_threshold(auc_mtx, regulon_name)
ylim = plt.ylim()
plt.plot([threshold]*2, ylim, "r:")
plt.ylim(ylim)
plt.xlabel("AUC")
plt.ylabel("//")
After Change
:param bins: The number of bins to use in the AUC histogram.
:param threshold: The threshold to use for binarization. If None then this will be derived automatically.
if ax is None:
ax=plt.gca()
auc_mtx[regulon_name].hist(bins=bins,ax=ax)
if threshold is None:
threshold = _derive_threshold(auc_mtx, regulon_name)