ChAMP包学习(1)

idat数据类型
该包包含两个测试数据集,一个是HumanMethylation450 data (.idat),另一个是模拟的EPIC data,可以用来测试ChAMP中可用的函数。这可以通过将目录指向testDataSet来加载,如下所示:

library("ChAMP")
testDir=system.file("extdata",package="ChAMPdata")
myLoad <- champ.load(testDir,arraytype="450K")

450k肺肿瘤数据集仅包含8个样本,4个肺肿瘤样本(T)和4个对照样本(C),我们使用这个数据集来展示ChAMP的功能。


对于EPIC仿真数据集,我们可以使用以下代码进行加载:

data(EPICSimData)

该仿真数据集包含16个样本,它们实际上都来自一个样本,但都被修改为DMP和DMR,并且单独存在一些误差方差。在这16个样本数据集中,我们模拟了8个样本作为对照,8个样本作为case,样本在EPICSimData对象的pd中标记。在此数据中,我们从bumphunter package的clusterMaker()函数中随机选取5000个区域。在每个区域,我们随机选择一些连续的CpGs作为DMR,然后增加或减少该DMR的beta值。因此,该数据集中的DMR(4700+)应该小于5000,因为一些模拟DMR只包含1-2个CpGs,所以在champ.DMR()函数中不被认为是DMR。用户可以使用这些数据测试ChAMP在EPIC数据上的功能。


450k肺肿瘤数据集仅包含8个样本,4个肺肿瘤样本(T)和4个对照样本(C),接下来,我们使用450K这个数据集来展示ChAMP的功能。

ChAMP包用于分析Illumina甲基化beadarray数据(EPIC和450k),并提供了一个集成分析450k和EPIC分析方法。这包括各种不同的数据导入方法(例如从.idat文件或beta值矩阵导入)和质量控制图。
Type-2 probe correction methods include SWAN1, Peak Based Correction (PBC)2 and BMIQ3 (the default choice).
minfi包也提供了流行的功能归一化方法。
奇异值分解(SVD)方法允许处理批次效应,纠正多个批次处理效果,采用Combat包实现。采用RefbaseEWAS9可以校正细胞类型的异质性。
ChAMP还包括一个函数,该函数可以从450k或EPIC数据中推断拷贝数变异。

为了识别差异甲基化区域(DMR), ChAMP除了提供Bumphunter和DMRcate这两种检测差异甲基化的方法,还提供了新的Probe Lasso Method。

对于那些需要寻找差异甲基化区块( Differentially Methylated Blocks)的人,新版ChAMP包含了一个检测这些区块的功能。

基因集富集分析(GSEA)也是可能的,ChAMP的新版本结合了一些方法,可以纠正由于基因中探针表示不均匀而造成的偏差。

此外,新版ChAMP还整合了FEM 包,该包可以在用户指定的基因网络(表现出不同表型间甲基化的差异)中推断基因模块。

虽然还有一些管道和包可以用于450k或EPIC阵列分析(比如IMA、minfi、methylumi、RnBeads和wateRmelon)。ChAMP提供了一个更加全面完整的分析管道,从读取原始数据文件到最终的三级分析结果,例如GSEA,它为研究人员简化了甲基化阵列分析。新版ChAMP还提供了一系列闪亮的、基于程序的web浏览器交互分析函数(GUI函数),以帮助科学家查看ChAMP的结果。这需要一个基于web浏览器的交互式框架,用于本地或远程调用图形系统。

安装

if (!requireNamespace("BiocManager", quietly=TRUE))
    install.packages("BiocManager")
BiocManager::install("ChAMP")

Test Data

该包包含两个测试数据集,一个是HumanMethylation450 data (.idat),另一个是模拟EPIC data,可以用来测试ChAMP中可用的函数。可以像下面这样将目录指向testDataSet来加载它

library("ChAMP")
testDir=system.file("extdata",package="ChAMPdata")
myLoad <- champ.load(testDir,arraytype="450K")
image.png
image.png
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,445评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,889评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,047评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,760评论 1 276
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,745评论 5 367
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,638评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,011评论 3 398
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,669评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,923评论 1 299
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,655评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,740评论 1 330
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,406评论 4 320
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,995评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,961评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,197评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,023评论 2 350
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,483评论 2 342

推荐阅读更多精彩内容