Interpretation of "less" and "greater" after running gage with `same.dir = TRUE`
0
0
Entering edit mode
4.8 years ago
tmms ▴ 10

Hello

I am having troubles interpreting the results after running gage, with same.dir = TRUE, on a log2 fold changes that were obtained after an analysis with DESeq2. I am not sure what the meaning of up- and downregualted pathways is. Are this pathways where all the genes are up- or downregulated? Additional, up- or downregulated compared to what?

I know that the sign of the log2 fold changes reflect in which condition the expression is up- or downregulated compared to chosen control. In my case, a positive log2 fold change corresponds to a gene that is more expressed in diseased cells compared to healthy cells, while a negative log2 fold change corresponds to a gene that is less expressed in diseased cells compared to healthy cells.

Thanks in advance.

RNA-Seq gage functional_analysis pathways kegg • 1.1k views
ADD COMMENT
0
Entering edit mode

After digging somewhat deeper, I found out that same.dir = TRUE performs a one-side two-sample t-test. The difference between "greater" (upregulated?) and "less" (downregulated?) is that first one calculates a p-value by using the right tail of the distribution while the second uses the left tail. I verified this by summing the p-values of the same pathway of greater and less. They sum to 1.

I am still not sure how to link such a result to the initial experiment.
Let's say I want to compare two groups: disease and control. I analyze the data in such a way that control is my reference group. Therefore, a positive log fold change corresponds to a gene that more expressed in the disease compared to the control and a negative log fold change corresponds to a gene that is less expressed in the diseases compared to the control.
Let's say I get a significant hit on a pathway in the "greater" category after running gage. This means that the mean log fold change in the gene set is larger than the overall mean log fold change.
So does this mean that the gene set is upregulated in the disease compared to the control?

ADD REPLY

Login before adding your answer.

Traffic: 1830 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6