run_BHq: Run Benjamini-Hochberg procedure
In svteichman/prelim.knockoffs: Implements the knockoff method by Barber & Candes

Runs the either the original Benjamini-Hochberg procedure, the BH procedure with a log correction, or the BH procedure with whitened noise.

1	run_BHq(X, y, fdr, method = c("base", "log", "white_noise"))

`X`	A nxp design matrix.
`y`	A vector of n responses.
`fdr`	The desired false discovery rate.
`method`	either "base", "log", or "white_noise" (default: "base").

A list containing:

`X`	n-by-p design matrix.
`y`	vector of observed responses.
`Z`	vector of test statistics.
`thres`	data-dependent threshold.
`selected`	list of selected variables.