FDR

请问在多重比较中我想用Benjamimi-Hochberg法计算FDR,那原始的p值如何快速计算,FDR求得后如何用于两两比较。想知道这些问题的R语言实现