是的,比对用clean data就可以,bwa需要参考基因组,比对之前要建立索引才行。
回答于 2019-11-14 07:16
输入的文件格式有问题,你的文件与示例对一下,看看格式是否一致。
回答于 2019-11-13 09:36
你看看两个参考基因组是否有差别,如果有差别会有这个现象,发文章,用的哪个基因组引用文献就行;
回答于 2019-11-12 20:23
1.你看下你的CDS序列ID在不在GFF里面,检查一下;如果不在说明两个文件有问题。 2.你的GFF文件第三列只有mRNA信息,缺失GENE信息建议添加一下:https://www.omicsclass.com/article/816 如果检查完毕,就可以继续往下分析了;
回答于 2019-11-11 16:47
通常只有的网络错误,都很诡异,毕竟长城在这里。 解决方式也简单,加上代码: options( 'download.file.method.GEOquery' = 'libcurl' ) 即可! 再次运行后如下: > options( 'download.file.method.GEOquery' = 'libcurl' )> ###获取数据> gset = getGEO(GSE, GSEMatrix =TRUE, AnnotGPL=TRUE,destdir=workdir)...
回答于 2019-11-11 16:42
数据下载,可根据TCGA网站的筛选条件添加GDCquery()参数,详情见: https://www.omicsclass.com/article/1060
回答于 2019-11-11 16:00
可以查看hmmsearch的帮助,自定义设置搜索参数:官方说明 http://eddylab.org/software/hmmer/Userguide.pdf
回答于 2019-11-11 15:57