RDS文件转成Scanpy输入

首先安装loomR包

# Install devtools from CRAN
install.packages("devtools")
# Use devtools to install hdf5r and loomR from GitHub
install.packages("hdf5r")
devtools::install_github(repo = "mojaveazure/loomR", ref = "develop")

首先安装loomR包

# 读入seurat处理后的rds文件
library(Seurat)
library(loomR)
library("SeuratDisk")
packageVersion("Seurat")
packageVersion("loomR")
packageVersion("hdf5r")
sdata <- readRDS(file = "/Users/liji/PycharmProjects/docker_R/total_24h.diff_PRO.rds")
# seurat对象转换为loop文件
sdata.loom <- as.loom(x = sdata, filename = "/Users/liji/PycharmProjects/docker_R/total_24h.diff_PRO.loom", verbose = FALSE)
# Always remember to close loom files when done
sdata.loom$close_all()