from wdqs import Client
client = Client()

query = '''
SELECT ?discoverer ?name (COUNT(?asteroid) AS ?count)
 WHERE {
   ?asteroid wdt:P31 wd:Q3863 .
   ?asteroid wdt:P61 ?discoverer .
  SERVICE wikibase:label {
    bd:serviceParam wikibase:language "en" .
    ?discoverer rdfs:label ?name
   }
 }
 GROUP BY ?discoverer ?name
 ORDER BY DESC(?count)
 LIMIT 1000
'''
result = client.query(query)

from pandas import DataFrame
df = DataFrame(result)
df

% matplotlib inline
import matplotlib.pyplot as plt

plt.hist([discoverer['count'] for discoverer in result])
plt.show()