Why am I getting a warning in "mafdr" when using my list of p-values?

I am getting the warning below when trying to calculate FDR for some p-values on a list. Even if I try with a smaller subset of the p-values on my list I get this warning. Could some one please explain why the mafdr function is complaining?
>> pvalues = [0.8708
0.0709
0.3546
0.4540
0.1767
0.0333
0.9174
0.6118
0.2853
0.8558
0.8419
0.9736
0.5124
0.4884
0.5656
0.9730]
>> [fdr qValues] = mafdr(pvalues);
Warning: The estimated PI0 is greater than 1. Please check the p-values are valid
or try a different lambda method
PI0 is set to 1.
> In mafdr>bootstrapchooser at 316
In mafdr>storeyFDR at 220
In mafdr at 190
>> qValues
qValues =
0.9736
0.5668
0.9736
0.9736
0.9426
0.5333
0.9736
0.9736
0.9736
0.9736
0.9736
0.9736
0.9736
0.9736
0.9736
0.9736
PS I am new to this kind of analysis and MATLAB.

 Accepted Answer

Wayne King
Wayne King on 30 Jul 2012
Edited: Wayne King on 30 Jul 2012
Hi, that is apparently a bug in the bootstrapchooser helper function. Do you have access to the R2012b prerelease version of MATLAB?
That behavior is fixed in the R2012b prerelease version.

2 Comments

Thank you very much for the fast answer. No, I don't have access to that version. Is it possible to replace the function file with the one in R2012b to get rid of the bug?
I don't know what other dependencies might exist, but you would have to have access to the R2012b version of mafdr.m -- which would mean you have access to the prerelease.

Sign in to comment.

More Answers (1)

Does this work for you?
[fdr,qvalues] = mafdr(pvalues,'lambda',0.15);

Categories

Asked:

JMA
on 30 Jul 2012

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!