-
Notifications
You must be signed in to change notification settings - Fork 64
/
spm_P_peakFDR.m
47 lines (42 loc) · 1.85 KB
/
spm_P_peakFDR.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
function [Q] = spm_P_peakFDR(Z,df,STAT,R,n,ui,Ps)
% Return the corrected peak FDR q-value
% FORMAT [Q] = spm_P_peakFDR(Z,df,STAT,R,n,ui,Ps)
%
% Z - height {minimum over n values}
% df - [df{interest} df{residuals}]
% STAT - Statistical field
% 'Z' - Gaussian field
% 'T' - T - field
% 'X' - Chi squared field
% 'F' - F - field
% R - RESEL Count {defining search volume}
% n - Conjunction number
% ui - feature-inducing threshold
% Ps - Vector of sorted (ascending) p-values
%
% Q - FDR q-value
%__________________________________________________________________________
%
% References
% J.R. Chumbley and K.J. Friston, "False discovery rate revisited: FDR and
% topological inference using Gaussian random fields". NeuroImage,
% 44(1):62-70, 2009.
%
% J.R. Chumbley, K.J. Worsley, G. Flandin and K.J. Friston, "Topological
% FDR for NeuroImaging". NeuroImage, 49(4):3057-3064, 2010.
%__________________________________________________________________________
% Copyright (C) 2009-2012 Wellcome Trust Centre for Neuroimaging
% Justin Chumbley & Guillaume Flandin
% $Id: spm_P_peakFDR.m 5160 2012-12-21 16:58:38Z guillaume $
% Expected Euler characteristic for level ui
%--------------------------------------------------------------------------
[P, p, Eu] = spm_P_RF(1, 0, ui, df, STAT, R, n);
% Expected Euler characteristic for level Z
%--------------------------------------------------------------------------
[P, p, Ez] = spm_P_RF(1, 0, Z, df, STAT, R, n);
% Uncorrected p-value for peaks using Random Field Theory
%--------------------------------------------------------------------------
Z = Ez / Eu;
% q value using the Benjamini & Hochberch False Discovery Rate procedure
%--------------------------------------------------------------------------
Q = spm_P_FDR(Z, df, 'P',n, Ps);