Local Spatial Analysis of the Crash Frequency of Food Delivery Motorcyclists vs. Nondelivery Motorcyclists in relation to Points of Interest

Te COVID-19 pandemic has increased the demand for online food delivery services (OFDS), leading to an increase in related crashes over the last few years. While recent studies have focused on nonmotorised vehicles (such as bicycles or e-bikes), few researchers have examined the role of motorcycles and the possible spatial relationships with various points of interest (POIs). In addition, most crash and POIs studies have utilized typical restaurant datasets instead of specifc restaurants partnered with OFDS, which might bias the impact of trafc safety estimation. To address these gaps, a geographically weighted negative binomial model (GWNBR) was used to determine the factors contributing to OFDS-related motorcycle accidents and account for spatial heterogeneity. Te results indicated that areas with more restaurants, intersections, and shopping malls (only signifcant on weekends) tended to have more OFDS motorcycle crashes. Te results should inspire more efective policies for delivery drivers, given the increasing popularity of OFDS.


Introduction
Anyone can order food from a variety of restaurants by using simple online apps.Te order is picked up from the restaurant by a courier and delivered to nearby customers.Tis industry has grown exponentially worldwide [1,2].For instance, in Australia, nearly four million people used online food delivery services (OFDS) in 2020, a number that has grown by 81.1% since 2015 [3].Mexico also experienced an upswing of more than 70% in visits to OFDS websites [4].It must be noted that the popularity of this service was further increased by the COVID-19 quarantines, during which many restaurants partnered with these companies to stay in business [5].
Te surging demand has led to intense competition among delivery service companies that are vying to ofer speedier and more cost-efective options [6].Motorcycles are the favored vehicle due to their capacity to circumnavigate through congested roads, which is crucial for maintaining fast delivery times.Tis ability is particularly relevant in Asia, where dense road networks and narrow alleys are common.Hence, motorcycles, known for their wide availability, popularity, and low cost, emerged as the primary choice of delivery companies.Notably, in Taiwan alone, these companies hired approximately 88,000 OFDS workers in 2020, which underscores this sector's substantial impact on employment rates [7].However, the proliferation of OFDS startups has brought about a noticeable increase in motorcycle-related crashes in recent times.A study conducted by Chen [8] clarifed this pattern, illustrating that the proportion of motorcycle riders deemed at fault rose from 31.0% in 2017 to 61.1% in 2020 in the context of all trafc accidents involving OFDS motorcycle riders in Taiwan.Strikingly, the numbers of Korean hospital patients revealed that the percentage of severe crash injuries among these individuals was much higher compared to nondelivery motorcyclists from 2014 to 2018 [9].A similar pattern manifested in South Korea, where the fatality rate attributed to delivery-related crashes constituted roughly 12% of all fatal trafc accidents, a stark contrast to the 5% attributable to nondelivery crashes [1].Tis concerning trend may be rooted in the demanding environment of delivery drivers, characterized by substantial workloads and stringent time constraints for each order.In their pursuit of maximizing earnings and averting customer complaints, delivery personnel resort to reckless driving, including speeding, weaving through trafc, texting while driving, and ignoring trafc signals [10,11], practices which difer signifcantly from nondelivery motorcyclists who tend to prioritize safety over speed.
Over the past several years, researchers have compared the driving behaviors of nondelivery and delivery motorcyclists as well as other crash-related factors, particularly those linked to trafc violations and driver characteristics [1,12].For example, Wang et al. [13] conducted an extensive observational study involving 600 delivery drivers, in which 480 participants were interviewed to assess their driving behaviors and identify factors that may contribute to injury crashes in China.Te fndings revealed that 91.3% of the respondents frequently exceeded speed limits, ran red lights, and disregarded trafc signals, which led to 76.5% of the interviewees being involved in at least one trafc accident.Li et al. [14] conducted a similar analysis which revealed that delivery-related vehicles in Shanghai, China, were 5.4% more likely to be involved in red-light running incidents compared to personal vehicles.
Surprisingly, there are few motorcycle trafc studies in which the built environmental factors are considered (e.g., points of interest).According to Qin et al. [15], OFDS motorcycle riders may resort to riskier driving behaviors and travel routes compared to nondelivery motorcyclists who are typically commuters.Furthermore, the partners of the OFDS are no longer limited to well-known restaurants; they also include local eateries, convenience stores, and supermarkets that deliver to hotels, schools, and residential areas.Tis expansion may raise the crash risk for OFDS drivers [16,17].In other words, these types of points of interest (POIs) may increase the use of OFDS, contributing to higher collision rates.
Te present study focuses on specifying and comparing the POIs associated with OFDS and nondelivery motorcycle crashes.Specifcally, OFDS-related crashes, defned as motorcycle crashes in this context, involve drivers employed by OFDS platforms to deliver food to customers.In contrast, nondelivery motorcycle crashes encompass accidents where motorcycles are not involved in food delivery activities.In order to establish a more robust relationship with the exposure specifc to OFDS, various risk factors were considered including partnerships with restaurants afliated with OFDS platforms (Uber Eats and Food Panda), trafc conditions, road geometry, and population density.To enhance the methodological approach, this study used the geographically weighted negative binomial regression (GWNBR) model, which is unique and ideal for the analysis of OFDS motorcycle crashes, as it efectively addresses the spatial heterogeneity and the overdispersion issues commonly observed in crash data [18,19].Tis approach serves to provide a local-level analysis of the factors infuencing crash patterns in the context of OFDS motorcycles and nondelivery motorcycles.

