老师我做存在缺失分析部分,awk -F"\t" '{if($0!~/^Fam_N/){for(i=2;i<=27;i++){print $1"\t"$i}}}' PAV_gene.list |sed "s/;$//"|awk -F"\t" '{gsub(";","\n"$1"\t",$0);print}'|awk -F"\t" '{if($2!=""){print}}' > Fam_gene.list这一步,得到的泛基因编号为什么不是连续的呢?PAV_gene.list中每行都对应的有基因啊,老师的课程里都是连续的
你有多少个样本啊,我在课程里面说了要把这个地方改成你对应的样本数+1,你是不是没改,你改一下:
awk -F"\t" '{if($0!~/^Fam_N/){for(i=2;i<=你的样本数+1;i++){print $1"\t"$i}}}' PAV_gene.list |sed "s/;$//"|awk -F"\t" '{gsub(";","\n"$1"\t",$0);print}'|awk -F"\t" '{if($2!=""){print}}' > Fam_gene.list