Analyzing Local Spatio-Temporal Patterns of Police Calls-for-Service Using Bayesian Integrated Nested Laplace Approximation

This research investigates spatio-temporal patterns of police calls-for-service in the Region of Waterloo, Canada, at a fine spatial and temporal resolution. Modeling was implemented via Bayesian Integrated Nested Laplace Approximation (INLA). Temporal patterns for two-hour time periods, spatial patterns at the small-area scale, and space-time interaction (i.e., unusual departures from overall spatial and temporal patterns) were estimated. Temporally, calls-for-service were found to be lowest in the early morning (02:00–03:59) and highest in the evening (20:00–21:59), while high levels of calls-for-service were spatially located in central business areas and in areas characterized by major roadways, universities, and shopping centres. Space-time interaction was observed to be geographically dispersed during daytime hours but concentrated in central business areas during evening hours. Interpreted through the routine activity theory, results are discussed with respect to law enforcement resource demand and allocation, and the advantages of modeling spatio-temporal datasets with Bayesian INLA methods are highlighted.


Introduction
In law enforcement, the implementation of computer-aided dispatch systems that automatically record location-and time-specific information for police calls-for-service has facilitated the storage, retrieval, and analysis of large volumes of spatio-temporal data [1][2][3][4][5].Paired with technological advances in geographic information systems, crime mapping, or the analysis and visualization of historical spatial and temporal police call-for-service patterns, has become an important part of police operations and strategy [6].Police calls-for-service encompass reported crimes, emergency response, traffic management, and other police-facilitated services [6], and compared to official crime data, provide a more comprehensive indicator of overall police resource demand.Increasingly, call-for-service data is available at fine spatial and temporal resolutions, enabling detailed analysis of local space-time variation.
Past crime mapping research at the small-area scale often analyzes spatial and temporal patterns separately.Purely spatial analyses identify hotspots, or groups of small areas that exhibit clustering or high-spatial autocorrelation [7].Assuming a null hypothesis of spatial randomness (i.e., no spatial autocorrelation) a variety of statistical tests have been applied to identify spatial clusters probability distributions via numerical integration rather than an MCMC-based iterative process [28,43].Past research comparing MCMC and INLA approaches has found that INLA reduces computational time and retains reliable parameter estimates [28].Computational efficiency is important when analyzing spatio-temporal data because various prior distributions and assumptions for spatial, temporal, and space-time random effects parameters can be tested and compared [44].
This research analyzes spatio-temporal patterns of police calls-for-service at the small-area level in the Region of Waterloo, Canada.Calls-for-service for one year (2011) were temporally aggregated to 12 time periods of 2 hours, representing a 24 hour day, and geographically aggregated to 755 census dissemination areas, resulting in a total of 9060 space-time units.Aggregated calls-for-service were analyzed to provide broad insight into the space-time patterns of law enforcement resources and because this provides sufficient count data to capture meaningful spatial, temporal, and space-time patterns at this resolution.Hourly patterns were analyzed because crime and police resources have been shown to vary by hour of day more than any other predictor [45] and anticipated spatial and temporal autocorrelation can be appropriately accounted for via CAR processes.
This paper first reviews the study region and call-for-service data.Next, the spatio-temporal model is detailed and four types of space-time interactions are described, each of which holds specific assumptions regarding the spatial and/or temporal structure of residual space-time patterns.Overall spatial and overall temporal patterns, as well as space-time interaction and space-time hotspots, are visualized, interpreted through the routine activity theory, and discussed in the contexts of understanding and informing police resource allocation.In conclusion, we reflect on the use of Bayesian spatio-temporal models and INLA to model spatio-temporal datasets.

Study Region
The Regional Municipality of Waterloo, Ontario, Canada, is composed of the cities of Waterloo, Kitchener, and Cambridge, as well as four rural townships.In 2011, the region had a population of 506,107 distributed across 755 census dissemination areas (DAs).For reference, DAs are the smallest areal units that cover the entirety of Canada and are delineated such that residential populations are between 400 and 700 [46].

