如何以行而非列的形式为Seurat对象添加MetaData？

问题场景 #

你用这段代码给Seurat对象加元数据：

Addmetadata <- AddMetaData(object=seurat,metadata = transposemeta,col.name = "cell_type__ontology_label")

结果报错：

Error in [[<-:
! None of the meta data requested was found in the data frame

你已经把元数据转置成了「细胞在列、属性在行」的结构，但AddMetaData()只认「细胞在行、属性在列」的格式，需要调整后再添加。

解决方法 #

单个属性添加

如果只加cell_type__ontology_label这一个属性，直接把该行转成带细胞名的向量即可：

# 提取目标属性对应的列（细胞为向量名）
cell_type_vec <- transposemeta["cell_type__ontology_label", ]
# 对齐Seurat对象的细胞顺序（避免匹配错误）
cell_type_vec <- cell_type_vec[colnames(seurat)]
# 添加到Seurat对象
seurat <- AddMetaData(object = seurat, metadata = cell_type_vec, col.name = "cell_type__ontology_label")

多个属性批量添加

如果要一次性加所有行里的属性，先把元数据转置回「细胞在行、属性在列」的格式，再传入函数：

# 转置元数据，让细胞成为行名
meta_df <- t(transposemeta)
# 对齐Seurat对象的细胞顺序
meta_df <- meta_df[colnames(seurat), ]
# 批量添加所有元数据列
seurat <- AddMetaData(object = seurat, metadata = meta_df)

原理说明 #

AddMetaData()要求传入的元数据满足：

• 数据框格式：行名是细胞ID，每一列对应一个属性
• 向量格式：向量的名字是细胞ID，值对应属性内容

你之前转置后的结构刚好反过来，所以函数找不到匹配的细胞ID，调整结构对齐后就能正常添加了。

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