I used fry() in limma package and it gave me this as results
NGenes Direction PValue FDR PValue.Mixed FDR.Mixed KEGG_APOPTOSIS 81 Down 2.853559e-24 5.307619e-22 3.188755e-62 2.281186e-61 KEGG_CYTOKINE_CYTOKINE_RECEPTOR_INTERACTION 161 Down 1.037278e-23 9.646683e-22 9.176848e-246 5.689646e-244 KEGG_CHEMOKINE_SIGNALING_PATHWAY 150 Down 1.723760e-23 1.068731e-21 5.378129e-200 2.000664e-198 KEGG_NATURAL_KILLER_CELL_MEDIATED_CYTOTOXICITY 103 Down 6.604034e-22 3.070876e-20 1.150026e-118 2.139049e-117 KEGG_EPITHELIAL_CELL_SIGNALING_IN_HELICOBACTER_PYLORI_INFECTION 61 Down 3.996045e-20 1.486529e-18 2.675739e-18 7.899801e-18 KEGG_PHOSPHATIDYLINOSITOL_SIGNALING_SYSTEM 67 Down 1.318473e-19 4.087265e-18 6.360791e-58 4.381878e-57 KEGG_OLFACTORY_TRANSDUCTION 42 Down 5.682502e-17 1.509922e-15 2.971345e-15 8.373791e-15 KEGG_NON_HOMOLOGOUS_END_JOINING 13 Up 3.624166e-16 8.426186e-15 4.728688e-01 6.925480e-01 KEGG_RNA_DEGRADATION 56 Up 3.536264e-14 7.135116e-13 1.445358e-18 4.336074e-18 KEGG_JAK_STAT_SIGNALING_PATHWAY 102 Down 3.883424e-14 7.135116e-13 6.883971e-95 9.849374e-94 KEGG_LEUKOCYTE_TRANSENDOTHELIAL_MIGRATION 85 Down 4.337482e-14 7.135116e-13 2.992323e-87 3.478576e-86
And so on and so forth. My questions are how I can get the genes involved in the gene sets and how I can interpret the p-value. I want to validate the gene by qPCR but I couldn't find which ones are meaningful in the gene sets. And also there are a lot of p < 0.05 which kind of indicates that it is overly sensitive. Is there an adjusted p-value or would it make sense that I just do my own p-value adjustment according to the number of gene sets are tested?
Hi Gordon,
Thank you for suggesting the use of camera. I am quite new to the field and saw a review saying that self-contained gene set analysis is more specific and sensitive, so I just followed without thinking much. I will have a look at the paper and vignette again.
Thanks a lot for your help.
P.S. Sorry I just missed part of answers. I have deleted the follow up questions. Many thanks for your help.