这是我的数据的快照:
structure(list(CPUBID = c(1000001L, 1000002L, 1000003L, 10001L,
1000201L, 1000203L, 10003L, 1000801L, 1000802L, 1000803L, 1001L,
1001101L, 1001102L, 1001601L, 1002401L, 1002402L, 1002403L, 1002601L,
1002602L, 1002604L), MPUBID = c(10000L, 10000L, 10000L, 100L,
10002L, 10002L, 100L, 10008L, 10008L, 10008L, 10L, 10011L, 10011L,
10016L, 10024L, 10024L, 10024L, 10026L, 10026L, 10026L), CYRB = c(1982L,
1984L, 1988L, 1985L, 1986L, 1992L, 1993L, 1984L, 1986L, 1988L,
1983L, 1987L, 1992L, 1977L, 1981L, 1984L, 1998L, 1980L, 1981L,
1984L), twinfam = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, …