DESeq2 adjusted p-values of 0 when p-value > 0.05

0

DESeq2 assigns adjusted p-values of 0 in some cases when the p-value is above 0.05:

name    baseMean          log2FoldChange   lfcSE             stat              pvalue              padj  gene_description
BCL6B   1.16290763279826  4.8304810831086  2.84617260952113  1.69718486747764  0.0896617151432926  0     BCL6B_transcription_repressor
C4B     1.16290763279826  4.8304810831086  2.84617260952113  1.69718486747764  0.0896617151432926  0     complement_C4B_(Chido_blood_group)
GATD3B  1.16290763279826  4.8304810831086  2.84617260952113  1.69718486747764  0.0896617151432926  0     glutamine_amidotransferase_like_class_1_domain_containing_3B
IL12A   1.16290763279826  4.8304810831086  2.84617260952113  1.69718486747764  0.0896617151432926  0     interleukin_12A
KCNJ13  1.16290763279826  4.8304810831086  2.84617260952113  1.69718486747764  0.0896617151432926  0     potassium_inwardly_rectifying_channel_subfamily_J_member_13
MRC1    1.16290763279826  4.8304810831086  2.84617260952113  1.69718486747764  0.0896617151432926  0     mannose_receptor_C-type_1
PNLDC1  1.16290763279826  4.8304810831086  2.84617260952113  1.69718486747764  0.0896617151432926  0     PARN_like,_ribonuclease_domain_containing_1
SGSM3   1.16290763279826  4.8304810831086  2.84617260952113  1.69718486747764  0.0896617151432926  0     small_G_protein_signaling_modulator_3
TRPV2   1.16290763279826  4.8304810831086  2.84617260952113  1.69718486747764  0.0896617151432926  0     transient_receptor_potential_cation_channel_subfamily_V_member_2

Coincidentally these genes aren't really expressed in my sample.

If I sort to filter based on both pvalue and padj then the results make sense.

Does anyone know why the weird behavior of padj being 0?

Thanks!


RNA-seq


DEseq2

• 198 views



Source link