Rna-seq 去除批次效应ComBat-seq
ComBat-seq
library(sva)
df=read.table("all_gene_count_raw.xls",head=T,check.names=F,comment.char="")
head(df)
batch=c("3D","3D","3D","7D","7D","7D","3D","3D","3D","7D","7D","7D")
#group=c(0,0,0,0,0,0,1,1,1,1,1,1)
df[,-c(1,13)] <- ComBat_seq(df[,-c(1,13)], batch=batch, group=NULL)
write.table(df,file = "all_gene_count.xls",sep="\t",quote = F,row.names = F)
参考：https://github.com/zhangyuqing/ComBat-seq
- 发表于 2024-01-09 10:54
- 阅读 ( 1402 )
- 分类：转录组
