|
|
@@ -0,0 +1,649 @@
|
|
|
+---
|
|
|
+category: literaturenote
|
|
|
+citekey: jigeographicallyweightedpoissonregression2021
|
|
|
+title: "Geographically weighted poisson regression under linear model of coregionalization assistance: Application to a bicycle crash study"
|
|
|
+authors: "Ji, Shujuan; Wang, Yuanqing; Wang, Yao"
|
|
|
+year: 2021
|
|
|
+date: 2021-00-00 2021
|
|
|
+doi: 10.1016/j.aap.2021.106230
|
|
|
+publication: "Accident Analysis & Prevention"
|
|
|
+url: "https://doi.org/10.1016/j.aap.2021.106230"
|
|
|
+zotero_key: RVGQ646R
|
|
|
+zotero_storage: ZI2S6HDF
|
|
|
+collections: transpordi_planeerimine
|
|
|
+folder: 001_artiklid/Road Safety Indicators
|
|
|
+firstAuthor: "Ji, Shujuan"
|
|
|
+status: converted
|
|
|
+---
|
|
|
+
|
|
|
+
|
|
|
+Contents lists available at [ScienceDirect](www.sciencedirect.com/science/journal/00014575)
|
|
|
+
|
|
|
+# Accident Analysis and Prevention
|
|
|
+
|
|
|
+journal homepage: [www.elsevier.com/locate/aap](https://www.elsevier.com/locate/aap)
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+# Geographically weighted poisson regression under linear model of coregionalization assistance: Application to a bicycle crash study
|
|
|
+
|
|
|
+Shujuan Ji , Yuanqing Wang \* , Yao Wang
|
|
|
+
|
|
|
+*Key Laboratory of Transport Industry of Management, Control and Cycle Repair Technology for Traffic Network Facilities in Ecological Security Barrier Area, Chang'an University, P.O. Box 487, Xi'an 710064, China*
|
|
|
+
|
|
|
+# ARTICLE INFO
|
|
|
+
|
|
|
+*Keywords:* Semi-parametric GWPR Spatial non-stationary Spatial correlation Bicycle lane
|
|
|
+
|
|
|
+#### ABSTRACT
|
|
|
+
|
|
|
+While cycling benefits individuals and society, cyclists are vulnerable road users, and their safety concerns arouse more macro-level spatial crash studies. Our study intends to investigate the spatial effects of population, land use, and bicycle lane infrastructures on bicycle crashes. This was done by developing a semi-parametric Geographically Weighted Poisson Regression (sGWPR) model which deals with the issue of spatial correlation and spatial non-stationarity simultaneously. It is a model that combines both constant and geographically varying parameters. To determine which parameter is fixed or non-stationary, previous studies suggest monitoring the Akaike Information Criterion (AICc). Yet, relying only on AICc might bury some spatial associations. So, in this study, we propose a Linear Model of Coregionalization (LMC) to assist the decision. Here, we use bicycle crash data across the metropolitan area of Greater Melbourne to establish sGWPR models suggested by AICc and LMC, respectively. Comparing the two sGWPR models, we found the sGWPR model under LMC results performs as well as sGWPR models suggested by AICc from the AICc perspective, and a 22.5% improvement in the mean squared error (MSE). It also uncovers more details about the spatial relationship between bicycle crashes and bicycle lane intersection density (BLID), an effect not suggested under AICc results. The parameters of BLID, a new measurement of bicycle lane facilities proposed by us, vary over space across analysis zones in Greater Melbourne.
|
|
|
+
|
|
|
+### **1. Introduction**
|
|
|
+
|
|
|
+Cycling can deliver health, environmental, economic, and social benefits. Bicycle use in Australia has steadily grown. For example, in Melbourne City and Greater Melbourne, about 41% and 37% of residents participated in cycling respectively, due to the implementation of a series of Bicycle Plan in 2012 ([Melbourne, 2013](#page-14-0)). Despite the benefits of cycling, cyclists are one of vulnerable road users who are exposed to a higher risk of crashes in traffic. In the world, approximately 1.35 million people die each year as a result of road traffic crashes, and 26% of all road traffic deaths are among cyclists, [\(WHO \(World Health Organiza](#page-14-0)[tion\), 2018](#page-14-0)). The high traffic death rate for cyclists impels micro-level and macro-level spatial crash studies. This study focuses on the latter.
|
|
|
+
|
|
|
+An important issue is spatial dependence which is the propensity for nearby locations to influence each other and to possess similar attributes. Bicycle crash occurrence is more likely to be similar at points two meters apart than at points two kilometers apart. Another possible issue is spatial non-stationarity, a condition where the relationships between the independent and dependent variables vary over space
|
|
|
+
|
|
|
+([Fotheringham et al., 2002](#page-13-0)). Ignorance of spatially non-stationarity causes biased parameter estimates and misleading significance tests ([Anselin and Griffith, 1988](#page-13-0)). The Geographically Weighted Poisson Regression (GWPR) technique has been used to fully capture spatial dependence in the residuals and spatial non-stationarity simultaneously in macro-level crash count analysis ([Huang et al., 2018, Ariannezhad](#page-14-0) [et al., 2020, Li et al., 2013, Xu and Huang, 2015, Zhang et al., 2015,](#page-14-0) [Gomes et al., 2019, Bao et al., 2017, Briz-Redon et al., 2019, Arvin et al.,](#page-14-0) [2019, Hezaveh et al., 2019, Yao et al., 2015](#page-14-0)). For example, [Huang et al.](#page-14-0) [\(2018\)](#page-14-0) used the GWPR model to find that both the strength and direction (either positive or negative relationship) of relationships between crashes and the built environment vary over space.
|
|
|
+
|
|
|
+Semi-parametric Geographically Weighted Poisson Regression (sGWPR) is one of the GWPR techniques. It is a mixed model that combines both geographical varying and constant parameters to accommodate situations in which some parameters may not vary geographically ([Nakaya et al., 2005\)](#page-14-0). Nakaya and his colleagues (2005) suggested the decision whether a parameter varies geographically or not should follow the stepwise regression model selection process by
|
|
|
+
|
|
|
+*E-mail address:* [wyqing@chd.edu.cn](mailto:wyqing@chd.edu.cn) (Y. Wang).
|
|
|
+
|
|
|
+<sup>\*</sup> Corresponding author.
|
|
|
+
|
|
|
+monitoring Akaike Information Criterion (AICc). The comparison is between a full GWPR with all varying parameters and a sGWPR model in which all but one is allowed to vary geographically. Then select the one with a smaller AICc value in each comparison only if the AICc difference between two models is greater than 2. However, model selection or variable selection (it is a variable selection of which the variable's parameter should be fixed or varied) exclusively depending on the AICc might ignore the important details in the spatial relationships ([Ribeiro](#page-14-0) [et al., 2016\)](#page-14-0).
|
|
|
+
|
|
|
+To overcome this, our study aims to establish an sGWPR technique under Linear Model of Coregionalization (LMC) to provide an insight in appropriate model frameworks for macro-level crash analysis. LMC is a geostatistical tool broadly used in applied geology, ecology, epidemiology and other geosciences [\(Wackernagel, 2003, Pelletier et al., 2008,](#page-14-0) [Ribeiro et al., 2013, Zawadzki and Badura, 2013\)](#page-14-0) but less used in crash studies. Since it is a convenient tool to estimate the spatial relationship between variables and to quantify the spatial dependence of variables ([Goulard and Voltz, 1992, Isaaks and Srivastava, 1989](#page-14-0)), a LMC estimation, instead of an AICc measure, is applied to the model selection. In addition, under the framework, we also tend to assess bicycle safety implications socio-economic, road network, and traffic characteristics. During the assessment, we develop a new variable, bicycle lane intersection density (BLID) which no previous studies have mentioned, to measure the effect of intersections between bicycle lanes and other types of lanes.
|
|
|
+
|
|
|
+The remainder of this article is organized as follows: A literature review for macro-level bicycle crashes is provided in the second part. The next section describes data and the study area. The GWPR, sGWPR and LMC used in our study are presented in the third part. This is followed by results and discussions. The last section is the conclusion.
|
|
|
+
|
|
|
+## **2. Literature review of cashes analyses at Macro level**
|
|
|
+
|
|
|
+# *2.1. Built environment factors contributing to bicycle crashes*
|
|
|
+
|
|
|
+One of the road characteristic factors is the attribute of bicycle facilities or road intersections, such as intersection density [\(Siddiqui et al.,](#page-14-0) [2012, Raihan et al., 2019, Osama and Sayed, 2016, Lee et al., 2015](#page-14-0)), bicycle lane length ([Wei and Lovegrove, 2013](#page-14-0)), bicycle lane types ([Morrison et al., 2019](#page-14-0)) and bike network connectivity [\(Zhang et al.,](#page-15-0) [2015, Osama and Sayed, 2016](#page-15-0)). Some studies found road intersection density and bicycle lane length are all positive to bicycle crashes. [Mor](#page-14-0)[rison et al. \(2019\)](#page-14-0) point out that bicycle lane effect on bicycle risks depends on the bicycle lane types and the roadway conditions or characteristics. For example, the crash odds on roadway sections with bicycle lanes increase (compared to the roadway sections without bicycle lanes) while bicycle lanes were associated with reduced crash odds where bus routes and tram stops were present.
|
|
|
+
|
|
|
+The pattern of land use is an important predictor of crashes. Land use is classified as commercial, residential, industrial etc. Different patterns of commercial use have different effects on bicycle crashes. In studies of [Narayanamoorthy et al. \(2013\) and Osama and Sayed \(2016\),](#page-14-0) the commercial use is insignificant in bicycle crashes. Retail land use is not always significantly associated with an increase in cyclist casualties ([Priyantha Wedagama et al., 2006, Cai et al., 2016](#page-14-0)). However, the commercial use [\(Zhang et al., 2015\)](#page-15-0), big box stores (a retail land use with a building area of 50,000 square feet or more) ([Dumbaugh and Li,](#page-13-0) [2010\)](#page-13-0) and "log of rooms of hotel, motel, and time shares per square mile" are positive in non-motorist crashes ([Lee et al., 2015](#page-14-0)).
|
|
|
+
|
|
|
+Residential use has a positive effect on bicycle crashes in most studies. The number of dwelling units ([Siddiqui et al., 2012, Osama and](#page-14-0) [Sayed, 2016, Guo et al., 2018](#page-14-0)) and percentage of residential use [\(Amoh-](#page-13-0)[Gyimah et al., 2016](#page-13-0)) are associated with an increase of cyclist crashes. Only in the study of [Zhang et al. \(2015\)](#page-15-0) was number of dwelling units an insignificant factor. Some studies explored the impact of mixed land use development on bicycle crashes. Different calculations of a land use mix index have different effects: the land use balance index [\(Amoh-Gyimah](#page-13-0) [et al., 2017\)](#page-13-0) was associated with an increase of bicycle crashes; land use mixture was negatively associated with the increase ([Chen and Shen,](#page-13-0) [2016\)](#page-13-0).
|
|
|
+
|
|
|
+# *2.2. Other factors contributing to bicycle crashes*
|
|
|
+
|
|
|
+Commonly, bicycle crashes are predicted by socio-economic and demographic factors such as percentage of households without motor vehicles, median household income, and age. [Amoh-Gyimah et al.](#page-13-0) [\(2016\)](#page-13-0) found that a young population (under 19 years old) is less likely to be involved in bicycle crashes while an elderly population (above 65 years old) has a positive relationship with total and minor-level crashes. [Saha et al. \(2018\)](#page-14-0) discovered that a population aged from 18 to 64 all has a positive association with cyclist crashes. The variable of a household's car ownership is positive with more cycling crahes ([Saha et al.,](#page-14-0) [2018, Siddiqui et al., 2012](#page-14-0)).
|
|
|
+
|
|
|
+Previous studies have explained the risk or exposure characteristics on macro-level bicycle crash studies. Traffic exposure measures of vehicle traffic volume ([Dumbaugh and Li, 2010, Guo et al., 2018, Zhang](#page-13-0) [et al., 2015](#page-13-0)), and bike traffic volume ([Saha et al., 2018, Guo et al., 2018](#page-14-0)), or measures of heavy vehicles ([Cai et al. 2016\)](#page-13-0) are significant factors in bicycle crashes. Other significant exposure variables include measurements of commuters with different transportation modes (cycling, transit riding or driving) ([Zhang et al., 2015, Lee et al., 2017, Jacobsen,](#page-15-0) [2003\)](#page-15-0) and population ([Amoh-Gyimah et al., 2016, Siddiqui et al., 2012](#page-13-0)).
|
|
|
+
|
|
|
+### *2.3. Crash prediction models*
|
|
|
+
|
|
|
+The appropriate analytical modelling approaches and effective explanatory variables are essential to the accuracy of crash prediction. Numerous crash spatial analysis methods have been explored in the previous literature and [Ziakopoulos and Yannis \(2020\)](#page-15-0) discussed those different modelling approaches. Some advanced approaches such as zero inflation [\(Cai et al., 2016, Raihan et al., 2019, Qin et al., 2004](#page-13-0)), autoregressive prior models ([Besag et al., 1991, Zhai et al., 2018b,](#page-13-0) [Huang and Darwiche, 2010](#page-13-0)), latent class [\(Yasmin and Eluru, 2016,](#page-15-0) [Heydari et al., 2017, Eluru et al., 2012](#page-15-0)), Bayesian ([Miaou. 2003, Li et al.,](#page-14-0) [2019, MacNab, 2004](#page-14-0)) and machine learning([Cai et al., 2019, Bao et al.,](#page-13-0) [2019\)](#page-13-0) have been widely applied to deal with the data problems of excess zeros, spatial autocorrelations, spatial heterogeneity, spatial spillover, and spatiotemporal effects.
|
|
|
+
|
|
|
+Apart from the GWR techniques mentioned in the Introduction, random parameter techniques (e.g. random parameter negative binomial (RPNB)) ([Xu and Huang, 2015, Amoh-Gyimah et al., 2016, Zhai](#page-14-0) [et al., 2018b, Wang et al., 2017](#page-14-0)) and Bayesian spatially varying coefficient (BSVC) models also can capture the variation relationships over space ([Congdon, 1997, Congdon, 2003, Xu et al., 2017\)](#page-13-0). The two approaches are different from GWR techniques: the former is statistical models drawn only from some mathematical distribution and structures, while the latter assumes local spatial regression models in each spatial unit. From this view, GWR techniques could provide powerful and intuitive descriptions and explanations for researchers and policy makers ([Pirdavani et al., 2014\)](#page-14-0). On the other hand, [Xu and Huang](#page-14-0) [\(2015\)](#page-14-0) employ RPNB and sGWPR approaches to account for the spatial variation's relationships between the number of crashes and traffic, road, socio-demographic characteristics, and found sGWPR is the more appropriate tool for zonal crash analysis than RPNB in terms of mean absolute and AIC values.
|
|
|
+
|
|
|
+Another issue in crash data may be overdispersion, so the Poissonbased models (e.g., negative binomial (NB)) have been used to deal with this issue. Some studies proved the GWPR has a better model performance than the NB model. Recently, Geographically Weighted Negative Binomial Regression Models (GWNBR) was proposed by [da](#page-13-0) [Silva and Rodrigues \(2014\)](#page-13-0) for capturing overdispersion. The GWNBR model incorporates the negative binomial distribution instead of the
|
|
|
+
|
|
|
+<span id="page-2-0"></span>Poisson distribution in the GWR framework. [Gomes et al. \(2017\), Tang](#page-14-0) [et al.\(2020\) and Obelheiro et al.\(2020\)](#page-14-0) applied the GWNBR using % gwnbr macro in SAS [\(da Silva and Rodrigues, 2016\)](#page-13-0) to estimate the spatial effect of traffic exposure, socio-economics, road networks, and land use on crash occurrence.
|
|
|
+
|
|
|
+[Gomes et al. \(2017\)](#page-14-0) compared the performance of the GWPR and GWNBR models. In terms of AICc value, the GWNBR has a better model performance than the GWPR and the NB, while based on the root mean squared error, the GWPR has a more powerful prediction. Further, as mentioned in the Introduction, the sGWPR can fit the situation where some parameters may not vary geographically, but the GWNBR cannot accommodate this kind of situation. Finally, SAS can estimate GWNBR rather than sGWPR, and GWR 4.0 cannot support the estimation of GWNBR but can be used for the estimation of sGWPR instead. Given all these considerations, we focus on the sGWPR results in our study and the GWNBR is ready for future research.
|
|
|
+
|
|
|
+Overall, the literature indicates the built environment is an important factor of bicycle crash prediction. Road networks (especially intersections), bicycle lane infrastructures and land use are the most frequent measurements of a built environment. However, there is a lack of consistent findings between bicycle crashes and land use, especially commercial use. [Siddiqui et al. \(2012\)](#page-14-0) also stated that research highlighting the effects of a built environment on bicycle crashes is limited. Therefore, more research is still needed to capture the interactions
|
|
|
+
|
|
|
+between built environment and bicycle crashes. We presented a new measurement—bicycle lane intersection density, for bicycle lane infrastructure characteristics. This study mainly focuses on the GWPR techniques to explore the spatial effect of commercial use, exposure characteristics and bicycle infrastructures on bicycle crashes.
|
|
|
+
|
|
|
+## **3. Data Description**
|
|
|
+
|
|
|
+# *3.1. Traffic analysis zones*
|
|
|
+
|
|
|
+Based on the 2016 Australian Census by the Australian Bureau of Statistics (ABS, 2017), Great Melbourne consists of 9 Statistical Area level 4 (SA4) zones, 309 Statistical Area level 2 (SA2) zones and 10,849 Statistical Area level 1 (SA1) zones. Comparing with using SA1 as units, SA2 zones are used as spatial analysis units in this study because more bicycle crash data can be aggregated in each unit, and several units without bicycle crashes can be avoided.
|
|
|
+
|
|
|
+# *3.2. Bicycle crash data*
|
|
|
+
|
|
|
+Bicycle crash data from three years (from 2015 to 2017) in Great Melbourne, Australian, were collected[\(VicRoads, 2019a](#page-14-0)). Our study integrates the number of bicycle crashes over three years into SA2s as a dependent variable to investigate the association with contributing
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+**Fig. 1.** Study area and spatial distribution of total bicycle crashes (2015–2017).
|
|
|
+
|
|
|
+<span id="page-3-0"></span>factors [\(Amoh-Gyimah et al., 2016, Gomes et al., 2019\)](#page-13-0). We found 51.6% of bicycle crashes occurred in an inner area. The spatial distribution of SA2 zones, an inner area (one of SA4s), and bicycle crashes are shown in [Fig. 1](#page-2-0).
|
|
|
+
|
|
|
+### *3.3. Bicycle lane intersection density (BLID) at macro level*
|
|
|
+
|
|
|
+We have created a new measurement, bicycle lane intersection density (BLID), combining the bicycle lane and road intersection characteristics simultaneously to counter the effect of bicycle lanes on bicycle crashes. Two reasons are provided to support the variable development. First, bicycle lanes are attractive to cyclists and benefit their safety [\(Schepers et al., 2017\)](#page-14-0). Roadway segments with bicycle lanes could decrease crash risks, compared to those without [\(Lusk et al.,](#page-14-0) [2011\)](#page-14-0). But the effectiveness of bicycle lanes differs according to other roadway characteristics or conditions. [Morrison et al. \(2019\)](#page-14-0) underline bicycle lanes with the roadway conditions or characteristics (i.e., speed limit, bus routes, tram routes, bridges, one-way flow, traffic lane width), which are associated with bicycle crash risks. They, however, did not consider the intersection characteristic, which is one of the important roadway characteristics in the impact of bicycle crashes. Second, the length of bicycle lanes [\(Wei and Lovegrove, 2013\)](#page-14-0) used in bicycle crash analysis did not consider roadway characteristics and could not measure the risks of bicycle lanes when the bicycle lane encounters intersections.
|
|
|
+
|
|
|
+A bicycle lane link is a road link with an on-road bicycle lane or offroad bicycle lane. Bicycle lane intersection (BLI) is the endpoint of a bicycle lane link that connects to other road links. BLID in SA2 *i* is calculated by:
|
|
|
+
|
|
|
+$$BLID_i = \frac{BLI_i}{tLBL_i \times area_i} \tag{1}$$
|
|
|
+
|
|
|
+Where*BLIi*is the number of bicycle lane intersections and*tLBLi*is the total length of bicycle lane links within SA2 *i* and*areai*area of SA2 *i*.
|
|
|
+
|
|
|
+The bicycle lane intersections were identified by three steps based on two shape files: principal bicycle network ([VicRoads, 2016](#page-14-0)) and road network (Open Street Map [\(OSM \(Open Street Map\), 2019\)](#page-14-0)) in Software ArcGIS 9.2:
|
|
|
+
|
|
|
+Step 1: The buffers of bicycle lane links with a 20-meter buffer distance were established using the Buffer tool. The bicycle lane links are
|
|
|
+
|
|
|
+usually drawn on the road central line (shown in Fig. 2) and a 20-meter radius is half of the maximum width of roads in Great Melbourne ([VicRoads, 2019b](#page-14-0)). Therefore, the buffer with a 20 m distance could cover the intersections connecting bicycle lanes and other roads;
|
|
|
+
|
|
|
+Step 2: The point file of road intersections was generated by applying the Feature Vertices To Points tool of GIS to the road network data;
|
|
|
+
|
|
|
+Step 3: Using the Intersect tool, a new point file of bicycle lane intersections was generated by intersecting the point file of road intersections created by Step 2 and the buffers created by Step 1. The bicycle lane intersections might repeat, since buffers were overlapped in the end of bicycle lane links. So, the points data of road intersections was coded with unique numbers. When summarizing the identified bicycle lane intersections within each spatial unit, the repeated bicycle lane intersections were excluded in terms of a unique number. Most SA2s have bicycle lanes. For the 32 SA2s without any bicycle lanes, the values of BLID are equal to zeros.
|
|
|
+
|
|
|
+### *3.4. Land use*
|
|
|
+
|
|
|
+Residential and commercial areas are important trip productions and attractions; the land use is considered as crash generation variable. In this study commercial use is grouped into large-scale commercial use and small-scale commercial use in this study. Land use data was collected from Victorian Resources Online ([VRO \(Victorian Resources](#page-14-0) [Online\), 2017\)](#page-14-0). A typical size of shops along the roadside or in a community is less than a threshold, 80,000 square feet (about 7432 square meters). About 96% of the commercial areas are less than the threshold and only account for about 33.4% of all commercial use. Thus, the area of the commercial use less than 80,000 square feet is defined as smallscale commercial use in our study.
|
|
|
+
|
|
|
+# *3.5. Other variables*
|
|
|
+
|
|
|
+Socio-economic and risk or exposure variables [\(Table 1\)](#page-4-0) suggested in most prior studies are also considered in this study. Socio-economic data come from the 2016 Census ([ABS \(Australian Bureau of Statistics\).,](#page-13-0) [2017\)](#page-13-0). Daily vehicle kilometers traveled (DVKT) is calculated by the Annual Average Daily Traffic (AADT) data and the roadway length within each SA2 ([VicRoads, 2019c\)](#page-14-0).
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+**Fig. 2.** Identification of bicycle lane intersections.
|
|
|
+
|
|
|
+<span id="page-4-0"></span>**Table 1** Descriptive statistics of spatial units (309 SA2s).
|
|
|
+
|
|
|
+| Variables | Description | Min | Max | Mean | S.D. |
|
|
|
+|---------------------|-----------------------------------------------------|-----|--------|-------|-------|
|
|
|
+| Bicycle | Total number of bicycle | 0 | 236.0 | 11.9 | 22.1 |
|
|
|
+| crashes | crashes | | | | |
|
|
|
+| BCR | Bicycle crash rate per | 0 | 1.05 | 0.06 | 0.12 |
|
|
|
+| | thousand vehicle<br>kilometers travelled | | | | |
|
|
|
+| Socio-economic | | | | | |
|
|
|
+| Per14 | Percentage of | 0.0 | 39.1 | 18.5 | 5.0 |
|
|
|
+| | population below 14 | | | | |
|
|
|
+| | years old % | | | | |
|
|
|
+| Per15_24 | Percentage of | 0.0 | 45.9 | 13.1 | 4.4 |
|
|
|
+| | population between 15<br>and 24 years old % | | | | |
|
|
|
+| Per25_34 | Percentage of | 0.0 | 40.7 | 15.7 | 6.6 |
|
|
|
+| | population between 25 | | | | |
|
|
|
+| | and 34 years old % | | | | |
|
|
|
+| Per35_44 | Percentage of | 0.0 | 33.3 | 14.3 | 3.0 |
|
|
|
+| | population between 35 | | | | |
|
|
|
+| | and 44 years old % | | | | |
|
|
|
+| Per45_64 | Percentage of<br>population between 45 | 0.0 | 35.3 | 23.9 | 5.3 |
|
|
|
+| | and 64 years old % | | | | |
|
|
|
+| Per64 | Percentage of | 0.0 | 32.5 | 13.5 | 5.2 |
|
|
|
+| | population above 64 | | | | |
|
|
|
+| | years old % | | | | |
|
|
|
+| Mincome | Mediate personal | 0.0 | 1.6 | 0.7 | 0.2 |
|
|
|
+| | income weekly /1000\$ | | | | |
|
|
|
+| PerNC | Percentage of household | 0.0 | 75.0 | 8.0 | 8.0 |
|
|
|
+| Risk or exposure | with no cars% | | | | |
|
|
|
+| Ln(pop) | Log of population(pop) | 0 | 10.5 | 9.4 | 1.2 |
|
|
|
+| Perbyk | Percentage of | 0.0 | 33.3 | 0.8 | 2.3 |
|
|
|
+| | commuters cycling to | | | | |
|
|
|
+| | work % | | | | |
|
|
|
+| DVKT | Daily vehicle kilometers | 3.0 | 1534.3 | 298.9 | 255.8 |
|
|
|
+| | traveled in thousands | | | | |
|
|
|
+| PerPT | Percentage of<br>commuters to work by | 0.0 | 36.9 | 11.5 | 8.0 |
|
|
|
+| | public transit (Tram, | | | | |
|
|
|
+| | transit, and light rail) % | | | | |
|
|
|
+| Road infrastructure | | | | | |
|
|
|
+| Road network | Density of road network | 0.5 | 41.2 | 9.0 | 4.0 |
|
|
|
+| density | [km/km2<br>] | | | | |
|
|
|
+| Sidewalk | Density of sidewalk | 0.0 | 60.9 | 1.7 | 4.2 |
|
|
|
+| density | [km/km2<br>] | | | | |
|
|
|
+| Residential<br>road | Density of residential<br>road [km/km2<br>] | 0.1 | 31.8 | 7.1 | 4.1 |
|
|
|
+| density | | | | | |
|
|
|
+| Arterial | Density of Arterial | 0.0 | 14.1 | 0.7 | 1.2 |
|
|
|
+| Density | mileage [km/km2<br>] | | | | |
|
|
|
+| Secondary | Density of Secondary | 0.0 | 15.8 | 0.4 | 1.0 |
|
|
|
+| density | road mileage [km/km2<br>] | | | | |
|
|
|
+| Tertiary | Density of Tertiary road | 0.0 | 14.1 | 1.3 | 1.4 |
|
|
|
+| Density | mileage [km/km2<br>] | | | | |
|
|
|
+| BLID | Bicycle lane intersection<br>density [/(km*km2<br>] | 0.0 | 19.2 | 1.5 | 2.1 |
|
|
|
+| Land use | | | | | |
|
|
|
+| PerLC | Percentage of large | 0.0 | 44.8 | 2.0 | 3.78 |
|
|
|
+| | scale commercial area | | | | |
|
|
|
+| | [/km2<br>] | | | | |
|
|
|
+| PerSC | Percentage of small | 0.0 | 28.7 | 2.0 | 3.33 |
|
|
|
+| | scale commercial use | | | | |
|
|
|
+| | area[/km2<br>] | | | | |
|
|
|
+| PerSR | Percentage of single | 0.0 | 74.6 | 37.2 | 15.2 |
|
|
|
+| | family residential area<br>[/km2<br>] | | | | |
|
|
|
+| PerMR | Percentage of multiple | 0.0 | 25.1 | 3.2 | 3.5 |
|
|
|
+| | family residential area | | | | |
|
|
|
+| | [/km2<br>] | | | | |
|
|
|
+| SA2 area | The area of SA2[km2<br>] | 1.3 | 729.8 | 32.3 | 84.1 |
|
|
|
+
|
|
|
+*Note:* All variables in this study are continuous except Crashes. S.D. is standard deviation.
|
|
|
+
|
|
|
+#### **4. Methodology**
|
|
|
+
|
|
|
+We intend to establish a Semi-parametric Geographically Weighted Poisson Regression model (sGWPR) to fit spatial variations in associations between contributing factors and bicycle crashes. The model and the variable selection procedure are shown in Fig. 3. The potential factors were selected after possible multicollinearity identification by Pearson correlation and variance inflation factor (VIF) analysis. Then we put the potential factors into full Geographically Weighted Passion Regression (full GWPR) and monitored the AICc to estimate which best owns the lowest AICc value and all its parameters vary over study area. After that, we tried two approaches to establish the best sGWPR model under the selection of a smaller AICc value and the results of the Linear Model of Coregionalization (LMC), respectively. Finally, we compared the goodness of fit between these two sGWPR models using the Bayesian Information Criterion (BIC) and mean squared error (MSE). GWPR and sGWPR models were performed by the GWR 4.0, and the LMC results are from R 4.0.3 "gstat" package [\(Pebesma, 2004](#page-14-0)).
|
|
|
+
|
|
|
+### *4.1. Linear model of coregionalization (LMC)*
|
|
|
+
|
|
|
+[Nakaya et al. \(2005\)](#page-14-0) point out that because the meaning of covariates may change depending on different geographical contexts, and some relationships might be intrinsically different across space, we need to consider geographical variations in the associations of analysis. In our study, we apply the LMC to explain changed covariates of variables based on an analysis of auto-variogram and cross-variogram between variable pairs.
|
|
|
+
|
|
|
+The general statistical framework of LMC is:
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+**Fig. 3.** sGWPR modeling selection procedure. *Note*s: Pearson correlation coefficient larger than 0.6 indicates high multicollinearity ([Evans, 1996](#page-13-0)) ; VIF larger than 5 indicates a significant multicollinearity between these variables [\(Neter](#page-14-0) [et al., 1996; O](#page-14-0)'brien, 2007; [Xu and Huang, 2015; Hezaveh et al., 2019](#page-14-0)).
|
|
|
+
|
|
|
+<span id="page-5-0"></span>
|
|
|
+$$\gamma(\mathbf{h}) = \begin{bmatrix} \gamma_{11}(h) & \cdots & \gamma_{1p}(h) \\ \vdots & \ddots & \vdots \\ \gamma_{p1}(h) & \cdots & \gamma_{pp}(h) \end{bmatrix} = \sum_{d=1}^{D} C^{d} * g^{d}(\mathbf{h}) \\
|
|
|
+= \sum_{d=1}^{D} \begin{bmatrix} c_{11}^{d} & \cdots & c_{1p}^{d} \\ \vdots & \ddots & \vdots \\ c_{p1}^{d} & \cdots & c_{pp}^{d} \end{bmatrix} * g^{d}(\mathbf{h}) \tag{2}$$
|
|
|
+
|
|
|
+*γ*(*h*)represents a simple variogram matrix of two sets of variables, and it can be calculated by (3):
|
|
|
+
|
|
|
+$$\gamma_{ij}(\mathbf{h}) = \frac{1}{2N(\mathbf{h})} \sum_{(i|\hat{\mathbf{h}})=1}^{N(\mathbf{h})} (x_i(s) - x_i(s+\mathbf{h}))(x_j(s) - x_j(s+\mathbf{h}))$$
|
|
|
+(3)
|
|
|
+
|
|
|
+*h* and *s* are both vectors. If*j* ∕= *i*, *γij*(*h*)is for cross variograms ofXiandXj; if*j* = *i*,*γij*(*h*) is for the auto-variogram of one variableXi. *xi*and*xj*are observations of the variablesXiandXj(i, j = 1,…,*N*(*h*)).*N*(*h*) means the number of SA2s separated by a distance*h*, *s*and*s* + *h*are two SA2s pairs separated by a distance*h*.
|
|
|
+
|
|
|
+The LMC model is defined by*C<sup>d</sup>*(coregionalization matrix) and*g<sup>d</sup>*(*h*)variogram function at spatial scale *d* (d = 1,2…D)*.C<sup>d</sup>*are unknown non-negative matrices and corresponding coefficient matrices of variogram functions at different spatial scale *d. cd ij* is known as partial sill (psill).
|
|
|
+
|
|
|
+In this study, we considered 3 spatial scales (*d* = 0,1,2)—the distance effect without the effect of anisotropy influence. We set a linear combination of three basic functions, nugget effect (4), periodic (5) and Gaussian (6) to fit LMC model and the same basic model must appear in each auto- and cross variogram ([Isaaks and Srivastava, 1989,](#page-14-0) p390).
|
|
|
+
|
|
|
+$$g^{0}(\mathbf{h}) = \begin{cases} 0, & \text{if } \mathbf{h} = 0 \\ 1, & \text{otherwise} \end{cases} \tag{4}$$
|
|
|
+
|
|
|
+$$g^{1}(\boldsymbol{h}) = 1 - \cos(-2\pi \boldsymbol{h}/a) \tag{5}$$
|
|
|
+
|
|
|
+$$g^2(\mathbf{h}) = 1 - \exp(-3\mathbf{h}^2/b^2)$$
|
|
|
+ (6)
|
|
|
+
|
|
|
+*a* and *b* are spatial ranges.
|
|
|
+
|
|
|
+Based on coregionalization matrices*C<sup>d</sup>*, we can calculate overall spatial variance or spatial variability*Qd*in each spatial scale *d* by (2) ([Isaaks and Srivastava, 1989](#page-14-0)). From (7), we could estimate the contribution of each scale to overall variation, which could provide the contribution details of each scale to overall variability.
|
|
|
+
|
|
|
+$$Q^{d} = \sum_{d=0}^{D=2} C^{d} \tag{7}$$
|
|
|
+
|
|
|
+The spatial covariance*Qd xij*between variablesXiandXj(*j* ∕= *i*) can be calculated by summing the partial sills of cross variograms:
|
|
|
+
|
|
|
+$$Q_{x_{ij}}^{d} = \sum_{i=0}^{D-2} c_{x_{ij}}^{d}$$
|
|
|
+ (8)
|
|
|
+
|
|
|
+If*j* = *i*, (8) presents the variance of univariate variableXi, which is the sum of partial sills of Xi. Based on coregionalization matrices, we can decompose spatial variability of each variable and estimate the relative contribution (percentage) to each spatial scale *d.* The ratio of the nugget effect can present the degree of spatial dependence [\(Isaaks and Srivas](#page-14-0)[tava, 1989](#page-14-0)). For example, a LMC model can be exclusively defined by a nugget effect variogram function entailing a complete no spatial dependence (no spatial correlation). The contribution of two variables at each scale *d* could provide the importance of geographical effect in the relationship between explanatory variables and bicycle crashes.
|
|
|
+
|
|
|
+# *4.2. GWPR and sGWPR model*
|
|
|
+
|
|
|
+GWPR, allowing each parameter to vary across study area, can
|
|
|
+
|
|
|
+clearly capture some important local variations in the relationship between bicycle crash counts and explanatory variables with spatial data about socio-economic, traffic, road infrastructure and land use. The model specification of the GWPR:
|
|
|
+
|
|
|
+$$O_i \sim P\widetilde{o}isson\left[E_i exp\left(\beta_0\left(u_i, v_i\right) + \sum_{k=2}^p \beta_k\left(u_i, v_i\right) x_{ik}\right)\right] \tag{9}$$
|
|
|
+
|
|
|
+Where*Oi*is the total crashes took place 2015–2017 in SA2 *i*, (*ui,vi*) is a center coordinate of SA2 *i*. *xik* are explanatory factors. *Ei* is the exposure variable (i.e., daily vehicle kilometers traveled (DVKT) in the current study) of SA2 *i*. *p* is the number of expanding factors; The parameters*βk*(*ui, vi*)at each SA2 can be estimated by solving a maximization problem of the geographically weighted log-likelihood ([Nakaya et al.,](#page-14-0) [2005\)](#page-14-0):
|
|
|
+
|
|
|
+$$maxL(u_i, v_i) = \sum_{i=1}^{N} (- \hat{O}_j(\widehat{\beta}(u_i, v_i)) + O_j log \hat{O}_j(\widehat{\beta}(u_i, v_i))) \times w_{ij} d_{ij}$$
|
|
|
+(10)
|
|
|
+
|
|
|
+where *Oj* is the number of crashes occurring in SA2 *j*, and *O* (̂*β*(*ui, vi*) is the predicted number of crashes in SA2 *j* at regression point SA2 *i*:
|
|
|
+
|
|
|
+$$\bigwedge_{j}^{\hat{\rho}}(\widehat{\beta}(u_{i}, v_{i}) = E_{j}exp(\beta_{0}(u_{i}, v_{i}) + \sum_{k=0}^{p} \beta_{k}(u_{i}, v_{i})x_{jk})$$
|
|
|
+(11)
|
|
|
+
|
|
|
+*wij*(*j* = 1,2,…,*n*) is the weight of SA2 *j* in calibration of the model for SA2 *i*.*wij*can be calculated by adaptive bi-square as suggested by ([Li](#page-14-0) [et al., 2013\)](#page-14-0):
|
|
|
+
|
|
|
+$$w_{ij} = \begin{cases} \left(1 - d_{ij}^2 / \theta_{i(n)}^2\right)^2 d_{ij} < \theta_{i(n)}, d_{ij} > \theta_{i(n)} \\ 0 \text{ otherwise} \end{cases}$$
|
|
|
+ (12)
|
|
|
+
|
|
|
+Here,*dij*represents the Euclidean distance between SA2 *i* and SA2 *j*; *θi*(*n*)is an adaptive bandwidth size defined as the longest distance between SA2 *i* and SA2 *n*. *n* is number of near SA2s incorporated into the model in calibration for SA2 *i*.
|
|
|
+
|
|
|
+An important extension of GWPR is semi-parametric Geographically Weighted Poisson Regression (sGWPR). The model specification is:
|
|
|
+
|
|
|
+$$O_i \sim Poisson \left[ \mathbb{E}_i exp(\beta_0(u_i, v_i) + \sum_{j=2}^l \beta_j x_{ij} + \sum_{k=l+1}^p \beta_k(u_i, v_i) x_{ik}) \right] \tag{13}$$
|
|
|
+
|
|
|
+where*βj*the *j*th global parameter, which indicates no significant geographical variation across the study area. We described the selection of variable whose parameter is global in the Introduction section and more details can be found in [Nakaya et al.\(2005\).](#page-14-0)
|
|
|
+
|
|
|
+The Akaike information criterion (AIC) can evaluate prediction error and thereby it also can be a model selection indicator. When the sample size is small, however, AIC would overfit: it could select models with too many parameters. To iron out such potential overfitting, the secondorder AIC—AICc is recommended in this study. More details can be found in [Nakaya et al. \(2005\).](#page-14-0)
|
|
|
+
|
|
|
+### *4.3. Model performance assessment*
|
|
|
+
|
|
|
+The goodness of fit of the GWR models was assessed through the Akaike information criterion (AICc) ([Fotheringham et al., 2002, Gomes](#page-13-0) [et al., 2019, Hezaveh et al., 2019](#page-13-0)). Unlike the AIC, the Bayesian Information Criterion (BIC) penalizes the model more for its complexity. The lower the BIC, the less complex the model. The BIC is calculated as follow:
|
|
|
+
|
|
|
+$$BIC = -2LL + log(N)K (14)$$
|
|
|
+
|
|
|
+where *LL* is the log-likelihood of the model, *N* is the number of SA2s, *K* is the number of parameters in the model.
|
|
|
+
|
|
|
+Mean Squared Error (MSE), assessing the estimation error, is typi-
|
|
|
+
|
|
|
+cally used in other traffic crash studies for model comparison [\(Arvin](#page-13-0) [et al., 2019, Amoh-Gyimah et al., 2016, Lu et al., 2014](#page-13-0)). We also used MSE to compare the performance between sGWPR models under monitoring AICc results and under LMC results, respectively. The smaller value of MSE indicates better performance of model estimation. MSE has the form:
|
|
|
+
|
|
|
+$$MSE = \frac{\sum_{i=1}^{N} |\hat{y}_{i} - y_{i}|^{2}}{N}$$
|
|
|
+ (15)
|
|
|
+
|
|
|
+where*yi*is the crash count of SA2 *i*, *yi* is the expected crash count of each SA2 *i* obtained by GWR models, and *N* is the number of SA2s.
|
|
|
+
|
|
|
+# **5. Results**
|
|
|
+
|
|
|
+# *5.1. Variable selection and model selection for a full GWPR*
|
|
|
+
|
|
|
+Through Pearson correlation analysis, if the correlation of two variables is higher than 0.6, then the two were maintained in separate models but not together. After the Pearson correlation analysis, we monitored the AICc value, we selected a full GWPR model with smallest AICc = 905.44 [\(Table 3\)](#page-7-0) and incorporated 10 explanatory variables (Table 2). The VIF values range from 1.21 to 2.02 and the Pearson correlations fall into the range of − 0.49 to 0.50, which indicates no significant multicollinearity in our models (shown in Table 2).
|
|
|
+
|
|
|
+# *5.2. Variable selection for sGWPR model with monitoring AICc values*
|
|
|
+
|
|
|
+When applying GWPR models, we tried adaptive bi-square to determine the best bandwidth. The best bandwidth from the adaptive bisquare method represents the number of SA2s which were incorporated to estimate the GWPR models. Now, we determine which parameter of variable should vary geographically by comparing the AICc values of two models: one is the full GWPR model (M0), another is a sGWPR model in which all parameters except one vary geographically (M1- M10). The comparison is shown in [Table 3](#page-7-0).
|
|
|
+
|
|
|
+By comparing the AICc value of the full GWPR model M0 with that of the sGWPR model (M1-M10), we can find a sGWPR model with only invariant Sidewalk density (Δ*AIC*cof M5 = -3.03), Bicycle lane intersection density (BLID,Δ*AIC*cof M7 = -3.81 or PerSC (Δ*AIC*cof M8 = -5.25) would get better than the full GWPR model (AICc of M0 = 905.44) ; the sGWPR model with only one variable constant (M1-M4, M6 and M9-M10) is a worse model, respectively (the AICc ranges from 916.80 to 1145.07). According to the AICc results and the specification of variable choice for the sGWPR model in [Nakaya et al. \(2005\),](#page-14-0) if the sGWPR model with only invariant variable (M1-M10) could reduce the AICc value at least 2 comparing with the AICc value of full GWPR M0, the invariant variable should be incorporated into the final optimal model M11. Thus, the optimal model (M11 with AICc = 896.65 in [Table 3](#page-7-0)) is a semi-parametric one, in which parameters associated with three variables (Sidewalk density, BLID and PerSC) are held constant and other parameters associated with the remaining explanatory variables are
|
|
|
+
|
|
|
+allowed to vary geographically.
|
|
|
+
|
|
|
+## *5.3. Variable selection for sGWPR model with LMC results*
|
|
|
+
|
|
|
+Due to different scale of variable measurements, the scale of psills in a coregionalization matrix from LMC can be dramatically different. That leads to no general metric to estimate the relative spatial variability contribution on overall variability, mentioned in section 4.1. Thus, when estimating a spatial variability contribution, all explanatory variables here are standardized to have zero mean and unit variance in the LMC estimations.
|
|
|
+
|
|
|
+89.1% of bicycle crashes involves vehicles in our dataset. In addition, [Fournier et al. \(2019\)](#page-13-0) claim the bicycle crash rate is related to vehicle and bicycle exposure simultaneously; but we lack bicycle exposure, so Bicycle crashes is dealt with the bicycle crash rate per vehicle per kilometer travelled (BCR) in the LMC estimation. BCR in SA2 *i* is calculated by:
|
|
|
+
|
|
|
+$$BCR_i = \frac{y_i}{E_i} \tag{16}$$
|
|
|
+
|
|
|
+where*yi*is an observed number of bicycle crashes in SA2 *i* and*Ei*is an exposure variable i.e., daily vehicle kilometers traveled (DVKT) of SA2 *i*.
|
|
|
+
|
|
|
+In this study, we applied sample variograms to describe how the spatial continuity changes (*γ*(*h*)) with distance, along with the accompanying details in [Fig. 4](#page-7-0) (dots). We can find the periodic trends are tense and dramatic in our data set, and these trends are different variogram shapes from the general one where the*γ*(*h*)increases slightly and erratically as distance [\(Isaaks and Srivastava, 1989](#page-14-0)).
|
|
|
+
|
|
|
+Due to the periodic behavior in the simple variograms, we introduce the periodic function to capture this behavior. The linear combination of LMC in this study is:
|
|
|
+
|
|
|
+$$\gamma(h) = Nug(0) + Per(40) + Gau(59)$$
|
|
|
+ (17)
|
|
|
+
|
|
|
+*Nug*(0)is the nugget effect,*Per*(40)is a smaller spatial scale where the LMC is fitted with periodic variogram function within a 40 km range, and*Gau*(59)is a larger spatial scale where the LMC is fitted with a Gaussian variogram function within the 40 km and 59 km range. The fitted LMC model is depicted in [Fig. 4](#page-7-0) (continuous solid line). Based on coregionalization matrices, the spatial variability of overall and each variable was decomposed at different scales and presented in the [Table 4](#page-9-0). The higher ratio of the nugget effect can present the lower spatial dependence ([Isaaks and Srivastava, 1989](#page-14-0)), thus a variable with high nugget effect should be invariant.
|
|
|
+
|
|
|
+The first row shows the nugget, periodic and gaussian effects capture at 48.7%, 12.0%, 39.3% of the overall variation respectively. The 51.3% overall variation is explained by a variance appearing in the nearby SA2s up to 59 km. The second row shows that 30.8% of BCR (bicycle crash rate) variability is explained by the nugget effect; the remaining variability, 10.0% is explained by periodic behavior within a smaller scale (up to 40 km), and 59.2% is interpreted by the variance in the range of 40–59 km. BCR variance represents greater spatial dependence, which implies that the GWPR is an appropriate method to capture the variable
|
|
|
+
|
|
|
+**Table 2** Pearson correlation matrix and VIF for the explanatory variables used in the study.
|
|
|
+
|
|
|
+| Type | Pearson correlation | | | | | | | | | | VIF |
|
|
|
+|-----------------------------|---------------------|------|--------|--------|--------|------|--------|--------|------|------|------|
|
|
|
+| Variables | A | B | C | D | E | F | G | H | I | J | |
|
|
|
+| Per14 (A) | 1.00 | | | | | | | | | | 1.92 |
|
|
|
+| Mincome (B) | − 0.02 | 1.00 | | | | | | | | | 1.21 |
|
|
|
+| Ln(pop) (C) | 0.31 | 0.12 | 1.00 | | | | | | | | 1.60 |
|
|
|
+| Perbyk (D) | − 0.44 | 0.43 | 0.05 | 1.00 | | | | | | | 1.45 |
|
|
|
+| Sidewalk density (E) | − 0.27 | 0.21 | 0.04 | 0.33 | 1.00 | | | | | | 1.31 |
|
|
|
+| Residential road density(F) | − 0.27 | 0.12 | 0.12 | 0.50 | 0.32 | 1.00 | | | | | 1.67 |
|
|
|
+| BLID (G) | − 0.39 | 0.31 | − 0.02 | 0.31 | 0.33 | 0.49 | 1.00 | | | | 1.91 |
|
|
|
+| PerSC (H) | − 0.49 | 0.28 | 0.12 | 0.26 | 0.32 | 0.36 | 0.50 | 1.00 | | | 2.02 |
|
|
|
+| PerSR (I) | 0.37 | 0.02 | 0.33 | − 0.18 | − 0.20 | 0.16 | − 0.10 | − 0.21 | 1.00 | | 1.42 |
|
|
|
+| PerMR (J) | − 0.22 | 0.17 | 0.18 | 0.21 | 0.06 | 0.40 | 0.31 | 0.26 | 0.10 | 1.00 | 1.32 |
|
|
|
+
|
|
|
+<span id="page-7-0"></span>**Table 3** Comparison of Negative binomial (NB), Full GWPR (M0) and sGWPR (M1-M12) models.
|
|
|
+
|
|
|
+| Model | Explanation | Best bandwidth size n/θ | AICc | ΔAICc | MSE | BIC |
|
|
|
+|---------|-------------------------------------|-------------------------|---------|---------|--------|---------|
|
|
|
+| Poisson | Poisson regression | None | 2202.27 | 1296.83 | 65.56 | 2246.60 |
|
|
|
+| NB | Negative binomial | None | 1790.10 | 884.66 | 159.55 | 1830.28 |
|
|
|
+| M0 | Full GWPR | 58/56.6 km | 905.44 | 0 | 36.17 | 1180.47 |
|
|
|
+| M1 | sGWPR_ Per14 | 110 | 1145.07 | 239.63 | / | / |
|
|
|
+| M2 | sGWPR_Mincome | 58 | 938.28 | 32.84 | / | / |
|
|
|
+| M3 | sGWPR_Ln(pop) | 47 | 1093.67 | 188.23 | / | / |
|
|
|
+| M4 | sGWPR_Perbyk | 54 | 916.80 | 11.36 | / | / |
|
|
|
+| M5 | sGWPR_Sidewalk density | 57 | 902.41 | − 3.03 | / | / |
|
|
|
+| M6 | sGWPR_Residential road density | 57 | 904.83 | − 0.61 | / | / |
|
|
|
+| M7 | sGWPR_BLID | 57 | 901.63 | − 3.81 | / | / |
|
|
|
+| M8 | sGWPR_PerSC | 58 | 900.19 | − 5.25 | / | / |
|
|
|
+| M9 | sGWPR_PerSR | 62 | 955.58 | 50.14 | / | / |
|
|
|
+| M10 | sGWPR_PerMR | 54 | 991.37 | 85.93 | / | / |
|
|
|
+| M11# | sGWPR_Sidewalk density, BLID, PerSC | 46/55.1 km | 896.65 | − 8.79 | 37.34 | 1168.45 |
|
|
|
+| M12# | sGWPR_Sidewalk density, PerSC | 54/55.7 km | 895.98 | − 9.46 | 28.94 | 1163.22 |
|
|
|
+
|
|
|
+*Notes:* sGWPR\_X indicates a semi-parametric geographically weighted Poisson regression model in which only the parameter of the variable X is invariant. Δ*AIC*cis the difference of AICc of a model from that of the full GWPR(M0). M11# is the model under AICc results, while M12# is the model under LMC results. the third column *n/ θ*represents the number of samples incorporated into models and an adaptive bandwidth size defined as the longest distance between two samples (in equation [\(12\)\)](#page-5-0).
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+**Fig. 4.** Plots of auto-variograms of variables(dots) and the LMC fitted line. *Note:* The horizonal axes represent distance *h* ranging from 0 to 59 km.
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+**Fig. 4.** (*continued*).
|
|
|
+
|
|
|
+variance.
|
|
|
+
|
|
|
+For the explanatory variables, we can see the Sidewalk density and PerSC are more monopolized by the nugget effect (higher ratio over 60%) while remaining variables are less dominated by the nugget effect (lower ratio less than 40%). For the variables with a lower nugget ratio, the spatial dependence of Per14, Mincome and Ln (pop) is more staged in the smaller spatial scale, ranging up to 40 km (over 40%). The spatial dependence of Residential road density, BLID, PerSR (Percentage of
|
|
|
+
|
|
|
+<span id="page-9-0"></span>**Table 4** Spatial variability contribution (%) of variable on overall variability at each scale.
|
|
|
+
|
|
|
+| Variable | Nugget | Periodic(0,40<br>km) | Gaussian(40 km,59<br>km) |
|
|
|
+|-----------------------------|----------------|----------------------|--------------------------|
|
|
|
+| Overall | 48.7 | 12.0 | 39.3 |
|
|
|
+| BCR (bicycle crash<br>rate) | 30.8 | 10.0 | 59.2 |
|
|
|
+| Per14 | 29.1(5.0) | 41.2(41.8) | 29.7(53.2) |
|
|
|
+| Mincome | 39.6<br>(12.3) | 49.4(65.9) | 11.0(21.8) |
|
|
|
+| Ln(pop) | 24.4<br>(36.9) | 74.9(29.2) | 0.7(33.9) |
|
|
|
+| Perbyk | 31.0<br>(34.5) | 37.1(14.9) | 31.9(50.6) |
|
|
|
+| Sidewalk density | 69.2<br>(23.0) | 30.2(65.8) | 0.6(11.2) |
|
|
|
+| Residential road<br>density | 33.5<br>(86.0) | 2.5(3.3) | 64.0(10.7) |
|
|
|
+| BLID | 31.6<br>(36.2) | 6.0(15.2) | 62.4(48.6) |
|
|
|
+| PerSC | 62.6<br>(12.2) | 9.8(9.2) | 27.6(78.6) |
|
|
|
+| PerSR | 38.7<br>(16.6) | 11.4(38.0) | 49.9(45.4) |
|
|
|
+| PerMR | 33.9(7.1) | 17.9(6.2) | 48.2(86.7) |
|
|
|
+
|
|
|
+*Note:* Values in brackets represent spatial variability contribution (%) between explanatory variables at each scale.
|
|
|
+
|
|
|
+single residential area) and PerMR (Percentage of multiple residential area) is mainly presented in the range of 40–59 km (over 48%), and the spatial dependence of Perbyk equally appears in these two spatial scales (37.1% in the smaller spatial scale and 31.9% in the larger spatial scale).
|
|
|
+
|
|
|
+The size of the non-nugget effect (ranging up to a 40 km and range of 40 km-59 km) represents spatial dependence. The lower non-nugget effect is in Sidewalk density and PerSC (30.8% and 37.4%, respectively); the higher non-nugget effect is in the remaining variables: Per14, Mincome, Perbyk, Ln (pop), Residential road density, BLID, PerSR and PerMR (from 61% to 75%). These variables with higher non-nugget effect should be allowed to vary and the ones with lower non-nugget effect should be held constant in the sGWPR model to explain the geographical variation in BCR. From the spatial variances between BCR and explanatory variables (the values in brackets in Table 4), we can find non-nugget effect ranges from 14% (variance with Residential road density) to 95% (variance with Per14).
|
|
|
+
|
|
|
+Thus, based on the nugget effect, we fitted a sGWPR model (M12 in [Table 3\)](#page-7-0) in which the parameters associated with Sidewalk density and PerSC are held constant, and the parameters associated with remaining variables vary spatially.
|
|
|
+
|
|
|
+# **6. Discussion**
|
|
|
+
|
|
|
+# *6.1. Model performance comparation and parameter estimation*
|
|
|
+
|
|
|
+Before the comparison of sGWPR under AICc values and under LMC results, we run the non-spatial Poisson model (Negative binomial model, NB and Poisson model) for comparison. From [Table 3](#page-7-0), the traditional NB and Poisson model perform worser than the GWR models since they have higher AICc, MSE, and BIC values. The worst performance of NB indicates the adaptability of spatial models (GWPR and sGWPR models), which is consistent with many previous studies, such as [Xu and Huang](#page-14-0) [\(2015\).](#page-14-0)
|
|
|
+
|
|
|
+The current research relied on the AICc value to identify the spatial association in the sGWPR model. Based on AICc results, it suggested a sGWPR in which Sidewalk density, BLID and PerSC are held constant and the remaining variables vary spatially in this study. No evidence suggests that the spatial relationship between BLID and bicycle crashes would vary over space. However, LMC results provide such evidence suggesting the spatial varying association between BLID and bicycle crashes, since 69.4% BLID variability and 63.8% variance between BLID and bicycle crashes can be explained by spatial dependence at two spatial scales (ranging up to 40 km and a range of 40 km-59 km). The difference of M11 and M12 in AICc values is less than 2 [\(Table 3\)](#page-7-0) and that means the M12 (sGWPR under LMC results) does not reduce the performance from AICc perspective. We also applied BIC and MSE to measure model performance of M11 and M12; however, their results show the M12 has lower BIC (1163.22, − 5.23) with less complexity, lower MSE (278.94, − 8.4) than M11 and its accuracy improves by 22.5% ((37.34–28.94/37.34) (see in [Table 3](#page-7-0)).
|
|
|
+
|
|
|
+Table 5 summarizes model estimation results of the Full GWPR(M0) and sGWPR (M11 and M12). In Full GWPR and sGWPR models, where different crash prediction models are developed for each individual SA2, the parameters are described by the mean, minimum (Min), lower quartile (Lower), median, upper quartile (Upr), and maximum (Max) of values.
|
|
|
+
|
|
|
+Two general observations are worthy of notice: [\(1\)](#page-3-0) With respect to the five statistics indicators (Mean, Min, Lower, Median, Upr and Max), the direction (either positive or negative relationship) of relationships between crashes and each explanatory variable across SA2s is consistent in the corresponding counterparts in M0, M11and M12. For example, the parameters of Per14 in M0, M11and M12 are all negative or negative in each SA2, although the signs of a few parameters changed from negative to positive or from positive to negative; [\(2\)](#page-5-0) The varying parameters of each variable in M0 always fall into the range of
|
|
|
+
|
|
|
+**Table 5** Model estimations of Full GWPR(M0) and sGWPR (M11 and M12).
|
|
|
+
|
|
|
+| M0-Full GWPR | | | | | | |
|
|
|
+|-------------------|--------------|----------|--------|--------|--------|--------|
|
|
|
+| Variable | Mean | Min | Lower | Median | Upr | Max |
|
|
|
+| Intercept | − 6.45 | − 20.03 | − 9.11 | − 7.26 | − 3.97 | 10.53 |
|
|
|
+| Per14 | − 1.70 | − 24.95 | − 8.34 | − 0.14 | 3.27 | 25.75 |
|
|
|
+| Mincome | 1.22 | − 2.91 | 0.18 | 0.97 | 1.97 | 11.54 |
|
|
|
+| Ln(pop) | 0.10 | − 1.43 | 0.02 | 0.15 | 0.35 | 1.20 |
|
|
|
+| Perbyk | 29.13 | − 34.87 | 10.25 | 18.87 | 37.97 | 317.08 |
|
|
|
+| Sidewalk density | 0.03 | − 1.05 | − 0.05 | 0.03 | 0.15 | 0.71 |
|
|
|
+| Residential road | 0.06 | − 0.32 | 0.02 | 0.06 | 0.12 | 0.28 |
|
|
|
+| density | | | | | | |
|
|
|
+| BLID | 0.01 | − 0.74 | − 0.06 | 0.02 | 0.07 | 0.76 |
|
|
|
+| PerSC | 7.24 | − 42.96 | 2.13 | 6.56 | 13.32 | 35.27 |
|
|
|
+| PerSR | 0.53 | − 3.09 | − 0.62 | 0.93 | 1.82 | 5.88 |
|
|
|
+| PerMR | 5.88 | − 15.23 | − 1.40 | 3.32 | 14.06 | 32.31 |
|
|
|
+| M11-sGWPR | | | | | | |
|
|
|
+| Variable | Mean | Min | Lower | Median | Upr | Max |
|
|
|
+| Intercept | − 6.39 | − 16.74 | − 9.95 | − 6.83 | − 3.50 | 12.36 |
|
|
|
+| Per14 | − 2.25 | − 28.45 | − 9.29 | − 1.28 | 3.23 | 33.98 |
|
|
|
+| Mincome | 1.28 | − 5.02 | − 0.04 | 1.11 | 2.28 | 9.79 |
|
|
|
+| Ln(pop) | 0.10 | − 1.81 | − 0.01 | 0.18 | 0.42 | 0.97 |
|
|
|
+| Perbyk | 22.48 | − 302.04 | 7.71 | 18.17 | 47.81 | 253.73 |
|
|
|
+| Sidewalk density# | 0.04* (5.83) | | | | | |
|
|
|
+| Residential road | 0.05 | − 0.35 | 0.01 | 0.05 | 0.11 | 0.27 |
|
|
|
+| density | | | | | | |
|
|
|
+| BLID# | 0.05* (5.48) | | | | | |
|
|
|
+| PerSC# | 5.08* (8.38) | | | | | |
|
|
|
+| PerSR | 0.56 | − 4.45 | − 0.66 | 0.67 | 1.98 | 6.12 |
|
|
|
+| PerMR | 5.95 | − 20.45 | − 3.30 | 4.70 | 15.22 | 65.14 |
|
|
|
+| M12- sGWPR | | | | | | |
|
|
|
+| Variable | Mean | Min | Lower | Median | Upr | Max |
|
|
|
+| Intercept | − 6.71 | − 21.07 | − 9.47 | − 7.09 | − 3.98 | 10.38 |
|
|
|
+| Per14 | − 1.86 | − 20.65 | − 7.67 | − 0.58 | 2.72 | 29.07 |
|
|
|
+| Mincome | 1.34 | − 3.33 | 0.10 | 0.98 | 2.16 | 10.58 |
|
|
|
+| Ln(pop) | 0.13 | − 1.65 | 0.03 | 0.17 | 0.39 | 1.28 |
|
|
|
+| Perbyk | 27.44 | − 102.55 | 9.75 | 17.95 | 43.63 | 233.29 |
|
|
|
+| Sidewalk density# | 0.02*(4.16) | | | | | |
|
|
|
+| Residential road | 0.06 | − 0.36 | 0.02 | 0.06 | 0.12 | 0.24 |
|
|
|
+| density | | | | | | |
|
|
|
+| BLID | 0.02 | − 0.64 | − 0.07 | 0.03 | 0.11 | 1.20 |
|
|
|
+| PerSC# | 5.75* (8.90) | | | | | |
|
|
|
+| PerSR | 0.47 | − 3.94 | − 0.75 | 0.65 | 1.80 | 5.61 |
|
|
|
+| PerMR | 5.94 | − 16.55 | − 1.97 | 3.64 | 13.41 | 36.34 |
|
|
|
+
|
|
|
+*Note*: #The variable holds a constant coefficient across the study area. \*The coefficient is significant at 10% level.
|
|
|
+
|
|
|
+<span id="page-10-0"></span>corresponding counterparts in M11 and M12; similarly, the varying parameters of each variable in M12 always fall into the range of corresponding counterparts in M11. This might indicate as the number of variables with fixed parameters increases, the effects of the variables with varying parameters on crash occurrence estimated by the GWPR technique would be enhanced.
|
|
|
+
|
|
|
+The reversal of relationship is a common counterintuitive signs problem in GWPR ([Chow et al., 2006, Li et al., 2013, Huang et al., 2018](#page-13-0)). One reason may be that counterintuitive signs of variables are insignificant in some zones due to the possible basis of calibrating sGWPR[\(Xu](#page-14-0) [and Huang, 2015, Munira and Sener, 2020, Hezaveh et al., 2019, Zhai](#page-14-0)
|
|
|
+
|
|
|
+[et al., 2018a\)](#page-14-0). This can be measured by mapping t-value with 1.65 at a 10% confidence level. For avoiding false positives or negatives and the Simpon's Paradox [\(Tu et al., 2008\)](#page-14-0), we corrected the t-value at a 10% confidence level as proposed by [da Silva and Fotheringham \(2016\)](#page-13-0). *α*can be calculated by equation [\(14\):](#page-5-0)
|
|
|
+
|
|
|
+$$\alpha = \frac{\zeta}{p_e} p \tag{18}$$
|
|
|
+
|
|
|
+*α*is the corrected*α*, *ζ*is the usual*α*, *pe*is the effective number of independent parameters, *p*is the is the number of parameters. For the results of M13, *ζ*,*pe*and*p*is 0.10, 81.73, and 9, so the corrected*α*is 0.011. So we
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+**Fig. 5.** Maps of local parameters and the significance of (a) Per14, (b) Mincome, (c) Ln (pop), (d) Perbyk, (e) Residential road density, (f) Bicycle lane intersection density (BLID), (g) Percentage of single-family residential area (PerSR) and (h) Percentage of multiple-family residential area (PerMR) in sGWPR M13.
|
|
|
+
|
|
|
+
|
|
|
+
|
|
|
+**Fig. 5.** (*continued*).
|
|
|
+
|
|
|
+mapped the t-value with 2.575 at the 10% level of confidence instead of the value of 1.65 in [Fig. 5.](#page-10-0) In this study, we only discuss the significant relationships.
|
|
|
+
|
|
|
+After the correction, some significant parameters still exist with reversal signs for a variable. [Hadayeghi et al. \(2010\)](#page-14-0) ascribed the issue to the local multicollinearity, however, effective diagnostic tools have not been developed to examine this issue up to now. Another explanation might be that unobserved factors which affect the traffic conditions are not considered. This explanation for each counterintuitive relationship with bicycle crashes is discussed in section 6.3.
|
|
|
+
|
|
|
+*6.2. Specification of spatial scale analysis in bicycle crash variance by LMC*
|
|
|
+
|
|
|
+In this study, we explained the bicycle crash variance based on an analysis of auto-variograms and cross-variograms between variable pairs ([Fig. 4\)](#page-7-0). For the spatial scale dependence analysis, we used the Mean Center tool in ArcGIS 10.2 to produce location of mean center of bicycle crashes occurring in each SA2 and considered the distance between mean centers in estimating auto-variograms and crossvariograms. We simplified the distance calculation, and it may make some bias in the estimation of variograms.
|
|
|
+
|
|
|
+We established the Linear Model of Coregionalization (LMC) to fit the variogram and estimated the coregionalization matrix or psill in R 4.0.3 with an R "gstat" package suggested by [Pebesma \(2004\)](#page-14-0). Here we do a succinct introduction for this estimation: first, we fitted each variogram model to a direct or cross variogram (there are 20 basic variogram models in "gstst" package); then each of the partial sill coefficient matrices was approached by its (in a least square sense) closest positive definite matrices; finally, the "best reasonable" LMC was chosen in this way.
|
|
|
+
|
|
|
+There are some practical difficulties in the choice of reasonable variogram models and spatial ranges to fit the spatial patterns. We found the periodic behavior is tense and dramatic in almost variable pairs ([Fig. 4\)](#page-7-0). We applied each variogram function to fit the periodic behavior, and found the periodic function was the reasonable one which could capture this behavior. This periodic behavior appearance may associate with the geographical distribution of population, road facilities and land use, such as urban centers. Applying the *fit. variogram* model in the "gstst" package, spatial range a (equation [\(6\)\)](#page-5-0) in the LMC estimation was determined, and we found the range of 40 km is more reasonable in the LMC estimation. Similarly, in Gaussian function, 59 km was chosen as spatial range b (equation [\(7\)\)](#page-5-0). Apart from this, an adaptive bandwidth size of three models M0, M11, M12 (56.6 km, 55.1 km, 55.7 km, respectively in [Table 4](#page-9-0)) is less than 59 km. Thus, we argued the range up to 59 km could reasonably capture the spatial dependence of variable pairs.
|
|
|
+
|
|
|
+## *6.3. Specification of the sGWPR model with LMC results*
|
|
|
+
|
|
|
+The constant parameters of Sidewalk density and Percentage of small-scale commercial use area (PerSC) are 0.02 (t-value = 4.16) and 2.75 (t-value = 8.90), respectively. We can find as these two variables increase, the number of bicycle crashes increase (significant at 10% level). Higher Sidewalk density potentially means the more walking spaces where there might be more cycling mixing with walking; thus, it is better to separate bicycles and pedestrians reducing the deterrence of real or perceived safety concerns since 60 percent of Victorians are curious about cycling (according to the report of Transport for [Victoria](#page-14-0) [\(2018\)](#page-14-0).The effect of PerSC on bicycle crashes is in accordance with conclusions of studies of [Dumbaugh and Li \(2010\) and Cai et al. \(2016\)](#page-13-0). Shopping strips or sidewalk alongside arterial roads has intense pedestrian activity, and it is a potential conflict area between bicycles and pedestrians, or even other vehicles.
|
|
|
+
|
|
|
+The results for the varied parameters associated with the remaining explanatory variables in sGWPR model under LMC results is shown in [Fig. 5.](#page-10-0) According to the Metropolitan Transport Plan, the Greater Melbourne is divided into three areas: the inner area (shown in [Fig. 5](#page-10-0)), middle area, and outer area. The significant parameters of Per14 are mainly concentrated in the inner city and partial outer areas (Melbourne West and Mornington Peninsula). The increase in percentage of the young (under 14 years old) decreases the number of bicycle crashes of the inner city or outer area. This is consistent with and supports the findings in which an increase in percentage of the younger would reduce traffic crashes ([Boufous et al., 2012, Ladron](#page-13-0) ´ de Guevara et al., 2004). It is mainly because those under 14 years old do not cycle for commuting as much as other age cohorts [\(Amoh-Gyimah et al., 2016\)](#page-13-0).
|
|
|
+
|
|
|
+The significant parameters of mediate personal weekly income (Mincome) are concentrated in the middle area and part of the outer area while parameters of Mincome are insignificant in the inner area ([Fig. 5](#page-10-0) (b)). From the Australian Bureau of Statistics 2013, we compared the median household income between the inner area and the rest of the area, and found the inner area has a higher median household income. lower income areas have a higher probability of bicycle crashes. Perhaps people with lower income are more likely to cycle. In a higher income area, there appears to be an insignificant relationship between income and bicycle crashes. [Siddiqui et al \(2012\)](#page-14-0). Estimated a Bayesian model of fixed global parameters to accounting for spatial correlation and found median household income has an insignificant effect on bicycle crashes. Maybe such a global method covers the spatially varying relationships. Since all the significant parameters are positive, we claim the Mincome would increase bicycle crash occurrence.
|
|
|
+
|
|
|
+The two exposure variables log of population (Ln (pop)) and percentage of commuters cycling to work (Perbyk) both have a positive and significant association with bicycle crashes in majority SA2s: 94% and 89% SA2s with significant effect are positive ([Fig. 5](#page-10-0) (c) and (d)). Our findings support the previous studies ([Amoh-Gyimah et al., 2016, Kim](#page-13-0) [et al., 2007\)](#page-13-0).
|
|
|
+
|
|
|
+Residential road density displays a significantly positive effect in more than half of the study area, and 96% of SA2s with significant relationships are positive [\(Fig. 5](#page-10-0) (e)). This result is consistent with the finding of [Dumbaugh and Li \(2010\)](#page-13-0), which states that residential street increases crash incidence because wider rights-of-way in residential streets lead to higher speeds ([Smith and Appleyard, 1981\)](#page-14-0). Research on residential street design reducing bicycle crash occurrence is further needed.
|
|
|
+
|
|
|
+The sGWPR model with AICc measurement does not suggest that the relationship between bicycle lane intersection density (BLID) and bicycle crashes should vary over space. According to the suggestion of our approach, however, the association of BLID and bicycle crashes should vary over space, because the variability of BLID can be explained by 68.4% of spatial variability patterns, and the variance between BLID and bicycle crashes can also be explained by 63.8% of spatial variability. Our approach applies a geostatistical tool capturing details in the spatial variations of associations in bicycle crash analysis, which cannot be captured by only measuring the AICc values. BLID, a newly defined measurement proposed by us, could better explain the bicycle crash occurrence. Intersection density (the ratio of the number of intersections and road length in each unit) is a common variable to explain crash occurrence. We incorporate intersection density, instead of BLID, into the original model M12 and find that it performs worse than the original model M12 (AICc of 913.42 vs 895.98).
|
|
|
+
|
|
|
+BLID is significantly associated with bicycle crashes in the Inner area and the east of Inner area (78% of SA2s with significant relationships) (Shown in [Fig. 5](#page-10-0) (f)). As the number of bicycle lane intersections increases, bicycle crashes are likely to increase. Bicycle lane intersections are potential conflict areas between bicycles and other vehicles. According to the Metropolitan Transport Plan, the government plans to continue building the Principal Bicycle Network (PBN) on arterial roads. We suggest the government should also establish a program to improve visibility using green paving or improve the standards for bicycle facilities by separating cyclists from traffic streams at the intersections of arterial roads where the PBN is built.
|
|
|
+
|
|
|
+In some SA2s which are located in the Mornington Peninsula and the south of the Greater Melbourne, however, BLID is significantly associated with a decrease in bicycle crashes. It seems counterintuitive but is also observed in [Huang et al \(2018\)](#page-14-0) and is reasonable. The reason might be that the data of unobserved factors that affect the traffic or infrastructure condition of intersections in each SA2s were not collected. For example, better traffic control signals may help ease traffic conflicts at intersections to reduce bicycle crashes. Other factors may include intersection configuration or intersection sight distance, etc. More detailed traffic control data and other related data are needed to help confirm these hypotheses. Another possibility is the effect of BLID on bicycle crashes might intrinsically vary across space. For instance, people's preferences and attitudes for crossing intersections might vary over space due to factors which are enormously complex and impossible to identify in the form of a single explanatory variable, such as cycling culture, habits and customs.
|
|
|
+
|
|
|
+The significant effect of PerSR or PerMR on bicycle crashes vary over space both only in strength but also in direction: the relationships are mainly positive. The PerSR (percentage of single-family residential areas) has a positive effect on bicycle crash occurrence in 60.8% of SA2s. 75% of SA2s with significant relationships between PerSR and bicycle <span id="page-13-0"></span>crashes are positive. We observe the distribution of SA2s with significantly positive effect. The majority of these SA2s are concentrated on the Inner area and part of Western Melbourne and Eastern Melbourne (shown in [Fig. 5](#page-10-0) (g)). The PerMR (percentage of multiple-family residential area) shows a positive relationship with bicycle crashes in 67% of SA2s. 72% of SA2s with significantly positive relationships between PerMR and bicycle crashes are positive. These SA2s are mainly located in the Southern Melbourne, shown in [Fig. 5](#page-10-0) (h). These findings support the [Kim and Yamashita \(2002\)](#page-14-0), who showed that crashes occur in areas with mixed residential and commercial land use. We also observe that the PerMR or PerSR has a significant influence on bicycle crashes in a small part of the study area. It might be that there is more occurrence of disconnected residential subdivisions in some SA2s. Disconnected residential subdivisions are reported fewer crashes than residential communities with an accessibility to heavy traffic ([Marks, 1957, Merlin](#page-14-0) [et al., 2020\)](#page-14-0).
|
|
|
+
|
|
|
+# **7. Conclusion**
|
|
|
+
|
|
|
+Our study of crashes aims to identify the social economic and built environment factors in the relationship with bicycle crashes over a Metropolitan area. We modeled spatial scale patterns by the Linear Model of Coregionalization (LMC) method. Its results provided evidence, combined with the sGWPR approach, to determine what spatial local (its parameters varying) or global (its parameters held constant) association is in the bicycle crash occurrence. We found that there is no performance degradation between GWPR models enhanced by LMC results and under AICc values, but there are yields some insights into spatial association and model specification in crash studies. This model uncovers spatial and significant associations between bicycle crashes and bicycle lane intersection density (BLID) as well as characteristics of traffic exposure, social economics, road networks and land use in Greater Melbourne. The effect of residential road density, percentage of single-family area, percentage of multiple-family area, and BLID positively varies across the space, and Sidewalk density and Percentage of small-scale commercial use area (PerSC) with constant parameters are positively associated with bicycle crashes.
|
|
|
+
|
|
|
+### **8. Limitations**
|
|
|
+
|
|
|
+Due to the limited methods and access to data sources, we must acknowledge two limitations. One is that the Traffic volume data (2017) does not fully cover the crash data (2015–2017). The possible change of traffic flow during the three years might influence the accuracy of estimation. Another limitation is that [Fig. 4](#page-7-0) shows different variogram shapes with tense and dramatic erratic in this study to the general trend with visible increase of *γ*(*h*)with distance, with slight periodic trends ([Isaaks and Srivastava, 1989](#page-14-0)). [Marcotte \(2012\)](#page-14-0) argued that more variogram function structures are needed to fit different variogram shapes in the LMC estimation. That would be our further research.
|
|
|
+
|
|
|
+### **CRediT authorship contribution statement**
|
|
|
+
|
|
|
+**Shujuan Ji:** Conceptualization, Methodology, Software, Writing original draft, Visualization. **Yuanqing Wang:** Supervision, Validation. **Yao Wang:** Visualization, Writing - review & editing.
|
|
|
+
|
|
|
+# **Declaration of Competing Interest**
|
|
|
+
|
|
|
+The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
|
|
|
+
|
|
|
+### **Acknowledgements**
|
|
|
+
|
|
|
+The authors are grateful to Australian Bureau of Statistics (ABS) for
|
|
|
+
|
|
|
+sharing data and we also thank to Runing Ye, PhD, for her contributions at the early stages of this research. The authors are also grateful for the support from National Natural Science Foundation of China (51878062, 51908462), the Natural Science Foundation of Shaanxi Province (2020JM-246), the 111 project of Sustainable Development of Transportation in Western Urban Agglomeration (No. B20035) and the Fundamental Research Funds for the Central Universities, Chang'an University (300102210214).
|
|
|
+
|
|
|
+# **References**
|
|
|
+
|
|
|
+- ABS (Australian Bureau of Statistics). (2017). 2016 census. https://datapacks. censusdata.abs.gov.au/datapacks/.
|
|
|
+- Amoh-Gyimah, R., Saberi, M., Sarvi, M., 2016. Macroscopic modeling of pedestrian and bicycle crashes: A cross-comparison of estimation methods. Accid. Anal. Prev. 93, 147–159. <https://doi.org/10.1016/j.aap.2016.05.001>.
|
|
|
+- Amoh-Gyimah, R., Saberi, M., Sarvi, M., 2017. The effect of variations in spatial units on unobserved heterogeneity in macroscopic crash models. Analytic Methods in Accident Research 13, 28–51. [https://doi.org/10.1016/j.amar.2016.11.001.](https://doi.org/10.1016/j.amar.2016.11.001)
|
|
|
+- [Anselin, L., Griffith, D.A., 1988. Do spatial effects really matter in regression analysis?](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0020) [Papers in Regional Science 65 \(1\), 11](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0020)–34.
|
|
|
+- Ariannezhad, A., Karimpour, A., Wu, Y.-J., 2020. Incorporating mode choices into safety analysis at the macroscopic level. Journal of Transportation Engineering, Part A: Systems 146. <https://doi.org/10.1061/jtepbs.0000337>.
|
|
|
+- Arvin, R., Kamrani, M., Khattak, A.J., 2019. How instantaneous driving behavior contributes to crashes at intersections: Extracting useful information from connected vehicle message data. Accid. Anal. Prev. 127, 118–133. [https://doi.org/10.1016/j.](https://doi.org/10.1016/j.aap.2019.01.014) [aap.2019.01.014](https://doi.org/10.1016/j.aap.2019.01.014).
|
|
|
+- Bao, J., Liu, P., Ukkusuri, S.V., 2019. A spatiotemporal deep learning approach for citywide short-term crash risk prediction with multi-source data. Accid. Anal. Prev. 122, 239–254. <https://doi.org/10.1016/j.aap.2018.10.015>.
|
|
|
+- Bao, J., Liu, P., Yu, H., Xu, C., 2017. Incorporating twitter-based human activity information in spatial analysis of crashes in urban areas. Accid. Anal. Prev. 106,
|
|
|
+- 358–369. <https://doi.org/10.1016/j.aap.2017.06.012>. Besag, J., York, J., Molli´[e, A., 1991. Bayesian image restoration, with two applications in](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0045) [spatial statistics. Ann. Inst. Stat. Math. 43 \(1\), 1](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0045)–20.
|
|
|
+- Boufous, S., de Rome, L., Senserrick, T., Ivers, R., 2012. Risk factors for severe injury in cyclists involved in traffic crashes in victoria, australia. Accid. Anal. Prev. 49, 404–409. <https://doi.org/10.1016/j.aap.2012.03.011>.
|
|
|
+- Briz-Redon, A., Martinez-Ruiz, F., Montes, F., 2019. Investigation of the consequences of the modifiable areal unit problem in macroscopic traffic safety analysis: A case study accounting for scale and zoning. Accid. Anal. Prev. 132, 105276 [https://doi.org/](https://doi.org/10.1016/j.aap.2019.105276) [10.1016/j.aap.2019.105276](https://doi.org/10.1016/j.aap.2019.105276).
|
|
|
+- Cai, Q., Abdel-Aty, M., Sun, Y., Lee, J., Yuan, J., 2019. Applying a deep learning approach for transportation safety planning by using high-resolution transportation and land use data. Transportation Research Part A: Policy and Practice 127, 71–85. [https://doi.org/10.1016/j.tra.2019.07.010.](https://doi.org/10.1016/j.tra.2019.07.010)
|
|
|
+- Cai, Q., Lee, J., Eluru, N., Abdel-Aty, M., 2016. Macro-level pedestrian and bicycle crash analysis: Incorporating spatial spillover effects in dual state count models. Accid.
|
|
|
+- Anal. Prev. 93, 14–22. [https://doi.org/10.1016/j.aap.2016.04.018.](https://doi.org/10.1016/j.aap.2016.04.018) Chen, P., Shen, Q., 2016. Built environment effects on cyclist injury severity in automobile-involved bicycle crashes. Accid. Anal. Prev. 86, 239–246. [https://doi.](https://doi.org/10.1016/j.aap.2015.11.002)
|
|
|
+- [org/10.1016/j.aap.2015.11.002.](https://doi.org/10.1016/j.aap.2015.11.002) [Chow, L.F., Zhao, F., Liu, X., Li, M., Ubaka, I., 2006. Transit ridership model based on](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0080) [geographically weighted regression. Transp. Res. Rec. 1972, 105](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0080)–114.
|
|
|
+- [Congdon, P., 1997. Bayesian models for spatial incidence: A case study of suicideusing](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0085) [the bugs program. Health Place 3 \(4\), 229](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0085)–247.
|
|
|
+- [Congdon, P., 2003. Modelling spatially varying impacts of socioeconomicpredictors on](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0090) [mortality outcomes. J. Geogr. Syst. 5, 161](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0090)–184.
|
|
|
+- da Silva, A.R., Fotheringham, A.S., 2016. The multiple testing issue in geographically weighted regression. Geographical Analysis 48, 233–247. [https://doi.org/10.1111/](https://doi.org/10.1111/gean.12084) [gean.12084.](https://doi.org/10.1111/gean.12084)
|
|
|
+- da Silva, A.R., Rodrigues, T.C.V., 2014. Geographically weighted negative binomial regression—incorporating overdispersion. Statistics and Computing 24 (5), 769–783. <https://doi.org/10.1007/s11222-013-9401-9>.
|
|
|
+- [da Silva, A.R., Rodrigues, T.C.V., 2016. 2016. A sas macro for geographically weighted](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0105) [negative binomial regression, AS Global Forum](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0105).
|
|
|
+- Dumbaugh E. and Li W. (2010). Designing for the safety of pedestrians, cyclists, and motorists in urban environments. Journal of the American Planning Association, 77, 69-88.10.1080/01944363.2011.536101.
|
|
|
+- Eluru, N., Bagheri, M., Miranda-Moreno, L.F., Fu, L., 2012. A latent class modeling approach for identifying vehicle driver injury severity factors at highway-railway crossings. Accid. Anal. Prev. 47, 119–127. [https://doi.org/10.1016/j.](https://doi.org/10.1016/j.aap.2012.01.027) [aap.2012.01.027](https://doi.org/10.1016/j.aap.2012.01.027).
|
|
|
+- [Evans, J.D., 1996. Straightforward statistics for the behavioral sciences. Brooks/Cole](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0120) [Publishing Company, Pacific Grove, California.](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0120)
|
|
|
+- [Fotheringham, A.S., Brunsdon, C., Charlton, M., 2002. Geographically weighted](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0125) [regression: The analysis of spatially varying relationships. John Wiley](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0125) & Sons Inc., [New York.](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0125)
|
|
|
+- N. Fournier E. Christofa M.A. Knodler Jr. A mixed methods investigation of bicycle exposure in crash rates Accident Analysis & Prevention 130 2019 54 61 10.1016/j. aap.2017.02.004.
|
|
|
+
|
|
|
+- <span id="page-14-0"></span>Gomes, M., Cunto, F., da Silva, A.R., 2017. Geographically weighted negative binomial regression applied to zonal level safety performance models. Accid. Anal. Prev. 106, 254–261. <https://doi.org/10.1016/j.aap.2017.06.011>.
|
|
|
+- [Gomes, M.M., Pirdavani, A., Brijs, T., Pitombo, C.S., 2019. Assessing the impacts of](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0140) [enriched information on crash prediction performance. Accid. Anal. Prev.](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0140)
|
|
|
+- Goulard, M., Voltz, M., 1992. Linear coregionalization model: Tools for estimation and choice of cross-variogram matrix. Math. Geol. 24, 269–286. [https://doi.org/](https://doi.org/10.1007/bf00893750) [10.1007/bf00893750](https://doi.org/10.1007/bf00893750).
|
|
|
+- Guo, Y., Osama, A., Sayed, T., 2018. A cross-comparison of different techniques for modeling macro-level cyclist crashes. Accid. Anal. Prev. 113, 38–46. [https://doi.org/](https://doi.org/10.1016/j.aap.2018.01.015) [10.1016/j.aap.2018.01.015](https://doi.org/10.1016/j.aap.2018.01.015).
|
|
|
+- Hadayeghi, A., Shalaby, A.S., Persaud, B.N., 2010. Development of planning level transportation safety tools using geographically weighted poisson regression. Accid. Anal. Prev. 42, 676–688. <https://doi.org/10.1016/j.aap.2009.10.016>.
|
|
|
+- Heydari, S., Fu, L., Miranda-Moreno, L.F., Jopseph, L., 2017. Using a flexible multivariate latent class approach to model correlated outcomes: A joint analysis of pedestrian and cyclist injuries. Analytic Methods in Accident Research 13, 16–27. [https://doi.](https://doi.org/10.1016/j.amar.2016.12.002) [org/10.1016/j.amar.2016.12.002](https://doi.org/10.1016/j.amar.2016.12.002).
|
|
|
+- Hezaveh, A.M., Arvin, R., Cherry, C.R., 2019. A geographically weighted regression to estimate the comprehensive cost of traffic crashes at a zonal level. Accid. Anal. Prev. 131, 15–24. [https://doi.org/10.1016/j.aap.2019.05.028.](https://doi.org/10.1016/j.aap.2019.05.028)
|
|
|
+- Huang, H.A.-A.M.A., Darwiche, A.L., 2010. County-level crash risk analysis in florida: Bayesian spatial modeling. Transp. Res. Rec. 2148 (1), 27–37. [https://doi.org/](https://doi.org/10.3141/2148-04) [10.3141/2148-04.](https://doi.org/10.3141/2148-04)
|
|
|
+- Huang, Y., Wang, X., Patton, D., 2018. Examining spatial relationships between crashes and the built environment: A geographically weighted regression approach.
|
|
|
+- J. Transp. Geogr. 69, 221–233. [https://doi.org/10.1016/j.jtrangeo.2018.04.027.](https://doi.org/10.1016/j.jtrangeo.2018.04.027) [Isaaks, E.H., Srivastava, R.M., 1989. An introduction to applied geostatistics. Oxford](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0180)
|
|
|
+- [University Press, New York](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0180). Jacobsen L. P. (2003). Safety in numbers: More walkers and bicyclists, safer walking and bicycling[j]. Injury prevention, 2003, 9(3):205-209. Injury Prevention, 9(3):205- 209.
|
|
|
+- [Kim, J.-K., Kim, S., Ulfarsson, G.F., Porrello, L.A., 2007. Bicyclist injury severities in](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0190) bicycle–[motor vehicle accidents. Accid. Anal. Prev. 39 \(2\), 238](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0190)–251.
|
|
|
+- [Kim, K., Yamashita, E., 2002. Motor vehicle crashes and land use: Empirical analysis](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0195) [from hawaii. Transp. Res. Rec. 1784, 73](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0195)–79.
|
|
|
+- Ladron ´ [de Guevara, F., Washington, S., Oh, J., 2004. Forecasting crashes at the planning](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0200) [level: Simultaneous negative binomial crash model applied in tucson, arizona.](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0200) [Transp. Res. Rec. 1897, 191](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0200)–199.
|
|
|
+- Lee, J., Abdel-Aty, M., Cai, Q., 2017. Intersection crash prediction modeling with macrolevel data from various geographic units. Accid. Anal. Prev. 102, 213–226. [https://](https://doi.org/10.1016/j.aap.2017.03.009) [doi.org/10.1016/j.aap.2017.03.009.](https://doi.org/10.1016/j.aap.2017.03.009)
|
|
|
+- Lee, J., Abdel-Aty, M., Jiang, X., 2015. Multivariate crash modeling for motor vehicle and non-motorized modes at the macroscopic level. Accid. Anal. Prev. 78, 146–154. [https://doi.org/10.1016/j.aap.2015.03.003.](https://doi.org/10.1016/j.aap.2015.03.003)
|
|
|
+- Li, Z., Chen, X., Ci, Y., Chen, C., Zhang, G., 2019. A hierarchical bayesian spatiotemporal random parameters approach for alcohol/drug impaired-driving crash frequency analysis. Analytic Methods in Accident Research 21, 44–61. [https://doi.org/](https://doi.org/10.1016/j.amar.2019.01.002) [10.1016/j.amar.2019.01.002.](https://doi.org/10.1016/j.amar.2019.01.002)
|
|
|
+- Li, Z., Wang, W., Liu, P., Bigham, J.M., Ragland, D.R., 2013. Using geographically weighted poisson regression for county-level crash modeling in california. Saf. Sci.
|
|
|
+- 58, 89–97. <https://doi.org/10.1016/j.ssci.2013.04.005>. Lu B., Charlton M., Harris P. and Fotheringham A. S. (2014). Geographically weighted regression with a non-euclidean distance metric: A case study using hedonic house price data. International Journal of Geographical Information Science, 28, 660-
|
|
|
+- 681.10.1080/13658816.2013.865739. [Lusk, A.C., Furth, P.G., Morency, P., Miranda-Moreno, L.F., Willett, W.C., et al., 2011.](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0230) [Risk of injury for bicycling on cycle tracks versus in the street. Injury Prevention. 17](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0230) [\(2\), 131](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0230)–135.
|
|
|
+- MacNab, Y.C., 2004. Bayesian spatial and ecological models for small-area accident and injury analysis. Accid. Anal. Prev. 36, 1019–1028. [https://doi.org/10.1016/j.](https://doi.org/10.1016/j.aap.2002.05.001) [aap.2002.05.001](https://doi.org/10.1016/j.aap.2002.05.001).
|
|
|
+- [Marcotte, D., 2012. Revisiting the linear model of coregionalization. Springer,](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0240) [Netherlands.](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0240)
|
|
|
+- [Marks, H., 1957. Subdividing for traffic safety.](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0245) *Traffic Quarterly* 11 (3), 308–325.
|
|
|
+- Melbourne C. o. (2013). Melbourne bicycle account: Cycling census 2013. Merlin, L.A., Cherry, C.R., Mohamadi-Hezaveh, A., Dumbaugh, E., 2020. Residential accessibility's relationships with crash rates per capita. Journal of Transport and Land Use 13, 113–128. <https://doi.org/10.5198/jtlu.2020.1626>.
|
|
|
+- Miaou, S.-P.L.D., 2003. Modeling traffic crash-flow relationships for intersections: Dispersion parameter, functional form, and bayes versus empirical bayes methods. Transp. Res. Rec. 1840 (1), 31–40. <https://doi.org/10.3141/1840-04>.
|
|
|
+- Morrison, C.N., Thompson, J., Kondo, M.C., Beck, B., 2019. On-road bicycle lane types, roadway characteristics, and risks for bicycle crashes. Accid. Anal. Prev. 123, 123–131. <https://doi.org/10.1016/j.aap.2018.11.017>.
|
|
|
+- Munira, S., Sener, I.N., 2020. A geographically weighted regression model to examine the spatial variation of the socioeconomic and land-use factors associated with strava bike activity in austin, texas. J. Transp. Geogr. 88 [https://doi.org/10.1016/j.](https://doi.org/10.1016/j.jtrangeo.2020.102865) [jtrangeo.2020.102865](https://doi.org/10.1016/j.jtrangeo.2020.102865).
|
|
|
+- Nakaya, T., Fotheringham, A.S., Brunsdon, C., Charlton, M., 2005. Geographically weighted poisson regression for disease association mapping. Stat. Med. 24, 2695–2717.<https://doi.org/10.1002/sim.2129>.
|
|
|
+- [Narayanamoorthy, S., Paleti, R., Bhat, C.R., 2013. On accommodating spatial](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0280) [dependence in bicycle and pedestrian injury counts by severity level. Transportation](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0280) [Research Part B: Methodological 55, 245](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0280)–264.
|
|
|
+
|
|
|
+- [Neter, J., Kutner, M.H., Nachtsheim, C.J., Wasserman, W., 1996. Applied linear statistical](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0285) [models. McGraw-Hill Education \(Irwin Chicago\).](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0285)
|
|
|
+- O'brien, R.M., 2007. A caution regarding rules of thumb for variance inflation factors. Qual. Quant. 41, 673–690. <https://doi.org/10.1007/s11135-006-9018-6>.
|
|
|
+- Obelheiro, M.R., da Silva, A.R., Nodari, C.T., Cybis, H.B.B., Lindau, L.A., 2020. A new zone system to analyze the spatial relationships between the built environment and traffic safety. J. Transp. Geogr. 84 [https://doi.org/10.1016/j.jtrangeo.2020.102699.](https://doi.org/10.1016/j.jtrangeo.2020.102699)
|
|
|
+- Osama, A., Sayed, T., 2016. Evaluating the impact of bike network indicators on cyclist safety using macro-level collision prediction models. Accid. Anal. Prev. 97, 28–37. [https://doi.org/10.1016/j.aap.2016.08.010.](https://doi.org/10.1016/j.aap.2016.08.010)
|
|
|
+- OSM (Open Street Map) (2019). Road network. https://www.openstreetmap.org/ #map=15/-37.6859/144.5712,.
|
|
|
+- [Pebesma, E.J., 2004. Multivariable geostatistics in s: The gstat package. Comput. Geosci.](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0310) [30, 683](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0310)–691.
|
|
|
+- B.D. Pelletier P. Larocque G. Fyles JW, Coregionalization analysis with a drift for multiscale assessment of spatial relationships be- tween ecological variables 2. Estimation of correlations and coeffi- cients of determination Environmental and Ecological Statistics 16 4 2008 467 94.
|
|
|
+- Pirdavani, A., Bellemans, T., Brijs, T., Kochan, B., Wets, G., 2014. Assessing the road safety impacts of a teleworking policy by means of geographically weighted regression method. J. Transp. Geogr. 39, 96–110. [https://doi.org/10.1016/j.](https://doi.org/10.1016/j.jtrangeo.2014.06.021) [jtrangeo.2014.06.021.](https://doi.org/10.1016/j.jtrangeo.2014.06.021)
|
|
|
+- Wedagama D. Priyantha M., Bird R. N. and Metcalfe A. V. The influence of urban landuse on non-motorised transport casualties Accident Analysis & Prevention 38 2006 1049 57 10.1016/j.aap.2006.01.006.
|
|
|
+- Qin, X., Ivan, J.N., Ravishanker, N., 2004. Selecting exposure measures in crash rate prediction for two-lane highway segments. Accid. Anal. Prev. 36, 183–191. [https://](https://doi.org/10.1016/s0001-4575(02)00148-3) [doi.org/10.1016/s0001-4575\(02\)00148-3](https://doi.org/10.1016/s0001-4575(02)00148-3).
|
|
|
+- Raihan, M.A., Alluri, P., Wu, W., Gan, A., 2019. Estimation of bicycle crash modification factors (cmfs) on urban facilities using zero inflated negative binomial models. Accid. Anal. Prev. 123, 303–313. <https://doi.org/10.1016/j.aap.2018.12.009>.
|
|
|
+- Ribeiro, M.C., Sousa, A.J., Pereira, M.J., 2016. A coregionalization model can assist specification of geographically weighted poisson regression: Application to an ecological study. Spat Spatiotemporal Epidemiol 17, 1–13. [https://doi.org/10.1016/](https://doi.org/10.1016/j.sste.2016.02.001) [j.sste.2016.02.001](https://doi.org/10.1016/j.sste.2016.02.001).
|
|
|
+- M.P. Ribeiro P. Llop E. Branquinho C. Sousa A.J. Pereira MJ, Multivariate geostatistical methods for analysis of relationships between ecological indicators and environmental factors at multiple spatial scales Ecological Indicators 29 2013 339 47.
|
|
|
+- Saha, D., Alluri, P., Gan, A., Wu, W., 2018. Spatial analysis of macro-level bicycle crashes using the class of conditional autoregressive models. Accid. Anal. Prev. 118, 166–177. <https://doi.org/10.1016/j.aap.2018.02.014>.
|
|
|
+- [Schepers, P., Twisk, D., Fishman, E., Fyhri, A., Jensen, A., 2017. The dutch road to a high](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0355) [level of cycling safety. Saf. Sci. 92, 264](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0355)–273.
|
|
|
+- Siddiqui, C., Abdel-Aty, M., Choi, K., 2012. Macroscopic spatial analysis of pedestrian and bicycle crashes. Accid. Anal. Prev. 45, 382–391. [https://doi.org/10.1016/j.](https://doi.org/10.1016/j.aap.2011.08.003)
|
|
|
+- Smith D. T. and Appleyard D. (1981). Improving the residential street environment. (report fhwa/rd-81/031). Washington DC: Federal Highway Administration, U.S. Department of Transportation.
|
|
|
+- Tang, J., Gao, F., Liu, F., Han, C., Lee, J., 2020. Spatial heterogeneity analysis of macrolevel crashes using geographically weighted poisson quantile regression. Accid. Anal. Prev. 148, 105833<https://doi.org/10.1016/j.aap.2020.105833>.
|
|
|
+- Y.-K. Tu D. Gunnell MS. G. Simpson's paradox, lord's paradox, and suppression effects are the same phenomenon – the reversal paradox Emerg Themes Epidemiol. 2008.
|
|
|
+- VicRoads (2016). 2016 principal bicycle network data. https://vicroadsopendatavicroadsmaps.opendata.arcgis.com/datasets/principal-bicycle-network-v2.
|
|
|
+- VicRoads (2019a). Crashes last five years. https://services2.arcgis.com/ 18ajPSI0b3ppsmMt/arcgis/rest/services/Crashes\_Last\_Five\_Years/FeatureServer/0.
|
|
|
+- VicRoads (2019b). Road width and number of lanes. https://vicroadsopendatavicroadsmaps.opendata.arcgis.com/datasets/24ccad5c745e4addabfcfb32c400ee83\_ 0/geoservice.
|
|
|
+- VicRoads (2019c). Traffic volume. https://vicdata.vicroads.vic.gov.au/server/rest/ services/Operations\_Traffic/FeatureServer/0.
|
|
|
+- [Victoria, T.f., 2018. Victorian cycling strategy 2018](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0400)–28: Increasing cycling for transport. [Transport for Victoria Department of Economic Development, Jobs, Transport and](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0400) [Resources, Melbourne Victoria http://www.economicdevelopment.vic.gov.au](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0400).
|
|
|
+- VRO (Victorian Resources Online) (2017). Victorian land use information system 2016. https://discover.data.vic.gov.au/dataset/victorian-land-cover-mapping-2016. DOI: 10.4226/92/590aba615a44d.
|
|
|
+- [Wackernagel, H., 2003. Multivariate geostatistics,3rd ed. Springer.](http://refhub.elsevier.com/S0001-4575(21)00261-X/h0410)
|
|
|
+- Wang, J., Huang, H., Zeng, Q., 2017. The effect of zonal factors in estimating crash risks by transportation modes: Motor vehicle, bicycle and pedestrian. Accid. Anal. Prev. 98, 223–231.<https://doi.org/10.1016/j.aap.2016.10.018>.
|
|
|
+- Wei, F., Lovegrove, G., 2013. An empirical tool to evaluate the safety of cyclists: Community based, macro-level collision prediction models using negative binomial regression. Accid. Anal. Prev. 61, 129–137. [https://doi.org/10.1016/j.](https://doi.org/10.1016/j.aap.2012.05.018) [aap.2012.05.018](https://doi.org/10.1016/j.aap.2012.05.018).
|
|
|
+- WHO (World Health Organization) (2018). Global status report on road safety 2018 summary.Https://www.Who.Int/publications/i/item/9789241565684.
|
|
|
+- Xu, P., Huang, H., 2015. Modeling crash spatial heterogeneity: Random parameter versus geographically weighting. Accid. Anal. Prev. 75, 16–25. [https://doi.org/10.1016/j.](https://doi.org/10.1016/j.aap.2014.10.020)
|
|
|
+- Xu, P., Huang, H., Dong, N., Wong, S.C., 2017. Revisiting crash spatial heterogeneity: A bayesian spatially varying coefficients approach. Accid. Anal. Prev. 98, 330–337. [https://doi.org/10.1016/j.aap.2016.10.015.](https://doi.org/10.1016/j.aap.2016.10.015)
|
|
|
+
|
|
|
+- <span id="page-15-0"></span>Yao, S., Loo, B.P., Lam, W.W., 2015. Measures of activity-based pedestrian exposure to the risk of vehicle-pedestrian collisions: Space-time path vs. Potential path tree methods. Accid. Anal. Prev. 75, 320–332. [https://doi.org/10.1016/j.](https://doi.org/10.1016/j.aap.2014.12.005) [aap.2014.12.005](https://doi.org/10.1016/j.aap.2014.12.005).
|
|
|
+- Yasmin, S., Eluru, N., 2016. Latent segmentation based count models: Analysis of bicycle safety in montreal and toronto. Accid. Anal. Prev. 95, 157–171. [https://doi.org/](https://doi.org/10.1016/j.aap.2016.07.015) [10.1016/j.aap.2016.07.015](https://doi.org/10.1016/j.aap.2016.07.015).
|
|
|
+- J.F. Zawadzki n., P Badura, H, Estimation of methane content in coal mines using supplementary physical measurements and multivariable geostatistics International Journal of Coal Geology 118 2013 33 44.
|
|
|
+- Zhai, L., Li, S., Zou, B., Sang, H., Fang, X., et al., 2018a. An improved geographically weighted regression model for pm 2.5 concentration estimation in large areas. Atmos. Environ. 181, 145–154. [https://doi.org/10.1016/j.atmosenv.2018.03.017.](https://doi.org/10.1016/j.atmosenv.2018.03.017)
|
|
|
+- Zhai, X., Huang, H., Gao, M., Dong, N., Sze, N.N., 2018b. Boundary crash data assignment in zonal safety analysis: An iterative approach based on data augmentation and bayesian spatial model. Accid. Anal. Prev. 121, 231–237. [https://](https://doi.org/10.1016/j.aap.2018.09.010) [doi.org/10.1016/j.aap.2018.09.010.](https://doi.org/10.1016/j.aap.2018.09.010)
|
|
|
+- Zhang, Y., Bigham, J., Ragland, D., Chen, X., 2015. Investigating the associations between road network structure and non-motorist accidents. J. Transp. Geogr. 42, 34–47. <https://doi.org/10.1016/j.jtrangeo.2014.10.010>.
|
|
|
+- Ziakopoulos, A., Yannis, G., 2020. A review of spatial approaches in road safety. Accid. Anal. Prev. 135, 105323<https://doi.org/10.1016/j.aap.2019.105323>.
|