Batch effect removal (Combat)
1
0
Entering edit mode
3.0 years ago
GSAENZDEPIP ▴ 20

Hi! I want to remove batch effect from a count data matrix (RNAseq experiment) using Combat function (sva package). I know that RemoveBatchEffect function (from Limma package) needs Log-transformed input, however Combat manual doesn't mention whether input data should be Log-transformed...

Thank you in advance!!

Combat manual: https://bioconductor.org/packages/release/bioc/manuals/sva/man/sva.pdf

RNA-Seq • 4.2k views
ADD COMMENT
1
Entering edit mode

Thank you! I have just found the answer in the tutorial: "For sequencing data, which are often represented as counts, a more suitable model may involve the use of a moderated log function"

ADD REPLY
1
Entering edit mode
3.0 years ago
h.mon 32k

The sva tutorial states one can extract the expression data from the expressionSet object:

The expression data can be obtained from the expression slot of the expression set.

edata = exprs(bladderEset)
ADD COMMENT

Login before adding your answer.

Traffic: 1079 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