甲基化芯片Beta值意义详解,以及minfi包使用
illumina甲基化芯片計(jì)算
illumina的軟件(Illumina GenomeStudio軟件)自動(dòng)就把450k芯片的IDAT原始文件轉(zhuǎn)換成了甲基化信號(hào)文件,一般是3~8列值每個(gè)樣本。
使用Illumina的Bead Studio軟件和甲基化模塊v3.2計(jì)算每個(gè)CpG靶標(biāo)的β值。
β值=來自甲基化珠粒類型的強(qiáng)度值/(來自甲基化的強(qiáng)度值+來自未甲基化珠粒類型的強(qiáng)度值+ 100)。
然后使用p值對(duì)數(shù)據(jù)進(jìn)行質(zhì)量過濾。
p值大于0.01的β值被認(rèn)為低于最小強(qiáng)度,閾值顯示為“NA”。
minfi包計(jì)算
minfi包的一個(gè)函數(shù)read.450k.exp也可以直接讀取.IDAT文件
平均β=信號(hào)B /(信號(hào)A +信號(hào)B + 100)
通過計(jì)算甲基化(信號(hào)A)和未甲基化(信號(hào)B)等位基因之間的強(qiáng)度比來確定DNA甲基化水平(β值)。
具體地,β值是由甲基化(M對(duì)應(yīng)于信號(hào)A)和未甲基化(U對(duì)應(yīng)于信號(hào)B)等位基因的強(qiáng)度計(jì)算的,熒光信號(hào)的比率β= Max(M,0)/ [Max( M,0)+ Max(U,0)+ 100]。
因此,β值的范圍從0(完全未甲基化)到1(完全甲基化)
具體的β值的意義是:
- 任何等于或大于0.6的β值都被認(rèn)為是完全甲基化的。
- 任何等于或小于0.2的β值被認(rèn)為是完全未甲基化的。
- β值在0.2和0.6之間被認(rèn)為是部分甲基化的。
總結(jié)
以上是生活随笔為你收集整理的甲基化芯片Beta值意义详解,以及minfi包使用的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: selenium结合sikuliX操作F
- 下一篇: DIY NAS服务器之OMV 5.6入坑