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