首页 > 其他 > 详细

统计文件某列各项或某项出现次数(VCF文件中每个染色体信息)

时间:2021-09-02 08:36:39      阅读:15      评论:0      收藏:0      [点我收藏+]
#!/bin/bash
grep -v ^# $i|awk ‘{sum[$1+=1]END{for(i in sum)print i"\t" sum[i]}‘ >chr_info
#按大小排序
sort -n chr_info
#某列求和
awk ‘{sum +=$2;END {print sum}‘ chr_info
#符合要求的某列求和
awk ‘/[0,9]/ {sum +=$2;END {print sum}‘ chr_info
#统计某项出现次数,第一列为chr1,第一列为chr2
awk ‘{if($1=="chr1") ++sum1;if($1=="chr2") ++sum2}END{print "00""\t"sum1"\n""01""\t"sum2}‘ $i

统计文件某列各项或某项出现次数(VCF文件中每个染色体信息)

原文:https://www.cnblogs.com/xiaosagege/p/15212160.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!