Factors Associated with OFDS Motorcycle Crashes.
During peak hours, online food orders, particularly takeout, prompt concentrated order aggregations in densely populated zones such as residential, commercial, and urban centers [20].Tis scenario necessitates that OFDS, often on motorcycles, navigate through the bustling streets and intersections of Asian cities, where they share the road with a diverse mix of vehicles, pedestrians, and cyclists [13].Consequently, a notable correlation was frequently observed between vehicle density and OFDS motorcycle-related crashes [21].Te increased width of these intersections is often accompanied by longer trafc signal intervals ranging from 60 to 120 seconds [15].In addition, main roads are frequently divided by physical barriers such as fences and green belts, forcing trafc to move to the nearest front or rear intersection to execute U-turns [22].Given this chaotic environment, instances of trafc violations among delivery couriers are prevalent, with intersections serving as hotspots for such behaviors [15,23].
Te popularity of these services continues to grow.Han et al. [24] and Li and Zhang [25] showed that in China, students and white-collar workers use OFDS because they are fast, convenient, and do not require leaving the workplace or school.Similarly, 15 Italian cities saw an increase of 137% in the use of OFDS by ofce employees at lunchtime [26].In other words, ofce and school areas become the main destination for OFDS which might increase the crash likelihood.Tis platform is also frequently used by tourists who have meals delivered to their hotels [27] in order to have more time for sightseeing.Most are unfamiliar with the area and others wish to avoid long lines at famous restaurants.Tus, in this study, universities, schools, ofces, and hotels were treated as POIs (destinations with a concentration of customers) that are associated with a high number of OFDS motorcycle crashes.

