hanatFDR: Compute the False Discovery Rate for an anatomical hierarchy

In Mouse-Imaging-Centre/RMINC: Statistical Tools for Medical Imaging NetCDF (MINC) Files

Compute the False Discovery Rate for an anatomical hierarchy

Description

Compute the False Discovery Rate for an anatomical hierarchy

Usage

hanatFDR(buffer)

Arguments

buffer

What to compute FDR on

Details

Uses the Benjamini & Yekutieli (2001) algorithm for computing FDR to account for dependence (which is obviously present in a hierarchy where parent values are directly derived from their children).

Stores the values using the same name as the incoming statistic but with a q prepended. I.e. tvalue.Sexmale becomes qtvalue.Sexmale. Note that the original, input, is modified with these new additions as well.

Value

the anatomical hierarchy with FDR added (see details)

Examples

## Not run: 
vols <- addVolumesToHierarchy(hdefs, allvols)
hLm <- hanatLm(~Sex, gfBasic, vols)
hLm <- hanatFDR(hLm)
thresholds(hLm)

## End(Not run)

Mouse-Imaging-Centre/RMINC documentation built on Nov. 12, 2022, 1:50 p.m.