最近在做甲基化分析,看到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同理

相关文章:

  • 2022-02-05
  • 2021-05-24
  • 2021-08-25
  • 2022-12-23
  • 2021-09-29
  • 2021-10-24
猜你喜欢
  • 2021-12-15
  • 2023-03-15
  • 2021-04-11
  • 2021-05-08
  • 2021-07-09
  • 2022-12-23
相关资源
相似解决方案