I want to perform p adjustment for Benjamini and Hochberg’s FDR using R. However when I use the build-in function the result is not what I would expect (almost exclusively the same values). However, when I extract the function (with minor modifications) and have it as its own it gives a different result. I am curious as to why that is.I want to perform p adjustment for Benjamini an