Jupyter Notebook加载AnnData对象时R内核崩溃求助

1. 检查H5AD文件完整性 #

• 用终端工具h5ls检查文件结构，确认无损坏或缺失关键组：

  h5ls ../adata/adata_celltypes_labelled_final.h5ad
  

  如果文件损坏，用Python重新导出修复：

  import anndata
  adata = anndata.read_h5ad("../adata/adata_celltypes_labelled_final.h5ad")
  adata.write("../adata/adata_repaired.h5ad", compression="gzip")
  

2. 调整readH5AD参数定位问题 #

• 关闭use_hdf5改用内存读取，同时开启verbose追踪崩溃节点：

  sce <- readH5AD(file, use_hdf5=FALSE, verbose=TRUE)
  

• 单独读取obs/var排查矩阵问题：

  sce_obs <- readH5AD(file, use_hdf5=TRUE, verbose=TRUE, only.obs=TRUE)
  sce_var <- readH5AD(file, use_hdf5=TRUE, verbose=TRUE, only.var=TRUE)
  

3. 调整R内存限制 #

• 针对Windows系统增大内存配额：

  memory.limit(size=64000) # 单位MB，按需调整
  

• Linux/macOS可在启动R时指定内存：

  R --max-mem-size=64G
  

4. 重装相关依赖包 #

• 即使之前环境正常，也可能存在包依赖损坏，重装rhdf5和SeuratDisk：

  remove.packages(c("rhdf5", "SeuratDisk"))
  BiocManager::install("rhdf5")
  install.packages("SeuratDisk")
  

5. Python中转转换格式 #

• 若上述方法无效，用Python读取后转存为R兼容的RDS格式：

  import anndata
  import rpy2.robjects as ro
  from rpy2.robjects import pandas2ri
  from rpy2.robjects.packages import importr
  
  pandas2ri.activate()
  base = importr('base')
  sce = importr('SingleCellExperiment')
  
  adata = anndata.read_h5ad("../adata/adata_celltypes_labelled_final.h5ad")
  ro.globalenv['counts'] = adata.X.todense()
  ro.globalenv['col_data'] = adata.obs
  ro.globalenv['row_data'] = adata.var
  
  ro.r('sce_obj <- SingleCellExperiment::SingleCellExperiment(assays=list(counts=counts), colData=col_data, rowData=row_data)')
  base.saveRDS(ro.r('sce_obj'), file="../adata/sce_converted.rds")
  

  再在R中读取：

  sce <- readRDS("../adata/sce_converted.rds")
  

内容的提问来源于stack exchange，提问作者Lauren