The meat and potatoes of this function's code was extracted from limma::kegga, originally written by Gordon Smyth and Yifang Hu.

plot_ora_bias(x, selected, feature.bias, title = "DE status vs bias", ...)