Factors Associated with Nondelivery Motorcycle Crashes.
Several studies of nondelivery motorcycle crashes have shown that POIs, such as commercial centers, schools, and public transportation stops, attract many pedestrians, which subjects these areas to heavy trafc, conficts, and collisions [28][29][30][31][32]. Ivan et al. [29] found that primary schools are a hotspot for pedestrian-motorcycle crashes due to unsafe crossing conditions.Increased pedestrian activity also occurs near commercial areas, shopping centers, bus stops, and subway stations.Tus, there is an increase in trafc accidents at these locations [33,34].
Chen et al. [28] found that few motorcycle crashes occur in touristy areas with many hotels because visitors are more likely to travel by tour bus and public transit than by 2 Journal of Advanced Transportation motorcycle, so trafc patterns in these areas are relatively simple.Similarly, open spaces, such as parks and green spaces, are not as populated as other POIs, so there is less trafc and fewer collisions [35].Areas with many POIs, pedestrians, and heavy trafc tend to experience high rates of motorcycle accidents [30,36,37].
2.3.Crash Prediction Models.Several models have been used in existing crash studies to defne the spatial dependency and heterogeneity of crash frequency and predictor factors [38].
In order to address the spatial heterogeneity of the crash data, researchers such as Flask et al. [39], Jia et al. [30], and Truong et al. [40] mainly focused on spatial global regression; however, the parameter estimates remained constant across the study area, so it was impossible to accurately determine the relationship between environmental conditions and crash frequency.For example, Cheng et al. [41] successfully applied spatial modeling to predict crash frequency in trafc analysis zones (TAZs) in California; however, their model only accounted for spatial dependence and disregarded spatial heterogeneity.GWR-based models were utilized in several later studies to allow the coefcients to change throughout the study area in order to more accurately predict crash frequency [42].For example, Hezaveh et al. [43] and Bao et al. [44] used a geographically weighted Poisson regression (GWPR) at various spatial units, which performs better than the traditional global Poisson model because it accounts for spatial heterogeneity.However, it must be noted that GWPR has a limited ability to predict high variability in crash data, as typically there are a higher number of crashes in urban than in rural areas, so GWPR might be not suitable for use in large study areas (whether urban or rural) or study units that include a variety of land uses.When DaSilva and Rodrigues [45], Amoh-Gyimah et al. [46], Gomes et al. [18], Obelheiro et al. [47], and Mathew et al. [48] used a geographically weighted negative binomial regression model (GWNBR) for crash frequency prediction; they found it to be more accurate than the GWPR and traditional global modeling, possibly because it addresses spatial heterogeneity and the overdispersion problem simultaneously, thus avoiding the underestimation of standard errors and misleading inferences regarding coefcients.Yet, despite its potential use in practical implications for local OFDS policy regulation and road safety improvement, it has not been utilized to assess the association between OFDS motorcycle crashes and POIs.
Te present study, therefore, specifed the POIs that may be associated with OFDS motorcycle crashes, through the novel approach of examining the restaurants that partnered with OFDS (such as Uber Eats and Food Panda).Furthermore, in order to explain this relationship more accurately and provide new evidence that OFDS workers are in more danger than nondelivery motorcyclists, several risk factors were considered, such as trafc conditions, road geometry, and population.Unlike any previously published study, the GWNBR model was utilized in this study because, as previously stated, it accounts for the spatial heterogeneity and overdispersion found within the crash data [19].

