MixfMRI (version 0.1-0)

False Discovery Rates for Spatial Signals: False Discovery Rates for Spatial Signals using Benjamini and Heller (2007)

Description

Compute q-values Benjamini and Heller's (2007) approach for controlling FDR for spatial signals.

Usage

fdr.bh.p1(p, w = rep(1, length(p)), q = 0.05)
  fdr.bh.p2(p, w = rep(1, length(p)), q = 0.05)

Arguments

p

a p-value vector. No NA is allowed and all values are in [0, 1].

w

a weight vector for p-values.

q

a desired cutoff for adjusting p-values.

Value

Return the number of rejected hypotheses and all corresponding q-values for the input p-values.

Details

These functions implement first two procedures in Benjamini and Heller (2007) for controlling FDR for spatial signals.

References

http://maitra.public.iastate.edu/

See Also

qvalue().

Examples

# NOT RUN {
library(MixfMRI, quietly = TRUE)
set.seed(1234)
da <- gendataset(phantom = shepp1fMRI, overlap = 0.01)
p <- da$pval[!is.na(da$pval)][1:100]
fdr.bh.p1(p)
fdr.bh.p2(p)
# }