The function fdp performs different multiple test procedures for
controlling the false discovery proportion (FDP).
function [O] = fdp(Input,n1,samp) returns the number of rejected hypotheses,
the rank (O(:,1)), the indices of the rejected hypotheses (O(:,2)) and the unadjusted
p-values (O(:,3)) for the procedure B (conservative) of Korn et al. (2004)
with the significance level alpha=0.05.
----------
INPUT
These input arguments are required:
Input: data matrix with the size [n,k]
n1: number of patients in group one (0 < n1 <= n ),
restricted by the kind of samp
samp: kind of sample
single sample 'single' (n1 = n)
paired sample 'paired' (n1 = n/2; n must be even)
independent sample 'indept' (n1 < n)
[...] = fdp(...,'PARAM1',VAL1,'PARAM2',VAL2,...) specifies additional
parameters and their values. Valid parameters are the following:
Parameter Value
'gamma' gamma = 0.1; (the default)
for an other value: 0< gamma <= 0.5
'B' number of permutations
default: 500
B must be in the intervall
500 <= B <= 2^n1 for single and paired sample
(for 2^n1 < 500 : B = min(B,2^n1))
500 <= B <= n! / n1!*(n-1) for independent sample
(for n! / n1!*(n-n1)! < 500 : B = min(B,n!/n1!*(n-n1)!))
'test' Value for single sample
'ttest' to compute the t-Test
assumption : normal(gaussian) distribution
'wilcox' to compute the Wilcoxen signed rank test
assumption : symmetrical distribution
'sign' (the default) to compute the sign-test
assumption : none
Value for paired sample
'ttest' to compute the t-Test
assumption : normal(gaussian) distribution
'wilcox' to compute the Wilcoxen signed rank test
assumption : symmetrical distribution
'sign' (the default) to compute the sign-test
assumption : none
Value for independent sample
'ttest' to compute the t-Test
assumption : normal(gaussian) distribution
'wilcox' (the default) to compute the Wilcoxen rank test (Wilcoxen-Man-Whitney-Test)
assumption : none
'tail' The alternative hypothesis against which to compute
p-values for testing the hypothesis of no differences.
Choices are:
tail Alternative Hypothesis
'~=' (the default) "there is a significant difference" (two-sided test)
'>' "the values of group 1 are higher than the values of group 2" (one-sided test)
'<' "the values of group 1 are smaller than the values of group 2" (one-sided test)
---
'proc' 'ProcBv' (the default) chooses the procedure B (conservative) of Korn et al. (2004)
'ProcBe' chooses the procedure B of Korn et al. (2004)
'TL' chooses the procedure of Troendle (1995) and the extention of van der Laan et al.
'LR1' chooses the procedure of Lehmann and Romano (2005) with some dependence
assumptions or asymtotic control (see Romano and Wolf
(2005) "Control of Generalized Error Rates in Multiple Tetsing")
'LR2' chooses the procedure of Lehmann and Romano (2005) without any dependence assumptions (conservative!)
'HL' chooses the procedure of Holm and the extention of van der Laan et al.
---
'alpha' 0.05 (the default) significance level
alpha must be a scalar and in the interval 0 < alpha <= 0.2
OUTPUT
[O] = fdp(Input,n1,samp) returns the rank (O(:,1)),
the indices of the rejected hypotheses (O(:,2)) and
the adjusted p-values (O(:,3)).
-----------
REFERENCES
[1] Hemmelmann, C., Horn, M., S�e, T., Vollandt, R., Weiss, S. (2005):
New concepts of multiple tests and their use for evaluating
high-dimensional EEG data, Vol 142/2 pp 209-217.
Copyright (C) 2006 by Claudia Hemmelmann <claudia.hemmelmann@mti.uni-jena.de>
Institute of Medical Statistics, Computer Sciences and Documantation
University of Jena
This work was supported by DFG Project VO 683/2-1
This is part of the BIOSIG-toolbox http://biosig.sf.net/
---
***
This library is free software; you can redistribute it and/or
modify it under the terms of the GNU Library General Public
License as published by the Free Software Foundation; either
Version 2 of the License, or (at your option) any later version.
This library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Library General Public License for more details.
You should have received a copy of the GNU Library General Public
License along with this library; if not, write to the
Free Software Foundation, Inc., 59 Temple Place - Suite 330,
Boston, MA 02111-1307, USA.
--------------------------------------------------------------------------