1
Spatial Analysis of Bureaucrats’ Attempts to Resist Political Capture in a Developing
Democracy: The Distribution of Solar Panels in Ghana

Appendix (Online): Technical Discussion of Data and Methods
Part A: Dataset Creation
District Shapefiles
For the analysis in “Spatial Analysis of Bureaucrats’ Attempts to Resist Political Capture
in a Developing Democracy: The Distribution of Solar Panels in Ghana,” we created a new
spatial dataset by merging several types of data. We use a shapefile for Ghana’s districts that was
created by the RS/GIS Lab at the University of Ghana-Legon. Existing shapefiles were of 2012
jurisdictions, requiring us to create new shapefiles that represent district boundaries in 2008,
before the solar panel commenced. This was necessary because some districts were createdbetween 2008 and 2012 by splitting jurisdictions. After identifying the new districts created in2012, we were able to merge them back together to form their 2008 counterparts.

Fitting Variables to Shapefiles
Most variables in the analysis required some form of transformation in order to fit into
these shapefiles. We took three approaches to variable transformation, depending on data type.
First, for variables that originated at the constituency level (NDC vote share 2008; voter turnout
2008; and ethnic fractionalization), we created district-level data by computing average values
among the constituencies in the district. These averages are weighted by population (although
population-weighted averages and averages that are not weighted by population are almost
identical). Out of 230 total constituencies, we implement this transformation for 97
constituencies that fall within 37 districts. Second, for geo-located data (solar panels and WorldBank projects) we spatially joined the existing data to the district shapefile. Third, for line
shapefiles (electric grid and roads), we used the Intersect feature in ArcMap. We then took theresultant file and used the Dissolve feature in ArcMap to reduce the road and electric grid files
down to one unit per district. The Calculate Geometry tool then allowed us to calculate line
lengths. For electric grid per capita, we divided the length of grid by population. For roaddensity, we divided the length of roads by area of the district. We display descriptive statistics inTable A1 below.

2
Table A1: Descriptive Statistics
Mean
Dependent Variable
Solar projects
Political Variables
Reward NDC loyalists(2008 NDC vote share)
Induce swing voters(NDC vote share volatility)
Mobilize voters where turnout level varies(turnout volatility)
Need Variables
Electric grid per capita (*1000)
Road network density
Control Variables
Ethnic fractionalization
Population density
World Bank projects
Percent nonvoters 2008
Poverty rate
Inequality level (Gini)
Number of health facilities
Literacy
Female ratio

Min

Max

SD

N

7.28

0

202

27.38

170

0.49

0.11

0.96

0.18

170

0.05

0

0.21

0.04

170

0.15

0

0.27

0.04

170

0.28
0.12

0
0.03

4.35
0.42

0.44
0.05

170
170

0.41
391.59
5.36
0.29
30.52
38.69
21.71
0.68
0.51

0.06
7.97
0
0.16
1.95
27.2
3
0.2
0.47

0.82
10124.2
52
0.46
92.4
64
276
0.93
0.55

0.2
1311.76
5.14
0.05
19.96
6.02
26.32
0.18
0.02

170
170
170
170
170
170
170
170
170

Part B: Additional Methodological Discussion and Analysis
Accounting for Spatial Autocorrelation: Selecting Spatially Lagged Dependent Variable(SLDV) Models for Use in the Analysis
In the body of our manuscript, we present SLDV models that account for spatial
clustering. It has been suggested that we might instead aggregate our district data into higherlevel units (regions) in order to employ multi-level models, because fixed effects for given
spatial units offer leverage to examine within-unit effects. Such an approach is not appropriate
for several reasons. First, the approach does not account for clustering. Second, region fixed
effects are not theoretically relevant for the questions asked. Third, given the way thatbureaucrats selected locations, a “nested” approach is not as well-matched, theoretically, as the
SLDV approach. We do present OLS models in the paper and additional models in thisappendix, and those models cluster standard errors by region.
We thus highlight modelling strategies that account for spatial clustering. To confirm that
this was appropriate, first, we use an LM test to examine whether an OLS model sufficiently fits
our data, or if the model needs a spatially lagged variable or spatial error term. Spatially lagged
dependent variable models account for the possibility that the values of the dependent variable in
one location are influenced by the values in nearby locations. Spatial error models, on the other
hand, are intended to account for omitted variables that are spatially autocorrelated. Such