Police Call-for-Service Data
Call-for-service data were provided by the Waterloo Regional Police Service for 1 January 2011 to 31 December 2011.Calls-for-service include violent crimes (count = 2696), property crimes (13,704), disorder (18,949), bylaw complaints (10,598), motor-vehicle-related calls (31,105), 911 calls (65,617), and a variety of other police-related services [47,48].Locations reflect the closest intersection to call location and calls were summed to the DA-scale.Descriptive statistics for total calls-for-service are shown in Table 1.Expected counts were calculated as the product of study region rate and DA population (= (total call-for-service count/total population) x DA population) to reflect underlying population variation.In general, areas with high total counts were located in the centre of the study region, closely corresponding to the central business areas of Waterloo, Kitchener, and Cambridge, as well as in areas towards the periphery of the study region in the east and south (Figure 1A).Compared to observed counts, expected counts were lower in central business areas but higher in areas adjacent to central business areas (i.e., areas that have larger residential populations) (Figure 1B).Temporally, calls-for-service were categorized into 12 time periods, each representing 2 hours of a 24 hour day (Table 2).Calls were highest during afternoon and evening hours and lowest in the early morning (Figure 2).In total, there were 290,027 calls-for-service distributed over 9,060 spacetime units (= 755 DAs × 12 time periods).The spatial patterns of calls-for-service are shown in Figure 3 for 04:00-05:59 (Figure 3A), the time period with the lowest observed count, and for 14:00-15:59 (Figure 3B), the time period with the highest observed count.Compared to 14:00-15:59, there were fewer areas with high counts (>50) during 04:00-05:59 and high-count areas were located closer to the main transportation corridor connecting Waterloo, Kitchener, and Cambridge.From 14:00 to 15:59, there were many more areas with high call-for-service counts and these areas were geographically dispersed throughout the study region.Temporally, calls-for-service were categorized into 12 time periods, each representing 2 hours of a 24 hour day (Table 2).Calls were highest during afternoon and evening hours and lowest in the early morning (Figure 2).In total, there were 290,027 calls-for-service distributed over 9060 space-time units (= 755 DAs × 12 time periods).The spatial patterns of calls-for-service are shown in Figure 3 for 04:00-05:59 (Figure 3A), the time period with the lowest observed count, and for 14:00-15:59 (Figure 3B), the time period with the highest observed count.Compared to 14:00-15:59, there were fewer areas with high counts (>50) during 04:00-05:59 and high-count areas were located closer to the main transportation corridor connecting Waterloo, Kitchener, and Cambridge.From 14:00 to 15:59, there were many more areas with high call-for-service counts and these areas were geographically dispersed throughout the study region.

Spatio-Temporal Modeling
Bayesian spatial and spatio-temporal models are often hierarchical with three levels [39].In this context, hierarchical is used to describe the structure of Bayesian models: the first level models the observed data with a likelihood function (probability distribution) and specifies a process model that includes spatial, temporal, and space-time random effects; the second level assigns prior distributions for model parameters; and the third level specifies hyperpriors for parameters of second-level prior distributions (e.g., variances).Prior distributions quantify researcher uncertainty in parameter estimates and impose spatial and/or temporal structure for random effects parameters [30,49].
At the first level, call-for-service count for small-area i (1, …, 755) and time period t (1, …, 12) is represented by Oit and sampled from a Poisson distribution (Equation ( 1)).The Poisson mean, μit, is the product of expected call-for-service count (Eit) and relative risk (rit) (Equation ( 2)).Poisson models are common in spatial and spatio-temporal analysis of count data at the small-area scale [50].

Oit~Poisson(μit)
(1) Equation ( 3) decomposes relative risk into overall risk for the study region (α) and spatial, temporal, and space-time random effects.This is the general nonparametric model structure used for space-time analysis [44].The parametric structure models space-time interaction by incorporating a

