Hello,
I am trying to compare the results from 2 different sequencing machines. We run the same samples/libraries on both machines, analyzed using cellranger count and are just trying to compare how different the results are. I would like to do a co-embedded UMAP for both datasets. Would I use seurat to merge the 2 datasets and then use the following commands to produce the UMAP?
pbmc.rna <- NormalizeData(pbmc.rna)
pbmc.rna <- FindVariableFeatures(pbmc.rna)
pbmc.rna <- ScaleData(pbmc.rna)
pbmc.rna <- RunPCA(pbmc.rna)
pbmc.rna <- RunUMAP(pbmc.rna, dims = 1:30)
DimPlot(rna, reduction = "umap")
Is there a way to do this without having to go through most of the seurat pipeline?
Thank you
1 answer
Yes; if you want to visualize the impact that sequencing machine has (i.e., retain all batch effects) you can use Seurat::merge to directly combine the two objects prior to any analysis. Seurat will (if I recall) add a suffix on the redundant barcodes and treat them as different cells.
Log in to answer this question.
-
CCA integration
written by kayah •
-
scRNA integration problem with cluster not being separated
written by mnx0723 •
-
sci-RNA-seq
written by kilcdincer •
-
SIngle cell analysis
written by JY •
-
comparision of umap single cell
written by synat.keam •
-
UMAP graph using DimPlot pre/post integration
written by ralbero •
-
RunCCA() crash
written by Beth •
-
scATAC-seq and scRNA-seq integration issue
written by inayat45shaikh •
-
Cannot convert the inetgarted seurat object to the cds for monocle 3
written by shweta.johari13 •
-
How to add Seurat Object name in DimPlot title?
written by bk11 •
the following may be helpful: https://satijalab.org/seurat/archive/v3.1/merge_vignette.html