Seurat依舊是單細胞分析的利器，也有越來越的的單細胞分析工具基于python開發(fā)。在數據分析時，有時需要將seurat的結果遷移到python下游分析，手動遷移seurat的cluster，umap，pca等結果很麻煩，現在和大家分享一款數據格式轉換工具，幫助大家完成R與python多平臺的多工具的數據分析。
工具來自這篇文章——Comparison of visualization tools for single-cell RNAseq data
這里是github上? ? README.md
Installation
sceasy is installable either as a bioconda package:
conda install -c bioconda r-sceasy
or as an R package:
devtools::install_github("cellgeni/sceasy")
which will require the biconductor packages BiocManager and LoomExperiment:
if (!requireNamespace("BiocManager", quietly = TRUE))
? ? install.packages("BiocManager")
BiocManager::install(c("LoomExperiment", "SingleCellExperiment"))
To use sceasy ensure the anndata package is installed:
conda install anndata -c bioconda
Optionally, if you plan to convert between loom and anndata, please also ensure that the?loompy?package is installed:
conda install loompy? -c bioconda
You will also need to install reticulate package:
install.packages('reticulate')
Usage
Before converting your data please load the following libraries in your R session:
library(sceasy)
library(reticulate)
use_condaenv('EnvironmentName')
loompy <- reticulate::import('loompy')
Seurat to AnnData
sceasy::convertFormat(seurat_object, from="seurat", to="anndata",
? ? ? ? ? ? ? ? ? ? ? outFile='filename.h5ad')
AnnData to Seurat
sceasy::convertFormat(h5ad_file, from="anndata", to="seurat",
? ? ? ? ? ? ? ? ? ? ? outFile='filename.rds')
Seurat to SingleCellExperiment
sceasy::convertFormat(seurat_object, from="seurat", to="sce",
? ? ? ? ? ? ? ? ? ? ? outFile='filename.rds')
SingleCellExperiment to AnnData
sceasy::convertFormat(sce_object, from="sce", to="anndata",
? ? ? ? ? ? ? ? ? ? ? outFile='filename.h5ad')
SingleCellExperiment to Loom
sceasy::convertFormat(sce_object, from="sce", to="loom",
? ? ? ? ? ? ? ? ? ? ? outFile='filename.loom')
Loom to AnnData
sceasy::convertFormat('filename.loom', from="loom", to="anndata",
? ? ? ? ? ? ? ? ? ? ? outFile='filename.h5ad')
Loom to SingleCellExperiment
sceasy::convertFormat('filename.loom', from="loom", to="sce",
? ? ? ? ? ? ? ? ? ? ? outFile='filename.rds')
