Learn R Programming

⚠️There's a newer version (3.2.4) of this package.Take me there.

pwrFDR (version 2.8.4)

FDR Power

Description

This is a package for calculating two differing notions of power, or deriving sample sizes for specified requisite power in multiple testing experiments under a variety of methods for control of the distribution of the False Discovery Proportion (FDP). More specifically, one can choose to control the FDP distribution according to control of its (i) mean, e.g. the usual BH-FDR procedure, or via the probability that it exceeds a given value, delta, via (ii) the Romano procedure, or via (iii) my procedure based upon asymptotic approximation. 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 package will compute power, sample size or any other missing parameter required for power based upon (i) the mean of the TPP which is the average power (ii) the probability that the TPP exceeds a given value, lambda, via my asymptotic approximation procedure. The theoretical results are described in Izmirlian, G. (2020) Stat. and Prob. letters, , and an applied paper describing the methodology with a simulation study is in preparation.

Copy Link

Version

Install

install.packages('pwrFDR')

Monthly Downloads

265

Version

2.8.4

License

GPL (>= 2)

Maintainer

Grant Izmirlian

Last Published

September 26th, 2020

Functions in pwrFDR (2.8.4)

cCDF.VoR

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

The Distribution family object
detail

The detail extraction function for simulated power objects
cCDF.Rom

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

Ensemble power or sample size under selected control of the FDP
CDF.Pval.HA

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

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

Evaluate pwrFDR on a grid.
paste

The paste operator
controlFDP

Helper function for the BHCLT FDP control method
CDF.Pval

CDF of pooled (H0 and HA) population p-values
gentempfilenm

Generate a tempfile name
if.y.z

A helper function -- substitute y's with a specified 'z'.
sd.rtm.ToM

Extractor function for asymptotic sd[T_m/M_m] under selected FDP control method
es.ROC

sd.rtm.Rom

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

sd.rtm.VoR

Extractor function for asymptotic sd[V_m/R_m] under selected FDP control method
logitInv

Computes the inverse logit transform
%over%

Division operator with divide by zero clobbering
CDF.Pval.ua.eq.u

Function which solves the implicit equation u = G( u alpha)
nna

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

A helper function -- substitute 'NA's with a specified 'x'.
CDF.Pval.ar.eq.u

Calculates the fixed point for the Romano procedure.
if.0.rm

A helper function-- remove if zero.
logit

Computes the logit transform
criterion

BH-FDR and Romano Criterion
arg.vals

Extracts the full argument list and call attribute.