Data and Methodology
3.1.Data.Tis study utilized the 2020 motorcycle crash reports from the Taipei Department of Transportation.Taipei City was selected as the study site because after 2019, the crash dataset was separated into OFDS-related accidents (platform partner restaurants are on-demand food applications) and nondelivery accidents.Te study area encompasses both the urban region in the central-western part of the city and the suburbs, located on the periphery of the downtown area, including the districts of Beitou, Shilin, Neihu, Wenshan, and Nangang (refer to Figure 1 for an illustration of the study area).
Te crash dataset is accessible to the public through the government website (https://data.gov.tw), which provides free and open access to the data.Te statistics show that out of the 64,795 total number of motorcycle crashes in 2020, 2,314 involved OFDS drivers.Te crash dataset for this study includes location data (latitude and longitude) and other nonspatial information, such as crash ID codes, causal factors, collision type, severity, number of fatalities and injuries, road surface, weather, lighting conditions, crash date, vehicle type maneuvers, and directions.Other unrelated variables are not included due to space constraints.
Four explanatory variables are used: POIs, trafc conditions (number of bus stops and vehicle density in each village), road geometry, and population in Taipei City.Te general POIs data (excluding restaurants) were downloaded from OpenStreetMap, including the number of supermarkets, shopping malls, schools, universities, hotels, and hospitals within the study area.
Unlike previous literature, this study utilized restaurants that partnered with OFDS such as Uber Eats and Food Panda.Tese OFDS are commonly used by citizens and have many partner restaurants.For restaurant data (limited to those partnered with OFDS), the Python "Selenium" library and "BeautifulSoup" library were used to collect information on the restaurants partnered with OFDS in the study area.
Because the current open datasets are not limited to online options, the Selenium web driver was used to mimic how the user would search for this information, while BeautifulSoup was used to flter the webpage content and scrape the data into csv fles.Tese tools can open the browser automatically and access OFDS websites in just a few steps, as shown in Figure 2. Te procedure begins with accessing OFDS websites (Food Panda and Uber Eats in this case).Selenium was used to simulate how a customer would enter a pseudogeolocation (based on WGS84) into the search bar and how this program would scan the websites for restaurants closely associated with OFDS.Te URLs of nearby restaurants were recorded based on query outcomes within the pseudolocation.Tis process was iterated across various locations until it covered the entire study area.Next, pertinent elements such as each restaurant's name, address, and coordinates were collected from the URLs and individually extracted using the Python "BeautifulSoup" library.Tis study also utilized this library to create a restaurant database in the csv format.Subsequently, duplicate entries (sharing identical names, coordinates, addresses, and foors) were removed.For more in-depth information on Python-based web browsing with Selenium, please see references from Raghavendra [50] and Yudin [51].
Te other explanatory variable is trafc fow, which is automatically generated by using trafc fow data from vehicle detectors on main streets in Taipei City.Road geometry data were collected using ArcGIS to calculate the key characteristics for the road network, including the number of intersections and total road length for each village.Tis study employed 456 villages as the study unit (Figure 1), as they constitute the smallest administrative entity and align with the practical implementation scale for local governments.

Negative Binomial (NB) Model.
As described in previous studies, crash frequency is assumed to be Poissondistributed where the mean is assumed to be gammadistributed [52,53].Tis gives rise to the Poisson-gamma or negative binomial (NB) distribution.Te NB (gamma Poison) distribution is particularly well-suited for count data when there is evidence of overdispersion, which occurs when the variance of the data is higher than what would be expected under a binomial distribution [52].In this case, the occurrence of motorcycle crashes may exhibit varying levels of dispersion due to factors such as a variety of road conditions, trafc density, and weather conditions.In addition, the NB model was used in this study because of the general relationship between each prediction variable and crash frequency [19,54,55] and shows the efects of not accounting for spatial heterogeneity.Te relationship is shown as follows: where NB stands for the negative binomial distribution, y i is the number of motorcycle crashes in the i th (i � 1, . . ., n) spatial unit, x ik represents the k th explanatory variable for a spatial unit i, β k (k � 0, 1, . . ., p) are the coefcients for variable k, which is the ofset variable, and α is the overdispersion parameter.

Geographically Weighted Negative Binomial Model.
A geographically weighted-based model is used to determine a nonstationary relationship between crash frequency and related factors.Tis model allows its parameters to vary over space, in order to determine the spatial heterogeneity between the dependent variable (motorcycle crash frequency) and independent variables.Most trafc crash prediction studies show that the GWR-based model performs better than other conventional nonspatial methods [56].GWR is derived from the frst law of geography, which states that occurrences near location i have a greater infuence on the estimation of β k (u i , v i ) than those farther away.A kernel function, such as the Bisquare, represents the magnitude of this infuence which is illustrated as follows: where d ij is the distance between spatial units i and j and b i(k) represents the adaptive bandwidth.GWNBR is the most commonly used model for crash frequency prediction.
In an early study, Nakaya et al. [57] developed a GWPR to construct a Poisson regression model with various explanatory variables.Te Poisson model cannot be used for analyzing crash data since the mean is assumed to be the same for similar characteristics, other than the spatial components, which is not true [19].Later, Xu and Huang [58] used the GWNBR, which accounts for overdispersed count data.Tis model can address the problems of spatial heterogeneity and overdispersion, as shown in the following equation: where β k (u i , v i ) is the coefcient for the explanatory variable x k in location (u i , v i ), k � 1, . .., n, y i is the number of crashes in the i th spatial unit, and α is the overdispersion parameter.
Terefore, this study assessed the performance of the GWNBR, which was used as a local crash prediction model; by comparing the results (in terms of AIC, log-likelihood, and R-square), it yielded to those of commonly used models such as NB and GWR.

Descriptive Statistics.
Te descriptive statistics and expected association for these candidate variables are shown in Table 1 based on an analysis of the literature [30,59,60].Te results of web scraping revealed that 13,318 restaurants partnered with online OFDS.Supermarkets and shopping malls totaled 233 and 54, respectively.Te numbers of notable customer-centric locations such as schools, universities, hotels, and hospitals were 276, 26, 20, and 32, respectively.Furthermore, the expected association is that all POIs would have a positive relationship with trafc crashes [15][16][17].For example, supermarkets, universities, and hospitals would have a positive association with trafc crashes, meaning that the areas with a high density of these POIs would tend to have a higher crash frequency.However, the association of schools and restaurants with trafc crashes is inconsistent across studies [59,61], which requires further exploration.
Tis study also uses visual results to confrm this spatial relationship among related factors and crashes.To establish an initial correlation between crash frequency and the number of POIs, the present study used the kernel density estimation method to better understand and visualize the relationship between crashes and POI hotspots.A predefned search bandwidth and cell size (100 meters and 20 meters) were selected based on Silverman's rule of thumb in order to better separate the hotspots (smaller bandwidth) and the smoothness of the fgure (smaller cell size) [62,63].Te red circle in Figure 3(a) shows that crash hotspots for OFDS personnel are concentrated in the central business district (CBD) of Taipei City, which is an area with many restaurants, supermarkets, shopping malls, and hospitals (Figure 4).Tis area also has a high population density and a large commercial zone (purple areas in Figure 5(b)).Te nondelivery crash hotspots (black circles in Figure 3(b)) are dispersed throughout the central, western, and southern regions of Taipei City in the commercial and school zones.
Figure 6 illustrates the temporal trends in OFDS and nondelivery motorcycle crashes.Both types of crashes exhibit bimodal distributions with two peak crash frequencies.However, the peak crash frequency for OFDS motorcycles is observed during lunch and dinner hours, specifcally at 12 PM and 6 PM.In contrast, nondelivery motorcycle crashes Filter the page content (Name, address, and coordinate) List of partner restaurants in csv format.

Data cleaning (Remove duplicate data)
Figure 2: Web scrapping (web scraping is a data collection technique that involves extracting information from websites [49]) partnered restaurant with OFDS (Selenium is a browser automation tool in Python which is used to control a web browser and interact with websites much as a human user would).

Model Comparison.
A nonspatial global model (NB) and spatial local model (GWNBR) were developed to determine the efect of overdispersion on crash analyses and spatial heterogeneity, using the methodology that is described.Tese NB-based models were used because the dependent variable, which is the frequency of crashes, has a skewed distribution.Te results for the coefcient are presented in Tables 2-5.Te nonspatial model (NB) only estimates a single coefcient, so the parameters of the NB model are constant over the entire study area.However, the coefcients for the spatial models (GWNBR) are spatially varied, so they are described using the minimum, median, and maximum values for the coefcients.Te variables andcoefcients highlighted in bold within the GWNBR model indicate that their coefcient values demonstrate statistical signifcance at a 5% level in over 50% of the total study units.Please refer to the percentage of signifcant coefcients column for the specifc proportion of study units withsignifcance levels of 5% for each variable.Four criteria were used to determine the performance of the models; recall that models were estimated for four dependent variables, as described in Table 1.In order to show the model complexity and model ftting, the Akaike information criterion (AIC) is used to assess the model performance.Te model with the lowest AIC value has a better ft with the observed crash data.Te model with the highest R-square and log-likelihood values and the lowest RMSE value also indicates a better model ft with higher likelihood and lower prediction errors.Te model-ftting results for OFDS-related crashes and nondelivery crashes are shown in Tables 2-5.Te GWNBR model performed better than the NB model, as expected.Te results show that the spatial models outperform the nonspatial model because the local model incorporates more spatial heterogeneity than the NB, for which coefcients can vary spatially across observations.
Tere is no residual clustering within the GWNBR models because its p value for Moran's I is higher than 0.05.Otherwise, for the NB model residuals, the Moran's I value is positive with a p value lower than 0.05.In other words, the residual clustering becomes insignifcant when the local model and overdispersion in the data are incorporated.
Terefore, the spatial autocorrelation between the model residuals (spatial heterogeneity) and the overdispersion of crash frequency is explained by the GWNBR model.
Furthermore, the variance infation factor (VIF) was used as a metric to evaluate multicollinearity.VIF gauges the degree to which multicollinearity contributes to an increase in the variance of estimated regression coefcients.A VIF value exceeding 1 indicates the presence of multicollinearity, with higher values denoting more pronounced multicollinearity [64,65].Typically, VIF values surpassing 5 or 10 are deemed indicative of substantial multicollinearity, as suggested by Akinwande et al. [66] and Tsagris and Pandis [67], although specifc thresholds may vary depending on the specifc context.However, it is noteworthy that all independent variables in the model yielded VIF values below 5, signifying the absence of multicollinearity issues in the crash risk model.

Discussion
Te results of this study show that the GWNBR model performs better because it can account for spatial heterogeneity.Te results show that areas with more restaurants, major intersections, and shopping malls (only on weekends) tend to experience more crashes involving OFDS drivers.Te nondelivery crash model uses the same predictor variables but schools are added.Each POI has a diferent spatial coefcient distribution in relation to OFDS-related crashes on weekends and weekdays.Te distributions for local coefcient estimations and their signifcance are plotted in Figures 7 and 8. Areas with a signifcance level of more than 50% were included in the parameter analysis [47].

Te Association of POI, Trafc Conditions, and Road
Geometry.Te results reveal a signifcant and positive correlation between the presence of restaurants as POIs and  Journal of Advanced Transportation both types of motorcycle crashes in 69.30%-90.57% of the study area, as illustrated in Tables 2-5.Tis fnding implies that areas with a higher density of restaurants tend to experience increased crash rates for both OFDS and nondelivery motorcycles.Tis result is consistent with that of Lin et al. [17], which shows that more crashes involved OFDS drivers (mean coefcient value of 0.009 on weekdays and 0.008 on weekends) in with a higher number of  0.181 * (0.0001) 0.110 (0.213) * Signifcance levels of 5%.Te variables and coefcients highlighted in bold within the GWNBR model indicate that their coefcient values demonstrate statistical signifcance at a 5% level in over 50% of the total study units.Please refer to the percentage of signifcant coefcients column for the specifc proportion of study units with signifcance levels of 5% for each variable.restaurants than nondelivery crashes (mean coefcient value of 0.007 on weekdays and 0.005 on weekends).Tis outcome is expected, given the simplicity of the majority of OFDS routes from the restaurant to the customer, while routes for nondelivery motorcyclists vary widely.Regarding nondelivery crashes, studies underscore that patrons often park near restaurants, disrupting trafc fow and increasing the likelihood of accidents [68][69][70].Furthermore, the results of this study show that commercial districts (shopping malls) have a signifcant positive relationship with crashes involving OFDS drivers in 67.54% of the study area (Table 3).Tis result aligns with previous studies [71,72], highlighting the fnding that food courts in shopping malls which are associated with OFDS contribute to an increased probability of crashes in the vicinity of these commercial centers.Such zones, which are characterized by high pedestrian and vehicle density, often witness frequent trafc accidents [73,74].
Te other POI in this study is schools, which show a signifcant positive association with nondelivery crashes in 62.94% of the study area (Table 4).Tis result is consistent with that of prior studies, such as the study by Ivan et al. [29], which identifed several vehicle crash hotspots near schools.Most motorcycle-pedestrian crashes occur when the driver does not give pedestrians enough time to cross the street, so intersections are prime areas for trafc accidents [75].Interestingly, this study reveals that schools and universities are not signifcantly correlated with OFDS motorcycle crashes, a fnding that diverges from a study by Li et al. [26], which stated that OFDS are commonly used by high school and university students.One possible reason is that the majority of schools in Taiwan ofer students several canteens and lunch box options, which reduces the likelihood that they will order food from an outside source, which also may be too expensive.
In terms of trafc conditions, the number of intersections is signifcantly and positively associated with an increase in both types of motorcycle crashes in 69.96%-91.01% of the study area (Tables 2-5).Tis result is consistent with that of earlier studies citing high incidence rates of drivers running red lights at intersections [1,12].However, the coefcient value shows that this behavior has more of an efect on OFDS drivers than on other motorcyclists (0.438 and 0.314, respectively).Li et al. [26] discovered that these two groups have diferent reasons and temporal patterns for running red lights, indicating that OFDS couriers are more likely to violate trafc signals during lunch and dinner hours due to time pressures, whereas other motorcyclists are more likely to commit trafc violations during the morning rush hour (8-9 AM).Tis temporal distinction results in an increased occurrence of right-angle collisions during the morning rush hour, leading to more severe injuries and fatalities, particularly for those involved in collisions with other vehicles.

Spatial Distribution Trend for the Estimated Parameter.
Te GWNBR analysis revealed a notable spatial distribution trend for the estimated parameters, particularly concerning the restaurant variable.Te coefcient values exhibited an increasing trend from urban to rural areas, spanning from west-central to the northeast, for both OFDS and nondelivery motorcycle crashes (refer to Figures 7(a), 7(b), 8(a), and 8(b)).Tis pattern aligns with expectations, as rural areas with national parks, such as northern Taipei, tend to have fewer restaurants.Tus, the addition of one restaurant in such an area increases the demand for OFDS, consequently elevating the likelihood of accidents involving delivery drivers.Tis observation is consistent with fndings from Vahedi and Efati [74] that commercial POIs, such as restaurants, are the key factors in crashes in rural areas.
However, it is crucial to note that each POI exhibited diferent spatial coefcient distributions with regard to OFDS-related activities depending on whether it is a weekday or a weekend (represented by coefcients "a" and "b").On weekdays, the coefcient values were higher in the northern, western, and south-western regions of Taipei (indicated by the blue circles B and C in Figure 7(a)), whereas in the east, they were lower (illustrated by the green circle A in Figure 7(b)) compared to weekends.
Te higher coefcients on weekdays can be attributed to a higher number of deliveries to commercial districts, schools or universities, ofces, and industrial complexes situated in the central-western and eastern parts of Taipei City.Tis aligns with the fndings of Li et al. [26], indicating that OFDS are particularly popular among white-collar workers.Tis fnding explains the higher coefcient values for restaurants near these areas during the week.On weekends, the delivery focus shifts to residential areas mainly in the northern and southern regions of Taipei (represented by the green circle).However, a reduction in the coefcient distribution of nondelivery crashes occurred in most areas, except in western Taipei (indicated by the blue circle D in Figure 8(b)).
Figures 7(c), 7(d), 8(c), and 8(d) illustrate that the coefcient values for the intersection variable in urban areas were higher than in suburban locals (central to southern) for both nondelivery and OFDS crashes.While the trend for coefcient distribution remained similar, the overall value decreased on weekends, potentially due to reduced demand for OFDS and less trafc.Consequently, the probability of motorcycle crashes decreased throughout Taipei City.It is essential to note that the magnitude of the coefcients for two diferent models (weekdays and weekends) cannot be used to determine the efect of the same variable because the datasets are diferent and weekends are much shorter than weekdays.
Regarding other POIs, shopping malls were found to be signifcantly associated with OFDS crashes only on weekends (refer to Figure 7(e)).Shopping malls have a similar coefcient distribution to other variables in that the coefcient value increases from urban to suburban and rural areas (central to north and south).A high number of shopping malls in rural areas heighten the probability of OFDS crashes as malls contribute to increased trafc, particularly on weekends [48,76].
Furthermore, the presence of schools had a signifcant positive association with nondelivery crashes on weekdays, with the coefcient value increasing from urban to rural areas of Taipei (see Figure 8(e)).Tis result is expected because, on weekdays, trafc density increases signifcantly when students are walking to school and parents are dropping their kids of or picking them up.In addition, during school hours, there is less OFDS motorcycle trafc because most people are working.

Conclusions
Tis study has documented the development of several statistical models to defne the association between built environmental (e.g., POIs) factors and OFDS motorcycle crashes.In addition, this study also investigated the difference in the signifcant variable related to OFDS and Journal of Advanced Transportation nondelivery motorcycle crashes.Te traditional NB model and GWR local models (GWNBR) were estimated.Te GWNBR outperforms the traditional model for all dependent variables because it captured spatial heterogeneity.Tus, the results suggest the use of the local model for better estimating crash frequency.Despite its better performance, the optimum bandwidth for the GWNBR model is large (84% of the study area), which indicated that a small spatial dependence between the frequency of OFDS and nondelivery crashes exists in the study area.
Te local model results are mostly nonstationary, but some built environmental variables have a more homogenous relationship with OFDS and nondelivery crashes.Te results of the study by Obelheiro et al. [47] show that a built environmental variable can signifcantly afect crashes if it has a homogenous direction (either positive or negative) and the signifcant coefcient is high (at least 50%).Te results show that the restaurant and intersection variables have a more stable and strong relationship with OFDS and nondelivery crashes during the weekends and weekdays.However, restaurants have a greater efect on food deliveryrelated crashes because they have a higher coefcient value.
In terms of other POIs, shopping malls have a signifcant efect on OFDS-related crashes only on weekends.Te school variable only has a signifcant efect on nondelivery crashes on weekdays.
In terms of policy implications, the results highlight the need for efective strategies to reduce OFDS-related crashes, with a specifc focus on areas with a high density of restaurants, schools, shopping malls, and numerous intersections.For instance, restaurants and local authorities could collaborate to establish dedicated drive-through or pick-up lanes exclusively for OFDS drivers.Tis could help streamline the process, reduce time spent in congested areas, and minimize interaction between OFDS motorcycles and pedestrians or other road users.In addition, for particularly popular drop-of locations, such as shopping malls and schools, deploying additional police ofcers to direct trafc during peak hours in city centers (areas with higher coefcients) could prove efective for reducing trafc violations and the overall number of crashes involving delivery drivers [77].Te higher crash rates at intersections are a major area of concern.To address this issue, retroftting these intersections with variable message signs (VMS) could infuence delivery drivers to choose routes [78] that would help them avoid trafc congestion and reduce delays.Tis would lower the risk associated with weaving through heavy trafc and enhance the safety of all road users.Lastly, collaborative eforts between local authorities and OFDS platforms to launch public awareness campaigns are suggested to educate both delivery drivers and pedestrians about safe driving and the importance of obeying trafc laws.
To the best of our knowledge, this study is one of the frst to defne the relationship between POI and OFDS motorcycle crashes that utilized the GWNBR model to analyze their association at a zonal level.According to our fndings, using the GWNBR model is an approach that could yield valuable results, which contribute to the expanding literature on the impacts of built environmental factors on trafc safety.Another contribution of this study shows that restaurants with OFDS have distinct relationships between POIs and OFDS motorcycle crashes.
However, this study contains some limitations.First, this study is not able to consider the size or area of the POIs due to the dearth of comprehensive data.Tis is especially notable for restaurants, hotels, and supermarkets.Furthermore, there is a lack of information about small restaurants, which are not involved with OFDS.Future studies should consider this information to better illustrate the diference between the impact of larger and smaller POIs in the regions under study.Second, future studies might use the number of students as a prediction variable which could better describe the impact on OFDS motorcycle crashes.Lastly, this study recommends that future investigations encompass extended periods of time to enable temporal analysis of OFDS motorcycle crashes, particularly during mealtimes [20,79].In addition, consideration of the operational hours of each POI category could enhance the temporal diversity of the data, consequently improving the model's performance when assessing temporal patterns.As a possible avenue for improvement, adopting alternative models such as geographical and temporal weighted regression (GTWR) might be benefcial.If spatial and temporal heterogeneity were included, this model could provide a more comprehensive framework for analyzing the complex dynamics of OFDS crash data.

Figure 1 :
Figure 1: Illustration of the study area.

Table 1 :
Summary of statistics for variables for each village.

Figure 5 :
Figure 5: (a) Population density distribution and (b) land use in Taipei City.

Figure 8 :
Figure 8: Te spatial distribution for parameter estimation for nondelivery motorcycle crashes.(a) Restaurants on weekdays.(b) Restaurants on weekends.(c) Intersections on weekdays.(d) Intersections on weekends.(e) Schools on weekdays.

Table 2 :
Comparison of the results for OFDS-related crashes on weekdays for diferent models.Signifcance levels of 5%.Te variables and coefcients highlighted in bold within the GWNBR model indicate that their coefcient values demonstrate statistical signifcance at a 5% level in over 50% of the total study units.Please refer to the percentage of signifcant coefcients column for the specifc proportion of study units with signifcance levels of 5% for each variable. *

Table 3 :
Comparison of results for OFDS-related crashes on weekends for diferent models.

Table 4 :
Comparison of results for nondelivery motorcycle crashes on weekdays for diferent models.Signifcance levels of 5%.Te variables and coefcients highlighted in bold within the GWNBR model indicate that their coefcient values demonstrate statistical signifcance at a 5% level in over 50% of the total study units.Please refer to the percentage of signifcant coefcients column for the specifc proportion of study units with signifcance levels of 5% for each variable. *

Table 5 :
Comparison of results for nondelivery motorcycle crashes on weekends for diferent models.