Abstract
Wyoming has one of the lowest seat belt law compliance rates in the United States. Understanding the factors associated with seat belt non-use is crucial for the development of appropriate mitigation measures that aid in increasing seat belt usage rates. In this study, Bayesian binary logistic regression models were developed to: i) investigate the factors associated with daytime seat belt non-use using observational survey data; and ii) assess the factors related to seat belt non-use in fatal crashes while accounting for possible intra-class correlation (effects of unobserved factors common to occupants in the same vehicle). The results indicated that males were less likely to use seat belts than females. Also, trips made in pick-up trucks, those made within urban areas, those made during weekdays, and those made by means of vehicles with Wyoming license plates were likely to involve seat belt non-use. Furthermore, counties with thriving oil and gas industries had lower seat belt use rates. It was also inferred that young occupants, rear seat occupants, impairment, weekdays, traveling around midnight, traveling on collector roads, and traveling on local roads were associated with seat belt non-use. A strong correlation in seat belt non-use among occupants within the same vehicle warrants the implementation of the hierarchical logistic model to investigate seat belt use among occupants involved in fatal crashes. The findings of this study highlight remarkable trends in seat belt use habits in Wyoming, which will aid legislatures in proposing road safety measures aimed at raising seat belt usage proportions.
Motor vehicle crashes are among the leading causes of death in the United States. The seat belt is known to be an effective safety device in preventing fatalities and severe injuries resulting from motor vehicle crashes. By restraining people in their seats during crashes, seat belts prevent the occupants from being ejected from the vehicle. According to Deutermann ( 1 ), ejected occupants involved in rollover crashes are more likely to be killed than restrained occupants, and seat belt use lowers ejection risk. Also, the National Highway Traffic Safety Administration (NHTSA) reported that seat belt usage lowers fatality and serious injury likelihoods for front-seat passengers by almost half ( 2 ). Buckling up prevented roughly 14,955 fatalities in 2017 ( 2 ).
Despite the effectiveness of seat belts in lowering mortality rates and lessening severe injury counts, seat belt non-use still remains a concern. Wyoming has one of the lowest seat belt law compliance rates in the United States according to NHTSA’s National Occupant Protection Use Survey (NOPUS) data ( 3 ). A comparison for Wyoming’s seat belt use rates with nationwide average is depicted in Figure 1.

