Seurat 对象添加任意低维嵌入数据

简介: 本文分享了一种将任意来源的低维嵌入数据添加到 Seurat reductions数据槽中的实现方法

Seurat自带的RunPCA (pca)和RunTSNE (tsne)等数据降维方法产生的低维嵌入信息会格式化为DimReduc类,存储在DefaultAssay(scRNA)的reductions数据槽。当使用这些自带的方法时,会自动创建DimReduc类并填充到数据槽。当然我们也可以将由其它方法产生的维度信息存储到reductions数据槽中。
reductions 数据槽

代码示例

创建一份自定义二维数据表

#random a cell position dataframe
mini.sc <- readRDS("seurat.rds")
cell_position <- data.frame(x=runif(n = ncol(mini.sc),min = 1,max = 2),y=runif(n = ncol(mini.sc),min = 1,max = 2))
rownames(cell_position ) <- colnames(mini.sc)
colnames(cell_position ) <- c('locus_1','locus_2')#注意维度名必须要以键名+数字的格式(key)
> cell_position
        locus_1  locus_2
cell.1 1.435976 1.189780
cell.2 1.497868 1.001837
cell.3 1.691928 1.877578
cell.4 1.760313 1.134111
cell.5 1.155401 1.022741
cell.6 1.849457 1.939137
cell.7 1.946818 1.292949
cell.8 1.588419 1.164327
cell.9 1.502251 1.399103

将cell_position插入到seurat对象的reduction数据槽中

mini.sc[["cell.locus"]] <- CreateDimReducObject(embeddings = as.matrix(cell_position),key = "locus_", assay = "RNA")
#通过seurat函数DimPlot绘图
DimPlot(mini.sc,reduction = "cell.locus",group.by = "class",raster=FALSE,pt.size = 8,shape.by = "class")+
    NoLegend()+
    scale_x_continuous(limits = c(1,2))+scale_y_continuous(limits = c(1,2))

目录
相关文章
|
6月前
|
存储 数据可视化 Python
单细胞RNA速率分析: scVelo入门教程
单细胞RNA速率分析: scVelo入门教程
单细胞RNA速率分析: scVelo入门教程
|
机器学习/深度学习 算法 网络架构
【YOLOv8改进 - Backbone主干】EfficientRep:一种旨在提高硬件效率的RepVGG风格卷积神经网络架构
【YOLOv8改进 - Backbone主干】EfficientRep:一种旨在提高硬件效率的RepVGG风格卷积神经网络架构
|
机器人
给 Mac 添加右键菜单「使用 VSCode 打开」
如何在 Mac 下右键文件或文件夹,直接通过菜单项「用 VSCode 打开」。
641 2
|
机器学习/深度学习 自然语言处理 算法
机器学习和深度学习之间的区别
机器学习和深度学习在实际应用中各有优势和局限性。机器学习适用于一些数据量较小、问题相对简单、对模型解释性要求较高的场景;而深度学习则在处理大规模、复杂的数据和任务时表现出色,但需要更多的计算资源和数据,并且模型的解释性较差。在实际应用中,需要根据具体的问题和需求,结合两者的优势,选择合适的方法来解决问题。
581 0
|
传感器 Web App开发 编解码
基于51单片机的智能热水器设计
基于51单片机的智能热水器设计
298 0
|
数据可视化 关系型数据库 数据挖掘
scRNA分析|一(尽)文(力)解决你的单细胞火山图问题
scRNA分析|一(尽)文(力)解决你的单细胞火山图问题
1484 0
|
Linux
edac-utils工具如何使用
【6月更文挑战第1天】edac-utils工具如何使用
903 3
|
移动开发 Linux
Linux 下基于R语言打开hdf5(.h5)文件异常的处理方法
本文分享了如何在 Linux 中关闭 hdf5 文件锁定,让rhdf5 能够调用 HDF5库解析 h5文件的处理方法,以供参考
1208 0
|
IDE Linux 开发工具
Linux本地RStudio工具安装指南及远程访问配置安装RStudio Server
Linux本地RStudio工具安装指南及远程访问配置安装RStudio Server
|
数据挖掘 索引
单细胞不同样本数据整合-解决AnnData合并时ValueError: cannot reindex from a duplicate axis问题
单细胞不同样本数据整合-解决AnnData合并时ValueError: cannot reindex from a duplicate axis问题