- 论坛徽章:
- 25
|
本帖最后由 wh7211 于 2021-05-27 16:06 编辑
回复 1# chenjiali1
从你发的结果看,应该是“第5列除以第5列第6列之和”吧?
- cat 1
- chr1 3003639 3003640 1.0 2 0
- chr1 3003720 3003721 0.75 3 1
- chr1 3003639 3003640 1.0 2 0
- chr1 3003720 3003721 0.0 0 2
- chr1 3003639 3003640 1.0 2 1
- chr1 3003720 3003721 0.0 0 2
- awk '{a=$1" "$2" "$3;if(!b[a]++){c[a]=++d};e[a]+=$5;f[a]+=$6}END{PROCINFO["sorted_in"]="@val_num_asc";for(i in c){printf("%s %.7f %d %d\n",i,e[i]/(e[i]+f[i]),e[i],f[i])}}' 1
- chr1 3003639 3003640 0.8571429 6 1
- chr1 3003720 3003721 0.3750000 3 5
复制代码 |
|