Spatially join these to the dataset on the PUMA, and count the number of felonies per square feet in that borough as the dependent variable, and the independent variable the different proportion of the five hierarchies location in that borough.