Speed Limit Changes and Driver Behaviour: A Spatial Lag Model

Findings Speed is a major factor in road safety and the interplay of the different factors affecting speed choice is not completely understood. This paper presents a study of operating speeds in Quebec, Canada, between 2000 and 2018 on provincial highways (except freeways) in response to speed limit changes. This paper shows that driving speed is spatially correlated. Two statistical models are then compared, with and without a spatial component. The spatial model provides a better fit and demonstrates that speed behaviour is spatially correlated. The resulting models confirm that drivers only partly adjust their speeds after a change and that several road features like side curbs, lateral and median strips are associated with driving speed.


Research Question and Hypothesis
The lack of compliance to speed limits is pervasive. Typically, 40 % to 50 % of motorists drive above the speed limit (OECD/ECMT 2006). Speed is an essential factor in road safety because it directly affects the frequency and severity of crashes (Aarts and Van Schagen 2006;FHWA 2017). Driver behaviour is influenced by many external factors such as road design, road environment, traffic control devices, the presence of other vehicles, and weather conditions (Gargoum, El-Basyouny, and Kim 2016;Goldenbeld and van Schagen 2007;Wilmot and Khanal 1999). Transport authorities try to change the drivers' behaviour through education, speeding controls, changes to the road design or road environment, and through technology such as driver assistance systems in vehicles (Wallén Warner and Åberg 2008). However, the effects of these measures are often limited in time and space (Comte, Varhelyi, and Santos 1997). This paper examines speed choice by drivers in response to a change in the speed limit without any known modification to the road design and the immediate environment. This paper tests whether driving speed is spatially correlated and compares two statistical models, a regular linear model, and a spatial lag model.

Methods and data
In total, 76 road segments (sites) with a decrease (n=12) or an increase (n=64) of their speed limit between 2006 and 2017 were analysed, over a period of six years before and after the speed limit changes. Those sites are located across the province (see Figure 1). The targeted road segments are provincial highways (national, regional, and collector roads), excluding local roads and freeways. First, the individual vehicle speed data were collected at the study sites by the Quebec Ministry of Transportation (MTQ) and the research team. For each site, speeds were recorded at different periods, usually for three to five consecutive days at a time, before and after the change (the records done by the research team were only three hours, following the MTQ guidelines). There are between 2 and 17 records at each site. Only free-flow speeds are used in this analysis. Figure 2 illustrates the differential of the 85 th percentile of the average speeds (after minus before) per site according to the change in speed limits. This figure shows the non-linear pattern of the changes in speed behaviour by the drivers according to the speed limit change: drivers do not all comply with the speed limit, nor do they adjust their speed according to the magnitude of the speed limit change.
Secondly, for each road segment and period during which speed was recorded, data about the road design and environment were extracted from public spatial databases. These explanatory (independent) variables are described in Table 1.
To select an appropriate model using these variables, we applied the stepwise model selection procedure, which automatically selects the best subset of variables from a pool of candidates using the AIC (Akaike Information Criterion).
The speed ratio, defined as the average speed for each speed record divided by the speed limit, is the dependent variable. It is used to compare the average operating speeds between sites that do not have the same speed limit (Morency Since the speed behaviour of drivers is influenced by other drivers and the road environment, the social interactions and culture, there may be some spatial correlations between the observed speeds. Unlike linear models, a spatial model considers the neighborhood effect of observations. The model (with the spatial lag) can be represented as where is the speed ratio at site , is the vector of variables for the site i, the vector of coefficients, the spatial weight for sites i and j, a spatial autoregression scalar parameter, the errors are supposed to be normally distributed with standard deviation and n is the number of sites. depicts a lack of spatial correlation, while implies a positive or negative spatial autocorrelation between neighboring sites. For a proper treatment of spatial linear model and their variations, see Diggle (2003) and Bivand, Pebesma, and Gómez-Rubio (2013), among others.
To examine the spatial autocorrelation in the data, i.e., test whether a variable at one site depends on neighboring observations, the spatial autocorrelation for a variable can be measured through Moran's I:  The hypothesis that Moran's I is 0 can be tested (Bivand (2019)): a positive (resp. negative) spatial autocorrelation indicates that similar (resp. dissimilar) values appear close (resp. far) to each other in space. There are different methods to compute the spatial weights (Bivand, Pebesma, and Gómez-Rubio (2013)): here, two sites are neighbour within a reasonable maximum distance of influence (a threshold of 50 km was chosen based on the distribution of site distance by some trial and error). Three sites have no neighbour, and other sites have 4.9 neighbours on average. If there is spatial correlation for the speed ratio and for the spatial correlation of the residuals of a regular linear model, a spatial model can be fitted to the data.

Findings
Moran's I for the speed ratio is -0.053 ( ). A spatial and regular model were therefore fitted to the selected explanatory variables, only for the 49 sites with reduced speed limits. Moran's I for the residuals of the regular model Speed Limit Changes and Driver Behaviour: A Spatial Lag Model Findings is -8.860 ( ). The estimated parameters of the proposed models are presented in Table 2. The spatial autoregression parameter, ( ), which implies a negative spatial autocorrelation between neighboring sites. Furthermore, comparing the spatial and regular models in terms of AIC shows the spatial regression provides a better fit to the data. The models are very similar, with the same coefficient signs and the same significant variables (except for the year 2011). Most of the variable associations with the speed ratio are expected. There is an increase in the speed ratio after a speed limit decrease, which means that driving speed increases relative to the speed limit and is consistent with Figure 2. This increase is larger in urban areas based on the area and distance to the nearest city center. Road features like side curbs, lateral and median strips are associated with lower speed ratios and may be used as traffic calming measures. There is no clear temporal trend.
Our results demonstrate that speeds are spatially correlated. It also confirms that drivers only partly adjust their speeds after a speed limit change. More research is needed to understand how to increase compliance and how social and cultural factors may play a role in driving speed.