甲基化分析中limma包中的logFC值是怎么来的

最近在做甲基化分析,看到logFC值,突然好奇这个值是怎么转化过来的。

测试了一下,limma的logFC值与glm函数的Estimate值是一样的。

以下是limma包的结果,可以看到 logFC值为-0.09222974:


                logFC    AveExpr          t      P.Value    adj.P.Val        B
cg123       -0.09222974 0.47475441  -9.964224 1.504900e-12 5.492876e-07 17.63665

以下是glm函数的结果:

Coefficients:
                         Estimate Std. Error t value Pr(>|t|)    
(Intercept)              0.527457   0.007100  74.285  < 2e-16 ***
pheno$group_listcontrol  -0.092230   0.009393 -9.819  3.26e-12 ***

logFC值用公式来表示就是mean(case)-mean(control)
case指的是case样本的甲基化值,control同理

posted @ 2020-09-02 17:01  橙子牛奶糖  阅读(3511)  评论(0编辑  收藏  举报