Seat belt use rates in Wyoming and nationwide. (3)
In 2019, the nationwide seat belt use rate was 90.7%. Yet, in Wyoming, it was 78.3%. In addition, although the overall seat belt use rate has gradually increased over the years in the U.S., it decreased from 86.3% in 2018 to 78.3% in 2019 in Wyoming as shown in the Figure 1 ( 3 ).
Among several factors, seat belt law non-compliance poses an obstacle to efforts to reduce the number of fatal and severe injury crashes in Wyoming. According to the Insurance Institute for Highway Safety (IIHS), the fatality rate per 100,000 residents in Wyoming was 19.2 compared with the nationwide rate of 11.2 in 2018 ( 4 ). Also, Wyoming exhibited a larger proportion of unrestrained occupants involved in fatal crashes (50.6%) than the national average, which was 43.1% ( 4 ).
Understanding the factors associated with seat belt non-use is crucial for the development of appropriate measures that aid in increasing seat belt law compliance rates. Boakye et al. ( 5 ) showed that proper seat belt use is influenced by several factors including the occupant’s demographic characteristics, interpersonal factors (interacting with passengers, seat location), seat belt laws, and environmental factors. Road safety researchers ( 6 – 9 ) examined the association between seat belt laws and seat belt use rates among different states. McCartt and Northrup ( 6 ) and Lee et al. ( 7 ) found that seat belt use rates were greatly dependent on the seat belt laws. Likewise, Boakye and Nambisan ( 8 ) found a positive association between seat belt laws and seat belt use rates. Harper ( 9 ) revealed that switching from secondary to primary laws increased the number of buckled fatally injured occupants. Several studies also stated that primary seat belt law enforcement is more effective at increasing seat belt compliance rates and lowering traffic fatalities compared with secondary law enforcement ( 10 ). In addition to types of seat belt law enforcement (primary versus secondary), increasing fines for seat belt non-compliance and types of coverage with regard to vehicle occupants (whether seat belts are required for all seating positions versus only front-seat occupants) have significant impacts on seat belt use rates ( 11 – 13 ). A detailed description of the studies is provided in the literature review below.
Other studies have demonstrated significant associations between drivers’ seat belt use choices and those of their passengers ( 5 , 14–19), indicating correlations among occupants in the same vehicle. Williams and Shabanova ( 16 ) and Cooper et al. ( 15 ) indicated that a driver’s seat belt use is positively associated with the presence of other passengers and their seat belt use. In particular, Williams and Shabanova ( 16 ) reported that seat belt use among teenage drivers increased when traveling with parents or older passengers and decreased when traveling with peers. Moreover, the seat belt use of both front-seat occupants has been found to be highly correlated. Recently, Chakrabarty et al. ( 19 ) concluded that driver seat belt use choice has a substantial influence on the front-seat passenger’s seat belt use choice and vice versa signifying a strong positive correlation among front-seat occupants’ seat belt use. The finding is consistent with the previous studies ( 17 , 18 ). Boakye et al. ( 5 ) also found a positive association between seat belt use choices among front-seat occupants at night. That is, front-seat passengers tend to not wear seat belts if their accompanying drivers do not buckle up.
Most of the seat belt studies were limited to either analysis of fatal crash data or direct observational study. Fatal crash data may not represent the general motorist population since they are restricted to the persons involved in the fatal crashes while observational surveys compile information from front-seat occupants during daytime conditions only. In seat belt studies, no study was found simultaneously investigating the fatal crash data and using direct observational data. To this end, this study is focused on developing a Bayesian model, which includes both types of data. Moreover, this paper contributes to the body of literature by drawing a comparison between the two types of data for their sample population and modeling results.
The objectives of this study were: i) to investigate the factors associated with daytime seat belt non-use by implementing observational survey data; and ii) to develop a statistical model to assess the factors related to seat belt use in fatal crashes using the observational survey data as prior information while accounting for possible intra-class correlation effects (those of the unobserved factors influencing occupants’ seat belt use choices within the same vehicle). The binary logistic regression structure with Bayesian inference was employed to develop two separate models using the observational survey data and the NHTSA’s Fatality Analysis Reporting System (FARS) data comprising crashes in the years 2010 to 2019
In this study, prior information of seat belt use choice interpreted from an observational survey conducted in Wyoming was used. Also, factors associated with seat belt non-use in fatal crashes considering unobserved vehicle-specific heterogeneity effects were investigated. The findings of this study might assist local legislatures in suggesting appropriate programs to increase seat belt usage rates which would reduce the injury severity of motor vehicle crashes in Wyoming.
Literature Review
The road safety literature is replete with studies investigating the impact of seat belt usage on crash severity and evaluating the association between seat belt laws and seat belt usage rates. In this section, a brief account of the recent traffic safety literature pertaining to the objectives of this work is provided.
Studies assessing the efficacies of the seat belt laws indicated that such laws are the most effective approach for increasing seat belt use rates ( 6 , 7 ). In particular, it was found that the primary seat belt laws (whereby drivers can be halted for unbuckled occupants) are more effective than the secondary seat belt use laws (whereby drivers can only be halted if another traffic infraction was committed). Boakye et al. ( 5 ) examined the association between the adult seat belt use rates and the seat belt laws in the U.S. by processing data from FARS. They implemented the Tukey multiple comparison test and concluded that the states with primary seat belt laws, applicable to each vehicle’s occupant, had larger seat belt use rates than those with secondary or no seat belt laws. Houston et al. also found that seat belt use rates were higher in the states with primary seat belt laws than those with secondary seat belt laws ( 20 ). In addition, imposing higher fines along with standard enforcement has demonstrated significant impacts on seat belt use rates in the states with primary seat belt laws. Harper ( 9 ) explored the influence of enhancing seat belt law enforcement policies on road fatalities by employing a Bayesian data augmentation approach. The author inferred that although switching from secondary to primary laws resulted in a marginal influence on fatality toll reductions, the number of buckled fatally injured occupants rose by 16%. McCartt and Northrup ( 6 ) explored the seat belt use rates of deceased teenage drivers, aged 16 to 19, by implementing a multivariate logistic regression structure. The results showed that seat belt use rates among deceased teenage drivers varied widely among multiple states and were strongly dependent on the seat belt laws. The latter finding was in line with that of Lee et al. ( 7 ).
Previous studies have also shown that demographic, behavioral, and physiological attributes of vehicle occupants are associated with seat belt use habits. Lerner et al. ( 21 ) investigated several demographic characteristics associated with seat belt use choice pertaining to adults injured in road crashes. They concluded that elderly drivers, females, and wealthy individuals were more likely to buckle up. In addition, drivers were found to be more willing to use the seat belt as opposed to their passengers. Beck et al. ( 22 ) compared the vehicle occupant fatality tolls of rural areas to those of urban areas in relation to the seat belt use habits of the victims, aged 18 or above, across the United States. The results indicated that rural areas had larger rates of unbelted occupant fatalities, regardless of the seat belt laws. In another study, Boakye et al. ( 5 ) investigated the factors that contributed to seat belt use choice of front-seat passengers during evenings by employing the generalized estimating equation technique. The results revealed that male passengers were less disposed to use the seat belt than female passengers and that trips made along local roads discouraged seat belt use. Also, front-seat passengers were less inclined to buckle up when riding in cars and pick-up trucks as opposed to when riding in sports utility vehicles (SUVs). Recently, Chakraborty et al. ( 19 ) explored the seat belt use habits of front-seat occupants by employing the bivariate probit modeling framework. The findings indicated that the seat belt use habits of front-seat occupants differed depending on their socio-demographic characteristics. Age and education levels of the occupants also influenced seat belt use choice. Briggs et al. ( 23 ) examined high school students’ seat belt use choices. The authors concluded that students, aged 16 or above, were unwilling to buckle up when riding as passengers as opposed to when driving. Shinar et al. ( 24 ) and Houston et al. ( 20 ) found that occupants with higher education and income levels were more likely to buckle up. Seat belt use rates were also found to increase with age ( 25 , 26 ). Nambisan et al. ( 27 ) interpreted that seat belt use also depended on vehicle characteristics. They investigated the impacts of vehicle age and type on occupants’ seat belt use habits using five years of crash data, 2014 through 2018, collected from FARS. The researchers concluded that occupants of more recent vehicles were likely to buckle up relative to those of dated vehicles.
Several studies demonstrated that occupants in the same vehicle influenced each other in deciding whether to buckle up. Afghari et al. ( 26 ) defined this as the “vehicle atmosphere”. Han ( 14 ) examined passenger seat belt use choice and found that it was substantially dependent on the drivers’ discretions on using the seat belt. The likelihood of passengers opting not to buckle up rose considerably when the driver was unrestrained relative to when the driver was restrained. Chakraborty et al. ( 19 ) also concluded that there was a strong correlation among front-seat occupants making decisions on seat belt use. Nukenine and Daniel ( 28 ) examined whether the seat belt use habits of rear seat passengers induced an effect on the buckled driver’s injury severity. The results revealed that seat beat use status and the count of rear seat passengers affected the driver’s injury severity. Also, the risk of front-seat occupants being severely injured rose when the rear seat occupants were unbuckled. Jermakian and Weast ( 29 ) found that misconceptions about the grave ramifications of not buckling up, lack of knowledge, or lack of comfort incurred as a result of proper restraint use were some of the reasons for not fastening the seat belt among rear seat passengers. These findings were consistent with those of Boyle and Lampkin ( 30 ) and Kidd and McCartt ( 31 ). Eluru and Bhat ( 32 ) demonstrated the need for capturing the mixed effects of the parameters on the outcome. They modeled the concurrent effects of seat belt use and the resulting crash severities as a function of multiple relevant factors. The authors incorporated random parameters and considered correlations among the parameters influencing the crash severity.
Recently, random-effects models and hierarchical models with Bayesian inference have become popular in road safety studies. This is a result of their ability to account for the correlations among the parameters and the hierarchical nature of the data ( 33 ). Huang and Abdel-Aty ( 34 ) proposed a multi-tiered data framework to highlight the need to select hierarchical Bayesian modeling approaches in crash data analysis. Haung et al. ( 35 ) examined the influential factors that contributed to intersection-related crashes by implementing the hierarchical Bayesian binomial logistic regression model. Haque et al. ( 36 ) applied several hierarchical models to model motorcycle crash frequencies at signalized intersections. Ahmed et al. ( 37 ) employed the hierarchical ordered logit model with the no U-turns Hamiltonian Markov chain Monte Carlo (MCMC) technique in their analysis of resulting severities of hazardous material truck crashes. Haq et al. ( 38 ) used the Bayesian binary logistic regression framework to model injury severities of crashes that occurred as a result of tire failures. Afghari et al. ( 39 ) implemented the Bayesian multivariate binary model with latent variables to assess whether vehicle occupants influenced each other’s seat belt use choices.
Seat belt use studies were limited to either analysis of fatal crash data or direct observational study. Fatal crash data may not represent the general motorist population since they are restricted to the persons involved in the fatal crashes while observational surveys compile information from front-seat occupants during daytime conditions only. The current study is focused on developing a Bayesian model, which combines both types of data. Seat belt use choice data were obtained from observational surveys and FARS. The approach accommodates potentially interrelated seat belt use habits of both drivers and passengers. It also permits parameters to vary randomly across vehicle occupants, accounting for unobserved heterogeneity effects, such as roadway characteristics, vehicle attributes, and driver behaviors. To the best of the authors’ knowledge, no previous studies on seat belt use involved two data sets concurrently to capture the variability in the data.
Data Description
Two data sets were used in this study. The first was collected from a seat belt observational survey conducted by trained observers during the first week of June 2019 while the second was collected from FARS. It is worth mentioning that observational seat belt surveys are conducted every year to explore seat belt use rates among front-seat occupants in vehicles across many states. Wyoming observational surveys have also been conducted every year since 2012 in compliance with the national standards ( 40 ). Therefore, the data set represents a wider time frame since it is collected once every year. To ensure that the status of seat belt use with regard to the front-seat occupants within the same vehicle was recorded, the Uniform Criteria for State Observational Surveys of Seat Belt Use 23 CFR (Code of Federal Regulations) § 1340 standards were followed when conducting the survey ( 40 ).
For the first data set, a total of 24,821 observations of drivers and front-seat passengers in 18,286 vehicles were collected from 289 sites in 17 counties of Wyoming. According to the NHTSA, sample counties are selected as the primary sampling unit of the survey based on the vehicle miles traveled or population ( 41 ). It is worth noting that site locations are selected in a way that results in residents being likely to be observed ( 40 ). Moreover, multiple observation sites were picked to capture seat belt usage from a variety of motorists and avoid biases in data ( 41 ). More than two-thirds of the data were collected from secondary and local roads to guarantee that a majority of the observed vehicle occupants were local residents of the respective counties ( 40 ). This survey was conducted to explore seat belt usage by county and the association between seat belt use and several variables. These included population density, weather, vehicle registration, occupant gender, vehicle type, and occupant type (driver or passenger). Standard protocols were followed from observer training to data analysis to ensure the reliability and accuracy of the observed data ( 40 ). Figure 2 shows trends in seat belt use among occupants riding in vehicles in Wyoming as per the survey results and Table 1 lists summary statistics of the survey’s variables.
Summary Statistics of the Observational Survey’s Variables
Note: SUV = Sports utility vehicle.