3variables are often omitted because they are either very difficult or impossible to effectively
measure.
The equation for spatial error models is as follows:
y=X(β) + ε, ε =λ(W) ε + u
In this equation, y is the dependent variable, X is a matrix of observations on the explanatory
variables, λ and β are parameters, ε is a vector of spatially autocorrelated error terms, W is thespatial weights matrix, and u is a vector of independent identically distributed (i.i.d.) errors.
To estimate spatially lagged dependent variable models, the following equation is used:
y=(ρ)Wy + X(β) + ε
In this equation, Wy is the spatially lagged dependent variable for weights matrix W, X is amatrix of observations on the explanatory variables, ε is a vector of error terms, and ρ and β areparameters.
While the technical components of SLDV and spatial error models are distinct, they
address closely related issues. Many models may exhibit spatially clustered error terms andspatial autocorrelation in the dependent variable. This produces a need to statistically diagnose
which type of model is more appropriate. Since our data include both violations of theassumptions of OLS, we employ another LM diagnostic test, presented in Table A2.
Use of the LM diagnostic was originally suggested by Anselin and Rey.1 This wasrefined by Anselin et al. to incorporate the robust forms of the statistics.2 In all of the simple LM
tests for error dependence and for a missing spatially lagged dependent variable, our test
statistics were statistically significant. This led us to look to robust versions of these tests. The
LM test for error dependence in the possible presence of a missing spatially lagged dependent
variable is not significant at any of the levels of analysis. Meanwhile, the LM test for a missing
spatially lagged dependent variable in the possible presence of error dependence is significant at
the 0.1 level or lower throughout all models. This indicates that the spatially lagged dependent
variable model is more appropriate than the spatial error model. The variables used in theseanalyses are not highly collinear (defined as Pearson’s r > 0.7), and the variance inflation factor(VIF) clearly indicates that global multicollinearity is not present. We estimate the spatial models
using the “spdep” package in R.3
Researchers often interpret the coefficients immediately estimated by SLDV models(Whitten, Guy D., Laron K. Williams & Cameron Wimpy 2019). Yet this approach conflates
spatial dependencies and the relationship between the independent and dependent variable.
Instead, we use the impacts function as a post-estimation step to isolate the direct effects of our
independent variables. Direct effects are the effects of independent variables at a location i on thedependent variable at location i. Indirect effects are the effects of independent variables at alocation i on the dependent variable on i’s neighboring locations j. Direct effect coefficients canbe interpreted like OLS coefficients. We obtain distributions of our direct effects by using a
Markov Chain Monte Carlo simulation process. We ran 1000 simulations and calculated zstatistics from these simulations.

1Anselin and Rey 1991.
Anselin et al. 1996.
3
See Bivand et al. 2005. Since the R spatial analysis community is transitioning from spdep to spatialreg, future
researchers may need to update their R script accordingly.
2

4
TABLE A2. Diagnostics for Spatially Lagged versus Spatial Error Model Fit
District Data
Statistic

p-value

Simple LM error test

5.43130.019779*

Simple LM spatial lag test

8.77180.003059*

Robust LM error test

1.06630.301779

Robust LM spatial lag test

4.40680.035795*

Portmanteau test: LMerr+RLMlag

9.83810.007306** p<0.05

