Valid Conjunction Inference with the Minimum Statistic

Quick links: Summary | Article | HBM Poster | Software | Related HBM Poster

Important! Please see below for information on bug fix to the SPM2 Conjunction modifications

Summary

Conjunction inference in SPM2 and SPM99 is based on the minimum statistic. The inference is based on the global null hypothesis that all effects conjoined are null. Thus a significant result only implies one or more of the effects are non-null.

To assert that all effects are non-zero, the null hypothesis that one or more of the effects are null must be tested. In the paper below we find the appropriate inference for this conjunction null based on the minimum statistic. It happens that the appropriate inference comes from simply comparing the minimum statistic to the usual (single image) thresholds, corrected or uncorrected; further, no assumption about independence is required, so effects like (A-C) and (B-C) can now be conjoined.

Put another way, valid conjunction inference can be performed by creating thresholded statistic images and then creating an intersection mask. Also equivalently, at any given voxel, the maximum P-value over a set of images to conjoin is exactly the conjunction P-value (corrected or uncorrected).

Below is an complete article, as submitted to NeuroImage (a slightly revised version is currently In Press, NeuroImage) and software modifications for SPM. Please contact me with any questions.

Article

Valid Conjunction Inference with the Minimum Statistic.
Thomas Nichols, Matthew Brett, Jesper Andersson, Tor Wager & Jean-Baptiste Poline.
NeuroImage, Volume 25, Issue 3, 15 April 2005, Pages 653-660
doi:10.1016/j.neuroimage.2004.12.005 Preprint

Abstract

In logic a conjunction is defined as an AND between truth statements. In neuroimaging, investigators may look for brain areas activated by task A AND by task B, or a conjunction of tasks (Price and Friston, 1997). Friston et al. 1999 introduced a minimum statistic test for conjunction. We refer to this method as the minimum statistic compared to the global null (MS/GN). The MS/GN is implemented in SPM2 and SPM99 software, and has been widely used as a test of conjunction. However, we assert that it does not have the correct null hypothesis for a test of logical AND, and further, this has led to confusion in the neuroimaging community.
In this paper, we define a conjunction and explain the problem with the MS/GN test as a conjunction method. We present a survey of recent practice in neuroimaging which reveals that the MS/GN test is very often misinterpreted as evidence of a logical AND. We show that a correct test for a logical AND requires that all the comparisons in the conjunction are individually significant. This result holds even if the comparisons are not independent. We suggest that the revised test proposed here is the appropriate means for conjunction inference in neuroimaging.

HBM 2004 Poster

When is a conjunction not a conjunction?
Matthew Brett, Thomas Nichols, Jesper Andersson, Tor Wager & Jean-Baptiste Poline.
Poster WE 137: Abstract | Poster

Software: Revised conjunction inference for SPM

To address the issues raised above, install the modifications below. These modifications produce a new question when performing a conjunction: Null hyp. to assess? Select 'Conjunction' to use the proposed fix, so that P-values are valid for rejecting the conjunction null hypothesis; or select 'Global' to use the old method, where P-values are valid for rejecting the global null of all effects being absent. (Thanks to Darren Gitelman & Su Watkins for reporting this problem and helping me debug it.)

The compressed tar file below contains the seven (7) m-files needed to implement these changes within SPM. The files will expand into a separate spm2_conj (or spm99_conj) directory. Include this directory before your spm2 directory, e.g.


	addpath '/my/path/spm2_conj'

SPM2 Modifications

Download: spm2_conj.tgz (21KB)

Files:

spm_getSPM.m
In this file the meaning of the parameter n has been changed. Previously n was simply the number of tests to conjoin. Now, abs(n) is the number of tests to conjoin, and n<0 flags the use of the global null; a positive n implies the default behavior, using the (new) conjunction null.

spm_uc_Bonf.m
spm_uc_FDR.m
spm_uc_RF.m
spm_uc_Bonf.m
spm_uc_FDR.m
spm_uc_RF.m
These functions have been changed to impliment the new defintion of the argument n.

Important Bug Fix!
On October 16, 2004, an important bug in the Conjunction modifications for SPM2 was found and fixed. The problem was only present for 'Global Null' tests, and resulted in FDR P-values being wrong (and greater than 1!).
On November 9, 2004, an another bug in the Conjunction modifications for SPM2 was found and fixed. Again, the problem was only present for 'Global Null' tests, but only had the impact of producing thresholds that were too high.
If you had previously downloaded the conjunction modifications, please grab the whole tar file again, or just the affected files (spm_getSPM.m for the Oct 16 fix, spm_uc_FDR.m for the Nov 9 fix.)

SPM99 Modifications

Imortant Note: These modifications require that you have installed the SPM99 FDR modifications. Put the spm99_conj directory ahead of your spm_fdr directory in your matlab path.

Download: spm99_conj.tgz (21KB)