写在前面没错,标题党就是我! 更新了逻辑,准确性提高了估计花了两个多小时,定位到了逻辑可调整的位置,一顿乱改。测试了下来, 增加了出图步骤,方便人工校准前述,海南道明寺和小庭子在夏老师的指导下,对数据库做了一次不错的更新,增加了 PHAS位点的可视化。 感兴趣的可以登录 http:/// 查看。 事实上,这个非常实用。人工校准,可以说是目前小RNA数据分析绝不可少的步骤。当然,在数据库上,重点就在于用户可以查看。图片是使用 ggplot2 绘制的,数据则是 python 脚本处理的。这个图的绘制方式,其实是文献中常见的PHAS位点可视化方式。项目这段时间不是我的工作重点,所以我前面主要只做了解。其中可以加速的步骤有很多,最简单的包括 py准备数据时可以注意内存损耗和频繁遍历,R绘图的调用要考虑终端调用的成本。 既然流程更新了,很明显,我认为这个可视化确实实用(Emmm,值得进一步加速),于是用 JIGplot 实现了类似的图。 是的,JIGplot 直接可以输出可交互的图片。然而在此处没啥用。不过 JIGplot 其实支持直接出图~~ 测试了下,大概2600+个 PHAS Loci。 单线程的情况下 10分钟搞定!换句话说,如果是开 60个线程,其实...加上开销,我估计 一分钟 搞定,没啥问题。 题外话,IGV-sRNA还是要提一下,IGV-sRNA,我自认为这个工具应是目前植物小RNA数据分析必备!当然,其人工校正 PHAS 位点的部分可以用通过上线取代。但是他的交互,和多组学数据整合,以及其他专门定制的植物小RNA数据人工探索功能,几乎是无法被取代~~~还是贴一下张图。 说实话,真的流弊,而且有用。之前我做的部分改进已经PR到IGV官方版并被接收啦。 写在最后这里贴上之前写了一半的推文,我应是不会完善他了。不过希望我的朋友们明白,我一直都会是积极乐观,也非常感谢各位的关注和支持。 Emmm... 最后还是要补充,流程属于sRNAminer的一部分,一个高效(超高速)且跨平台(超顺手)的 植物小RNA 数据分析工具,不过我估计他不会太快面世,毕竟....还有很多其他东西要整。所以感兴趣的,直接给PI: rxia@scau.edu.cn 丢邮件。 |
|