Compute Q values {MixfMRI}

Q-values using Benjamini and Hochberg (1995)

Description

Compute q-values Benjamini and Hochberg's (1995) approach for controlling FDR.

Usage

  qvalue(p, method = c("BH1995", "BY2001"))

Arguments

p

a p-value vector.

method

using method by either BH1995 or BY2001

Details

This function compute q-values using Benjamini and Hochberg's (1995) approach for controlling FDR. The function bh.fdr is originally written by Dr. Dan Nettleton.

The Benjamini and Yeekutieli's (2001) approach for controlling FDR using the function by.fdr is coded by Wei-Chen Chen.

Value

Return corresponding q-values for the input p-values.

Author(s)

Dan Nettleton.

Modified by Wei-Chen Chen.

References

Chen, W.-C. and Maitra, R. (2021) “A Practical Model-based Segmentation Approach for Accurate Activation Detection in Single-Subject functional Magnetic Resonance Imaging Studies”, arXiv:2102.03639.

See Also

dpval(), dmixpval().

Examples

library(MixfMRI, quietly = TRUE)
set.seed(1234)
da <- gendataset(phantom = shepp1fMRI, overlap = 0.01)
p <- da$pval[!is.na(da$pval)][1:100]
qvalue(p)