Occupants’ seat belt use proportions in Wyoming for 2019.
The seat belt usage rate varied from a low of 63.5% belted occupants in Sweetwater County to a high of 97.8% in Niobrara County (Figure 2). In this survey, areas with more than 5,000 residents were identified as “urban,” while areas with less than 5,000 residents were considered “rural.” Also, 76.1% of the observations were witnessed in rural areas. According to the survey’s results, seat belt use rates were higher in rural than in urban areas. Similarly, seat belt usage rates were found to be higher on primary roads than on secondary and local roads. Furthermore, occupants were more likely to wear their seat belts during weekends. Although males comprised more than 80% of all vehicle occupants, they were 13.7% less likely to use their seat belts than female occupants.
Seat belt use was observed for a variety of vehicles, but they were broadly categorized as passenger cars, SUVs, pick-up trucks, and passenger minivans. With regard to vehicle types, pick-up truck occupants were likely to have the lowest seat belt use rate. Other than the vehicle type, weather conditions were condensed into two categories: clear and not clear. Observers also recorded the occupants’ vehicle license plates, whether they were Wyoming plates or not. With that, the survey results indicated that occupants of out-of-state vehicles exhibited higher rates of seat belt use than their Wyoming counterparts.
The second data set was collected from NHTSA’s FARS database. The database was queried to select records of Wyoming’s fatal crashes occurring from 2010 to 2019. This data set included records of 2,705 vehicle occupants who were involved in 1,135 fatal crashes, regardless of whether they suffered a fatal injury or not. The crash records included variables related to occupant demographics, vehicle characteristics, roadway characteristics, crash type, and environmental conditions at the times of the crashes. Categorical and quantitative variables related to seat belt use were extracted from the crash records to ascertain the factors that contributed to seat belt use habits. Roadway types were classified as interstates, primary arterials, and collector/local roads. The times of crashes, included in this analysis, were categorized into three classes: daytime (5 a.m.–5 p.m.), nighttime (5 p.m.–midnight) and midnight (midnight–5 a.m.). Although crashes occurred on all seven days of the week, the days were merged into weekdays and weekends for reasonable interpretation of the results. Likewise, the months, during which the crashes occurred, were conjoined into quarters: January to March as the first quarter, April to June as the second quarter, July to September as the third quarter, and October to December as the fourth quarter. Occupant characteristics, such as gender, age, seat location, and type of occupant (driver or passenger) were also included in this analysis. According to NHTSA, occupants 65 years and older are considered as elderly residents ( 42 ) while the Wyoming child restraint law states that children aged 9 years and younger are required to use child restraint systems ( 43 ). For maintaining consistency with NHTSA and Wyoming laws, the occupant ages were classified into seven categories: children (9 years and younger), teenage/young occupants (10–18 years), 19–30 years, 31–40 years, 41–50 years, 51–65 years, and elderly occupants (66 years and older). Note that the intermediate groups were classified such that the sample size of each age group was adequate and roughly equal to those of the others. The occupants’ seating locations were dichotomized as front and rear seats. In total, 1,104 occupants (41% of the total occupants) were unbuckled at the times of the crashes. Detailed analysis on seat belt non-use is provided in the following paragraph. Table 2 presents the summary statistics of the data set extracted from FARS while Figure 3 depicts the seat belt use rates among occupants riding in vehicles in Wyoming.
2010 to 2019 Fatal Crash Data Variables’ Summary Statistics

