Interoperability between single-cell object formats
In this vignette, we demonstrate the ability to convert between Seurat objects, SingleCellExperiment objects, and anndata objects.
# install scater https://bioconductor.org/packages/release/bioc/html/scater.html library(scater) # install loomR from GitHub using the remotes package remotes::install_github(repo = # 'mojaveazure/loomR', ref = 'develop') library(loomR) library(Seurat) library(patchwork)
SingleCellExperiment is a class for storing single-cell experiment data, created by Davide Risso, Aaron Lun, and Keegan Korthauer, and is used by many Bioconductor analysis packages. Here we demonstrate converting the Seurat object produced in our 3k PBMC tutorial to SingleCellExperiment for use with Davis McCarthy’s scater package.
# download from satija lab https://www.dropbox.com/s/kwd3kcxkmpzqg6w/pbmc3k_final.rds?dl=0 pbmc <- readRDS(file = "../data/pbmc3k_final.rds") pbmc.sce <- as.SingleCellExperiment(pbmc) p1 <- plotExpression(pbmc.sce, features = "MS4A1", x = "ident") + theme(axis.text.x = element_text(angle = 45, hjust = 1)) p2 <- plotPCA(pbmc.sce, colour_by = "ident") p1 + p2