ANALYZING SPATIALLY VARYING RELATIONSHIPS

 

Several techniques for analyzing more complex instances of spatial inequality.  Often times, inequality results from a combination of factors and cannot be explained by mapping one variable alone.  In these cases, we can map how multiple factors overlap in space and determine places where relationships between factors are strong or weak.

EXAMPLE: HOUSEHOLDS WITHOUT CARS AND POVERTY

In this example we are interested in a program that encourages alternatives to owning a car.  The program is more successful in some neighborhoods with better transit and walking opportunities but in other neighborhoods households’ lack of a car may be a function of poverty.  We consider car ownership in each neighborhood compared to the level of poverty in the area with the assumption that neighborhoods with more car-free households than poverty are more successful at creating alternatives to car ownership.

Using the GeoDa program, we begin with a standard deviation map of the percent of households without a car.  (See previous tutorial for more information on this step.)
http://sites.google.com/site/spatialjustice/_/rsrc/1257550544047/spatialjustice/figure-inferential-statistics/No%20Car%20Rate%20Map.png?height=370&width=400
There are some noticeable areas where the  “No Car Rate” is high, but we would like to know how much this is due to higher poverty as opposed to households having alternatives to driving.  We next create a standard deviation map of poverty by Census Tract.  It appears to have some similarity to the map above.
http://sites.google.com/site/spatialjustice/_/rsrc/1257550835758/spatialjustice/figure-inferential-statistics/Pov%20Rate.png?height=375&width=400

We can eyeball some similarities between the No-car rate and the poverty and a more precise way to consider this relationship is with a scatter plot.  In GeoDa we specify a scatter plot of the percent of households without car versus percent poverty.
http://sites.google.com/site/spatialjustice/_/rsrc/1257550981625/spatialjustice/figure-inferential-statistics/Scatter%20-Car%20-%20Pov.png?height=378&width=400
The scatter plot suggests a linear relationship between poverty and the lack of a car, because most data points lie fairly close to a straight line.  The cluster of points well above the line are neighborhoods where the no-car rate is higher than the poverty rate.  We will investigate these neighborhoods further.  A more formal investigation of the relationship is done with a regression model (found under the “Methods” menu in GeoDa.
http://sites.google.com/site/spatialjustice/_/rsrc/1257551065721/spatialjustice/figure-inferential-statistics/Regression%20Set%20up.png?height=224&width=320

We specify the title and files name and check the option for the Moran’s I z-score (this option can increase the computation time). We set up a spatial lag model that considers poverty rate as an independent variable to explain car-less rate, the dependent variable.  The use of a spatial lag model takes spatial autocorrelation into account, and will help us pinpoint areas where the no-car rate is higher or lower than expected.  The spatial lag model requires that a weight file be selected to determine the spatial relationships.  (See the previous tutorial for more on how to create a spatial weights file). Once the model is set up as desired we hit “Run.”
http://sites.google.com/site/spatialjustice/_/rsrc/1257551201856/spatialjustice/figure-inferential-statistics/Lag%20Regression.png?height=400&width=297

Once the model has run, use the ‘Save’ button to save the predicted values, residuals and errors of the model.  These saved values are added to the attributes of the shapefile, allowing them to be mapped in later steps.
http://sites.google.com/site/spatialjustice/_/rsrc/1257551725444/spatialjustice/figure-inferential-statistics/Save%20Residuals.png?height=400&width=294

The regression result listing will appear after clicking “OK.”  This text contains the information for the overall model fit.

http://sites.google.com/site/spatialjustice/_/rsrc/1257551792249/spatialjustice/figure-inferential-statistics/OLS%20Results.png?height=260&width=400

 

Although we selected a spatial age model, results for an ordinary least square (OLS) regression appear first. These results show that a model that only considers the poverty rate has a fairly good fit in explaining the no-car rate.  This model serves as a baseline for comparing the spatial lag presented below.
http://sites.google.com/site/spatialjustice/_/rsrc/1257551892437/spatialjustice/figure-inferential-statistics/Spatial%20Lag%20Result.png?height=238&width=400
The spatial lag model includes an additional variable “W_NoCarRate” which is the spatial lag term.  The spatial lag term accounts for the effect of adjacent neighborhoods.  The fact that this term is statistically significant indicates that spatial factors are important in explaining the no-car rate.  We can then map the model results to better identify where spatial factors are important.
http://sites.google.com/site/spatialjustice/_/rsrc/1257552012403/spatialjustice/figure-inferential-statistics/Lag%20Error%20Map.png?height=370&width=400

A standard deviation map of the model prediction error, “PRDERR,”  highlights where the poverty rate was least effective in predicting the no-car rate.  This term does not account for the spatial lag term, only the estimate based on the poverty rate.  The fact there is clustering of error values indicates that spatial factors are important in this model (which is confirmed in the model results previously discussed).  The negative values in southern neighborhoods indicate that the no-car rate is lower than expected given the poverty rate.  Red and orange clusters are areas where the no-car rate is higher than expected.

http://sites.google.com/site/spatialjustice/_/rsrc/1257552332334/spatialjustice/figure-inferential-statistics/Lag%20Residual%20Map.png?height=382&width=400
Lastly, we map the model residual, which accounts for the spatial lag term.  Spatial clusters are less apparent on this map because the influence of neighboring tracts has been incorporated into the model.  The residual, or remaining error, appears random indicating that the spatial lag model has largely accounted for spatial differences.

Returning to our original research questions, we can argue that spatial factors are important to whether households have alternatives to car ownership.  In the southern part of the city, the absence of a car seems to be due primarily to poverty.  In the central part of the city households are more likely to do without a car regardless of income.  It follows that neighborhood improvements may be needed in the southern parts of the city if more households are to be encouraged to do without a car.