Occupants’ seat belt use rates in fatal crashes from 2010 to 2019 in Wyoming.
Seat belt non-usage rates varied across counties from a high of 61% unbuckled occupants in Crook County to a low of 21% unbuckled occupants in Teton County (Figure 3). Table 2 shows that seat belt non-use rates ranged from 37% in the second quarter to 45% in the first quarter of the crash years. Furthermore, seat belt non-use rates were higher during weekends than weekdays. They were also observed to be the highest during midnight conditions.
The lowest passenger seat belt non-use rate was observed on principal arterials and the highest was observed on collector/local roads. In regard to gender, males had larger proportions of seat belt non-use rates than females. Seat belt non-use rates were lower in children and elderly occupants. By contrast, they were found to be higher among teenagers and young adults (aged 10–18 and 19–30, respectively). In regard to occupant types, seat belt non-use rates were higher among passengers than among drivers. Likewise, rear seat occupants had a larger proportion of seat belt non-use than front-seat occupants. Finally, alcohol or drug involvement substantially reduced seat belt use rates.
Research Methodology
In this study, Bayesian binary logit models with both fixed and random effects were applied to investigate the factors associated with seat belt use. Bayesian models treat the parameters as random variables, unlike frequentist models, which consider the parameters as fixed. In Bayesian analysis, the observed data are used to update the prior beliefs about the behavior of the parameters when estimating their distributional properties ( 44 ). One of the major advantages of Bayesian frameworks is their ability to select parametric families for prior probability distributions. Three types of prior distributions are mainly used in Bayesian approaches for estimating parameters including:
Informative priors based on the literature or previous relevant studies;
Weakly informative priors, which provide more realistic inferences being selected from trends (means and standard deviations); and
Non-informative priors, which allow the information from the likelihood to be interpreted.
Past studies ( 37 , 45 ) have shown that the incorporation of weakly informative priors is recommended since they enhance the understanding of the parameters’ distributional properties by confining the resulting posterior probabilities obtained within a reasonable range. In this study, prior distributions were drawn from the preliminary data analysis of the observational survey data and assumed to be normally distributed. Weekly informative priors were selected in a way that parameter estimates can be calculated with minimum bias and at the same time model convergence can be attained in substantially less time.
Various MCMC sampling techniques have been widely used, such as the Gibbs sampler, Metropolis Hastings algorithm, or the combination of both using the freeware, WinBugs, to obtain samples from the posterior distribution. Recently, the no U-turns Hamiltonian MCMC method, an extension of the Hamiltonian MCMC sampling technique, has been gaining momentum given its effectiveness and robustness over the other sampling techniques ( 37 , 46 ). In this study, the Bayesian Regression Models (brms) package of the R® statistical software is implemented to estimate the models ( 47 ).
Bayesian logistic regression models were used to model the relationship between the dichotomous response variable (seat belt non-use versus use) and the explanatory variables including occupant characteristics, environmental attributes, and other factors while considering within vehicle correlations. The model structure is described as follows. Let Yij denote the dependent variable, which represents the seat belt use choice of the ith occupant within the jth vehicle. The outcome of this distribution is binary, either 0 (buckled) or 1 (unbuckled). The probability of not using the seat belt, Yij = 1, is denoted as P(Yij) and conforms to the logistic regression equation, expressed as:
Equation 1 represents the fixed effect binary logistic regression while Equation 2 presents the random effect binary logistic regression model. The term, α, is the model’s constant and uj represents the random effects of the vehicle. It should be noted that uj is assumed to be normally distributed with a mean of 0 and a standard deviation of σ, uj∼N(0, σ 2 ). Also, Xp denotes the explanatory variables and βp is their respective regression coefficients. It is worth mentioning that the fixed effect Bayesian binary logit model was used to analyze the observational survey data while the random effect Bayesian binary logit model was employed to investigate the FARS data.
In Bayesian analysis, the resulting posterior distribution, P(θ|Y), is equivalent to the product of the likelihood of the observed data, P(Y|θ), and the prior information of the parameters, P(θ), divided by the average likelihood, P(Y), as shown in Equation 3.
The average likelihood, P(Y), is also known as the marginal distribution of Y. It is computed as follows:
Since P(Y) is independent of θ, the posterior distribution is only proportional to the product of the likelihood and the prior as shown in Equation 5.
An intra-class correlation (ICC) coefficient was incorporated in this study to quantify the unobserved variability of seat belt non-use among occupants within the same vehicle. The following equation was employed to compute the ICC.
In Equation 6, σ 2 represents the variability of seat belt non-use rates across vehicles, and σO 2 denotes the variability within the vehicle, which is assumed to be fixed across vehicles and equal to π2/3 ≈ 3.293 for a hierarchical logistic distribution ( 37 ). An ICC value near 0 indicates a small variation that can be explained by the group-level differences and, thus, a non-random effect ordinary logit model may be adequate to fit the data. On the other hand, an ICC value near 1 indicates substantial variations among groups and the use of the hierarchical model is, therefore, justified ( 48 ).
In the analyses, Bayesian credible intervals (CIs) describe the uncertainty related to the covariates by assigning a range of values obtained from their posterior probability distributions. In this study, the 95th percentile Bayesian CI, which consists of the central portion (2.5%–97.5%) of the posterior distribution of the covariates, was used. Covariates, of which 95th percentile CIs did not contain 0, were identified as significant predictors ( 49 ).
Odds ratios (ORs) were estimated to interpret the effect of the significant variables on seat belt non-use. The OR is the odds of observing an unbuckled occupant in the presence of a particular parameter relative to those in the absence of the parameter, provided that all else was unchanged.
Before model fitting, multicollinearity was assessed among the explanatory variables using the correlation test, tolerance and variance inflation factor (VIF). Absolute correlation coefficients near 1 represent strong correlations between variable pairs and VIF values, while greater than 10 indicate the presence of strong multicollinearity ( 50 ). Also, tolerance levels less than 0.1 denote the existence of multicollinearity ( 51 ). In this study, the three measures’ results ensured the absence of multicollinearity among the explanatory variables.
Among several model goodness of fit measures, the Watanabe–Akaike information criterion (WAIC) ( 52 ), has been gaining popularity as a performance measure of Bayesian models. As opposed to the Akaike information criterion (AIC) and the deviance information criterion (DIC), WAIC is capable of averaging the likelihoods, P(Y|θ)s, over the posterior distribution ( 49 ). The WAIC includes the log-pointwise predictive density (LPPD) in its formulation, which is calculated as follows ( 49 ):
The term, S, is the number of MCMC samples drawn from the posterior distribution. A correction factor that accounts for the effective number of parameters, pWAIC, is included in the WAIC estimation. It is computed as ( 49 ):
Finally, the WAIC measure is computed by the following formula ( 49 ):
Posterior predictive assessments may also be undertaken to verify whether the model adequately explains the observed variability in the data. The rationale of those assessments is that they indicate whether the model can be applied to generate data that are similar to the observed data, provided that it exhibits a good fit ( 53 ). Both the WAIC metric and the posterior predictive assessments were employed in this study to evaluate model performances.
Results and Discussions
Two Bayesian models were developed to aid in understanding seat belt use habits in Wyoming. The first was a binary logit model with seat belt use as the response. The 2019 observational seat belt survey data for Wyoming was used for its estimation. The second was a random-effects model with the seat belt use as the outcome as well. It was estimated using the FARS data pertaining to Wyoming’s conditions from 2010 to 2019. Note that the means of the seat belt use choices, obtained from the observational survey, were incorporated as priors in the constant term of the latter model. As a method of statistical regularization, weakly informative priors shrank the parameter estimates toward 0 and reduced the Type I error rates ( 45 ). Similar to previous studies ( 37 , 38 , 54 , 55 ), incorporating such priors in this study helped to constrain posterior distributions obtained to a reasonable range, which in turn enabled to estimate the parameters with minimum bias. Weakly informative priors also facilitated attaining model convergence in substantially less time. Three Monte Carlo Markov chains were prepared for each model. With the aid of R’s “brms” package, the Hamiltonian MCMC sampling method was conducted to obtain the posterior distributions of the samples. For each chain, 3,000 iterations were simulated with 1,000 warm-up samples after rigorous fine-tuning of the model specifications to achieve convergence. In total, 6,000 posterior samples of the three chains were obtained. In addition, a target acceptance rate of 90% and a maximum tree depth of 10 were considered in the models. The maximum tree depth indicates the maximum number of allowable steps of no U-turn sampling. It is used to prevent premature termination ( 48 ). The target acceptance rate is used to control step size. The higher the target acceptance rate, the lower the step size resulting in less computation time.
Model convergence in this study was ensured by generating trace plots and results of the Gelman–Rubin convergence diagnostic, Rhat ( 56 ). The generated trace plots showed that the chains mixed well, indicating identical distributions of the simulations between and within chains. The Rhat measure, also known as the “potential scale reduction factor,” is the square root of the ratio of the variance of all the chains to the average within-chain variance ( 46 ). In the final model, variables were identified as significant based on the 95th percentile Bayesian CI. The effective sample size (ESS) metric gauges the level of autocorrelation of the samples within the chains. The closer the ESS to the total number of samples, the more independent samples are within the chains. In this study, the ESS results were found to be close to the total number of samples, which indicated that the samples were independent. The results of both models are presented in Tables 3 and 4. In the tables, the standard errors (SEs) are the approximate standard deviations of the sample populations. Figure 4 presents sample trace plots of several significant variables from the random-effects model.
Observational Survey Data Model Results
Note: SE = standard error; OR = odds ratio; WAIC = Watanabe–Akaike information criterion.
Fatality Analysis Reporting System Data Model Results model
Variables significant at the 90% credible interval.
Note: SE = standard error; OR = odds ratio; ICC = intra-class correlation.

