Apply a batch effect with "combatseq" in wgcna
0
0
Entering edit mode
2.6 years ago
Maryam • 0

I have 336 samples, with two conditions of fibrosis and normal , and male and female genders , and Chinese, Indian, Malay , Caucasian , Bruneian, and not reported nationalities, I want to apply the batch effect by combat to this study before doing WGCNA. and I use DESeq2 for normalization data

How do I use the following function? I'm very confused to use scripts and I can't set scripts for my data For example, what should I define as a batch? and how set counts in this script

ComBat_seq(
    counts,
    batch,
    group = NULL,
    covar_mod = NULL,
    full_mod = TRUE,
    shrink = FALSE,
    shrink.disp = FALSE,
    gene.subset.n = NULL
 )

 count_matrix <- matrix(rnbinom(400, size=10, prob=0.1), nrow=50, ncol=8)
 batch <- c(rep(1, 4), rep(2, 4))

 adjusted <- ComBat_seq(count_matrix, batch=batch, group=NULL)

Thank you for your help

deseq2 R batch combatseq wgcna • 1.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 1974 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6