Second, we use an additional LM test to examine whether a model with only a spatially
lagged dependent variable (SLDV) can effectively fit our data. If this LM test calculates astatistically significant parameter, then there is justification to use a spatial Durbin model(Elhorst 2010). In our case, this test fails to support the need for inclusion of spatially lagged
independent variables, so we have confidence that either an SLDV or spatial error model is mostappropriate.4
Robustness and Sensitivity Checks on Spatially Lagged Dependent Variable Models
Presentation of Additional SLDV Models
To confirm the validity of our results, we report in Table A3 a variety of model
specifications. Models 2 and 3 include a variable for whether a jurisdiction borders Lake Volta,
which is a dichotomous variable valued at 1 if the district touches the lake. (This variable is not avariable for the district being part of Volta Region, which is a formal administrative unit; see
Figure 3 in the main text.) We include the variable for bordering the lake both because we weretold that the program initially targeting island communities, and because there is clearlyclustering in that area. Additionally, Model 4 includes a dummy variable for whether turnout in a
district increased between 2004 and 2008. We include this variable to test whether the directionof turnout matters – whether more or fewer people voted in 2008 than in 2004 in a jurisdiction.
In some models, we dropped population density, since it is clear from the maps that solar PV
panels were not being targeted to major urban areas. In Model 5, we estimated an SLDV modelwith population instead of population density. Due to relatively high collinearity between
population and road density (0.65), we had to manually increase the tolerance for this model toconverge.

4

Another option has been suggested, but it is not currently possible to do. Specifically, spatial econometrics haslong been interested in the possibility of including both a spatially lagged dependent variable and a spatial error term
(Fischer and Getis 2009, p. 381-382). However, this is not possible if the spatial weights matrix is the same for the
spatially lagged dependent variable and the spatial error term (Anselin and Bera 1998).

5
Table A3: Direct Effects of Spatially Lagged Total Solar Projects
Political Variables
H1A: Reward NDC loyalists(2008 NDC vote share)
H1B: Induce swing voters(NDC vote share volatility)
H1C: Mobilize nonvoting incumbent supporters(nonvoters * NDC vote share)

Model 1Model 2Model 3Model 4Model 513.00
(1.04)

6.58(0.55)

7.52(0.66)

12.11(0.93)

15.11(1.30)

47.35(0.98)

48.59(1.05)

51.89(1.18)

40.68(0.85)

45.38(1.05)

123.30*(2.12)

109.20*(2.06)

102.80*(1.95)

103.90(1.80)

109.96*(1.92)

H1C: Mobilize voters where turnout level varies(turnout volatility)
Need Variables
H2A: Electric grid per capita
H2B: Road network density
Control Variables
Ethnic fractionalization-18.23(-1.32)-4952(-1.07)-100.20*(-2.14)-4315(-1.03)-74.32(-1.69)-3926(-0.95)-111.20*(-2.23)-4906(-1.17)-153.10*(-2.83)-4224(-1.02)-114.09*(-2.58)5.97(0.56)-3.11(-0.31)-5.82(-0.54)0.003(1.30)

2.36(0.20)
0.004(1.50)

4.69(0.42)

Population density
Population size
World Bank projects
Percent voters 2008
Poverty rate
Inequality level (Gini)
Number of health facilities
Literacy
Female ratio
Borders Lake Volta-0.27(-0.63)-41.08(-0.81)-0.1(-0.71)0.04(0.14)
0.05(0.59)-16.93(-1.05)-347.90*(-2.38)-0.54(-1.37)-33.05(-0.72)-0.1(-0.75)0.08(0.27)
0.06(0.73)-10.66(-0.69)-303.00*(-2.26)33.29*(4.65)-0.54(-1.39)-35.19(-0.75)-0.09(-0.69)0.13(0.41)-0.01(-0.10)-6.99(-0.45)-316.70*(-2.36)32.91*(4.35)

N170
170* p<0.05
Estimates are mean direct effects; simulated z-statistics are in parentheses

170

-0.14
(-0.32)-63.23(-1.26)-0.08(-0.56)0.03(0.07)-0.05(-0.50)-13.14(-0.79)-356.40*(-2.48)0.00004(1.78)-0.45(-1.09)-45.28(-0.90)-0.09(-0.61)
0.11(0.33)-0.17(-1.14)-13.20(-0.73)-359.24*(-2.60)

170

170

