u3_fdr_ss: False Discovery Rates for Spatial Signals using Benjamini and...

Description Usage Arguments Details Value Author(s) References See Also Examples

Description

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

Usage

1
2
  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.

Details

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

Value

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

Author(s)

Wei-Chen Chen.

References

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

See Also

qvalue().

Examples

1
2
3
4
5
6
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)

MixfMRI documentation built on April 26, 2018, 5:03 p.m.