sva_combat校正批次效应

安装包

options(BioC_mirror="https://mirrors.ustc.edu.cn/bioc/") 
options("repos" = c(CRAN="https://mirrors.tuna.tsinghua.edu.cn/CRAN/")) 

options()$repos 
options()$BioC_mirror
BiocManager::install('bladderbatch') # 把数据准备好,这是包里内置的数据集
BiocManager::install('sva')
library(sva)
library(bladderbatch)
data(bladderdata)
dat <- bladderEset[1:50,]

使用内置数据

pheno = pData(dat)
edata = exprs(dat)
batch = pheno$batch

使用Hierarchical clustering的方法去看一下聚类的情况

dist_mat <- dist(t(edata))
clustering <- hclust(dist_mat, method = "complete")
plot(clustering, labels = pheno$batch)
plot(clustering, labels = pheno$cancer)
image.png

校正批次效应,model可以有也可以没有,如果有,也就是告诉combat,有些分组本来就有差别,不要给我矫枉过正!

mod = model.matrix(~as.factor(cancer), data=pheno)
combat_edata <- ComBat(dat = edata, batch = pheno$batch, mod = mod)

dist_mat_combat <- dist(t(combat_edata))
clustering_combat <- hclust(dist_mat_combat, method = "complete")
plot(clustering_combat, labels = pheno$batch)
plot(clustering_combat, labels = pheno$cancer)
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容