Spatio-Temporal Modeling
Bayesian spatial and spatio-temporal models are often hierarchical with three levels [39].In this context, hierarchical is used to describe the structure of Bayesian models: the first level models the observed data with a likelihood function (probability distribution) and specifies a process model that includes spatial, temporal, and space-time random effects; the second level assigns prior distributions for model parameters; and the third level specifies hyperpriors for parameters of second-level prior distributions (e.g., variances).Prior distributions quantify researcher uncertainty in parameter estimates and impose spatial and/or temporal structure for random effects parameters [30,49].
At the first level, call-for-service count for small-area i (1, …, 755) and time period t (1, …, 12) is represented by Oit and sampled from a Poisson distribution (Equation ( 1)).The Poisson mean, μit, is the product of expected call-for-service count (Eit) and relative risk (rit) (Equation ( 2)).Poisson models are common in spatial and spatio-temporal analysis of count data at the small-area scale [50].

Oit~Poisson(μit)
(1) Equation ( 3) decomposes relative risk into overall risk for the study region (α) and spatial, temporal, and space-time random effects.This is the general nonparametric model structure used for space-time analysis [44].The parametric structure models space-time interaction by incorporating a

Spatio-Temporal Modeling
Bayesian spatial and spatio-temporal models are often hierarchical with three levels [39].In this context, hierarchical is used to describe the structure of Bayesian models: the first level models the observed data with a likelihood function (probability distribution) and specifies a process model that includes spatial, temporal, and space-time random effects; the second level assigns prior distributions for model parameters; and the third level specifies hyperpriors for parameters of second-level prior distributions (e.g., variances).Prior distributions quantify researcher uncertainty in parameter estimates and impose spatial and/or temporal structure for random effects parameters [30,49].
At the first level, call-for-service count for small-area i (1, . . ., 755) and time period t (1, . . ., 12) is represented by O it and sampled from a Poisson distribution (Equation ( 1)).The Poisson mean, µ it , is the product of expected call-for-service count (E it ) and relative risk (r it ) (Equation ( 2)).Poisson models are common in spatial and spatio-temporal analysis of count data at the small-area scale [50].
Equation ( 3) decomposes relative risk into overall risk for the study region (α) and spatial, temporal, and space-time random effects.This is the general nonparametric model structure used for space-time analysis [44].The parametric structure models space-time interaction by incorporating a time covariate that assumes linearity [44,51].Main spatial effects (u i + s i ) measure the spatial pattern of calls-for-service and account for spatial autocorrelation and overdispersion [52].Main temporal effects (γ t + φ t ) capture overall time trend for the study region.Space-time interaction (Ψ it ) measures departures from main spatial and main temporal effects.

Prior Distributions
A vague prior of a normal distribution with mean 0 and variance 1000 was specified for α.Vague priors provide little information and let the observed data dominate posterior estimates.The prior for unstructured spatial effects (u i ) was a normal distribution with mean zero and variance σ 2 u (Equation (4a)).Spatially structured random effects (s i ) detect spatial autocorrelation and were assigned an intrinsic conditional autoregressive (ICAR) prior distribution, with spatial adjacency matrix W and variance σ 2 s (Equation (4b)).For reference, W is a symmetrical matrix where adjacency for a given small area includes all small areas sharing a vertex.ICAR is a common prior distribution for spatial random effects parameters, where the expected mean of s i is equal to the mean of adjacent s i 's [29,52,53].The variance of s i is controlled by σ 2 s and is inversely proportional to the number of neighbors of area i.Note that the ICAR prior imposes spatial structure.
Unstructured temporal effects (γ t ) were assigned a prior of a normal distribution with mean zero and variance σ 2 γ (Equation (5a)).Similar to s i , the prior distribution for structured temporal effects (φ t ) was an ICAR process with temporal weight matrix P and variance σ 2 φ , where φ t borrows strength from adjacent time periods (Equation (5b)).This is analogous to the spatial ICAR model (Equation (4a,b)).Alternative priors where temporal adjacency is conditional only on previous time periods are available [54,55], however ICAR was considered suitable in this context because calls-for-service at time t should be correlated with incidents at both t − 1 and t + 1 [18].
Without prior information regarding residual space-time structure, we tested four prior distributions for Ψ it .Each imposes different assumptions regarding spatial and/or temporal structure [44].Testing many types of space-time interaction would generally be prohibitively time consuming using MCMC-based algorithms, however INLA enables efficient analysis and model checking of multiple prior distributions.
Type I space-time interaction assumes that Ψ it for each space-time unit is exchangeable and independently and identically distributed (i.e., u i and γ t interact, as shown in Equation ( 3)).This type of space-time interaction is suitable in the case that residual space-time patterns in calls-for-service for a given space-time unit of analysis are not correlated with adjacent areas or times.Type II interaction is composed of spatially unstructured but temporally structured effects, assuming that small-area incident risk is similar over time but independent in space such that adjacent areas do not have similar patterns of calls-for-service (i.e., u i and φ t interact).Type III interaction is composed of spatially structured but temporally unstructured effects and assumes that small-area incident risk exhibits spatial autocorrelation for each time period but are independent in time (i.e., s i and γ t interact).Type IV space-time interaction captures both spatial and temporal structure and assumes that calls-for-service for a space-time unit of analysis are both spatially and temporally correlated (i.e., s i and φ t interact) [40,44].This specification is reasonable when the trend of calls-for-service for one area is similar over time and geographically adjacent areas exhibit similar trends.The variance for all types of Ψ it is controlled by σ 2 Ψ .Further details regarding the assumptions and modeling of inseparable space-time interaction can be found in Knorr-Held [44].

