1.Python读取.h5ad文件(比R方便)

import anndata
import pandas as pd
adata=anndata.read("/home/R/R_data/Seurat/PBMC10/output/adata.h5ad")
#adata.X.todense()#将稀疏矩阵转成普通矩阵
#X=pd.DataFrame(adata.X.todense())
#cell_name=adata.obs.index
#chr_name=adata.var.index
#X.index=cell_name
#X.columns=chr_name
#X=X.T #行为peak,列为cell
#X.to_csv('/home/R/R_data/Seurat/PBMC10/output/SCALE_ATAC.tsv',sep='\t')

2.R读取.h5ad文件(没有python方便 略写)

library(reticulate)
reticulate::py_config()
Sys.which('python')  # 该python 下要安装了anndata
SCALE_atac <- "output/adata.h5ad"
ad <- import("anndata", convert = FALSE)
final_ad <- ad$read_h5ad(SCALE_atac)#获取anndata数据成功
Logo

为开发者提供学习成长、分享交流、生态实践、资源工具等服务,帮助开发者快速成长。

更多推荐