omicsgene
omicsgene - 生物信息
实名认证 专业认证

性别: 北京 - 北京市 注册于 2018-04-20

擅长:重测序,遗传进化,转录组,GWAS

生物信息

向TA求助
14308金币数
80900 经验值
442个粉丝
主页被访问 86219 次

4089 个回答

0 赞同

哥,您好,如何检查GFF文件里面的mRNA ID是否对应一致,或者如...

脚本认为ID=与;之间的为gene或mRNA的ID,建议你把GFF文件里面的gene:  CDS:  transcript: 删除; 然后再看看fasta里面的ID与GFF里面的ID是否一致。 这样保持所有地方的ID一致,就不会因为找不到对应的ID而出错。 如果上面按照我说的都改了,脚本这里27行就可以删除gene:    获取ID信息了:

回答于 2019-11-08 10:22

0 赞同

ggplot,超过12种颜色,scale_color_brewer(type="seq",palette=...

这是ggplot2里面内置的色板,如果颜色不够可以使用手动设置色板scale_colour_manual():https://ggplot2.tidyverse.org/reference/scale_manual.html 颜色的获取可以参考:https://www.omicsclass.com/article/746 学习R语言基础与绘图:、R语言画图、R语言快速入门与提高

回答于 2019-11-08 09:36

0 赞同

为什么预测基因家族顺势原件位置信息出来是空表

检查GFF文件里面的mRNA   ID是否对应一致。 参考这个:https://www.omicsclass.com/question/1941

回答于 2019-11-07 10:04

0 赞同

转录组问题

有基因组了为啥还要做无参转录组呢? 如果,是后来出来的基因组,建议重新按照有参转录组分析一下。而不是用unigene比基因组;

回答于 2019-11-07 10:03

0 赞同

R语言怎么写循环,对多个元素进行重复操作?

可以用列的索引解决; 建议学习一下R语言基础,这些问题不需要循环就可以解决; R语言画图、R语言快速入门与提高

回答于 2019-11-07 10:01

0 赞同

基因 定位到染色体时,获取染色体的位置信息报错[fai_build_core...

可能你的fasta文件不是标准的格式,你看看Chr2序列折叠之后长度是否一致;

回答于 2019-11-06 17:46

0 赞同

老师你好,关于基因家族成员的鉴定

不同的人鉴定出来的基因数量可能会有些许差异,正常的。只要分析的步骤严禁的,一般都能发文章的;

回答于 2019-11-06 14:47

0 赞同

用 hisat2 建索引时出现 “is not reverse-deterministic, so rev...

建立索引这一步很消耗内存,几百G以上内存的节点运行一下试试;

回答于 2019-11-06 14:46

0 赞同

为什么第一次搜索结构域搜出来是空的?

如果输入的文件都没有错误的话,可能真的就是没有;

回答于 2019-11-06 14:44

0 赞同

基因家族加倍复制时,blastall所有基因比对所有基因,程序运行了...

运行时间和电脑的配置与基因的数量有关。如果电脑性能差,基因数量大,运行会慢一些,耐心等待。 2018年买的5000元的笔记本电脑,40000 vs 40000基因比对,大概运行一天时间,你可以参考一下;

回答于 2019-11-04 17:09