File:Filtered R mean.png: Revision history

Jump to navigation Jump to search

Diff selection: Mark the radio buttons of the revisions to compare and hit enter or the button at the bottom.
Legend: (cur) = difference with latest revision, (prev) = difference with preceding revision, m = minor edit.

12 March 2024

  • curprev 20:1720:17, 12 March 2024Brb talk contribs 528 bytes +528 Use sample mean instead of variance for each gene as the filter statistic. <syntaxhighlight lang='r'> # Follow the previous code chunks M2 <- rowMeans(exprs(ALL_bcrneg)) theta <- seq(0, .80, .01) R_BH <- filtered_R(alpha=.10, M2, p2, theta, method="BH") which.max(R_BH) # 10% <---- so theta=0.1 is the optimal; only 10% genes are removed # 11 max(R_BH) # [1] 270 plot(theta, R_BH, type="l", xlab=expression(theta), ylab="Rejections", main="BH cutoff = 0.1") abline(v=.1, lty=2) <...