美文网首页
R调用reticulate处理数据

R调用reticulate处理数据

作者: 倪桦 | 来源:发表于2022-04-03 22:13 被阅读0次
library(reticulate)
######加载python模块######
sc <- import("scanpy")
pd <- import("pandas")

adata = sc$read('./matrix.h5ad') ###载入scanpy输出的h5ad文件

#######导出基因名和样本信息################
meta = adata$obs
gene <- adata$var

#############导出矩阵并转置,scanpy和Seurat的行列是反的#############
adata2 = adata$X
adata2 = adata2$T

adata2@Dimnames[[1]] = rownames(gene)
adata2@Dimnames[[2]] = rownames(meta)

merge <- CreateSeuratObject(adata2)
merge <- AddMetaData(merge, meta)

相关文章

网友评论

      本文标题:R调用reticulate处理数据

      本文链接:https://www.haomeiwen.com/subject/zidqsrtx.html