故事起因是一篇文献的monocle分析结果(如下图),这篇文章之前我们公众号还投稿解读过((读文献)10分单细胞文章就是这么发的)。这个图有几个问题需要实现,首先肯定不是monocle自带函数完成的,其次我们要解决的就是提取monocle数据,在ggplot中实现轨迹和散点的完成,以及添加饼图。这个图其实很有意义,这里可以明确的看出来不同state拟时区间不同样本的比例,对于解释生物学意义很有帮助。
![](https://img-blog.csdnimg.cn/img_convert/0623873eb0fac68edbcac121fbd43fd6.png)
![](https://img-blog.csdnimg.cn/img_convert/17401b55484f2f8f156ead98718f0d06.png)
部分内容参考了原作者的函数,一般我们拟时做的几个轨迹图如下:
library(monocle)
#monocle中的可视化函数
plot_cell_trajectory(mouse_monocle, color_by = "Pseudotime")
plot_cell_trajectory(mouse_monocle, color_by = "State")
plot_cell_trajectory(mouse_monocle, color_by = "celltype")
![](https://img-blog.csdnimg.cn/img_convert/11212f148292581032a3ab5dc68e2355.png)
接着,我们将需要的数据全部提取出来,在ggplot中重构这个图:需要做的就是提取数据和ggplot作图,更多精彩请至我的公众号-KS科研分享与服务
![](https://img-blog.csdnimg.cn/img_convert/97d40b31a649a8c482c0cec549fd1cde.png)