- Open Access
Spatial analysis of cluster randomised trials: a systematic review of analysis methods
Emerging Themes in Epidemiologyvolume 14, Article number: 12 (2017)
Cluster randomised trials (CRTs) often use geographical areas as the unit of randomisation, however explicit consideration of the location and spatial distribution of observations is rare. In many trials, the location of participants will have little importance, however in some, especially against infectious diseases, spillover effects due to participants being located close together may affect trial results. This review aims to identify spatial analysis methods used in CRTs and improve understanding of the impact of spatial effects on trial results.
A systematic review of CRTs containing spatial methods, defined as a method that accounts for the structure, location, or relative distances between observations. We searched three sources: Ovid/Medline, Pubmed, and Web of Science databases. Spatial methods were categorised and details of the impact of spatial effects on trial results recorded.
We identified ten papers which met the inclusion criteria, comprising thirteen trials. We found that existing approaches fell into two categories; spatial variables and spatial modelling. The spatial variable approach was most common and involved standard statistical analysis of distance measurements. Spatial modelling is a more sophisticated approach which incorporates the spatial structure of the data within a random effects model. Studies tended to demonstrate the importance of accounting for location and distribution of observations in estimating unbiased effects.
There have been a few attempts to control and estimate spatial effects within the context of human CRTs, but our overall understanding is limited. Although spatial effects may bias trial results, their consideration was usually a supplementary, rather than primary analysis. Further work is required to evaluate and develop the spatial methodologies relevant to a range of CRTs.
Randomised controlled trials assess the efficacy and safety of interventions [1, 2]. When it is difficult to allocate interventions at the individual level, for example due to logistical or financial restrictions, randomisation and allocation of interventions at a group level may be preferred, this is a cluster randomised trial (CRT) . CRTs also allow for estimation of spillover and herd effects; the apparent treatment effect on individuals who do not receive the intervention . Failure to account for spillover effects can result in biases that reduce the quality of trials and mean that absence of bias is no longer guaranteed by the randomisation, especially when the relationship between intervention and outcome is complex [5, 6].
Spatial effects are effects stemming from locational variation in the distribution of phenomena of interest or in the intensity of interaction between phenomena of interest. Such effects manifest as local variation in the estimated treatment effect over a study area. The existence of spatial effects suggests that global effects estimates are uncertain, and may under- or over-estimate the true effect dependent on location. When values of a single variable are related to nearby values of the same variable this is called spatial dependence, the existence of which is usually captured using spatial autocorrelation measures . Spatial dependence is a fundamental concept in spatial statistics  and stems from Tobler’s  1st law of geography that “everything is related to everything else, but near things are more related than distant things.” In agricultural field trials it is long established that the location of the data can impact on trial results . Incorporation of spatial methodology in agricultural trials is common,  but the impact of spatial effects in human CRTs have not been researched extensively.
Clusters in CRTs are often defined geographically and valid inference relies on the assumption that the clusters are independent irrespective of their nearness to one another . There is frequently an assumption of an absence of spillover; that movement of people and diseases occurs freely within a cluster but movement between clusters is negligible, non-existent, or not relevant . This assumption can be violated when there is movement of people or diseases across borders, such as mosquitoes flying between control and intervention households. If an intervention such as insecticide-treated bed nets provide a protective effect to nearby control households then ignoring mosquito mobility will result in underestimating the intervention effect of the trial. This could result in trials discarding effective interventions because the control and intervention are both receiving the benefit of the treatment.
Spillover may also be due to connections in social networks,  also violating the assumption of independence. In this paper, we consider spillover that can be estimated using GPS data which is often collected as part of trials and therefore do not consider social networks. Spillovers are more likely in trials that have spatially close clusters, and the effect is especially important when an individual’s outcome is affected by their proximity to other individuals with different exposure statuses.
One way to minimise the potential for spillover is to design a trial with well-separated clusters. In practice this may not be logistically or financially feasible as spatial effects can be present over distances of several kilometres . Furthermore, greater distances between clusters removes our ability to measure spillover spatial effects. To be able to measure spatial effects we need to have data on people nearby to one another and by separating the clusters we may no longer have such information. If proximity to an intervention affects non-treated individuals this is of usually scientific interest and something we should measure. There is a need to control for and estimate spatial effects in the analysis of a trial without adding extra complexity to the design.
We therefore focus only on spatial analysis methods used in CRTs in this review and do not consider alternative trial designs. A further reason for focusing on analysis methods is that this may enable analysis of existing and previous CRTs where redesign is not possible. This review is a diagnostic review of spatial analysis methods that have been used in CRTs. As such, it does not attempt to pose statistical solutions or determine the best way to account for spatial effects within CRTs. We will describe the state of the literature and aim to improve understanding and help inform further research into spatial effects within CRTs by (1) Identifying spatial analysis methods used in CRTs. (2) Summarising and grouping spatial methods (3) Assessing the impact of spatial effects.
Search terms and review process
The PRISMA guidelines  for systematic reviews and meta-analyses were followed for this review. It was conducted between January 2016 and September 2016. Ovid/Medline, Pubmed, and Web of Science databases were electronically searched and Mendeley was used to store articles. Search terms for CRTs and spatial effects are detailed in Table 1. Studies up to end of 2015 were included and only English language articles and search terms were considered.
Papers from each database were combined into a single spreadsheet containing the title, authors, journal, and year. Duplicates were removed automatically within the software and then manually during the title screen. The titles were screened to remove irrelevant papers such as individually randomised trials. Following this, abstracts were screened and the full texts of potentially relevant papers were independently reviewed by two reviewers and disagreements resolved. After selecting relevant papers the references of the articles were screened.
Inclusion and exclusion criteria
The inclusion criteria for studies were: (1) the study is a cluster randomised trial. (2) Spatial methods are used in the analysis of the study. We categorised a spatial method as one which accounts for the structure, location, or relative distances of the data. This includes direct estimation of an effect such as, the change in risk for those within 100 m of an intervention household or the use of a spatial model which account for spatial structure.
The exclusion criteria were (1) non-randomised studies (2) individual randomised trials and hybrid CRTs such as Double or pseudo-randomised studies as they were considered not to be cluster randomised trials (3) grey literature (4) studies where spillover effects are measured in a non-spatial way, for example comparison of vaccinated and non-vaccinated individuals within an intervention cluster. (5) Studies that account for spatial effects at the design stage only; for instance, using buffer zones or well-separated clusters (6) articles which were study protocols and therefore had not applied their methods yet.
The following variables were collected on each paper: title, year, journal, author, intervention, outcome, whether a map was presented, spatial analysis method.
A flow chart of the search process can be seen in Fig. 1 and the search terms in Table 1. The search terms returned 6997 records, reducing to 571 records after the title screen and duplicate removal. Of the 571 records, 40 abstracts were considered relevant for full text review by the reviewers. One study  was a replication analysis and it was decided to include the original study in the review instead of the replication. There are ten papers and thirteen trials in this review as some papers include multiple trials.
Whilst this review was being conducted a systematic review on health-related spillover in impact evaluations was released . It has the more general aim of attempting to summarise methods to estimate health related spillover in low and middle income countries. Our review is different as it only includes spatial methods used within cluster randomised trials and does not restrict by type of country. The results from both reviews were compared and did not lead to additional records being included.
This review contains ten papers published between 1998 and 2015, they relate to thirteen trials as some papers contained more than one trial. The trials took place around the globe with three taking place in Kenya, three in the United Kingdom, and the rest in Mexico, Venezuela, Ghana, Papua New Guinea, Vietnam, Haiti, and India. There is one stepped wedged trial  and all others are parallel cluster randomised trials. Six of the papers were a spatial reanalysis of a previously reported trial.
Seven trials focused on pathogens carried by mosquitoes. The intervention for six of these was insecticide treated bed nets or curtains and the other intervention was a drug. Two of the seven mosquito trials considered all cause child mortality as an endpoint, the other five looked at entomological endpoints.
There were two vaccine trials, the first looked at vaccine uptake in response to a mass campaign and the second evaluated vaccine effectiveness for a typhoid vaccine. One paper considered primary care and community based trials, within this they applied spatial methods to three different trials, one simulated and two real trials. The final trial looked at the impact of deworming on education and health within schools. Further details of the trails in the review are in Table 2 and in Additional file 1: Table S1.
Five trials used a two stage method in their primary analysis, the rest a one stage method . Seven papers included a map of the study location. The software used included Fortran, SAS, ArcGIS, HLM, Stata, MLWin, R, and WinBugs, four studies did not report software used. There were two methodology papers (four trials) that focused on how to adjust for spatial effects, and eight applied papers (nine trials) where analysis of the spatial effect was not the primary aim of their trial.
The studies took two approaches for analysing spatial effects, referred to in this paper as spatial variables and spatial models, expanded upon in Fig. 2. Nine of the trials analysed the effect of a spatial variable which is a measurement that relates to where the observations are located. The two types of spatial variables found in this review were straight line distance such as distance between participants, and density, for instance the number of treated participants within a 100 m radius. Four trials used spatial models by including the spatial structure of the participants using random effects statistical models. The spatial models are specified by measuring how participants are connected to one another, for example recording participants who are neighbours, described further in Fig. 2. Incorporating spatial structure into a model this way treats a spatial effect as an underlying unobserved process which may not be directly measurable. The two approaches make different assumptions about the type of spatial relationship and a variety of methods were used for each approach.
Straight line distance
Five trials [18,19,20,21,22] estimate spatial effects by measuring the distance between participants and a location of interest. In these studies, the location is either another participant or a feature which may affect the outcome, such as a health facility. Several studies analysed the effect of distance to more than one type of location.
The distance between each control participant to their nearest intervention participant was analysed in three trials [18, 19, 22], termed distance to nearest discordant observation. They also analysed proximity to nearest reservoir or health facility. Distances were categorised and the effect measured for each category. For example, Binka et al.  calculated a standardised mortality rate at five separate distance categories.
Kroeger et al.  considered whether distance to a participant with the outcome at the beginning of the study affects the odds of having the outcome at the end of the study. They tested at four separate distances and corrected for multiple testing. Ali et al.  included distance to school and nearest hospital in a model which assessed the intervention effect and accounted for cluster effects. This was the only trial that included straight-line distance in the primary analysis of their trial.
Four trials [21, 23,24,25] analysed the effect of density in the area surrounding the participants. They analysed the density of factors that may affect the outcome, for example the number of people vaccinated within 100 m. The methods differ by whether they used a count or a proportion and whether they focused on the treatment density or the risk of infection from surrounding individuals. Including density as a spatial variable assumes that number of objects within a certain distance is important as well as the distance to the nearest object.
Lenhart et al.  measured intervention density as the number of households with bed nets within 100 m of an observation. The study assessed spatial spillover through the correlation of change in baseline of outcome measure with number of bed net households within 100 m. In contrast, Miguel and Kremer  measured density as the proportion of children treated within 6 km of a school as well as the total number and accounted for this in their primary analysis. Ali et al.  included a proportion and count density in their primary analysis model by accounting for typhoid prevalence and number of private practitioners for the neighbourhoods of each participant.
Chao et al.  differs from the previous applied papers because they develop and test a new method to deal with spatial effects in CRTs. They define a variable called ‘potential exposure’ which is ‘the sum of the relative risks of all who live within 100 m of each person.’ . The ‘potential exposure’ controls for the spatial variation in risk surrounding an observation. They demonstrate that this variable can be used to account for spatially heterogeneous risk factors in the primary analysis of a trial.
The methods presented in the previous section assume that the underlying spatial process can be measured and spatial effects can be estimated from this measure. The remaining two papers include four trials that model spatial effects using a spatially structured random effect. This approach makes fewer assumptions about the mechanism of spatial process and allows for a range of local and global dependency structures .
Alexander et al.  investigated the spatial pattern of mosquito borne vectors. Adapted from a previous paper  they incorporated a distance parameter in the covariance matrix of a random effect within a negative binomial model . This distance parameter allows participants who are closer together to be more similar than participants that are further apart as shown in Fig. 2. They also estimate a distance decaying parameter in the covariance structure of the random effect, which they use to estimate the “half distance” which is the distance at which spatial correlation halves.
Silcocks and Kendrick  applied several types of spatial models to two primary care trials and one simulated trial. The model was a variation of a Besag, York and Mollie model  which contains a spatially structured random effect and a random effect for cluster. The spatial structure of the participants was represented using a spatial weights matrix and is included in the covariance of the random effect further described in Fig. 2. In primary care or community based trials people may reside in one area and receive treatment in another . Having a spatial and non-spatial random effect allowed participants to have membership to multiple clusters which they called a multiple membership model. They also consider a fixed north/south and east/west gradient covariate in their model evaluations.
Impact of spatial effects
All thirteen trials found evidence of a spatial effect within their studies. Seven trials report a protective spillover effect for participants who live close to an intervention. There is evidence that adjusting for spatial effects affects the precision and value of the estimated intervention effect. Chao et al.  saw that adjustment for spatial effects lowered the effect estimate of the intervention. The precision and intervention estimate changed in the three trials analysed by Silcocks and Kendrick . The study demonstrated that spatial models fitted better than a standard CRT random effect model by comparing the Akaike information criteria of the models. They are explicit that this is just illustrative but both studies conclude that spatial effects may need to be adjusted for in CRTs and that further research into methods is required. Despite this only two of the nine applied trials adjusted for spatial effects in the primary analysis of their trial.
This review has found multiple approaches to incorporating or measuring spatial effects in the context of CRTs, however these stem from only a few examples in the literature. Further, no conventional or standard approach was found. The approaches differ by whether they directly analyse a spatial variable or model the spatial structure. Spatial variables were either straight-line distance from a participant to a place of interest or a measure of the density surrounding a participant. For instance, distance between a control participant and an intervention participant. Spatial models included spatial structure in the covariance of the random effects model using a distance parameter or spatial weights matrix. Accounting for spatial structure affects both the precision and point estimates of treatment effects and failure to do so could give inaccurate results [25, 29].
The papers in this review are only a small proportion of the total number of CRTs that have been published. That only ten records were found suggests that spatial effects are not often considered in this area. Furthermore, despite evidence of spatial effects, they were rarely adjusted for in the primary analysis of the trial. It appears that the impact of location on analytical approach is at best an afterthought and in most cases ignored. The trials come from a variety of domains and although predominantly focused in infectious diseases, there may be implications for a broader range of trials particularly in trials of health services organisation  which are becoming more common . Therefore, a wide range of trials may not be accounting for spatial effects that bias results, however it is presently unclear as to what extent this may be an issue.
Further research is required to determine how much spatial effects impact trial results. Simulation studies may allow exploration of how the magnitude and extent of spatial autocorrelation may bias trial results. There have been some attempts to quantify how important spatial effects may be in trials more generally [32, 33]. Methods that allow for estimation of treatment effects whilst accounting for spatial effects could be investigated and also testing under simulation. However, there are several challenges to overcome such as whether we can estimate the true randomised intervention effect using a spatial model, and how such an effect estimate should be interpreted.
An alternative use of spatial data is to conduct additional analyses to complement analysis of the main trial. These analyses explore spatial component of the trial, and could improve understanding of the mechanism of the intervention effect. Spatial methods could be applied to previous trial data and a toolbox of standard approaches defined allowing future trials to predefine spatial analyses. This would allow for quantification of the distance over which spatial effects are present for different disease areas. A further area of research is in alternative CRT designs, although they are not the focus of this review this might include double randomisation or pseudo randomised trials where clusters are first randomised and subsequently individuals within clusters are randomised to allow for measurement of spillover effects [34, 35].
Multiple terms in the literature refer to spillover effects and many of them could refer to spatial effects with differing terminology between fields and researchers in the same fields. To add to confusion, these terms can have dual meanings for instance, an indirect effect can be the effect on an individual who does not receive the intervention or the effect of an intervention through a mediating variable. The search strategy attempted to include a broad range of terms for CRTs and spatial effects but at present, there is no established standard for citing the use of spatial data in the analysis of CRTs. Consequently, it is possible that trials have been missed. Although this is a weakness, comparison with a larger more general systematic review on health-related spillover in impact evaluations  did not result in the addition of any further trials. They had searched 19 databases and screened more than 34,000 records. Due to this the authors conclude the omission of further studies is likely to be minimal.
This review has focused only on the analysis stage of a CRT and it could be argued that adjusting for spatial effects is not necessary in a well-designed trial as clusters should be well-separated to minimise spillover effects . Trial designs such as the fried egg design  which incorporates a buffer around clusters could be used to attempt to eliminate or measure spatial spillover. However, in cases where spatial correlation is present over large distances  this may not be possible and could lead to the inability to detect the difference between no effect and everyone having an effect . Additionally, as spatial effects are rarely considered in trials, it may not be until after the design stage that the problem becomes apparent, if at all.
On the other hand, having clusters relatively close together could have advantages, because the measurement of spatial spillover effects is of scientific interest. Knowledge that an intervention provides indirect benefit based on proximity is useful to differentiate between interventions and to plan how to benefit the largest number of people. Furthermore, subjects who live further apart are more likely to be heterogeneous than those living close together due to cultural, geographical, and social differences which could make treatment differences harder to distinguish due to imbalance between clusters.
In conclusion although there have been a few attempts to control and estimate spatial effects within the context of human CRTs our understanding is limited. Although there are commonalities between approaches there is no consensus on how to account for spatial effects with in CRTs and more work needs to be done to evaluate and develop spatial methodology within the context of a range of CRTs.
cluster randomised trial
Pocock SJ. Clinical trials: a practical approach. West Sussex: Wiley; 1983.
Smith PG, Morrow RH, Ross DA, editors. Field trials of health interventions: a toolbox. 3rd ed. Oxford: Oxford University Press; 2015.
Hayes RJ, Moulton LH. Cluster randomised trials. Boca Raton: Chapman & Hall/CRC; 2009.
Hayes RJ, Alexander ND, Bennett S, Cousens SN. Design and analysis issues in cluster-randomized trials of interventions against infectious diseases. Stat Methods Med Res. 2000;9:95–116.
Hahn S, Puffer S, Torgerson DJ, Watson J. Methodological bias in cluster randomised trials. BMC Med Res Methodol. 2005;5:10.
English M, Schellenberg J, Todd J. Assessing health system interventions: key points when considering the value of randomization. Bull World Health Organ. 2011;89:907–12.
Getis A. A history of the concept of spatial autocorrelation: a geographer’s perspective. Geogr Anal. 2008;40:297–309.
Cressie N. Statistics for spatial data. Revised ed. New York: Wiley; 1991.
Tobler AWR. A computer movie simulation urban growth in Detroit region. Econ Geogr. 1970;46:234–40.
Wiebe GA. Variation and correlation in grain yield among 1,500 wheat nursery plots. J Agric Res. 1935;50:331–57.
Besag J, Kempton R. Statistical analysis of field experiments using neighbouring plots. Biometrics. 1986;42:231–51.
Lines J, Kleinschmidt I. Combining malaria vector control interventions: some trial design issues. Pathog Glob Health. 2013;107:1–4.
Benjamin-chung J, Berger D, Clark A, Falcao L, Tran D, Hubbard A, et al. The identification and measurement of health-related spillovers in impact evaluations : a systematic review. 3ie Syst Rev 22; 2015.
Howard SC, Omumbo J, Nevill C, Some ES, Donnelly CA, Snow RW. Evidence for a mass community effect of insecticide-treated bednets on the incidence of malaria on the Kenyan coast. Trans R Soc Trop Med Hyg. 2000;94:357–60.
Moher D, Liberati A, Tetzlaff J, Altman DG, Grp P. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. Phys Ther. 2009;89:873–80 (Reprinted from Annals of Internal Medicine).
Aiken AM, Davey C, Hargreaves JR, Hayes RJ. Re-analysis of health and educational impacts of a school-based deworming programme in western Kenya: a pure replication. Int J Epidemiol. 2015;44:1572–80.
Hemming K, Haines TP, Chilton PJ, Girling AJ, Lilford RJ. The stepped wedge cluster randomised trial: rationale, design, analysis, and reporting. BMJ. 2015;350:h391.
Hawley WA, Phillips-Howard PA, Ter Kuile FO, Terlouw DJ, Vulule JM, Ombok M, et al. Community-wide effects of permethrin-treated bed nets on child mortality and malaria morbidity in western Kenya. Am J Trop Med Hyg. 2003;68:121–7.
Gimnig JE, Kolczak MS, Hightower AW, Vulule JM, Schoute E, Kamau L, et al. Effect of permethrin-treated bed nets on the spatial distribution of malaria vectors in western Kenya. Am J Trop Med Hyg. 2003;68:115–20.
Kroeger A, Lenhart A, Ochoa M, Villegas E, Levy M, Alexander N, et al. Effective control of dengue vectors with curtains and water container covers treated with insecticide in Mexico and Venezuela: cluster randomised trials. BMJ. 2006;332:1247–52.
Ali M, Thiem VD, Park JK, Ochiai RL, Canh DG, Danovaro-Holliday MC, et al. Geographic analysis of vaccine uptake in a cluster-randomized controlled trial in Hue, Vietnam. Heal Place. 2007;13:577–87.
Binka FN, Indome F, Smith T. Impact of spatial distribution of permethrin-impregnated bed nets on child mortality in rural northern Ghana. Am J Trop Med Hyg. 1998;59:80–5.
Miguel E, Kremer M. Worms: identifying impacts on education and health in the presence of treatment externalities. Econometrica. 2004;72:159–217.
Lenhart A, Orelus N, Maskill R, Alexander N, Streit T, McCall PJ. Insecticide-treated bednets to control dengue vectors: preliminary evidence from a controlled trial in Haiti. Trop Med Int Heal. 2008;13:56–67.
Chao DL, Park JK, Marks F, Ochiai RL, Longini IM, Halloran ME. The contribution of neighbours to an individual’s risk of typhoid outcome. Epidemiol Infect. 2015;143:3520–7.
Banerjee S, Carlin BP, Gelfrand AE. Hierarchical modeling and analysis for spatial data. Boca Raton: Chapman & Hall/CRC; 2015.
Alexander ND, Moyeed RA, Hyun PJ, Dimber ZB, Bockarie MJ, Stander J, et al. Spatial variation of Anopheles-transmitted Wuchereria bancrofti and Plasmodium falciparum infection densities in Papua New Guinea. Filaria J. 2003;2:14.
Alexander N, Moyeed R, Stander J. Spatial modelling of individual-level parasite counts using the negative binomial distribution. Biostatistics. 2000;1:453–63.
Silcocks P, Kendrick D. Spatial effects should be allowed for in primary care and other community-based cluster RCTS. Trials. 2010;11:55.
Besag J, York J, Mollié A. Bayesian image restoration, with two applications in spatial statistics. Ann Inst Stat Math. 1991;43:1–20.
Eldridge SM, Ashby D, Feder GS, Rudnicka AR, Ukoumunne OC. Lessons for cluster randomized trials in the twenty-first century: a systematic review of trials in primary care. Clin Trials. 2004;1:80–90.
Diamond R, Imbens GW, Diamond R. Clustering, spatial correlations and randomization inference. 2010.
Baylis K, Ham A. How important is spatial correlation in randomized controlled trials? How important is spatial correlation in randomized controlled trials? 2015.
Sinclair B, Mcconnell M, Green DP, Sinclair B, Green DP. Detecting spillover effects: design multilevel experiments and analysis of Chicago. Am J Polit Sci. 2014;56:1055–69.
Baird S, Bohren JA, Mcintosh C, Ozler B, Bank W, Org B. Designing experiments to measure spillover effects. Policy Res. Work. Pap. No. 6824. 2016;.
Rosenbaum PR. Interference between units in randomized experiments. J Am Stat Assoc. 2007;102:191–200.
CJ, GLDT, and DL were involved in defining the search strategy. CJ and GLDT performed the search and screening of articles. CJ drafted the article with feedback, input and guidance from GLDT, DL, NA, and WJE. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
Consent for publication
Ethics approval and consent to participate
This work was supported by the Medical Research Council Trials Methodology Hub [Grant Number MR/L004933/1- Q42]. N Alexander receives salary support from the MRC Tropical Epidemiology Group [Award Reference Number MR/K012126/1]. This award is jointly funded by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement and is also part of the EDCTP2 programme supported by the European Union.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.