mFDP.adjust {mFDP}

compute mFDP-adjusted p-values

Description

Provides mFDP-adjusted p-values, given a vector of p-values.

Usage

mFDP.adjust(P, c="1/(2m)", s1=0, s2=0.1)

Arguments

P

A vector of (raw, i.e. unadjusted) p-values.

c

The tuning parameter, which influences the intercept and slope of the envelope. Should either be numeric or "1/(2m)" or "1/m".

s1

The smallest p-value threshold of interest.

s2

The largest p-value threshold of interest.

Value

A vector of mFDP-adjusted p-values. Some can be infinity - which can be interpreted as 1.

Examples

set.seed(5193)

### make some p-values
m=500        #the nr of hypotheses
nrfalse=100  #the nr of false hypotheses

tstats = rnorm(n=m)   #m test statistics
tstats[1:nrfalse] = tstats[1:nrfalse] + 3    #add some signal
P = 1 - pnorm(tstats)                        #compute p-values

P.adjusted = mFDP.adjust(P=P)     #mFDP-adjusted p-values. Be careful with interpretation.

min(P.adjusted) #0.0208