False Discovery Rates for Spatial Signals | R Documentation |
Compute q-values Benjamini and Heller's (2007) approach for controlling FDR for spatial signals.
fdr.bh.p1(p, w = rep(1, length(p)), q = 0.05)
fdr.bh.p2(p, w = rep(1, length(p)), q = 0.05)
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. |
These functions implement first two procedures in Benjamini and Heller (2007) for controlling FDR for spatial signals.
Return the number of rejected hypotheses and all corresponding q-values for the input p-values.
Wei-Chen Chen.
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.
qvalue()
.
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.