6
In all models, turnout volatility remains statistically significant and positive, and roaddensity is statistically significant and negative, even when the Borders Lake Volta dummy isincluded. This Lake Volta variable is significant when included in models, but the dummy
variable marking when voter turnout increased in a district is not.
OLS Models
Table A4 presents OLS results. Models 3 and 4 are the same as Models 1 and 2 in Table
2 of the main text. Models 1 and 2 in Table A4 exclude control variables. We cluster standard
errors by region. Our OLS results are consistent with our SLDV results, except more variables
are statistically significant. This is to be expected when spatial clustering is not fully taken into
account.

7
Table A4: OLS Analysis of Spatially Lagged Total Solar Projects
Political Variables
H1A: Reward NDC loyalists(2008 NDC vote share)
H1B: Induce swing voters(NDC vote share volatility)

Model 1Model 2Model 3Model 419.55*(9.16)-51.67(75.86)

19.69(12.47)-56.75(70.11)54.30(65.14)67.01(68.84)48.68(73.04)60.39(76.22)

H1C: Mobilize nonvoting incumbent supporters(nonvoters * NDC vote share)
H1C: Mobilize voters where turnout level varies(turnout volatility)

247.28(270.65)

262.97(253.41)133.47*(57.87)115.90*(52.26)145.22*(59.38)128.93*(55.95)-5351.00*(2703.17)-198.85**(72.80)-5297.54*(2671.00)-193.22**(69.69)-5059.93+(2669.37)-171.33*(72.61)-4939.84+(2613.32)-162.41*(69.07)21.92+(13.20)0.003+(0.001)-0.32(0.23)
29.04(43.02)23.93+(14.20)0.003+(0.001)-0.37(0.25)-88.91(114.51)

7.00(12.24)0.004+(0.002)-0.22(0.25)
36.97(50.49)-0.06(0.12)
0.14(0.30)-0.04(0.05)-8.70(18.24)-377.37(237.45)

9.27(12.25)0.004+(0.002)-0.28(0.27)-88.89(106.52)-0.07(0.12)
0.17(0.29)-0.04(0.05)-10.49(18.02)-372.23(231.38)

N
170
170
170
+ p<0.10 ; * p<0.05 ; ** p<0.01
Estimates are mean effects; standard errors clustered by region are in parentheses

170

Need Variables
H2A: Electric grid per capita
H2B: Road network density
Control Variables
Ethnic fractionalization
Population density
World Bank projects
Percent nonvoters 2008
Poverty rate
Inequality level (Gini)
Number of health facilities
Literacy
Female ratio

8Spatial Analysis Using Count Models
Ideal analyses for the data would use models that reflect the fact that our data moreclosely resembles a Poisson distribution than a normal distribution; a count model would be mostappropriate. Given limitations of existing statistical packages, however, it is not possible tocombine global analyses that account for spatial autocorrelation with local analysis using
geographically weighted regression, if count models are used. Specifically, the only available
statistical packages for count models with spatially autocorrelated data employ Bayesian analysis
techniques, whereas GWR models employ frequentist maximum likelihood approaches. Because
these rely on fundamentally different procedures for inference, we do not believe they should be
used together. Moreover, there does not yet exist a reliable statistical package for GWR withcount models (discussed more below).
As a result, the body of our paper presents models that assume normally distributed data
in Table 2, but we present spatial count models here to show that results are not highly sensitive
to model specification. To do so, we employ the relatively recent CARBayes package in R.5
These count models address spatial autocorrelation through the incorporation of vectors of
random effects. This approach is different than that of SLDV models, in which the dependent
variable is lagged to incorporate values of neighboring units on the right-hand side of theequation. These count model approaches tend to rely upon simulation-based techniques for
inference, as opposed to maximum likelihood.
Table A5 shows the results of estimating Poisson models that include a vector of random
effects. Specifically, we fit a multivariate spatial generalized linear mixed model, where therandom effects are modelled by a multivariate conditional autoregressive model.6 The
conditional autoregressive prior was proposed by Leroux et al., so we refer to these as Leroux
models.7 Between variable correlation is captured by a between variable covariance matrix withno fixed structure. Further details are provided in the package vignette.8 For our models, we ran
150,000 Markov Chain Monte Carlo (MCMC) simulations with 50,000 burn-in. The highnumber of simulations was necessary for us to obtain a reasonable acceptance rate for theMarkov chain, although it may also make these models more likely to detect statisticalsignificance.

