当我使用R的p.adjust函数来计算错误发现率时,我似乎得到了不一致的结果.基于所引用的纸张文档 经调整的p值应该这样来计算:
adjusted_p_at_index_i= p_at_index_i*(total_number_of_tests/i).
现在,当我跑步时,p.adjust(c(0.0001, 0.0004, 0.0019),"fdr")我得到了预期的结果
c(0.0003, 0.0006, 0.0019)
但是当我跑步时,p.adjust(c(0.517479039, 0.003657195, 0.006080152),"fdr")我得到了这个
c(0.517479039, 0.009120228, 0.009120228)
而不是我计算的结果:
c(0.517479039, 0.010971585, 0.009120228)
R对可以解释这两种结果的数据做了什么？