老师,我将ID.txt文件按照基因家族列表前缀进行了修改,并排序,之后进行perl /work/scripts/PAV_list.pl -pangenel foxtail_millet.GRASgene.pan-genes.list -spid new_ID2.txt -unigenel foxtail_millet.GRASgene.uniq.list -fam GRAS得到的PAV_draw.list和PAV_gene.list,为
,
得到的PAV_gene.list就有1970条,但GRAS133一会就有没有数据了,都是空白的,
之后作图,Rscript $scriptsdir/PAV_heatmap.r -i PAV_draw.list -sp 112,得到的图
请问,问题出在哪里,