[R] p.adjust documentation suggestion: simple statement that new "BH" method is an alias for old "fdr" method

Earl F. Glynn efg at stowers-institute.org
Mon Jun 6 23:31:45 CEST 2005


In R 2.0.1

> p.adjust.methods
[1] "holm"       "hochberg"   "hommel"     "bonferroni" "fdr"        "none"

In R 2.1.0

> p.adjust.methods
[1] "holm"       "hochberg"   "hommel"     "bonferroni" "BH"         "BY"
"fdr"
[8] "none"

One  might conclude that two new methods "BH" and "BY" were added.  BUT,
there's a clue in one of the comments of the examples for ?p.adjust:

       ## or all of them at once (dropping the "fdr" alias):
          p.adjust.M <- p.adjust.methods[p.adjust.methods != "fdr"]

Apparently, the old "fdr" that meant "Benjamini & Hochberg" now shares an
alias with the new "BH", which is "Benjamini & Hochberg".

Wouldn't a simple statement in the online documentation be useful that
explained that "fdr" and "BH" are aliases?

Is "fdr" soon to be deprecated and eventually shouldn't be used at all and
that "BH" is the way of the future?

http://cran.r-project.org/src/base/NEWS mentions "p.adjust() has a new
method '"BY'" but was silent on the apparently new alias "BH".

efg




More information about the R-help mailing list