https://wiki.taichimd.us/index.php?action=history&feed=atom&title=File%3AFiltered_R.png File:Filtered R.png - Revision history 2026-04-09T09:32:22Z Revision history for this page on the wiki MediaWiki 1.43.6 https://wiki.taichimd.us/index.php?title=File:Filtered_R.png&diff=38977&oldid=prev Brb: <syntaxhighlight lang='r'> theta <- seq(0, .80, .01) R_BH <- filtered_R(alpha=.10, S2, p2, theta, method="BH") which.max(R_BH) # 60% <---- so theta=0.6 is the optimal filtering threshold # 61 max(R_BH) # [1] 380 plot(theta, R_BH, type="l", xlab=expression(theta), ylab="Rejections", main="BH cutoff = 0.1") abline(v=.6, lty=2) </syntaxhighlight> 2024-03-12T01:41:40Z <p>&lt;syntaxhighlight lang=&#039;r&#039;&gt; theta &lt;- seq(0, .80, .01) R_BH &lt;- filtered_R(alpha=.10, S2, p2, theta, method=&quot;BH&quot;) which.max(R_BH) # 60% &lt;---- so theta=0.6 is the optimal filtering threshold # 61 max(R_BH) # [1] 380 plot(theta, R_BH, type=&quot;l&quot;, xlab=expression(theta), ylab=&quot;Rejections&quot;, main=&quot;BH cutoff = 0.1&quot;) abline(v=.6, lty=2) &lt;/syntaxhighlight&gt;</p> <p><b>New page</b></p><div>== Summary ==<br /> &lt;syntaxhighlight lang=&#039;r&#039;&gt;<br /> theta &lt;- seq(0, .80, .01)<br /> R_BH &lt;- filtered_R(alpha=.10, S2, p2, theta, method=&quot;BH&quot;)<br /> which.max(R_BH)<br /> # 60% &lt;---- so theta=0.6 is the optimal filtering threshold<br /> # 61 <br /> max(R_BH)<br /> # [1] 380<br /> <br /> plot(theta, R_BH, type=&quot;l&quot;,<br /> xlab=expression(theta), ylab=&quot;Rejections&quot;,<br /> main=&quot;BH cutoff = 0.1&quot;)<br /> abline(v=.6, lty=2)<br /> &lt;/syntaxhighlight&gt;</div> Brb