Hi friends,
I have a table of P-values like:
| Feature |
X |
type |
LDA |
P-value |
| Prevotella_copri |
5.17933651645 |
test |
4.60356342912 |
0.0922867588261 |
| Butyrivibrio_crossotus |
4.69683479993 |
control |
4.07248710826 |
0.735615882525 |
| Bacteroides_dorei |
4.60167732782 |
test |
3.96437649555 |
0.843672782494 |
| Akkermansia_muciniphila |
4.28836733193 |
control |
3.90356518613 |
2.98756798706E-05 |
| Bacteroides_cellulosilyticus |
4.27495706056 |
control |
3.89732714054 |
0.389612546234 |
| Bacteroides_intestinalis |
4.19463599386 |
control |
3.8851611951 |
0.0312193961172 |
| Bacteroides_eggerthii |
4.27603800675 |
test |
3.87752454514 |
0.086688300821 |
| Alistipes_putredinis |
4.6767463495 |
control |
3.86667535783 |
0.0419820771958 |
| Ruminococcus_bromii |
4.46456513989 |
test |
3.808952528 |
0.225163911724 |
| Bacteroides_uniformis |
4.85724578241 |
control |
3.80355808508 |
0.274845078066 |
| Eubacterium_sp_CAG_180 |
4.28515167506 |
test |
3.73437627148 |
0.0137754984604 |
| Bacteroides_vulgatus |
4.80795432151 |
control |
3.72498343619 |
0.978424756919 |
| Paraprevotella_xylaniphila |
4.1622177234 |
test |
3.71112597761 |
0.239742113249 |
| Eubacterium_rectale |
4.41796600305 |
test |
3.70791598112 |
0.306888374573 |
I want to calculate the adjusted P-value (False Discovery Rate, Benjamini-Hochberg) and make an additional column in the output. I also expect one additional file where only q-value <0.25 will be featured.
Can anyone please tell me how can I do that?
Thanks