% of change and log2 fold change (L2FC)?
2
2
Entering edit mode
2.2 years ago
Wox ▴ 470

I am trying to understand Gene expression data analysis. Could someone the percentage of change of the GE calculating from log2 fold change (L2FC)?

FoldChange = μ_Treatment / μ_Control

log2FC=Log2(μ_Treatment)-Log2(μ_Control)

RNA-Seq gene • 2.2k views
ADD COMMENT
7
Entering edit mode
2.2 years ago

seems correct to me. log2FC is more informative as you will see the direction (sign) and the effect.

Treatment = 10
Control = 100


FoldChange = 10/100 = 0.1 which is not informative

log2FC = log2(10)-log2(100) = -3.3 which tells the effect and direction.

ADD COMMENT
0
Entering edit mode

Thank you @geek_y.

After getting log2FC; e.g. -3.3, how can we calculate the % of change? is that %change = (Change/Original) x 100 ?

ADD REPLY
2
Entering edit mode

If you do 2^log2FC, which tells how many folds the change is. Here 2^3.3=9.8 folds. The sign gives information about the up/down regulation.

ADD REPLY
3
Entering edit mode
2.2 years ago
Wox ▴ 470

calculate fold change (FC) When comparing these log transformed values, we use the quotient rule of logarithms:

log(A/B) = log(A) - log(B)

log(A) = 4

log(B) = 1

Therefore:

log(A/B) = 4 - 1

log(A/B) = 3 This gives a 3-fold change.

Please note that in this case we are reporting the log(fold change). Biologists often use the log(fold change) because without taking the log, down regulated genes would have values between 0 and 1, whereas up regulated genes would have any value between 1 and infinity. This distribution makes graphing and further statistical analysis difficult. Taking the log typically makes the resulting values more normally distributed, which is better for further analysis.

ADD COMMENT

Login before adding your answer.

Traffic: 999 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.