单细胞|M3-4. 细胞聚类与轨迹推断

embedded/2024/11/23 5:32:13/
library(ggplot2)
library(RColorBrewer)
library(scales)
# 聚类细胞
cds <- cluster_cells(cds)# 基于UMAP进行分区
cds@colData$part <- partitions(cds, reduction_method = "UMAP")# 绘制分区图,展示不同的细胞群体
plot_cells(cds, color_cells_by = "partition", label_cell_groups = FALSE)

数据结构角度讲解

要解读这段代码的数据结构和内容,需要先明确每个数据对象和函数的功能。以下是各部分的数据结构解读:


1. cds 对象

cds 是一个 CellDataSet 对象,来自 Monocle3。它是单细胞分析的核心数据结构,包含以下主要内容:

  • 表达数据:

    • 存储在 cds@assays$data,通常为基因表达矩阵(如 counts 数据)。
  • 细胞元数据:

    • 存储在 cds@colData,是一个 DataFrame,包含每个细胞的相关信息,比如细胞类型、分区信息等。
  • 基因元数据:

    • 存储在 cds@rowData,是一个 DataFrame,包含基因的注释信息(如基因名、基因 ID)。
  • 降维结果:

    • 通常存储在 cds@reducedDims 中,包含 UMAP 或 PCA 的降维坐标。
  • 轨迹图数据:

    • 存储在 cds@principal_graph,包含轨迹推断生成的细胞间连接信息。

2. cluster_cells(cds)

这一步执行了 细胞聚类,结果被存储到 cds@colData 中的新列:

  • 新增列名为 cluster 或类似名称。
  • 每个细胞被分配到一个聚类标签(如 1, 2, 3 等整数)。

3. partitions(cds, reduction_method = "UMAP")

这一步执行了基于 UMAP 的 细胞分区

  • 输入:

    • cds: 包含细胞降维信息的 CellDataSet。
    • reduction_method = "UMAP": 指定基于哪种降维方法(通常为 UMAP)。
  • 输出:

    • 返回一个向量,包含每个细胞的分区标签(通常是整数,如 1, 2, 3)。
  • 存储:

    • 分区结果被存储到 cds@colData$part,每个细胞对应一个分区标签。

4. plot_cells 可视化

功能:

plot_cells 函数用于绘制细胞分布图,支持多种颜色映射方式,比如按分区、聚类或拟时序。

color_cells_by = "partition":
  • 输入:
    • color_cells_by = "partition": 将细胞按分区标签着色。
  • 绘图内容:
    • 细胞按降维(UMAP)坐标分布。
    • 不同分区用不同颜色表示。
label_cell_groups = FALSE:
  • 不为细胞群体标注标签,仅展示颜色分区。
数据结构:

绘图用到的数据来源于 cds 的以下部分:

  • 降维坐标: 存储在 cds@reducedDims$UMAP
  • 分区标签: 存储在 cds@colData$part

总结数据流

  1. 输入数据: 包括表达矩阵(counts)、细胞元数据和基因元数据。
  2. 中间过程:
    • 聚类结果存储在 cds@colData$cluster
    • 分区结果存储在 cds@colData$part
  3. 可视化:
    • UMAP 降维结果来自 cds@reducedDims$UMAP
    • 分区标签用于按颜色分组。

通过这些步骤,可以有效地展示不同细胞群体的分布及其分区关系。


http://www.ppmy.cn/embedded/139791.html

相关文章

Vue3+SpringBoot3+Sa-Token+Redis+mysql8通用权限系统

sa-token支持分布式token 前后端代码&#xff0c;地球号: bright12389

服务机器人三甲坎德拉:用智能化开启售后服务新篇章

随着市场需求快速增长&#xff0c;专注自动驾驶服务机器人研发的坎德拉近年来实现了快速发展&#xff0c;并且通过技术创新与优质服务收获了良好口碑。为了不断提升售后服务质量和客户服务水平&#xff0c;坎德拉与新一代智能客户服务解决方案提供商售后宝共同打造智能化的售后…

基于知识图谱的群体事件知识问答系统

在信息爆炸的时代&#xff0c;如何快速而准确地获取知识成为我们面临的一大挑战。为此&#xff0c;我们推出了基于知识图谱的群体事件知识问答系统&#xff0c;为用户提供了一种全新的知识探索体验。 系统简介 我们的系统借助Neo4j图数据库的强大能力&#xff0c;结合Django后…

【网络云计算】2024第47周-每日【2024/11/21】周考-实操题-RAID6实操解析3

文章目录 一、标记故障磁盘二、移除故障磁盘三、添加备用磁盘四、监控数据重建过程五、验证RAID配置和数据完整性注意事项 使用 mdadm命令将故障磁盘替换为备用磁盘&#xff0c;并恢复数据的过程大致如下&#xff1a; 一、标记故障磁盘 检查RAID状态&#xff1a; 首先&#xf…

docker desktop上下载的镜像传到虚拟机上使用

在 Windows 上使用 Docker Desktop&#xff0c;下载 prom/node_exporter 镜像 C:\Windows\System32>docker images REPOSITORY TAG IMAGE ID CREATED SIZE prom/node-exporter latest 71dc9668b154 4 months ago 23.3MB步骤 1: 保存…

Dockerhub镜像加速

一、背景 dockerhub由于被封锁和站点处于国外的原因&#xff0c;docker pull拉取镜像非常慢&#xff0c;有时候直接都无法拉取。严重妨碍了我们的学习进度以及日常使用。 总结了一些proxy代理的镜像站点&#xff0c;配置之后速度会有明显提升&#xff0c;大家可以参考使用。 二…

小程序中模拟发信息输入框,让textarea可以设置最大宽以及根据输入的内容自动变高的方式

<textarea show-confirm-bar"{{false}}" value"{{item.aValue}}" maxlength"301" placeholder"请输入" auto-height"{{true}}" bind:blur"onBlurTextarea" focus"{{true}}" bindinput"…

汽车资讯新高度:Spring Boot技术飞跃

6系统测试 6.1概念和意义 测试的定义&#xff1a;程序测试是为了发现错误而执行程序的过程。测试(Testing)的任务与目的可以描述为&#xff1a; 目的&#xff1a;发现程序的错误&#xff1b; 任务&#xff1a;通过在计算机上执行程序&#xff0c;暴露程序中潜在的错误。 另一个…