Questions tagged [batch-effects]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
1 vote
1 answer
26 views

What is the best way to perform batch ajustment between two experiments?

I have duplicate drug efficacy experiments performed six months apart. I would like to know what is the right way to perform statistical analysis to understand fold changes after treatment with a drug....
user avatar
  • 11
2 votes
0 answers
11 views

Subsetting with harmony

Currently merged two Seurat objects together and then ran Harmony for batch correction. Now I want to subset out a cell type of interest do I re run harmony or just do the standard PCA? ...
user avatar
  • 317
0 votes
1 answer
29 views

What is the correct order of flooring-normalization-batch correction for microarrays?

This question was also asked on Biostars I am trying to learn and understand the correct order of data processing steps for microarrays. I have data which already was analyzed by a researcher using ...
user avatar
  • 167
0 votes
1 answer
42 views

Comparing differential expression across samples - is batch effect correction needed?

I have a bulk RNA-seq dataset made up of control and treatment conditions for a range of cell lines. This dataset was generated in two batches, such that the cell lines are split between batches but ...
user avatar
0 votes
1 answer
25 views

References request: deciding what duplicate samples to include to avoid batch effects with previous studies

Suppose there is cohort with many samples to assay. In phase I of the project, subjects with a certain disease incidence and healthy controls were assayed. Now, in phase II of the project, we are ...
user avatar
  • 101
1 vote
1 answer
40 views

(scRNA-seq) I have a dataset with 4 conditions (experiments) should I be anchoring aka integrating data for each comparison when looking for DEGs?

I have an experiment with 4 conditions in it two of them are controls. When I am correcting for batch effects (either with PCA or CCA) should I restrict my comparisons between particular groups? So ...
user avatar
1 vote
0 answers
20 views

Batch Effects that are Close to Collinear with Main Effect

Two years ago, I started working on a project uses both RNAseq and ATACseq. It's supposed to be a simple Healthy Control (HC) vs disease study. The sample collection was done in 2 phases, and it was ...
user avatar
  • 177
-1 votes
1 answer
38 views

Assumptions of batch effect removal

What does removing a batch effect (e.g. with limma::removeBatchEffect) assume about the batch effect? Does it assume simply a constant batch effect for each level? ...
user avatar
  • 167
0 votes
0 answers
74 views

Correcting for batch effect in bulk RNA seq datasets

I have samples from many different bulk RNA seq studies that exhibit an evident batch effect (samples from different labs cluster together regardless of treatment). I've found a lot of papers (here, ...
user avatar
2 votes
1 answer
92 views

Batch effect correction using a subset of samples - using DE genes

I have some RNA-seq data with two very obvious batches as you can see in the PCA: The samples of interest (A - H) are from tumor tissues. In addition, there are data for two cell lines I and J. The ...
user avatar
  • 788
0 votes
0 answers
79 views

How to use combat in order to remove batch effects?

I have RNA seq data and I need to use combat to remove the batch effects. Somehow when I run it, it isnt actually doing anything. The code: ...
user avatar
  • 1
2 votes
1 answer
38 views

Should one be especially conscious of removing low counts, when applying batch correction?

When (as a mistake) I did not remove low counts at all (beside those that equal zero for all samples), I got the following ma plot (using Glimma) : On the right you see the individual counts. The ...
user avatar
  • 167
0 votes
1 answer
530 views

Details of DESeq2 modeling a batch effect

When correcting my data for a batch effect using removeBatchEffect, some of the gene expression values become negative. When searching for differentially expressed genes, I do not use the data above, ...
user avatar
  • 167
3 votes
1 answer
971 views

Removing Batch Effect in Heatmaps after Differential Gene Expression Analysis

I'm working on a dataset in which the first replicate of each group is one batch and the second replicate is in a second batch. After checking the PCA plot and ...
user avatar
0 votes
2 answers
3k views

how to solve “dim(X) must have a positive length” at running ComBat function in R

I used ComBat() for batch effect correction in my expression data. basically, that function inputs are expression data, Batch covariate, and Model matrix for the ...
user avatar
2 votes
2 answers
1k views

TPM or rlog(CPM) for comparing expression?

I want to see the expression of a gene in a group of patient amongst the entire cohort using my RNA-Seq data. While I can do a differential expression analysis with limma or DESeq2, I want to see how ...
user avatar
  • 115
0 votes
2 answers
414 views

Interpreting this PCA plot for RNA-seq

I have RNA-seq from two sequencing batches; Lab technician says that he has run the RNA expression quantification two times in bathes 1 and 2 for example ...
user avatar
  • 1,513
1 vote
0 answers
137 views

ComBat batch correction: understanding the model

A collaborator of mine is using ComBat for some RNA-seq data. I would like to understand what it's doing, and I have a specific question about the structure of the model. Equation 2.1 of the paper ...
user avatar
2 votes
1 answer
112 views

Combat for multiplatform batch correction

I have been looking ways to use data from various micro-array platform such as agilient, rosetta/merk. Some of the established method I came across is combat, but I'm not sure if it can be done for ...
user avatar
  • 1,538
3 votes
1 answer
178 views

sva for RNA-Seq data without known phenotype

I have been working on RNA-Seq data from two different cohorts, and they show very strong batch effect (~35% variance explained by 1st component in PCA). Since I am trying to do a class discovery from ...
user avatar
  • 115
4 votes
1 answer
224 views

How should I address batch effects in my experiment?

Let's say I have an RNA-Seq experiment, where I'm interested in the significantly differentiated genes between pre-treatment and post-treatment conditions. "rep" == biological replicate. ...
user avatar
  • 2,639
3 votes
1 answer
984 views

How to get the corrected matrix after SVA batch effect correction

I ran SVA to remove batch effects for my bulk RNAseq experiments, but now I need to somehow correct my data matrix in order to ...
user avatar
0 votes
1 answer
1k views

Batch Effects in RNA Seq Sample

I am using the R (using EdgeR) for the RNA Seq analysis, I had few batch effect samples like Control vs treatment. Could anyone tell me the best way to remove the batch effects. I have looked into ...
user avatar
  • 13
5 votes
3 answers
3k views

Is it possible to merge scRNAseq data from experiments with different design?

I have 4 different single-cell RNAseq experiments, each one representing a different sample of cell type population. I'd like to merge them to a single dataset. However, different cell types are ...
user avatar
  • 1,713
2 votes
2 answers
2k views

ComBat for batch effects removal on scRNA-seq data

Is it possible to use sva's ComBat for batch effects removal on scRNA-seq data? Is there any difference between RNA-seq and scRNA-seq data which doesn't allow to use ComBat on single-cell data? (I am ...
user avatar
  • 1,713
2 votes
1 answer
112 views

How hard is it to clean and QC gene expression microarray data?

I am a PhD student working on developing ideas for my dissertation papers. One of my planned papers will be working with some (human) gene expression data. I have had a class on working with gene ...
user avatar
  • 195
6 votes
1 answer
1k views

Order of batch effects removal, data imputation and library size normalization in scRNA-seq data

I am preprocessing scRNA-seq data. What is the best practice in use to run both ComBat for batch effects removal, data imputation (to mitigate dropout) and library size normalization? I thought that ...
user avatar
  • 1,713