Question about the usage of Package, qvalue
1
1
Entering edit mode
manozawa ▴ 10
@manozawa-22650
Last seen 21 months ago

Hi,

I am a very beginner of R and Bioconductor. I would like to calculate q-values from a list of p-values.

## So, following the instruction on https://www.bioconductor.org/packages/release/bioc/html/qvalue.html, I installed the package 'qvalue' to RStudio Version 1.1.463 as follows.

if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages("BiocManager")

## tmp <- qvalue::qvalue(test2)

Then, the results somehow gives a lot of q-values that are smaller than the corresponding p-values. In addition, many q-values show exactly the same q-value, 4.670181e-02.

You can download the input list of p-values (removed). Sorry for the very long one (4571 p-values).

I would greatly appreciate if you have any suggestions to solve this issue.

Best regards, Fumi Nozawa

1
Entering edit mode

0
Entering edit mode

You could consider the more conservative Benjamini-Hochberg method of FDR adjustment. This will never produce smaller adjusted p-values than the input. (It will still sometimes produce identical adjusted p-values even if all the input values are distinct, this is expected behaviour.)

p.adjust(test2, method="BH")

2
Entering edit mode
@gordon-smyth
Last seen 8 hours ago
WEHI, Melbourne, Australia

It is known that qvalue can give q-values smaller than the corresponding p-values and there have been several posts about this topic over the years, for example:

The following might also be of interest: https://support.bioconductor.org/p/49864/

0
Entering edit mode

Thank you very much for your information! That helps me a lot to understand the concept of p-value and q-value. I may ask further questions in the near future.