5

Lee 2013.
Within CARBayes, this is done by the MVS.CARleroux command.
7
Leroux et al. 2000.
8
Lee 2017.
6

9
TABLE A5. Leroux Model Estimated Effects on Total Solar Projects
Districts
Political Variables
H1A: Reward NDC loyalists(2008 NDC vote share)

0.5112(-3.6493 – 8.8148)

H1B: Induce swing voters(NDC vote share volatility)

8.3577(-11.2945 – 32.0298)

H1C: Mobilize voters where turnout level varies(turnout volatility)

24.2897(-16.7662 – 46.6628)

Need Variables
H2A: Electric grid per capita
H2B: Road network density
Control Variables
Ethnic fractionalization
Population density
World Bank projects
Percent voters 2008
Poverty rate
Inequality level (Gini)
Number of health facilities
Literacy
Female Ratio

1013.2836(-1493.8302 – 2523.9297)-143.8257*(-162.7773 – -0.8783)-9.6067(-14.6354 – 5.7254)-0.0036(-0.0073 – 0.0018)-1.2172*(-1.4687 – -0.0487)-35.2699(-52.0530 – 8.3183)
0.2187*(0.0599 – 0.2692)-0.4224*(-0.6377 – -0.0543)
0.2068(-0.0374 – 0.3287)
12.3708(-2.3480 – 20.8424)
86.3596(-26.4562 – 142.3604)

N*p<0.05
Range between 2.5th percentile and 97.5th percentile in parentheses.

170

These results are largely consistent with those in Table A3 in showing the interaction
between when need-based and political factors influence the distribution of solar projects.
Namely, need-based factors consistently are correlated with solar panel distribution. Meanwhile,these models suggest different stories about political factors. In the SLDV models (Table A3),
our measures for need and politics are significant. In the Leroux model (Table A5), however,need variables are significant, but political ones are not. Before discussing these results in greater
detail, we emphasize that the Leroux and SLDV models employ different procedures for
inference (Bayesian MCMC simulations vs. Frequentist maximum likelihood). The substantial

10
differences between these models make it extremely difficult to obtain perfectly consistent
results.
The Leroux model presented in Table A5 finds statistical significance for need variables.
This result corresponds with the stated goals of the program implementers, and provides
suggestive evidence, unlike the SLDV models, that they were able to avoid political capture.
Additional Information on and Results from Geographically Weighted Regression (GWR)
While spatially lagged dependent variable models offer many valuable insights for ouranalysis, they average effects across all units into a single coefficient estimate for each variable.
To test whether differences exist in variable relationships over space, we employ geographicallyweighted regression (GWR) as a sensitivity analysis. GWR has been used for a wide variety of
topics to explore the spatial heterogeneity in relationships, such as voter behavior during the US
Democratic realignment of the 1920s and 1930s,9 campaign contributing and volunteerism,10
regional effects of terrorism on economic growth,11 and even afforestation in Vietnam.12 Here,
we provide additional explanation of GWR analyses, as well as model results.
GWR estimates distinct coefficients for each unit in an analysis, allowing the researcher
to examine the logic of distribution employed in individual jurisdictions. In our case, it allows us
to check whether there appears to be cross-district variation in the logic of public goods
distribution. GWR is critical for our analysis because distributed goods are visibly clustered in aparticular area – the islands and remote villages in and around Lake Volta. GWR allows us todetermine whether the logics at play in this region are consistent in the rest of the country. In
doing so, we also illustrate the utility of using advanced spatial methods to explore and compare
average national phenomena with local-level variation.
Using GWR, locations closer to a given location i receive greater weight than locations
that are further away from location i. Such weight is assigned through a spatial weights matrix.
The resulting equation, and its spatial weights matrix, can be displayed in two equations asfollows:
Equation 1: yi 0(ui, vi) kk(ui, vi)xik iEquation 2: wij exp[–1/2(dij/b)2]
where “ui, vi” denotes the coordinates of the ith point in space and βk(ui, vi) is a realization of thecontinuous function βk(u, v) at point i.13 In the spatial weights matrix, dij is the distance between
locations i and j and b is the bandwidth. Bandwidth reflects the distance-decay of the weighting
function and affects the spatial smoothing of the estimates. For our analysis, we use the corrected
AIC statistic to calculate the optimal distance for our spatial weights.
Estimates from GWR could reveal unit-level estimates that vary substantially. This could
suggest that different logics exist across space. Alternatively, GWR results could reveal that thesame logic holds across space, but that the magnitude of the relationship varies across units.
Because GWR results do not report significance, interpreting them requires knowledge andanalysis of local context.14
9Darmofal 2008.
Cho and Gimpel 2010.
11
Ocal & Yildirim 2010.
12
Clement et al. 2009.
13
Fotheringham, Brunsdon, and Charlton 2002.
14
Clement et al. 2009.
10

