看你的文件应该没有问题,再检查一下文件结尾是否有多余的空行,有的话删除试一下; 如果还不行,建议使用我们提供的biolinux 里面安装了python 的mcscan,怀疑你的软件安装问题; https://www.omicsclass.com/article/526 基因家族分析课程:学习链接:基因家族分析实操课程、基因家族文献思路解读
回答于 2020-03-02 20:02
电脑打开网易云课堂,课程的左下角有参考资料下载,所有的课程都检查一遍都下载下来; 有个打包文件解压之后里面有所有的脚步文件;
回答于 2020-03-02 19:55
电脑打开网易云课堂,课程的左下角有参考资料下载,所有的课程都检查一遍都下载下来; 有个打包文件解压之后里面有所有的脚步文件;
回答于 2020-03-02 19:55
这个是由于一个基因有多个 entrez ID 导致的: 下面的代码是把多个entreID取了并集,建议取一个作为代表序列;可以提前用Excel处理好基因信息文件表格; DEG_list_kegg <- c()for(i in degs$entrezid){ DEG_list_kegg<-c(DEG_list_kegg,eval(parse(text = i)))}
回答于 2020-03-02 19:53
这个涉及到基因组序列的批量注释,你可以把基因组上所有基因的序列到kegg数据库里面同源比对一下; 看看哪些基因比对到你感兴趣通路中就找到你想要的基因列表了;
回答于 2020-03-02 19:43