[R] p.adjust problem

[m1ch43l at gmx-topmail.de]

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