run_BHq_log: Run Benjamini-Hochberg procedure with log correction

In svteichman/prelim.knockoffs: Implements the knockoff method by Barber & Candes

Description

Runs the adjusted Benjamini-Hochberg procedure with q/S(p) replacing q where S(p) = 1 + 1/2 + ... + 1/p.

Usage

run_BHq_log(X, y, fdr)

Arguments

X	A nxp design matrix.
y	A vector of n responses.
fdr	The desired false discovery rate.

Value

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.

svteichman/prelim.knockoffs documentation built on May 28, 2020, 5:14 p.m.