Sample variable trace plots.
Remarks on the Observational Survey Data Model Results
The results of the binary logit model, presented in Table 3, reveal remarkable insights into seat belt use choice. It should be noted that incorporating the vehicle identification as a variable to capture correlations among occupants in the same vehicle was worth the attempt. However, this variable was not available in the data set. Nevertheless, all estimated parameters were identified as influential since their Bayesian 95th percentile CIs did not contain 0. The impact of each parameter on the response was interpreted assuming all else was controlled. The survey results demonstrated that vehicle occupants in Campbell, Carbon, Converse, Johnson, Park, Sheridan, and Sweetwater counties were likely to be unrestrained by 2.97, 3.90, 2.34, 1.30, 2.05, 1.20, and 3.94 times, respectively compared with those of which data were sampled from Albany County. On the contrary, occupants observed in Crook were less likely to be unrestrained by 0.34 times. These findings were consistent with those of previous studies ( 5 , 22 ). Furthermore, the oil and gas industry is the single largest economic driver in Wyoming, while Sweetwater, Campbell, Carbon, Park, Converse, Sheridan, Johnson, and other counties are characterized by mining activities. From Table 3, it is evident that occupants observed in these counties were more likely to be unrestrained, which led to the inference that oil and gas industry workers might have been unwilling to use the seat belts. According to the Wyoming Department of Workforce Services, 38% of fatally injured workers were unbuckled throughout the years 2012 to 2018 ( 57 ). Statistics also suggested that oil and gas industry workers experience six times the rate of motor vehicle-related death compared with workers in other industries ( 58 ). Another study also found a substantial increase in crash severity and frequency given the expanding oil industry in North Dakota ( 59 ). However, observational survey data do not have information with regard to occupants’ occupations. Future studies should explore additional evidence if there is an association between low seat belt use and oil/gas industry workers’ trips. In this study, oil and gas production data from 2019 were collected from the Wyoming Oil and Gas Conservation Commission ( 44 ), and analyses were conducted to assess the correlations between seat belt non-use and oil/gas production (Table 5). As shown in the table, such correlations existed.
Correlation Matrix of Seat Belt Non-Use Proportions and Oil and Gas Production
The difference in seat belt use habits among the counties may also be a result of the level of urbanization, differences in occupant attitudes toward enforcement and seat belt use awareness programs across the different counties. Appropriate interventions for each county ought to be considered to increase seat belt usage rates.
In the case of the occupant gender, the results indicated that the estimated odds of not using the seat belt increased by 1.6 times for males, on average. This finding was consistent with several previous studies ( 5 , 19 , 60 ).
Significant differences in seat belt use choice were also found among occupants of different vehicle types. According to the results, pick-up truck occupants were more likely to refrain from buckling up compared with passenger car occupants (OR = 1.28). On the other hand, when riding in vans and SUVs, the estimated odds of seat belt non-use decreased by 0.27 and 0.39 times, respectively. In multiple past studies, lower seat belt usage rates among occupants of pick-up trucks were observed compared with those of other vehicle types including sedans, SUVs, and minivans ( 19 , 60 , 61 ). Combining media campaigns with high visibility programs, such as “Buckle up in your truck,” were demonstrated to be effective in prompting pick-up truck occupants to fasten their seat belts as per ( 5 ).
As shown in Table 3, vehicle occupants traversing within urban areas were more likely to be unbuckled by 1.54 times compared with those traversing within rural areas. The finding is in line with the research carried out by Lipovac et al. ( 62 ) and Huang et al. ( 63 ). They inferred that drivers were unwilling to wear seat belts when traveling short distances on low-speed roads in urban areas ( 62 , 63 ). Although these results, low seat belt use in urban areas, contradicted most of the previous studies ( 50 , 51 ), it might be a result of a compensation effect in which travelers perceived less risk of not using the seat belt when traveling shorter distances in urban areas. On the other hand, given geographically dispersed/sparsely located amenities in rural areas ( 64 ), trip-takers traveled long distances to reach their destinations ( 65 ) and, therefore, would be inclined to buckle up. On a different note, the notion of rural and urban is very different in this least populated state relative to the more populated and geographically smaller states which may also be a reason for the opposing finding ( 66 ).
When it comes to the influence of the day of the week, the estimated odds of not using the seat belt were 1.30 times higher on weekdays than on weekends. A similar result was found in a previous study ( 60 ). This finding might also lead to the interpretation that oil and gas industry workers were likely to be unrestrained while commuting short distances to their job sites during weekdays.
Other than the day of the week’s influence on seat belt use habits, the findings revealed that occupants of vehicles registered in Wyoming had a lower chance of buckling up than those of out-of-state vehicles (OR = 1.42). Plausibly, this might be because Wyoming’s residents perceive lower risks of not using seat belts when traveling on familiar roads which consequently makes them more susceptible to crashes. Historically, Wyoming residents are two to three times more likely to be unrestrained and killed in a crash than nonresidents ( 67 ). According to the Wyoming DOT, two-thirds of all drivers involved in crashes had a Wyoming driver’s license in 2020 ( 68 ). Further studies should be carried out to investigate the reasons behind the unwillingness to use seat belts among Wyoming residents.
When it comes to the road type, the likelihood of not using seat belts increased on local and collector roads compared with interstates. As a result of the low speed, short distances, and road facilities, travelers perhaps perceived lower risks on local and collector roads, discouraging the use of seat belts. The number of lanes per direction also induced effects on seat belt use choice. As shown in Table 3, the odds ratio was estimated as 1.13 when having two lanes relative to one, a finding in line with that of Rezapour and Ksaibati ( 69 ).
Remarks on the Fatality Analysis Reporting System Data Model Results
Table 4 presents the comparison between the results of the random-effects Bayesian model and of the fixed-effects model, which were developed to investigate the factors associated with seat belt non-use in Wyoming using the FARS data from 2010 to 2019. The WAIC of the random-effects model was found to be lower than that of the fixed-effects model, indicating that the random-effects model demonstrated better performance than the fixed-effects model. Descriptive statistics of the observational survey were incorporated as weakly informative priors into the model. As previously stated, studies have demonstrated the existence of correlations among vehicle occupants affecting seat belt use habits ( 5 , 14 ). In this study, the ICC value was computed to measure such correlations. The variance of the occupants’ seat belt use choices (estimated as 4.16) was used to calculate the ICC as follows:
The ICC value indicated that 55.8% of unexplained variations were a result of the presence of vehicle-specific unobserved factors associated with seat belt non-use ( 70 ). McElreath ( 48 ) recommended the use of hierarchical models when unobserved heterogeneity effects were observed in the model. Therefore, given the existence of vehicle-specific unobserved factors associated with seat belt non-use, the hierarchical logit model was employed to account for this unobserved heterogeneity.
The model parameters were interpreted provided that all else was unchanged. As shown in Table 4, vehicle occupants in Big Horn, Crook, Sheridan, and Sweetwater counties had 2.77, 4.06, 2.72 and 1.97 higher odds of being unrestrained, respectively, compared with those observed in Albany County. These findings were consistent with those of previous studies ( 8 , 22 ).
Other than the seat belt use rates by county, it could be inferred from the model results that occupants were more likely to be unrestrained during midnight (midnight–5 a.m.) conditions compared with when riding during daytime conditions. This result was in line with that of ( 5 ) who interpreted that refraining from using the seat belt during the evening might be a result of the misconception of the improbability of receiving a citation at night among other possible reasons. Several other studies have also shown that the time of day has a substantial impact on road users’ seat belt usage ( 61 , 71–73). Another study found that single-vehicle occupants are less likely to be restrained in the morning while multi-vehicle occupants are less likely to be restrained in the afternoon ( 60 ).
In the case of occupant ages, it was found that the age groups, 10 to 18 years and 19 to 30 years, had 1.62 and 1.52 higher odds, respectively, of being unrestrained compared with occupants of the age group, 51 to 60 years. By contrast, the estimated odds of being unrestrained decreased by 0.66 and 0.51 times for children (9 years and younger) and elderly occupants (66 years and older), respectively, compared with occupants aged 51 to 65. These inferences were consistent with other studies ( 14 , 21 ). Furthermore, NHTSA reported that, in 2017, 54% of unbuckled car occupants involved in crashes, aged between 13 and 15, were fatally injured in road crashes ( 2 ).
The seat location was another factor that was found to have an impact on seat belt use habits. The odds of rear seat passengers not using the seat belt were 2.41 times higher than those of front-seat passengers. The false perception that rear seats would be safer than front seats might make rear seat occupants less likely to buckle up. According to Findley et al. ( 74 ), buckled rear seat occupants not only protected themselves but also encouraged front-seat occupants to use their seat belts. As such, appropriate mitigation measures ought to be implemented to increase seat belt use rates among rear seat occupants. Alcohol involvement was also found to play a significant role in seat belt use choice. From the analysis results, the estimated odds of not using the seat belt increased by 9.39 times as a consequence of alcohol intoxication. Han ( 14 ) concluded that passengers were likely to refrain from using their seat belts when the drivers were driving under the influence of alcohol and Afghari et al. ( 26 ) reported that passengers would be less willing to buckle up when they themselves were riding under the influence.
Other than alcohol involvement, the roadway functional classification was also found to influence seat belt use choice. The likelihood of not using the seat belt was larger for collectors and local roads compared with interstates. This finding was in line with that of the observational survey model results and of previous studies ( 5 , 60 ). Occupant gender and crash dates were other variables that were investigated in this analysis. However, those parameters were not found to have an impact on seat belt use habits (Bayesian CI included 0).
Incorporating random parameters into the model to capture unobserved heterogeneity effects is another approach worth selecting ( 75 ). The authors attempted the model with the age groups of the occupants as random parameters. However, the WAIC value of the random parameters model indicated that the random-effects model exhibited a better fit than the random parameters model (3,305 and 2,854 for the random parameters model and the random-effects model, respectively).
Conclusions
Previous seat belt use studies have focused either on the analysis of fatal crash data or on the use of observational survey data. This study offers a contribution to the occupant protection literature in that prior information about seat belt use choice, obtained from an observational survey, was introduced in the investigation of the factors associated with seat belt use in fatal crashes. Also, unobserved heterogeneity effects among vehicles were considered. Two Bayesian binary logistic regression models were developed. One was estimated using records of daytime seat belt use data collected from a survey disseminated in Wyoming in 2019. The second was developed for fatal crash records, collected from NHTSA’s FARS data on crashes from the years 2010 to 2019 in the state.
The results of the observational survey model indicated that males were less likely to be restrained than females and occupants traveling in pick-up trucks had a lower chance of using their seat belt compared with occupants of passenger cars. Furthermore, urban area trips were characterized by a lower probability of seat belt use relative to rural area trips. Similarly, the results indicated that weekday trips were less likely to be characterized by restrained occupants as opposed to weekend trips. In addition, occupants of vehicles having Wyoming license plates were less inclined to use their seat belts than occupants of out-of-state vehicles. The results also showed that counties with vibrant oil and gas industries such as Campbell, Carbon, Converse, Johnson, Park, Sheridan, and Sweetwater had higher seat belt non-use rates. It is thus inferred that workers/employees of the oil and gas industries may refrain from using seat belts. This finding is crucial for legislatures since they may target this specific group.
The results of the second model revealed insights into the seat belt use choices of occupants involved in fatal crashes based on a variety of factors and unobserved heterogeneity effects among vehicles. The ICC value signaled a strong correlation in seat belt non-use among occupants of the same vehicle. It was also inferred that vehicle occupants, on whom records were sampled from Big Horn, Crook, Sheridan, and Sweetwater, were more likely to be unrestrained compared with those on whom records were sampled from Albany County. In addition, the results indicated that the occupants, representing the fatal crash data, were more likely to refrain from using the seat belts during the late night hours of midnight to 5 a.m. compared with those traveling during the daytime (5 a.m. to 5 p.m.). It was also interpreted that young occupants (10 to 18 years and 19 to 30 years) had a higher chance of being unrestrained relative to occupants, aged 51 to 65. On the contrary, children (9 years and younger) and elderly occupants (66 years and older) were less likely to be unrestrained in relation to occupants aged 51 to 65. Also, rear seat occupants were found to be less disposed to buckle up compared with front-seat occupants. In addition, it was inferred from the findings that alcohol impairment considerably raised the probability of seat belt non-use. The roadway functional classification was another factor that was found to influence seat belt use habits. That is, occupants were less likely to wear their seat belts when traveling along collector and local roads.
Other than the inferences of the individual explanatory variables from the models’ results, seat belts use choices may also vary across the occupant population. As previously stated, observation survey results are those of front-seat occupants traveling during the day in 2019 while the NHTSA’s FARS data represent occupants involved in fatal crashes during 2010 to 2019. Given the variability among occupant populations and time frame, inferences drawn from the observational survey data modeling results differed from those of the FARS data modeling results. For instance, in the results of the daytime observational survey data model, vehicle occupants observed in Campbell, Carbon, Converse, Johnson, Park, Sheridan, and Sweetwater counties were more likely to be unrestrained. Yet, in the results of the fatal crash data model, vehicle occupants observed in Big Horn, Crook, Carbon, Sheridan, and Sweetwater counties were more likely to be unrestrained. Specifically, we can see inverse results of seat belt use rate in the crook county in two models. In addition, the former model’s results indicated that male occupants were more likely to be unbuckled whereas the occupant gender variable was not found to have an impact on seat belt use choice in the results of the latter model. Likewise, the day of the week variable (weekday versus weekend) was identified as an influential parameter only in the former model. Previous studies also confirmed the variability of occupant seat belt use policy compliance rates across the population (5, 76 , 77 ). Therefore, considering the variability of occupants across the population is a crucial factor when proposing mitigation strategies and seat belt use campaigns.
Wyoming has one of the lowest seat belt use rates in the United State. The findings of this study highlighted noteworthy trends in seat belt use habits in the state. This will aid policymakers in their decision-making processes with regard to the development and implementation of road safety programs aimed at increasing seat belt use rates. Several management measures including upgrading to primary seat belt laws, raising penalties for violating seat belt laws, and nighttime enforcement programs might be effective in reducing the seat belt non-compliance rate in Wyoming ( 78 ). Past studies have demonstrated strong evidence to support upgrading seat belt laws from secondary to primary laws by showing that primary seat belt law is more effective in lowering seat belt non-compliance rates compared with secondary laws ( 8 , 9 , 20 ). In addition, imposing higher fines along with standard enforcement has demonstrated significant impacts on seat belt use rates. Since the FARS data model revealed that vehicle occupants were more likely to be unrestrained during midnight, it could be helpful to implement nighttime seat belt law enforcement programs to increase seat belt use. Also, counties with thriving oil and gas industries exhibit higher seat belt non-use rates in the observational survey data model. Effective outreach and education interventions for those counties such as targeted educational programs to point out the benefit of seat belt use, and the risks associated with not using seat belts ought to be considered. It was also inferred that young occupants and impairment were positively associated with seat belt non-use. Combining law enforcement and media coverage is proved to be effective in increasing seat belt use rates for target populations such as men, teens, and young adults ( 78 ).
Finally, there were several limitations to this study worth discussing. First, the unobserved heterogeneity effects of the factors that contribute to seat belt use choice from the observational survey data were not examined. Incorporating random parameters into the model to capture those unobserved heterogeneity effects is an approach worth selecting ( 75 ). Furthermore, various spatial autocorrelation techniques, such as geographical weighted regression (GWR), spatial econometric modeling, and the Besag–York–Mollié structure have been employed to model traffic crashes ( 79 , 80 ). The spatial autocorrelation ought to be considered in investigating geographical factors that relate to seat belt use choice in future studies.
Footnotes
Author Contributions
The authors confirm contribution to the paper as follows: study conception and design: S. Nazneen, A. Farid, K. Ksaibati; data collection: S. Nazneen, A. Farid, K. Ksaibati; data analysis and interpretation of results: S. Nazneen, A. Farid, K. Ksaibati; draft manuscript preparation: S. Nazneen. All authors reviewed the results and approved the final version of the manuscript.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The authors would like to thank the Wyoming Department of Transportation (DOT) for sponsoring this study and for providing the seat belt use survey data. This work is part of Project #RS02220 funded by the Wyoming DOT. The authors would also appreciate the effort of the Mountain Plains Consortium for providing matching funds under grant 69A3551747108 (FAST Act).
All opinions, stated in this manuscript, are solely those of the authors. The subject matter, all figures, tables, and equations, not previously copyrighted by outside sources, are copyrighted by Wyoming DOT, the State of Wyoming, and the University of Wyoming. All rights reserved copyrighting in 2021.
