分享

Admixture群体结构可视化 + Group!

 生信药丸 2024-05-16 发布于贵州

写在前面

前述,因项目需要,我快快写了一个基于「Admixture」软件的群体结构分析结果可视化的小功能。个中原因?那是有一个似乎已经发表的R包,实在太难用。明明简单的东西,总要搞得很复杂,似乎一切就是为了:(1)把事情搞复杂;(2)把论文发出去。我觉得完全没必要。反正用的上我就写出来,合适就分享出来,至于发表不发表,主要不是看我想不想,而是工作有需要,那写出来发表,如果没需要,就真的懒得写论文。毕竟这实在花时间。

Admixture Q matrix Viz 增加分组标签功能

回归主题,前述,TBtools 的「Admixture Qmatrix Viz」功能可以输出这样的图

说实话,确实还挺好看的。用起来也方便,于是我不时收到用户反馈,说能不能这样或者那样优化。

自然,我都说好好好,但是确实没下定决定去做哪怕一点点优化。但有一个用户朋友非常热情,他真的花了时间:

  1. 写了word文档,给我描述了具体需求和预期效果

  2. 提供了 demo 数据

为此,我含泪花了大概30min优化了这个功能,于是他现在支持了「对样品进行分组 Label」。其实群体一多,这些样品标签确实就会堆叠在一起,如果样品几百个,图片长度又有限制。

那么怎么做出一些分组信息:

  1. 自己对样品已经有分组预期,那么就按照 样品名\t分组名 的形式准备好分组信息即可

  2. 如果没有预期,那么可以按照下述的神操作来折腾。

首先,直接把上图导出为 PDF,随后复制样品名

放到 Excel 里面,增加一列分组

随后复制黏贴到输入就行了

随后点击可视化就可以了

当然,从上图可以看出,K=5 正好符合我们的预期(当然符合...因为这个“预期”我们本来就按照 K=5 来做的)。看看 K=4,可以看到标签交互出现。为什么会这样?因为他就是交互出现...我们按照 Qratio 排序,当K=4的时候,排序上会交叉。不过,似乎「Guangxi」和 「Hainan」的分组比较稳定。

写在最后

当然,估计还是要个体足够多,这个功能才会好看,或者,我可以把图片宽度缩小一下?这样是不是就好看多了?

如果你知道怎么用 TBtools 图片交互?当然也没必要了。你导出用矢量图编辑工具也挺好

附上示例数据

https://tbtools./s/9587378dabd94b

    转藏 分享 献花(0

    0条评论

    发表

    请遵守用户 评论公约

    类似文章 更多