根据第一列重复选择行

紫风8824

文件123.txt 格式如下：
#CHROM        POS        ID        REF        ALT        QUAL
1        130        .        T        C        31.4
1        145        .        C        T        3.52
1        155        .        T        C        3.22
1        159        .        C        T        3.22
1        176        .        G        A        7.37
2        186        .        C        T        39
2        208        .        T        C        4.64
2        1067        .        G        A        36.9
2        1120        .        G        A        10.2
2        1130        .        C        T        11.3
3        1132        .        T        C        10.3
3        1148        .        T        C        10.9
3        1150        .        T        C        20.7
3        1178        .        G        C        68.4
3        1179        .        G        A        77.3
4        1292        .        A        T        26
4        1298        .        C        A        25
4        1395        .        T        C        38.6
4        1399        .        T        C        10.5
4        1395        .        T        C        38.6
。。。
根据第一列重复的，随机选取3行，如第一列是1的有5行，随机取3行；第一行是2的也是这样。
得到的文件456.txt 可能是这样的
#CHROM        POS        ID        REF        ALT        QUAL
1        130        .        T        C        31.4
1        145        .        C        T        3.52
1        155        .        T        C        3.22
2        186        .        C        T        39
2        208        .        T        C        4.64
2        1130        .        C        T        11.3
3        1132        .        T        C        10.3
3        1148        .        T        C        10.9
3        1150        .        T        C        20.7
4        1292        .        A        T        26
4        1399        .        T        C        10.5
4        1395        .        T        C        38.6
请大家指点？谢谢！

MMMIX

回复 1# 紫风8824


    awk:

NR > 1 {
    if ($1 != last_one) {
        last_one = $1
        cnt = 1
        print
    } else if (cnt < 3) {
        cnt++
        print
    }
}

ly5066113

回复 1# 紫风8824

1. awk 'a[$1]++<3' 123.txt > 456.txt

复制代码

tuyajie

awk '{a[$1]++}(a[$1]<=3)' test.dat

hjfeng1988

你们都不随机