seurat icon indicating copy to clipboard operation
seurat copied to clipboard

Batch-Corrected Counts from RPCA with log-normalized counts

Open andrewdchen opened this issue 11 months ago • 1 comments

Similar question to the one raised in 8146 but for log normalized data.

But running integration with IntegrateLayers as demonstrated in the vignette, doesn't return a batch-corrected counts layer if you already log-normalized your data. I'm running code like the below:

obj <- IntegrateLayers(
  object = obj, method = RPCAIntegration,
  orig.reduction = "pca", new.reduction = "integrated.rpca",
  verbose = FALSE
)

and there's a layer called 'integrated.rpca' with a (Nx50) (50 for the #PCA dims in my 'pca' reduction). But there's no corrected counts that is (Nxp) (p being the number of genes), that would have been available if I used the v3 function IntegrateData.

Is there a way to use the v5 function whilst still outputting batch corrected counts?

andrewdchen avatar Feb 28 '24 20:02 andrewdchen