I am using the Bioconductor package gage
in combination with the gene sets from package msigdbr
to perform gene set enrichment analysis. The result object (I called it res
) of gage()
contains the data frames res$greater
and res$less
containing up or downregulated pathways. Are the gene sets up/downregulated in the experimental group (parameter samp
in function gage()
) or the control group (parameter ref
in function gage()
)? I am pretty sure the data frames are related to experimental group but a gene set which I suspected to be significant in greater
is significant in less
.
Thanks for your help.