Hyperprior Distributions
The third level of this Bayesian hierarchical model specifies prior distributions for variance parameters of second-level prior distributions.These are referred to as hyperprior distributions.All precision parameters (the reciprocal of variance) were assigned a vague Gamma(0.5, 0.0005) distribution [56].All models were tested with hyperprior distributions of Gamma(0.001,0.001) to ensure that results were not sensitive to hyperprior specification.

Model Implementation and Goodness-of-Fit
Modeling was completed in R using the R-INLA package [43] and results were exported and mapped in Quantum GIS 2.14.The Deviance Information Criterion (DIC) was used to assess model fit [57].Broadly, DIC is a Bayesian analogue of the Akaike Information Criterion and balances model goodness-of-fit and complexity [58].Smaller DIC values by at least five units indicate superior model fit [57].DIC is shown in Equation ( 6), where D is the mean deviance and p D is effective number of parameters, or the mean deviance (D) minus the deviance evaluated at the posterior estimated values ( D) (Equation ( 7)).Of note, highly correlated parameters (i.e., space, time or space-time random effects) are often counted as less than one "effective parameter" in DIC calculations [59].

Results
DIC values and the number of effective parameters (p D ) for the four spatio-temporal models testing space-time interaction are shown in Table 3.Note that DIC values from the models testing Type I, II, III, and IV space-time interaction do not exhibit extreme differences.This can be attributed to all models taking the form shown in Equation ( 3), but differing with respect to the prior distributions of the space-time interaction parameters.The model analyzing Type IV space-time interaction had the lowest DIC and is identified as the best fitting model.The degree to which spatial, temporal, and spatio-temporal structures influence posterior parameter estimates can be observed by comparing the effective number of parameters (p D ) and the number of data points analyzed (i.e., 755 DAs × 12 time periods = 9060).In this research, p D is small relative to the number of data points, indicating that there is considerable spatial and temporal autocorrelation of calls-for-service in the dataset [39].
Each parameter in Equation ( 3) can be individually visualized using model decomposition [60].Because this research applies a Poisson log-linear model, we take the natural logarithm of model parameters to estimate relative risk due to main temporal effects (= exp(γ t + ]φ t )), main spatial effects (= exp(u i + s i )), and space-time interaction (= exp(Ψ it )).For interpretation, main effects and space-time interaction values greater than 1 can be considered to have higher than average spatial/temporal or space-time risk.Posterior means and associated 95% credible intervals (CI) are reported.95% CIs are the interval in which the posterior mean has a 95% probability of occurring.

