送分题!已知拟南芥 XYZs 基因功能关键,一分钟求当前物种的 XYZs 是哪些?
昨天推了一个新功能,这个功能实在是表现得太让我吃惊,因为不仅仅已知符合我的预期,而且帮我找到课题分析中错漏的地方。
感兴趣地可以翻前天晚上的推文,快速且准确获取家族进化分支成员(以分钟计)。[ 最近推文少,阅读量似乎也达不到关注人数的十分之一了,Emmm,赶紧先转发朋友圈,看完了别忘了点再看。不然,感觉继续写起来就没意思了 ]
昨天我做了一些分析,发现了确实流弊。当然也存在一点点小的问题。软件逻辑和结果是没问题的。因为同时输出了文字结果和图片信息(当然,我还有 阳性集合)。具体可以如下
这种情况想面,按照前天晚上凌晨的版本,除非调参数次,否则这两个成员可能不会同时被很好地鉴定出来(注意,这是阳性集)。
我做了不少调整:
纳入合适的“外类”
分支如上,其实上述三个基因是一个分支,但是在无根树(按照有根树)的显示方式,看起来不是一个分支(按照找基友的逻辑,鉴定补全);于是,我又写了一个置根的功能,不管你的树有多少个节点,反正我把每个节点都依次当根来处理,看看这个分支是否一直稳定。
Emmm,可能有我朋友圈才看得到。
那么大体结果就如上。
功能推出之后,师妹就问了一个问题,能不能搞一个分支?!而不是一个节点。很多时候,我们知道,通路上某个节点对应的可能就是某个基因的几个家族成员(或者特定进化分支),既然如此,在我们感兴趣的物种上,需要找的也是一个分支。Emmm,花了一点时间.....于是现在支持分支抓取。
稍微验证了之前我做了几天的家族成员鉴定,然后定分支的结果....
比如拿拟南芥三个成员,找找荔枝的ACS成员分支,结果一致
比如ACO分支,一共是六个,结果一致
说实话,我被自己感动了!
赞 (0)