<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><pre class="bz_comment_text" id="comment_text_0"><font face="Helvetica">Hi,
I realized that the fdr.m file gives a wrong output for the case that the smallest p values are above the FDR threshold. This means that the smallest p-values might be rejected, whereas lower p-values could pass the test. I have filed it as a bug (#3008).</font></pre><pre class="bz_comment_text" id="comment_text_0"><font face="Helvetica">The current script only declares the p-values below or equal to the threshold as significant, whereas in fact all p-values that are smaller or equal to the largest p-value below or equal to this threshold should be declared significant according to the original paper ("find the largest P value that is below the line [threshold]. All voxels with P values less than or equal to this are declared active.”).</font></pre><pre class="bz_comment_text" id="comment_text_0"><font face="Helvetica">I assume that the following line is incorrect and should be replaced by the one below:</font></pre><pre class="bz_comment_text" id="comment_text_0"><font face="Helvetica">Incorrect: h = (ps<=pi);</font></pre><pre class="bz_comment_text" id="comment_text_0"><font face="Helvetica">Correct: h = ps<=(max(ps(ps<=pi)));
Cheers,
Kathrin</font></pre><div>
<div>--<br><blockquote type="cite" style="color: rgb(0, 0, 0);"><blockquote type="cite"></blockquote></blockquote>Kathrin Müsch, Ph.D.<br><blockquote type="cite" style="color: rgb(0, 0, 0);"><blockquote type="cite"></blockquote></blockquote>Department of Psychology<br>University of Toronto<br>Toronto, Canada<br><a href="http://www.honeylab.org">www.honeylab.org</a></div>
</div>
<br></body></html>