Main Temporal Effects
Main temporal effects are shown in Figure 4 and were lowest in the early morning between 02:00 and 03:59 (0.918, 95% CI: 0.894, 0.942) and highest in the evening between 20:00 and 21:59 (1.047, 95% CI: 1.025, 1.071).There was an increasing trend from early morning to the middle of the day, peaking at 10:00-11:59 (1.045, 95% CI: 1.024, 1.068), followed by a decrease at 12:00-13:59 and a gradual increasing trend through the evening.This temporal pattern is noticeably different from the observed data, where the highest number of calls-for-service occurred between 14:00 and 15:59.This can be attributed to the ICAR smoothing imposed on φ t as well as the simultaneous analysis of spatial variation of calls-for-service.
ISPRS Int.J. Geo-Inf.2016, 5, 162 8 of 15 time interaction values greater than 1 can be considered to have higher than average spatial/temporal or space-time risk.Posterior means and associated 95% credible intervals (CI) are reported.95% CIs are the interval in which the posterior mean has a 95% probability of occurring.

Main Temporal Effects
Main temporal effects are shown in Figure 4 and were lowest in the early morning between 02:00 and 03:59 (0.918, 95% CI: 0.894, 0.942) and highest in the evening between 20:00 and 21:59 (1.047, 95% CI: 1.025, 1.071).There was an increasing trend from early morning to the middle of the day, peaking at 10:00-11:59 (1.045, 95% CI: 1.024, 1.068), followed by a decrease at 12:00-13:59 and a gradual increasing trend through the evening.This temporal pattern is noticeably different from the observed data, where the highest number of calls-for-service occurred between 14:00 and 15:59.This can be attributed to the ICAR smoothing imposed on φt as well as the simultaneous analysis of spatial variation of calls-for-service.

Main Spatial Effects
Main spatial effects of calls-for-service are shown in Figure 5.In general, main spatial effects were highest in central areas of the study region and lowest in peripheral and rural areas, particularly in the northwest.The highest spatial risk was found in Area A (1113.06, 95% CI: 781.77, 1551.25) which is likely due to a residential population considerably smaller than all other areas, and consequently, a low expected calls-for-service count (Figure 5).This is a limitation of using residential population to calculate expected counts.Area B had the second highest spatial risk of 284.76 (95% CI: 204.93, 388.51).The characteristics of Areas B, C, and D are relevant to understanding the geographic distribution of police calls-for-service and are discussed below.

Space-Time Interaction
Space-time interaction and space-time hotspots for 10:00-11:59 and 20:00-21:59 are shown in Figures 6 and 7, respectively.These time periods were chosen because they correspond to the highest main temporal effects (Figure 4) and demonstrate the geographic variability in space-time interaction throughout the course of a day.Space-time hotspots were identified by monitoring the posterior probability (PP) that exp(Ψit) is greater than 1.PP can be interpreted as the Bayesian equivalent of a p-value [61].The closer PP is to 1, the stronger the evidence that areas are space-time hotspots.Compared with the point estimate (i.e., posterior mean) of space-time interaction (Ψit), PP is a more robust indicator for detecting hotspots since it accounts for the variance/uncertainty of Ψit and makes

Main Spatial Effects
Main spatial effects of calls-for-service are shown in Figure 5.In general, main spatial effects were highest in central areas of the study region and lowest in peripheral and rural areas, particularly in the northwest.The highest spatial risk was found in Area A (1113.06, 95% CI: 781.77, 1551.25) which is likely due to a residential population considerably smaller than all other areas, and consequently, a low expected calls-for-service count (Figure 5).This is a limitation of using residential population to calculate expected counts.Area B had the second highest spatial risk of 284.76 (95% CI: 204.93, 388.51).The characteristics of Areas B, C, and D are relevant to understanding the geographic distribution of police calls-for-service and are discussed below.
use of the full posterior probability distribution of Ψit, potentially minimizing the probability of detecting false hotspots [53].

