regions function

Aggregate regional data