Which one? Log2foldchange or t statistic
1
0
Entering edit mode
2.8 years ago
Fatemeh ▴ 40

Hello , When I prepare my gene_list for gseGo() I tried Log2foldChange or t statistic ,separately , to order the list. However, I got different results for each. Which one do you recommend?

I used this code after deseq2 to get the differentially expressed gene list:

res <- results(dds, contrast=c("group", "KO", "wt")).

Afterwards, I just created the list of genes using either "Lof2foldChange" or "t statistic" and the gene names as the row names.

Thanks for your help in advance.

gseGo • 1.6k views
ADD COMMENT
1
Entering edit mode

Definitely, you should not use logFC to rank the genes for gene set enrichment analysis. So the answer is 't-statistics' (or "Wald-statistics" in DESeq2 output ). If you want to use logFC consider this line of code to convert logFC to a metric that can take into account both statistical significance and effect size:

metric = -log10(pvalue)/sign(log2FC)
ADD REPLY
0
Entering edit mode

Thanks very much Hamid for your reply. It was a big help.

ADD REPLY
1
Entering edit mode
2.8 years ago
nn ▴ 30

https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-017-1674-0

I usually go with t statistics

ADD COMMENT
3
Entering edit mode

The t statistic does NOT blend significance and effect size. It's significance only.

ADD REPLY
0
Entering edit mode

You're right, I removed my mistake.

ADD REPLY
0
Entering edit mode

Thanks a lot for your answer and provided link.

ADD REPLY

Login before adding your answer.

Traffic: 2364 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