Space-Time Interaction
Space-time interaction and space-time hotspots for 10:00-11:59 and 20:00-21:59 are shown in Figures 6 and 7, respectively.These time periods were chosen because they correspond to the highest main temporal effects (Figure 4) and demonstrate the geographic variability in space-time interaction throughout the course of a day.Space-time hotspots were identified by monitoring the posterior probability (PP) that exp(Ψ it ) is greater than 1.PP can be interpreted as the Bayesian equivalent of a p-value [61].The closer PP is to 1, the stronger the evidence that areas are space-time hotspots.Compared with the point estimate (i.e., posterior mean) of space-time interaction (Ψ it ), PP is a more robust indicator for detecting hotspots since it accounts for the variance/uncertainty of Ψ it and makes use of the full posterior probability distribution of Ψ it , potentially minimizing the probability of detecting false hotspots [53].

Discussion
Results of this analysis help to understand spatial, temporal, and spatio-temporal law enforcement resource demand and provide a starting point for more detailed inquiry focused on possible explanations for space-time patterns.Providing context for the interpretation of results is the routine activity theory, which hypothesizes that crime offenses result from the convergence of motivated offenders, suitable targets, and a lack of capable guardianship in time and space [62].While the routine activity theory was originally proposed to explain crime offenses, it has been used extensively to analyze call-for-service data because it provides a theoretical lens to interpret spatio-temporal fluctuations in police data as a function of population-level movement patterns as well as small-area land-use characteristics [63].It should be noted that the types and quantities of police resources required when responding to different classes of calls-for-service varies and that this is one limitation of analyzing aggregated call-for-service data.

Interpreting Results of Spatio-Temporal Analysis
Main temporal effects of police calls-for-service are representative of the general shifts in resource demand throughout the course of a day and can, for example, be used to interpret overall staffing requirements.As shown in Figure 4, the main temporal effect peaks at two time periods, between 10:00 and 11:59 and between 20:00 and 21:59.Despite similar main temporal effects, the compositions of calls-for-service for these time periods were distinct.For example, there were a total of 1843 reported motor vehicle collisions during 10:00-11:59 compared to 767 during 20:00-21:59.Similarly, during 20:00-21:59 there were 1447 bylaw complaints compared to 495 during 10:00-11:59.Applied to police resource allocation, dispatch staffing should peak at these two time periods, however, higher levels of police units capable of responding to motor vehicle collisions should be staffed during 10:00-11:59 and higher levels of bylaw enforcement should be staffed during 20:00-21:59.
Main spatial effects highlight the general spatial pattern of police calls-for-service.Areas with high spatial effects have consistently high numbers of calls-for-service and are areas where police resources should be directed for all time periods.Exploring the characteristics of Areas B, C, and D provides insight into reasons these areas may have high spatial effects.The study region's largest shopping centre and a major transit node are located in Area B and frequently occurring incidents in this area were motor vehicle stops (1108), minor motor vehicle collisions (187), and traffic enforcement (169).Area C contains a large university and student housing and incidents were largely composed of lost or damaged property (1257), general emergency calls (220), and small theft (372).Area D contains a major highway that transects the study region and calls-for-service in this area were dominated by emergency calls (38,173), incidents recorded as part of routine detail (750), and vehicle stops (459).
Space-time interaction estimates area-specific departures from main spatial and main temporal effects [44] and indicates if additional resources are needed for specific areas apart from baseline spatial and temporal resource allocations.For 10:00-11:59, there were no areas with exp(Ψ it ) > 1.5 and areas with moderate space-time interaction were dispersed throughout the study region (Figure 6A).In comparison, during 20:00-21:59, there were 63 areas with exp(Ψ it ) > 1.5 and these areas were concentrated in and around central business areas (Figure 7A).This may be explained by the shifting patterns of routine activities of residents [62]: daytime employment and leisure activities are geographically dispersed in urban, suburban, and rural areas of the study region, while central business areas are where evening leisure activities are concentrated [64,65].
Space-time hotspots, identified by evaluating the posterior probability of exp(Ψ it ) > 1.0, are specific areas that require additional law enforcement resources to handle high numbers of calls-for-service.Posterior probability ranges from 0 to 1 and enables hotspot areas to be ranked based on the strength that exp(Ψ it ) > 1.0.The three areas with greatest PP for 10:00-11:59 and 20:00-21:59 are shown in Figures 6B and 7B (labelled as Areas 1, 2, and 3), respectively.During 10:00-11:59, space-time hotspots were found to be concentrated in the central business area of Kitchener, but located close to a local university during 20:00-21:59.
To further explore geographic variations in space-time interaction, we visualize space-time interaction and space-time hotspots for 02:00-03:59 (Figure 8A,B, respectively).This time period had the lowest main temporal effect (Figure 4).During this time period, space-time interaction was clustered around central business areas and extended along the major commercial and transportation corridor in the study region.While overall staffing may be lowest during this time period, calls-for-service were geographically clustered so that resources can be targeted to specific high-demand areas.To further explore geographic variations in space-time interaction, we visualize space-time interaction and space-time hotspots for 02:00-03:59 (Figure 8A,B, respectively).This time period had the lowest main temporal effect (Figure 4).During this time period, space-time interaction was clustered around central business areas and extended along the major commercial and transportation corridor in the study region.While overall staffing may be lowest during this time period, calls-forservice were geographically clustered so that resources can be targeted to specific high-demand areas.