11
To interpret GWR results, it is common practice to compare GWR results for each unit ofanalysis (districts in our case) to average marginal effects reported in OLS regressions of thesame variables.15 One can then evaluate the range of variation in a coefficient by simply
comparing the minimum, maximum and various percentile coefficient values for each variable ofinterest, which are estimated in GWR. GWR allows us to inspect whether the average marginal
effect from the OLS regression falls close to the marginal effect of the median data point in
GWR results. OLS average marginal effect values that are far above or below the median effectindicate the possibility of a strong spatial component to the variation in average marginal effect.
In simple terms, if the average marginal effect from OLS regression falls below the 25th
percentile or above the 75th percentile calculated in GWR estimation, some areas of the country
are likely driving the results, and we can examine which areas those are. Table A6 thus shows
OLS results for models with the variables that are in the GWR analysis, alongside minimum,
maximum and several percentile GWR results.
We start with the need-based variable, road density, which was significant in our global
models. The sign of the coefficients is consistent across the OLS result and the interquartile
range, suggesting that the basic logic of distributing projects where existing infrastructure islimited is consistent countrywide. At the same time, however, the OLS result for road density
falls just below the interquartile range, suggesting that the statistical significance in the global
models may be driven by particularly strong relationships in a subset of the units. We explore
this subnational variation qualitatively in the body of the text.
TABLE A6. GWR Results: Estimated Effects on Number of Solar Panels
Electric grid per capita
Road density
NDC vote share 2008
NDC vote share volatility
Turnout volatility

OLS-6477.44-162.81
19.55
36.85
177.21

District Data
Min
-22042.67-782.99-7.72-158.27-7.34

25%-10673.30-153.1213.33-36.4038.98

Median-5111.39-74.22
26.55-14.4277.65

75%-1129.57-23.1429.98
44.45
150.45

Max
921.3662.78
52.27
641.92463.97

Results are largely the same for the measure of political influence that is statisticallysignificant in the SLDV results, turnout volatility. The interquartile range results in Table A6
show that turnout volatility is positively associated with panel distribution in most of Ghana.
However, the mean coefficient on turnout volatility is larger than the coefficients in theinterquartile range. This means that the coefficient estimates for more than 75 percent of districts
falls below the average marginal effect. As with the effect of road density, we can infer thatwhile there is generally a positive relationship between turnout volatility and panel distribution,
particularly strong relationships in a subset of jurisdictions may be driving the statisticalsignificance in the global models.

15

Although it would be ideal to compare GWR results with the results from our global SLDV model, they are notstrictly comparable, because SLDV models account for spatial clustering while GWR models do not. As a result, it
is standard to compare GWR and OLS results. While OLS model results may produce biased coefficients due toviolations of the independent observation assumption, here our OLS results are similar in coefficient sign andmagnitude when compared to our SLDV results.

