Emulate R aggregate R: aggregate(employees$salary, list(job = employees$job), mean) agate: jobs = employees.group_by('job') aggregates = jobs.aggregate([ ('mean', agate.Mean('salary')) ])