- 综合排序
- 最热优先
- 最新优先
- 来自专栏生信技能树
拟时序的正确姿势（拟时序一本通02）
也就是说，我们做拟时序之前通常是要细分亚群到足够深入，需要确定被做拟时序分析的对象是有比较大生物学变化的可能性。 但是也有很多文章在使用拟时序分析的时候其实忽略这个前提，就显得很“可笑”，接下来我们就用3个案例来说明一下拟时序的正确姿势，错误示范，还有创新型的应用场景。 但是前面我们在 为什么做拟时序 提到了其实可以把拟时序分析简化成为了主要是为了展示差异细节，这个文章里面也是如此，就有了图(F)的monocle2的拟时序结果，因为我们关心的是GZMH高表达的效应（毒性 而且可以很清晰的看到正常的乳腺组织跟两个分子分型的肿瘤样品的差异细节： 差异细节 值得一提的是拟时序一本通专辑的目录是： 为什么做拟时序 （展示差异细节） 拟时序的正确姿势，错误示范，创新型拟时序 拟时序的多种算法大比拼 monocle2拟时序实战 降维聚类分群 拟时序 多种基础可视化 正向特殊可视化（目标基因表达量，基因集打分） 反方向特殊可视化（5种可视化） monocle3拟时序实战 SCORPIUS实战 Slingshot
51701编辑于 2024-03-25 - 来自专栏生信技能树
为什么做拟时序（拟时序一本通01）
所以如下所示，我们增加了两个有挑战性的分析，就是拟时序和转录因子分析。 上面的案例，其实就回答了我们今天推文的问题：为什么做拟时序，因为它已经是非常常规而且几乎是必备的单细胞分析条目，你很难跳过去。 而且，加入拟时序分析有时候确实是能更加说明问题。 尝试一下monocle2的拟时序分析 拟时序有非常多多种算法，后面我们再慢慢介绍，包括但不限于： monocle2拟时序实战 monocle3拟时序实战 SCORPIUS Slingshot 基于其它编程语言的拟时序分析 值得一提的是拟时序一本通专辑的目录是： 为什么做拟时序（展示差异细节） 拟时序的正确姿势，错误示范，创新型拟时序 拟时序的多种算法大比拼 monocle2拟时序实战 降维聚类分群 拟时序 多种基础可视化 正向特殊可视化（目标基因表达量，基因集打分） 反方向特殊可视化（5种可视化） monocle3拟时序实战 SCORPIUS实战 Slingshot实战 基于其它编程语言的拟时序分析 然而，本人的理解肯定是片面的
1K10编辑于 2024-03-25 - 来自专栏单细胞学习小组
day 8 拟时序分析
单样本输入数据输入数据是降维聚类分群注释的数据做拟时序分析通常不是拿全部的细胞，而是拿感兴趣的一部分。用subset提取子集即可。因为要使用差异基因来排序，所以要两类及以上细胞。 ()library(patchwork)p2+p1/p3#多样本的添加这个可以看去批次效应plot_cell_trajectory(sc_cds, color_by = 'orig.ident')经典拟时序热图这里图展示基因随着时间的渐变 )(100))基因轨迹图gs = head(gene_to_cluster)plot_cell_trajectory(sc_cds,markers=gs,use_color_gradient=T)基因拟时序点图
51810编辑于 2024-07-01 - 来自专栏生物信息学、python、R、linux
monocle 2拟时序分析
monocle做拟时序分析首先要构建CDS需要3个矩阵：expr.matrix、pd、fd，其次将Seurat中的对象转换为monocle识别的对象。 然后选择想要做拟时序依据的基因就可以了,如果已知开始和结束的细胞，将过程开始时收集的细胞与结束时收集的细胞简单地进行比较，并找到差异表达的基因，做拟时序依据的基因，根据时间点的差异分析选择基因通常非常有效 个矩阵：expr.matrix、pd、fd # 将Seurat中的对象转换为monocle识别的对象 #cds <- importCDS(GetAssayData(seurat.object)) #选择做拟时序的亚群
2.3K10发布于 2020-04-01 - 来自专栏生信技能树
monocle单样本拟时序分析
0.背景知识 做拟时序分析是为了探索自己感兴趣的几种细胞之间的发育关系，一般不是用全部类型的细胞来做的。例如本例中选择了CD14和CD16单核细胞。 如果让ai来说拟时序的目的那就比我说的多: 拟时序分析（Pseudo-time analysis）是一种用于理解细胞状态和细胞命运变化的计算生物学方法。 以下是进行拟时序分析的几个主要目的： 细胞状态推断：通过分析单个细胞的基因表达模式，推断细胞在生物学过程中所处的状态。 基因调控网络推断：通过分析基因表达随“拟时间”的变化，推断基因调控网络和信号传导途径。 疾病机理探索：在疾病研究中，拟时序分析有助于揭示疾病发生和发展的分子机制。 经典的拟时序热图 展示了一些基因是如何随着时间轨迹的变化而渐变的，这个渐变不同于findmarkers，是体现变化过程的，而不是直接给出差异表达的基因。
68310编辑于 2024-07-05 - 来自专栏生信技能树
monocle2拟时序实战细节剖析（拟时序一本通04）
我们也简单的展示了目前的可以做拟时序分析的软件的测评，详见：拟时序的多种算法大比拼（拟时序一本通03） 。 创造一个函数包装monocle2的拟时序分析 前面在为什么做拟时序 （展示差异细节）我们简单的演示了如何做一个极简的使用monocle2包的拟时序分析 ，就是构建好对象后的挑选基因，降维，排序这3个步骤即可 所以如果我们仅仅是拿0和1的编号亚群CD14单核细胞去做拟时序，这个时候我们就算是抹去个体差异后做拟时序，结果也会很诡异，像一个蜈蚣一样的，如下所示： 抹除个体差异后做拟时序 其实之前我们就摸索过，详见 ：听说你的拟时序图跑的像蜈蚣，那个时候是因为没有抹去个体差异后做拟时序所以结果很多分叉，但是这个时候是因为抹去个体差异后做拟时序结果很多分叉，所以我也尝试一下不抹除个体差异后做拟时序，确实是有改变，如下所示 提高降维聚类分群的分辨率 我们关心什么差异，就应该是做什么样的拟时序分析，需要牢记为什么做拟时序 （展示差异细节）。
3.2K20编辑于 2024-03-26 - 来自专栏生信技能树
拟时序的多种算法大比拼（拟时序一本通03）
这个纯粹就是生物信息学领域的“马太效应”，大家都用monocle2做拟时序，所以后来者就简单的追随即可，而且绝大部分人其实并不关心算法细节，仅仅是为了做拟时序而做，那么就无所谓选择哪个软件了。 ，R包SCORPIUS推荐 把基因表达量画在拟时序结果图上 拟时序分析就是差异分析的细节剖析 拟时序分析的10个步骤 另外，值得一提的是这个2019的综述其实还把他们的45个软件的测评过程和结果做成了一个易于使用的 而且可以超脱表达量矩阵本身，引入全新信息来开发拟时序工具。 拟时序的正确姿势，错误示范，创新型拟时序 拟时序的多种算法大比拼 monocle2拟时序实战 降维聚类分群 拟时序 多种基础可视化 正向特殊可视化（目标基因表达量，基因集打分） 反方向特殊可视化（5种可视化 ） monocle3拟时序实战 SCORPIUS实战 Slingshot实战 基于其它编程语言的拟时序分析 然而，本人的理解肯定是片面的，希望通过这10个推文的指引能让大家对单细胞转录组的明星分析方法拟时序有一个基础的认知
86510编辑于 2024-03-25 - 来自专栏生信学习111
单细胞5 拟时序分析
1 拟时序分析拟时序分析是为了探索自己感兴趣的几种细胞之间的发育关系，一般不是用全部类型的细胞来做的。 实在不行问问ai，回答可详细1.1 单样本拟时序分析#rm(list = ls()) #单样本library(Seurat)library(monocle)library(dplyr)load("sce.Rdata celltype是细胞类型注释，用以下代码添加> scRNA$celltype = Idents(scRNA) #做拟时序分析通常不是拿全部的细胞，而是拿感兴趣的一部分。用subset提取子集即可。 在做拟时序分析的时候，因为是采用差异基因进行排序的，所以要求是两类细胞或者两类以上（要选择的细胞亲缘关系要近一点，有分化的可能性，完全不挨着的细胞不太行）。 基因可以换成你感兴趣的啊plot_cell_trajectory(sc_cds,markers=gs, use_color_gradient=T)有一些变化不明显1.4.3 基因拟时序点图基因表达水平随拟时序变化的可视化方法
67810编辑于 2024-06-26 - 来自专栏生信技能树
monocle多样本拟时序分析
不同样本中的细胞基本是均匀分布在轨迹上的，说明前面的代码很好的去除了样本间的批次效应 plot_cell_trajectory(sc_cds, color_by = 'orig.ident') 经典的拟时序热图
72810编辑于 2024-07-05 - 来自专栏单细胞测序
单细胞测序—拟时序分析综合
单细胞测序—拟时序分析综合拟时序分析（Pseudotime Analysis）在单细胞测序（Single-cell RNA-seq）中是一个重要的分析步骤，主要用于研究细胞在发育过程或其他生物学过程中所经历的状态变化 拟时序分析可以帮助研究者识别出这些分化路径，并且理解在这些路径中哪些基因起了关键作用。理解细胞状态的动态变化：拟时序分析允许研究者推断出细胞状态的连续变化，而不仅仅是静态的分类。 为后续的拟时序分析创建必要的文件结构和工作目录。 这种方法适用于拟时序分析，因为它可以揭示细胞状态如何随时间或分化进程变化。orderCells()：这个函数是 Monocle2 中拟时序分析的关键步骤之一。 这在拟时序分析中是关键的一步，因为根状态通常代表最早的或未分化的细胞状态，而其他细胞状态将在拟时序轨迹中相对于这个根状态进行排序。
3K14编辑于 2024-08-30 - 来自专栏天意生信俱乐部
单细胞转录组 | 拟时序分析
拟时序分析概述 拟时序分析（Pseudotime analysis）是一种在单细胞RNA测序（scRNA-seq）数据中重建细胞发育轨迹的计算方法。 为什么需要拟时序分析？ 传统的scRNA-seq分析通常将细胞分为离散的类型，但这忽略了细胞状态是连续变化的事实。 拟时序分析有以下优势： 揭示细胞分化路径和发育轨迹 识别驱动细胞状态转变的关键基因 在没有时间序列实验的情况下推断细胞发育动态 发现新的细胞亚群和中间过渡状态 怎么做拟时序分析？ 拟时序分析的基本流程包括： 从预处理好的数据中创建拟时序对象 降维和聚类 学习细胞轨迹 排序细胞并计算伪时间 分析与伪时间相关的基因表达变化 目前常用的拟时序分析工具包括Monocle3、Slingshot 实际操作：基于Seurat对象的拟时序分析流程 以下是从Seurat对象开始进行拟时序分析的完整流程。前提是已经有一个处理好的Seurat对象，包含标准化数据、降维结果和细胞类型注释。
1.5K20编辑于 2025-03-13 - 来自专栏生信技能树
拟时序分析的10个步骤
降维聚类分群 然后你需要有背景知识才能做 拟时序分析 ：文章是这样描述：focusing on endodermal cells, which have a known trajectory from ，已经针对哪部分细胞亚群进行拟时序分析哦！ i 可以看到，有两个发育分支，接下来就可以展现一下重要的基因，在拟时序的两个分支上面： ? 然后是拟时序的的各自特征基因的生物学功能数据库注释： 这一套图表的代码也是在「生信技能树」使用Smart-seq2单细胞转录组数据探索小鼠性腺发育，链接是：https://www.bilibili.com ，这里略 拟时序分析是否可以对全部单细胞亚群呢？
2.6K20发布于 2021-07-06 - 来自专栏生物信息学、python、R、linux
单细胞数据拟时序分析-destiny
单细胞数据分析常用到建立trajectory和pseudoTime，拟时序分析可以用 Diffusion（ Destiny R package） #Diffusion PseudoTime Analysis
4.1K20发布于 2020-04-01 - 来自专栏用户7627119的专栏
GeneSwitches：拟时序分化轨迹开关基因分析
通过拟时序分化轨迹分析可以推断出生物过程中细胞的分化轨迹或细胞亚型的演化过程，用拟时间来表征，拟时间越早代表该细胞位于分化转化过程的早期，反之则位于分化转化过程的晚期。 然而，目前的拟时序轨迹分析软件缺乏对轨迹关键驱动基因的下游分析，2020年发表在Bioinformatics杂志上的分析软件GeneSwitches基于此现状，重点开发了拟时序分化轨迹开关基因分析方法， R2值，其中表达激活的开关基因与拟时序正相关(R2>0)，被定义为上调型开关基因(Up-regulation)，而表达沉默的开关基因与拟时序负相关(R2<0)，被定义为下调型开关基因(Down-regulation 拟时序相关性越高代表基因与该轨迹进程的关系越密切。得到了每个潜在开关基因的开关时间和相关性R2值后，将Top开关基因按其开关时间在拟时序上排序可视化，可以更直观地展示分化轨迹进程中的关键基因作用。 开关基因富集分析结果表明与细胞周期(Cell Cycle)相关的通路在拟时序早期下调，而与心脏功能相关的通路在拟时序后期上调。
2.8K53编辑于 2022-09-21 - 来自专栏青青青山的学习笔记
单细胞转录组之拟时序分析
plot_genes_in_pseudotime(cds[cg,],color_by = "Cluster") 图片图片前面根据差异基因，推断好了拟时序，也就是说把差异基因动态化了，后面就可以具体推断哪些基因随着拟时序如何的变化 my_cds_subset=cds# 拟时序数据和细胞位置在pData 中head(pData(my_cds_subset))# 这个differentialGeneTest会比较耗费时间，测试每个基因的拟时序表达 gene_short_name) -> my_pseudotime_genemy_pseudotime_gene=my_pseudotime_gene[,1]my_pseudotime_gene#绘制一个或多个基因的拟时序 ", nrow= 3,ncol = NULL )+ scale_color_nejm()ggsave('monocle_top6_pseudotime_by_cluster.pdf')图片将前50个随拟时序变化的基因做聚类热图 gene_short_name %in% head(my_row[my_row$cluster == 2,'gene'])))my_gene# 绘制分支处的基因拟时序轨迹
4.6K22编辑于 2022-07-10 - 来自专栏用户7627119的专栏
单细胞测序数据拟时序分析
今天我们就跟随王老师一起来看一下BD SeqGeq™之单细胞测序数据拟时序分析。 ? 什么是拟时序分析？ 不同状态的细胞连起来就是一个“轨迹”，类似根据时间关系判断出的发育轨迹变化，但这并不是真正的时间，而是通过基因表达来模拟细胞的发育演化过程，所以称为“拟”时序分析（Pseudotime analysis 目前拟时序分析最常用的方法 Monocle是由Trapnell实验室开发的，采用了无监督算法，将单个细胞按照拟时间排列在对应的轨迹上。 BD SeqGeq™ 支持拟时序分析 BD SeqGeq™目前将Monocle v2.0整合为插件。 下面就为大家详细展示如何在SeqGeq™中获取Monocle以及使用它进行拟时序分析。
5K20发布于 2020-08-06 - 来自专栏生信技能树
把基因表达量画在拟时序结果图上
学员在练习单细胞时遇到了一个拟时序难题，如下所示的右图，因为它在monocle这样的拟时序R包的官方文档并没有介绍，而我以前没有这样的需求，所以就没有分享这方面笔记。 ，我们产生了 output_of_phe2_monocle.Rdata 文件，就是拟时序分析的结果， 就可以进行各种各样的可视化啦！ 而且我们演示了3种可视化选择，包括：Cluster，Pseudotime，State ，请务必自行看前面的教程：拟时序分析就是差异分析的细节剖析，可视化的代码如下所示： library(ggsci) p1 ) + scale_color_npg() p3 ggsave('trajectory_by_State.pdf') library(patchwork) p1+p2/p3 两个单核细胞的亚群拟时序 表达量就跟拟时序的结果结合起来啦！
3.1K30编辑于 2022-01-10 - 来自专栏单细胞天地
拟时序分析的热图提取基因问题
昨天我在单细胞天地讲解了使用monocle2进行拟时序分析的方法，基本上跟着我的代码走一波就可以学会了，当然具体参数理解需要自行发力哦，见：使用monocle做拟时序分析（单细胞谱系发育） 用法只是最基础的知识而已 ，更多的时候，我们需要活学活用，比如课程学员提到的问题，就是因为做不到活学活用，他想知道下面的拟时序分析的热图提取基因，学员把基因按照发育顺序绘制了热图，而这些基因被他分成了3组，想拿基因去做GO/KEGG 我这里不能拿学员真实项目数据来演示，所以还是用我们的老朋友，拿scRNAseq包的表达矩阵测试，见：使用monocle做拟时序分析（单细胞谱系发育） 首先根据细胞发育谱系来绘制热图 因为前面的教程 使用 monocle做拟时序分析（单细胞谱系发育） 我们已经把细胞发育情况做出来了，就是NPC细胞跟另外3种细胞从生理上就不一样，所以是单独的发育轨迹，而 “GW16” and “GW21” ，“GW21+3 ” 这种孕期细胞，就可以很清晰的看到时间被反映在我们的拟时序分析结果了。
3.2K30发布于 2020-03-30 - 来自专栏生信技能树
拟时序分析就是差异分析的细节剖析
副标题： 所有的大样本量差异分析都可以转为拟时序分析 两个分组的差异分析仅仅是上下调吗？ 很多小伙伴在后台表示对单细胞数据分析里面的拟时序分析不理解，恰好最近看到了一个超级清晰明了的展现拟时序分析的作用的文献，分享给大家。 它完美的展现了差异分析为什么不够，为什么拟时序分析就是差异分析的细节剖析。 拟时序的结果 拟时序虽然就展现了上面的一个图，但是具体的代码步骤还是有点难度哦。 拟时序的差异基因表达量图 最后这个图，看起来有技术含量！
3.3K30编辑于 2021-12-27 - 来自专栏生信星球
做拟时序的你，苦igraph已久了吧
时间已过半年，他们还是没有解决这个问题。受苦最多的应该是我吧！我的学生搞不定就会问我呀！他们就不能更新一下让我少费口舌吗！唉~他们不更，让我来吧。
3K02编辑于 2025-07-14
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号