12
Maps of GWR results—which are presented in Figure 4 in the main text—help us tofurther interpret these results.16 Specifically, we see that across nearly the entire country, therelationship between solar panel distribution and road density is negative, and the relationshipbetween solar panel distribution and turnout volatility is positive. But we also see that there isconsiderable variation in the magnitude of these coefficients across the country for bothvariables.
Here, we also recognize the limitations of GWR. Local collinearities, the lack of a basemodel to explain coefficient variation, repeated use of data to estimate model parameters, andsensitivity to kernel bandwidth can lead to unreliable estimates of coefficient signs, coefficientmagnitudes, and standard error sizes.17 While these concerns make GWR unreliable for
explanatory assessments of statistical significance, GWR is suitable as it is used in our paper,
where the goal is exploratory analysis of variation in coefficient signs and magnitudes.
A Note on GWR Analyses Using Count Models
As explained earlier, local analysis with geographically weighted regression is mostfrequently performed using OLS as the base model. Shifting to Poisson models to reflect the
non-normal distribution of our dependent variable, GWR becomes more complex. The spgwr
package in R is capable of estimating GWR Poisson models, but even the authors of the package
are uncertain about whether scholars should use it. They state, “The use of GWR on GLM isonly at the initial proof of concept stage, nothing should be treated as an accepted method at thisstage”.18 Given this uncertainty, we chose not to use GWR Poisson models.
Hot Spot Analysis
In the main text of the paper, we overlay solar power projects, represented as dots, onmaps of the electric grid and NDC vote share in Figure 1. We also overlay the dots on choropleth
maps of road density and voter turnout volatility in Figure 5. These maps yield valuable
opportunities for descriptive analysis. However, misperceptions from visually examining thesevariables are possible. Because of this, we leverage spatial statistics as a robustness check, which
can be done with hot spot analysis.
Whereas GWR explores spatial variation in the relationship between variables at a
specific unit of analysis, hot spot analysis is a statistical technique that assesses where a variableof interest is particularly clustered or dispersed. This method goes beyond a simple “eyeball test”
of where individual solar projects are located because it is not susceptible to visualmisperceptions. For our discussion, hot spot analysis also cuts through the complexity ofsearching for patterns in the spatial distribution of 1242 solar power projects.
Hot spot analysis has been used to help explain a variety of political phenomena, such as
terrorism.19 A “hot spot” is a geographic area in which the variable of interest – in our case, solar
panel distribution – is clustered. A “cold spot” is therefore where there is a high degree ofdispersion at the local level. We use the Getis-Ord Gi* statistic for our hot spot analysis, as

16

For our maps, we display five gradients. The intervals for these gradients are calculated by ArcGIS using the
Natural Breaks (Jenks) method. As Cho and Gimpel (2010) state, the natural breaks method maximizes within-classhomogeneity and between-class heterogeneity.
17
Wheeler and Tiefelsdorf 2005, Wheeler 2014.
18
Bivand et al. 2017.
19
Nemeth, Mauslein, and Stapley 2014.

13
calculated by ArcGIS.20 This statistic is calculated using the equations in the figure below.21
Since the Gi* statistic is a z-score, we can intuitively interpret the positive and significant values
as hot spots and the negative and significant values as cold spots.
Figure A1: Getis-Ord Statistic Equations

We create the “weighted” points that are necessary for hot spot analysis with the Integratecommand in ArcGIS. With this command, we grouped points that are within 10 kilometers ofeach other. In our analysis, we display hot spots (local clustering, shown as flames) and cold
spots (local dispersion, shown as diamonds) at the 90%, 95%, and 99% confidence levels (Figure
A2). In Figure A2, the highest levels of clustering are marked with dark flames, and the highest
levels of dispersion are marked with dark diamonds. We can obtain additional insights byoverlaying our hot spot analysis on maps of road density and voter turnout volatility.22 These
maps reveal that the main cluster around Lake Volta is the only area in Ghana with both low roaddensity and high voter turnout volatility. This yields a combined motivation to place solar power
projects in this area: high need and high voter mobilization potential.

20

Nemeth, Mauslein, and Stapley 2014.
The figure comes from ArcGIS documentation here: http://pro.arcgis.com/en/pro-app/tool-reference/spatialstatistics/h-how-hot-spot-analysis-getis-ord-gi-spatial-stati.htm
22
It is worth noting that hotspot analysis does not recognize sub-national boundaries; districts are shown in Figure
A2 in order to display descriptive statistics of road density or voter turnout volatility only.
21

