Clipping income data

For a project I am working on I would like to intersect fire station response districts with income data to determine the median income per response district, however the borders of the layers obviously do not line up. What is the best method to go about manipulating the data to have the output show the median income based on the response districts.

I am using ArcMap 10.2.2 for this analysis. I have a polygon shapefile with median income information and a polygon shapefile of fire stations first in response areas.