如何在R中加载.h5mu数据创建Seurat对象？及SeuratDisk安装问题

一、SeuratDisk包的安装 #

SeuratDisk未收录在CRAN官方源中，因此无法通过install.packages()直接安装，需从GitHub仓库安装：

• 先安装依赖的devtools包（若未安装）：

install.packages("devtools")
library(devtools)

• 安装SeuratDisk：

install_github("mojaveazure/seurat-disk")

注意：安装前确保R版本≥4.0，若提示缺少依赖（如Seurat、hdf5r），直接用install.packages()安装对应依赖即可。

二、.h5mu数据加载到R的方案 #

.h5mu是MuData格式，用于存储多模态单细胞数据（包含多个AnnData子集），以下是几种可行方法：

方法1：通过SeuratDisk转换加载 #

利用SeuratDisk的Convert函数将.h5mu转为Seurat支持的格式，或直接读取指定模态：

library(SeuratDisk)
# 将.h5mu转换为.h5seurat格式后加载
Convert("your_data.h5mu", dest = "h5seurat")
seurat_obj <- LoadH5Seurat("your_data.h5seurat")

# 直接读取.h5mu中的单个模态（如RNA）
rna_anndata <- ReadH5AD("your_data.h5mu", name = "rna")

方法2：使用muonR直接读取MuData #

muonR是专门处理MuData的R包，支持直接加载.h5mu：

• 安装muonR：

install_github("PMBio/muonR")

• 加载数据：

library(muonR)
mdata <- read_mudata("your_data.h5mu")
# 查看所有模态
mdata$mod
# 提取指定模态（如蛋白质组数据）
prot_data <- mdata$prot

方法3：用hdf5r手动读取HDF5结构 #

若上述方法受限，可直接操作底层HDF5文件提取数据：

install.packages("hdf5r")
library(hdf5r)
# 打开.h5mu文件
h5_file <- H5File$new("your_data.h5mu", mode = "r")
# 查看文件内部结构
h5_file$ls(recursive = TRUE)
# 提取RNA模态的表达矩阵示例
rna_counts <- h5_file[["rna/X"]]$read()
# 操作完成后关闭文件
h5_file$close()

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