14
FIGURE A2. Hotspot Analyses

15
REFERENCES
Anselin, Luc, and Anil K Bera. 1998. "Spatial dependence in linear regression models with an
introduction to spatial econometrics." Statistics Textbooks and Monographs 155:237-290.
Anselin, Luc, Anil K Bera, Raymond Florax, and Mann J Yoon. 1996. "Simple diagnostic tests
for spatial dependence." Regional science and urban economics 26 (1):77-104.
Anselin, Luc, and Serge Rey. 1991. "Properties of tests for spatial dependence in linear
regression models." Geographical analysis 23 (2):112-131.
Ayee, Joseph RA. 2012. "The political economy of the creation of districts in Ghana." Journal of
Asian and African Studies:0021909612464334.
Bivand, Roger, Andrew Bernat, Marilia Carvalho, Yongwan Chun, Carsten Dormann, Stéphane
Dray, Rein Halbersma, Nicholas Lewin-Koh, Jielai Ma, and Giovanni Millo. 2005. "The
spdep package." Comprehensive R Archive Network, Version:05-83.
Bivand, Roger, Danlin Yu, Tomoki Nakaya, Miguel-Angel Garcia-Lopez, and Roger Bivand.
2017. "Package ‘spgwr’." R software package.
Cho, Wendy K. Tam, and James G. Gimpel. 2010. "Rough Terrain: Spatial Variation in
Campaign Contributing and Volunteerism." American Journal of Political Science 54(1):74-89. doi: 10.1111/j.1540-5907.2009.00419.x.
Clement, Floriane, Didier Orange, Meredith Williams, Corinne Mulley, and Michael Epprecht.
2009. "Drivers of afforestation in Northern Vietnam: Assessing local variations using
geographically weighted regression." Applied Geography 29 (4):561-576. doi:
http://dx.doi.org/10.1016/j.apgeog.2009.01.003.
Darmofal, David. 2008. "The political geography of the new deal realignment." American
Politics Research. 36 (6): 934-961.
Elhorst, J Paul (2010) Applied spatial econometrics: raising the bar. Spatial economic analysis
5(1): 9-28.
Fischer, Manfred M, and Arthur Getis. 2009. Handbook of applied spatial analysis: softwaretools, methods and applications: Springer Science & Business Media.
Fotheringham, AS, C Brunsdon, and Martin Charlton. 2002. Geographically weighted
regression: the analysis of spatially varying relationships. West Sussex, UK: John Wiley.
Lee, Duncan. 2013. "CARBayes: an R package for Bayesian spatial modeling with conditional
autoregressive priors." Journal of Statistical Software 55 (13):1-24.
Lee, Duncan. 2017. "CARBayes version 5.0: An R Package for Spatial Areal Unit Modellingwith Conditional Autoregressive Priors." https://cran.rproject.org/web/packages/CARBayes/vignettes/CARBayes.pdf.
Leroux, Brian G., Xingye Lei, Norman Breslow, M Halloran, and Berry Donald Elizabeth. 2000.
"Estimation of disease rates in small areas: a new mixed model for spatial dependence."
Statistical models in epidemiology, the environment, and clinical trials:179-191.
Nemeth, Stephen C, Jacob A Mauslein, and Craig Stapley. 2014. "The primacy of the local:
Identifying terrorist hot spots using geographic information systems." The Journal of
Politics 76 (02):304-317.
Öcal, Nadir, and Jülide Yildirim. 2010. “Regional effects of terrorism on economic growth inTurkey: A geographically weighted regression approach.” Journal of Peace Research 47(4): 477-489.
Wheeler, David C. 2014. "Geographically weighted regression." In Handbook of Regional
Science, 1435-1459. Springer.

16
Wheeler, David, and Michael Tiefelsdorf. 2005. "Multicollinearity and correlation among local
regression coefficients in geographically weighted regression." Journal of Geographical
Systems 7 (2):161-187.