[BioC] removeBatchEffect options: design and covariates
Ryan C. Thompson
rct at thompsonclan.org
Wed Aug 6 00:42:10 CEST 2014
The covariates argument is just a more general way to specify batch
effects to be corrected for. If you had more than 2 batch factors, or
if you had some continuous numerical covariates to correct for, you
would manually construct your own "batch design matrix" and pass it as
the covariates argument. If all you have is one or two batch factors,
then all you need to do is pass them for batch and batch2, and you
don't need to worry about the covariates argument.
On Tue 05 Aug 2014 03:28:51 PM PDT, Rao,Xiayu wrote:
> Hi, Ryan
> That's very clear. Thank you very much for your instructions. Then when should I include covariates, or what kind of covariates I should adjusted for ?? Covariates in the model should already be in the design matrix, right?
> -----Original Message-----
> From: Ryan C. Thompson [mailto:rct at thompsonclan.org]
> Sent: Tuesday, August 05, 2014 5:18 PM
> To: Rao,Xiayu
> Cc: bioconductor at r-project.org
> Subject: Re: [BioC] removeBatchEffect options: design and covariates
> When calling removeBatchEffect, you should use the same design that you used for limma, but with with batch effect term removed from the design. Then you would pass the batch effect factor as the batch argument instead. So, if the design matrix that you used for limma was constructed as:
> model.matrix(~Condition + Batch),
> then for removeBatchEffect, you would use design=model.matrix(~Condition), and batch=Batch. In other words, you take the batch effect out of your model design and pass it as the batch argument instead.
> On Tue 05 Aug 2014 03:12:26 PM PDT, Rao,Xiayu wrote:
>> I want to use removeBatchEffect() on the expression data (Elist) prior to drawing a heatmap based on the expression of sig diff genes. Those sig diff genes were generated from limma linear modelling, with the batch factor already included in the linear model.
>> I saw people use removeBatchEffect(y, batch=batch) and removeBatchEffect(y, batch=batch, design=design). I would very much like to know in what condition I should include the design matrix, and when to also include covariates ??? Any comments would be very appreciated. Thank you in advance!
>> removeBatchEffect(x, batch=NULL, covariates=NULL,
>> design=matrix(1,ncol(x),1), ...)
>> [[alternative HTML version deleted]]
>> Bioconductor mailing list
>> Bioconductor at r-project.org
>> Search the archives:
More information about the Bioconductor