> table(table(npdb$ID)) 1 2 3 4 5 6 8 11 15 22 73 6105 235 12 4 7 1 1 1 1 1 1 > tmp = table(table(npdb$ID)) > sum(tmp[-1])/sum(tmp) # percent 2 or more [1] 0.04145