甲基化分析中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同理


免责声明!

本站转载的文章为个人学习借鉴使用,本站对版权不负任何法律责任。如果侵犯了您的隐私权益,请联系本站邮箱yoyou2525@163.com删除。



 
粤ICP备18138465号  © 2018-2025 CODEPRJ.COM