Hi, When I use the following code I get two different results, even though n is length(p) by default.
--------------------------------------------------------- pvals <- c(5.722385e-02, NA, 1.018087e-01) p.adjust(pvals, method="bonferroni") p.adjust(pvals, method="bonferroni", n=length(pvals)) > p.adjust(pvals, method="bonferroni") [1] 0.1144477 NA 0.2036174 > p.adjust(pvals, method="bonferroni", n=length(pvals)) [1] 0.1716716 NA 0.3054261 --------------------------------------------------------- Does anybody know what the reason for the different results is and which version is correct? Thanks, Michael -- NEU: FreePhone - 0ct/min Handyspartarif mit Geld-zurück-Garantie! Jetzt informieren: http://www.gmx.net/de/go/freephone ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