Modeling Bigger Spatio-Temporal Datasets with INLA
Bayesian analysis of spatio-temporal datasets via MCMC algorithms is computationally expensive because they gradually alter posterior distributions and require researchers to monitor posterior estimates and check convergence diagnostics [42].INLA, in contrast, calculates posterior distributions through numerical integration and does not require consistent monitoring of convergence.We compared the time required to estimate Equation ( 3  While this research analyzes 9,060 space-time units, the advantages of INLA modeling may be more fully realized on bigger spatio-temporal datasets.When data are available for large spatial and temporal extents, INLA spatio-temporal models may be applied to analyze small-area data for a country (e.g., Canada was composed of 56,204 DAs in 2011) over many years to identify generalizable, regional, and neighbourhood-specific patterns in calls-for-service.When spatiotemporal point data are available, INLA can be extended via stochastic differential partial equation (SPDE) modeling [40].

Conclusions
The implementation of computer-aided-dispatch systems in law enforcement has led to the increasing availability of call-for-service datasets with both spatial and temporal information.These datasets have facilitated the growth of crime mapping, however past research has infrequently considered spatial and temporal patterns simultaneously.Space-time analysis requires methods that can identify spatial and temporal patterns, as well as space-time interaction, and account for both local spatial autocorrelation and temporal autocorrelation.
This research applied a Bayesian spatio-temporal modeling approach implemented via INLA to identify spatial, temporal, and space-time patterns of calls-for-service for two-hour time periods at the small-area level.Results were interpreted through the routine activity theory, which focuses on

Modeling Bigger Spatio-Temporal Datasets with INLA
Bayesian analysis of spatio-temporal datasets via MCMC algorithms is computationally expensive because they gradually alter posterior distributions and require researchers to monitor posterior estimates and check convergence diagnostics [42].INLA, in contrast, calculates posterior distributions through numerical integration and does not require consistent monitoring of convergence.We compared the time required to estimate Equation (3) with the simplest Type I space-time interaction on an IBM ThinkPad with a 2.4 GHz processor and 12 gigabytes RAM.Using MCMC in WinBUGS, model convergence required approximately 13 hours whereas INLA required approximately seven minutes.
While this research analyzes 9060 space-time units, the advantages of INLA modeling may be more fully realized on bigger spatio-temporal datasets.When data are available for large spatial and temporal extents, INLA spatio-temporal models may be applied to analyze small-area data for a country (e.g., Canada was composed of 56,204 DAs in 2011) over many years to identify generalizable, regional, and neighbourhood-specific patterns in calls-for-service.When spatio-temporal point data are available, INLA can be extended via stochastic differential partial equation (SPDE) modeling [40].

Conclusions
The implementation of computer-aided-dispatch systems in law enforcement has led to the increasing availability of call-for-service datasets with both spatial and temporal information.These datasets have facilitated the growth of crime mapping, however past research has infrequently considered spatial and temporal patterns simultaneously.Space-time analysis requires methods that can identify spatial and temporal patterns, as well as space-time interaction, and account for both local spatial autocorrelation and temporal autocorrelation.
This research applied a Bayesian spatio-temporal modeling approach implemented via INLA to identify spatial, temporal, and space-time patterns of calls-for-service for two-hour time periods at the small-area level.Results were interpreted through the routine activity theory, which focuses on population-level activity patterns in both space and time.It was shown that main temporal effects of calls-for-service peaked between 20:00 and 21:59 and were lowest between 02:00 and 03:59.Purely spatial call-for-service patterns showed high risk in central business areas and in areas with highways, a local university, and a regional shopping centre.Space-time interaction estimates departures from main spatial and temporal patterns and enables the identification of space-time hotspots, or areas with high posterior probability of exhibiting space-time interaction greater than 1.During daytime hours (10:00-11:59), space-time hotspots were dispersed throughout the study region but during evening (20:00-21:59) and early morning time periods (02:00-03:59) space-time hotspots were clustered in downtown areas and close to a local university.
One limitation of this research is that expected call-for-service counts were calculated using residential population.Recent research has shown that alternative measures of population (i.e., the denominator in crime rate calculations) may influence the results and interpretation of spatial analyses of crime [66].Future research may consider including alternative measures to population, such as ambient populations, workday or commuting populations, or populations estimated via social media data [67].A second limitation of this research is that we only consider the routine activity theory and that there are a variety of additional theoretical perspectives and covariates not included.Future studies analyzing a specific call-for-service type should include covariates such as neighbourhood socioeconomic status, local land-use composition, and police perception, particularly if analyzing specific types of calls-for-service [6,68].Analysis of specific call types may also provide more details regarding the types and quantities of law enforcement resources required.
Future research analyzing spatio-temporal datasets should consider analyzing service time (i.e., hours of law enforcement resources) to better understand overall resource allocation and match call-for-service data to service time to quantify the resources required by different call types.While this research described the spatial, temporal, and residual spatio-temporal patterns of crime, small-area risk factors were not included.Future research should explore associations between neighbourhood characteristics and calls-for-service by extending this spatio-temporal model to include covariates and consider analyzing point data to highlight specific addresses or intersections where calls-for-service exhibit high clustering [62].It would also be advantageous to analyze hourly patterns in context of daily, monthly, and seasonal patterns, which may be suited for a multilevel spatio-temporal modeling approach.Finally, future research may adapt the presented spatio-temporal Bayesian model for probabilistic forecasts [69,70].

15 Figure 1 .
Figure 1.Geographic distribution of (A) count of all police calls-for-service in 2011 and (B) expected count of calls-for-service.

Figure 1 .
Figure 1.Geographic distribution of (A) count of all police calls-for-service in 2011 and (B) expected count of calls-for-service.

Figure 2 .
Figure 2. Total calls-for-service by 2 hour time period.

Figure 5 .
Figure 5. Main spatial effects (= exp(u i + s i )) of police calls-for-service.
ISPRS Int.J. Geo-Inf.2016, 5, 162 11 of 15 ) with the simplest Type I spacetime interaction on an IBM ThinkPad with a 2.4 GHz processor and 12 gigabytes RAM.Using MCMC in WinBUGS, model convergence required approximately 13 hours whereas INLA required approximately seven minutes.

Table 1 .
Descriptive statistics for call-for-service counts in Waterloo Region, 2011.

Table 2 .
Descriptive statistics for the temporal variation of calls-for-service.Note that 00:00 refers to midnight and 12:00 refers to noon.

Table 2 .
Descriptive statistics for the temporal variation of calls-for-service.Note that 00:00 refers to midnight and 12:00 refers to noon.

Table 3 .
Deviance Information Criterion (DIC) values for models with Type I, II, III, and IV space-time interaction.