Learn R Programming

pwrFDR (version 3.2.4)

FDR Power

Description

Computing Average and TPX Power under various BHFDR type sequential procedures. All of these procedures involve control of some summary of the distribution of the FDP, e.g. the proportion of discoveries which are false in a given experiment. The most widely known of these, the BH-FDR procedure, controls the FDR which is the mean of the FDP. A lesser known procedure, due to Lehmann and Romano, controls the FDX, or probability that the FDP exceeds a user provided threshold. This is less conservative than FWE control procedures but much more conservative than the BH-FDR proceudre. This package and the references supporting it introduce a new procedure for controlling the FDX which we call the BH-FDX procedure. This procedure iteratively identifies, given alpha and lower threshold delta, an alpha* less than alpha at which BH-FDR guarantees FDX control. This uses asymptotic approximation and is only slightly more conservative than the BH-FDR procedure. Likewise, we can think of the power in multiple testing experiments in terms of a summary of the distribution of the True Positive Proportion (TPP), the portion of tests truly non-null distributed that are called significant. The package will compute power, sample size or any other missing parameter required for power defined as (i) the mean of the TPP (average power) or (ii) the probability that the TPP exceeds a given value, lambda, (TPX power) via asymptotic approximation. All supplied theoretical results are also obtainable via simulation. The suggested approach is to narrow in on a design via the theoretical approaches and then make final adjustments/verify the results by simulation. The theoretical results are described in Izmirlian, G (2020) Statistics and Probability letters, "", and an applied paper describing the methodology with a simulation study is in preparation. See citation("pwrFDR").

Copy Link

Version

Install

install.packages('pwrFDR')

Monthly Downloads

292

Version

3.2.4

License

GPL (>= 2)

Maintainer

Grant Izmirlian

Last Published

January 14th, 2025

Functions in pwrFDR (3.2.4)

cCDF.VoR

Computes the complimentary CDF for the false discovery proportion, V_m/R_m.
CDF.Pval.HA

CDF of p-values for test statistics distribted under HA.
detail

The detail extraction function for simulated power objects
dists

The Distribution family object
cCDF.Rom

Computes the complimentary CDF for the significant call proportion, R_m/m.
pwrFDR.grid

Evaluate pwrFDR on a grid.
CDF.Pval

CDF of pooled (H0 and HA) population p-values
cCDF.ToM

Computes the complimentary CDF for the true positive proportion, T_m/M_m.
controlFDP

Helper function for the BHFDX FDP control method
pwrFDR

Ensemble power or sample size under selected control of the FDP
logitInv

Computes the inverse logit transform
%over%

Division operator with divide by zero clobbering
paste

The paste operator
gentempfilenm

Generate a tempfile name
logit

Computes the logit transform
CDF.Pval.apsi.eq.u

Calculates the fixed point for the Romano procedure.
sd.rtm.VoR

Extractor function for asymptotic sd[V_m/R_m] under selected FDP control method
sd.rtm.Rom

Extractor function for asymptotic sd[R_m/m] under selected FDP control method
sd.rtm.ToM

Extractor function for asymptotic sd[T_m/M_m] under selected FDP control method
CDF.Pval.au.eq.u

Function which solves the implicit equation u = G( u alpha)
cc.ROC

Computes the optimal number of controls per case in hypothesis tests involving the ROC. Included here with the intent that it can be used in conjunction with pwrFDR to allow power/sample size calculation for multiple tests of ROC curve based hypothesis. See details.
if.0.rm

A helper function-- remove if zero.
criterion

BH-FDR and Romano Criterion
arg.vals

Extracts the full argument list and call attribute.
es.ROC

Computes the equivalent Z-test effect size in hypothesis tests involving the ROC. Included here with the intent that it can be used in conjunction with pwrFDR to allow power/sample size calculation for multiple tests of ROC curve based hypothesis. See details.
backsolve.seFDPoalpha

Find missing argument giving required se[FDP]/alpha (or se[TPP]/average.power)
basic.tmPrint

Wrapper to Print a Basic Nicely Formatted Table
if.na.x

A helper function -- substitute 'NA's with a specified 'x'.
if.y.z

A helper function -- substitute y's with a specified 'z'.
join.tbl

Combine pwrFDR Results
nna

A helper function-- turns a missing column into 'NA's inside of a with statement