Young drivers are overrepresented in crashes, and newly licensed drivers are at high risk, particularly in the months immediately post-licensure. Using a virtual driving assessment (VDA) implemented in the licensing workflow in Ohio, this study examined how driving skills measured at the time of licensure contribute to crash risk post-licensure in newly licensed young drivers.
This study examined 16 914 young drivers (<25 years of age) in Ohio who completed the VDA at the time of licensure and their subsequent police-reported crash records. By using the outcome of time to first crash, a Cox proportional hazard model was used to estimate the risk of a crash during the follow-up period as a function of VDA Driving Class (and Skill Cluster) membership.
The best performing No Issues Driving Class had a crash risk 10% lower than average (95% confidence interval [CI] 13% to 6%), whereas the Major Issues with Dangerous Behavior Class had a crash risk 11% higher than average (95% CI 1% to 22%). These results withstood adjusting for covariates (age, sex, and tract-level socioeconomic status indicators). At the same time, drivers licensed at age 18 had a crash risk 16% higher than average (95% CI 6% to 27%).
This population-level study reveals that driving skills measured at the time of licensure are a predictor of crashes early in licensure, paving the way for better prediction models and targeted, personalized interventions. The authors of future studies should explore time- and exposure-varying risks.
Newly licensed drivers are at high risk of crashes, particularly in the months immediately post-licensure. Inadequate skill at the time of licensure is one critical risk factor.
This population-level study reveals that particular driving skills measured at the time of licensure are a predictor of crashes early in licensure, paving the way for targeted interventions.
Young drivers (aged 15–20) are overrepresented in crashes; despite accounting for only 5.1% of licensed drivers in 2020, young drivers made up 12% of all drivers involved in police-reported crashes and 8.5% of those in fatal crashes.1 In addition, newly licensed young drivers are most at risk; average crash rates peak in the months immediately post-licensure and then decline to the level of experienced adults over the first few years of licensure.2–4 Despite the initial success of the Graduated Driver Licensing policy that restricts driving to lower-risk conditions,5,6 young driver crash rates remain high.1 However, not all young new drivers crash; only 20% crash in the first year of licensure.3 In addition, previous studies have revealed variability in risk among young drivers.7,8 Thus, identifying individual factors that contribute to increased crash risk affords the opportunity to direct resources to those most at risk.9
The authors of previous epidemiologic, naturalistic, and simulated driving studies have identified inadequate skills (eg, in situation awareness and speed management) as critical reasons for most novice young driver crashes.10–15 Thus, driver training is one obvious way to reduce crashes. However, most states have no on-road training requirements for young drivers, which can largely be attributed to the negative conclusions drawn from a single randomized controlled trial of pre-licensure training in DeKalb in 1983.16,17 More recent observational studies indicate that comprehensive young driver licensing policies (Graduated Driver Licensing, driver education, and on-road training) are associated with lower crash rates.18–21 However, these interventions are broad, time intensive, costly, and typically only for drivers <18 years of age. Recent studies of freely available online training targeting situational awareness and hazard detection and mitigation (common skill deficits in young drivers15,22 ) reveal promise for reducing young driver crashes.23–25 However, the findings are, again, in general populations of young drivers. Given limited resources, it is important to target skill training interventions toward those most at risk via individualized assessments of young driver skills.
The authors of a previous study described the development and successful implementation of a virtual driving assessment (VDA; Diagnostic Driving, Inc.) administered at the time of license examination in the state of Ohio26 to assess license applicants’ safety-critical driving skills (by exposing them to common serious crash scenarios). Previous studies have revealed that this tool can predict performance on the Ohio on-road license examination, which consists of a maneuverability test (part 1) assessing basic vehicle control via steering the car around markers, and a driving skills test (part 2) assessing handling turns, starts, stops, reverses, signal use, lane choice, and maintaining safe following distances. One study revealed that the VDA can classify drivers likely to fail the on-road examination with reasonably high accuracy,27,28 but with limited insight into the specific skill deficits predicting failure. A second study distilled VDA performance metrics (capturing things like speed control, lane position, and following distance) into 20 “Skill Clusters” that were then grouped into 4 major summary “Driving Classes.” Compared with the average odds of failing the licensing examination, both the “No Issues” and “Minor Issues” Driving Classes had lower odds of failing, whereas the “Major Issues” and “Major Issues with Dangerous Behavior” Driving Classes had higher odds of failing.
Building on this previous work, in our current study, we seek to determine if time-of-licensure performance on the VDA (ie, Skill Cluster or Driving Class membership) can predict the time to crash post-licensure. If results reveal that the VDA-measured skills at the time of licensure are predictive of crash risk early in licensure, this will advance the development of better prediction models and delivery of targeted and personalized interventions in the future (acknowledging that there are crash-associated skills, behaviors, and attitudes that the VDA does not measure in this study).
Methods
Dataset Preparation
This study used a deidentified analytic dataset of linked VDA, licensing, and crash record data (see Data Sources section below for more detail) prepared by a data operations team at the Children’s Hospital of Philadelphia (CHOP), in accordance with data privacy agreements between CHOP and the State of Ohio and Diagnostic Driving Inc. Thus, this study was considered exempt from Institutional Review Board oversight by CHOP. Our analytical sample consists of first-time license applicants in Ohio who were <25 years of age and took the VDA immediately before passing the state license examination, for which they were awarded their license, between July 2017 and December 2019. Ohio license and crash database records (up to March 15, 2020, and before the coronavirus pandemic lockdown) were exact-matched to VDA records by using a unique driver ID, which was then deidentified through an Honest Broker at CHOP. Through the deidentification process, all dates were converted to age in days from birthdate, including age at permit issuance, license examination attempt, license issuance, and any police-reported crash. From these, we derived the time since licensure and the time to first crash. The crash follow-up time for individuals was also rounded to 15-day intervals to protect the privacy of the data. The driver’s license addresses were geocoded into Federal Information Processing codes by using Esri geocoder and then matched to American Community Survey data28 to obtain Census tract-level sociodemographic variables.
Analytical Sample
Our analytical sample consisted of 16 914 newly licensed drivers who were <25 years of age, completed the VDA immediately before passing the license examination in Ohio, and were awarded their first license within a month of their license examination. See Fig 1 for sample derivation. These individuals were located at 5 licensing sites in Ohio in which the VDA was implemented (3 sites in the urban/suburban Columbus area, 1 other suburban site, and 1 rural site).
Data Sources
Ohio Licensing and Crash Record Databases
The Ohio licensing database contains detailed information on each driver’s interactions with the Ohio Bureau of Motor Vehicles (BMV), including driver demographics (date of birth, gender, address) and licensing history (permit dates, license examination dates, license examination outcomes, license issue dates, and license status). Of note, in this sample, the “gender” variable only includes “male” and “female” and, thus, hereafter is referred to as “sex.” Ohio crash records are collected from law enforcement agencies via a statewide uniform crash report, in compliance with current Model Minimum Uniform Crash Criteria standards.29 In Ohio, a crash must be reported if any personal injury or fatality occurs, and/or there is at least $1000 worth of property damage.
VDA Data
The VDA data were acquired from Diagnostic Driving Inc., which collected it via the VDA software (Ready-AssessTM) that was implemented in Ohio BMVs. The Ready-Assess VDA is a ∼15-minute-long self-directed workflow including a driving route that incorporates common serious crash risk scenarios,30,31 such as rear-end events, intersections, curved roads, merges, and hazard zones. These routes also have varied settings (urban and suburban), physical road features, and other road potential hazards (eg, crosswalks, merges, construction zones, vehicles, and pedestrians). The Ohio VDA dataset had 69 variables that capture operational and tactical driving skills across the drive, including aberrant or hazardous behaviors like simulated crashes and traffic light violations, as well as performance in several known domains of driving (see Table 1). Our team previously used a subset of these VDA metrics (see Table 1 and Walshe et al 2022)32 to derive the VDA Skill Clusters and 4 major Driving Classes: (1) No Issues (ie, careful and skilled drivers showing no issues on the VDA), (2) Minor Issues (ie, an average new driver with minor vehicle control skill deficits), (3) Major Issues (ie, drivers with more control issues and who take more risks), and (4) Major Issues with Dangerous Behavior (ie, drivers with even more control issues with “jackrabbit,” sudden, jerky movements, and more reckless and risk-taking behavior). See Table 2 for Skill Cluster and Driving Class labels and their weight on domains of VDA skill performance that differentiate them (for more specific details on Skill Cluster weightings on each variable (see Supplemental Appendix A plots in Walsche et al32 ). Of note, the original paper labeled the worst performing Driving Class as “Major Issues with Aggression,” which has since been relabeled as “Major Issues with Dangerous Behavior” to be more descriptive and avoid implying underlying personality traits.
VDA Skill Domain . | VDA Variables . | Definition (Unit) . |
---|---|---|
Speed-related | Max Speed | Maximum speed value (mph) |
Mean Speed | Mean speed value (mph) | |
SD Speed | SD of speed (mph) | |
Max Throttle | Maximum accelerator depression (%) | |
Mean Throttle | Mean accelerator depression (%) | |
SD Throttle | SD of accelerator depression (%) | |
Drive Duration | Time to complete assessment drive (seconds) | |
Throttle control | Max Acceleration | Maximum vehicle acceleration (mph/s) |
Mean Acceleration | Mean vehicle acceleration (mph/s) | |
SD Acceleration | SD of vehicle acceleration (mph/s) | |
Max Jerk | Maximum vehicle jerk (mph/s2) | |
Mean Jerk | Mean vehicle jerk (mph/s2) | |
SD Jerk | SD of vehicle jerk (mph/s2) | |
Braking control | Max Brake | Maximum brake depression (%) |
Mean Brake | Mean brake depression (%) | |
SD Brake | SD of brake depression (%) | |
Lane position | Max Heading-Error | Maximum angle between vehicle’s heading vector and road-following vector (degrees [0, 180]) |
Mean Heading-Error | Mean angle between vehicle’s heading vector and road-following vector (degrees [0, 180]) | |
SD Heading-Error | SD of angle between vehicle’s heading vector and road-following vector (degrees [0, 180]) | |
Max Lane Deviation | Maximum vehicle lateral displacement from the center of lane (meters) | |
Mean Lane Deviation | Mean vehicle lateral displacement from the center of lane (meters) | |
SD Lane Deviation | SD of vehicle lateral displacement from the center of lane (meters) | |
Max Road-Center Deviation | Maximum distance from vehicle to center of the road (meters) | |
Mean Road-Center Deviation | Mean distance from vehicle to center of the road (meters) | |
SD Road-Center Deviation | SD of distance from vehicle to center of the road (meters) | |
Route following | Off-Route | Driving off-road or off-route (count of incidences) |
Car following | Time to Collision <3 s | Time spent driving at <3 s to crash (seconds) |
Miles Driven TTC <3 s | Distance driven while <3 s to crash (miles) | |
Time to Collision 3–5 s | Time spent driving at 5–3 s to crash (seconds) | |
Miles Driven TTC 3–5 s | Distance driven at 5–3 s to crash (miles) | |
Rule following | Failures to Stop | Failures to stop at stop signs/red lights (count) |
Crash avoidance | Crashes | Crashes with vehicle, pedestrian or object (count) |
VDA Skill Domain . | VDA Variables . | Definition (Unit) . |
---|---|---|
Speed-related | Max Speed | Maximum speed value (mph) |
Mean Speed | Mean speed value (mph) | |
SD Speed | SD of speed (mph) | |
Max Throttle | Maximum accelerator depression (%) | |
Mean Throttle | Mean accelerator depression (%) | |
SD Throttle | SD of accelerator depression (%) | |
Drive Duration | Time to complete assessment drive (seconds) | |
Throttle control | Max Acceleration | Maximum vehicle acceleration (mph/s) |
Mean Acceleration | Mean vehicle acceleration (mph/s) | |
SD Acceleration | SD of vehicle acceleration (mph/s) | |
Max Jerk | Maximum vehicle jerk (mph/s2) | |
Mean Jerk | Mean vehicle jerk (mph/s2) | |
SD Jerk | SD of vehicle jerk (mph/s2) | |
Braking control | Max Brake | Maximum brake depression (%) |
Mean Brake | Mean brake depression (%) | |
SD Brake | SD of brake depression (%) | |
Lane position | Max Heading-Error | Maximum angle between vehicle’s heading vector and road-following vector (degrees [0, 180]) |
Mean Heading-Error | Mean angle between vehicle’s heading vector and road-following vector (degrees [0, 180]) | |
SD Heading-Error | SD of angle between vehicle’s heading vector and road-following vector (degrees [0, 180]) | |
Max Lane Deviation | Maximum vehicle lateral displacement from the center of lane (meters) | |
Mean Lane Deviation | Mean vehicle lateral displacement from the center of lane (meters) | |
SD Lane Deviation | SD of vehicle lateral displacement from the center of lane (meters) | |
Max Road-Center Deviation | Maximum distance from vehicle to center of the road (meters) | |
Mean Road-Center Deviation | Mean distance from vehicle to center of the road (meters) | |
SD Road-Center Deviation | SD of distance from vehicle to center of the road (meters) | |
Route following | Off-Route | Driving off-road or off-route (count of incidences) |
Car following | Time to Collision <3 s | Time spent driving at <3 s to crash (seconds) |
Miles Driven TTC <3 s | Distance driven while <3 s to crash (miles) | |
Time to Collision 3–5 s | Time spent driving at 5–3 s to crash (seconds) | |
Miles Driven TTC 3–5 s | Distance driven at 5–3 s to crash (miles) | |
Rule following | Failures to Stop | Failures to stop at stop signs/red lights (count) |
Crash avoidance | Crashes | Crashes with vehicle, pedestrian or object (count) |
SD, standard deviation.
. | . | Virtual Driving Assessment Skill Domain . | |||||||
---|---|---|---|---|---|---|---|---|---|
Driving classes, skill clusters . | Sample % . | Speed control . | Throttle control . | Brake control . | Lane position . | Route following . | Car following . | Rule following . | Crash avoidance . |
No issues | 28.3% | + | + | + | + | + | |||
1. Good steering control | 9.7% | + | |||||||
2. Cautious | 9.2% | + | + | + | + | ||||
3. Good brake & steering | 9.3% | + | + | ||||||
Minor issues | 39.0% | − | − | − | + | ||||
4. Skilled with hard throttle | 10.1% | − | + | ||||||
5. Jerky braking | 4.9% | − | − | + | |||||
6. Speeder, tailgater, rule breaker | 6.4% | − | − | − | |||||
7. Skilled average | 5.3% | + | |||||||
8. Below average control | 3.8% | − | − | ||||||
9. Quick with controlled braking | 4.4% | − | − | − | |||||
10. Mild control issues | 4.2% | + | − | ||||||
Major issues | 15.2% | + | + | + | − | − | + | − | |
11. Skilled rule breakers | 8.1% | + | + | − | |||||
12. Slow, poor control, rule breakers | 2.1% | + | + | − | + | − | |||
13. Less control, rule breakers | 3.4% | − | + | − | − | ||||
14. Extremely slow, poor control, rule breaker | 1.6% | + | + | + | − | − | + | − | |
Major issues with dangerous behavior | 17.5% | − | − | − | − | − | − | − | − |
15. Aggressive tailgaters | 4.1% | − | − | − | − | ||||
16. Extremely aggressive, reckless | 4.1% | − | − | − | − | − | − | ||
17. Controlled jackrabbit | 4.1% | − | − | − | |||||
18. Less controlled jackrabbit | 2.5% | − | − | − | |||||
19. Risky, poor control, jackrabbit | 2.0% | − | − | − | − | − | − | − | |
20. Risky, no control, jackrabbit* | 0.8% | − | − | − | − | + | − | − |
. | . | Virtual Driving Assessment Skill Domain . | |||||||
---|---|---|---|---|---|---|---|---|---|
Driving classes, skill clusters . | Sample % . | Speed control . | Throttle control . | Brake control . | Lane position . | Route following . | Car following . | Rule following . | Crash avoidance . |
No issues | 28.3% | + | + | + | + | + | |||
1. Good steering control | 9.7% | + | |||||||
2. Cautious | 9.2% | + | + | + | + | ||||
3. Good brake & steering | 9.3% | + | + | ||||||
Minor issues | 39.0% | − | − | − | + | ||||
4. Skilled with hard throttle | 10.1% | − | + | ||||||
5. Jerky braking | 4.9% | − | − | + | |||||
6. Speeder, tailgater, rule breaker | 6.4% | − | − | − | |||||
7. Skilled average | 5.3% | + | |||||||
8. Below average control | 3.8% | − | − | ||||||
9. Quick with controlled braking | 4.4% | − | − | − | |||||
10. Mild control issues | 4.2% | + | − | ||||||
Major issues | 15.2% | + | + | + | − | − | + | − | |
11. Skilled rule breakers | 8.1% | + | + | − | |||||
12. Slow, poor control, rule breakers | 2.1% | + | + | − | + | − | |||
13. Less control, rule breakers | 3.4% | − | + | − | − | ||||
14. Extremely slow, poor control, rule breaker | 1.6% | + | + | + | − | − | + | − | |
Major issues with dangerous behavior | 17.5% | − | − | − | − | − | − | − | − |
15. Aggressive tailgaters | 4.1% | − | − | − | − | ||||
16. Extremely aggressive, reckless | 4.1% | − | − | − | − | − | − | ||
17. Controlled jackrabbit | 4.1% | − | − | − | |||||
18. Less controlled jackrabbit | 2.5% | − | − | − | |||||
19. Risky, poor control, jackrabbit | 2.0% | − | − | − | − | − | − | − | |
20. Risky, no control, jackrabbit* | 0.8% | − | − | − | − | + | − | − |
The + and – indicate positive or negative loading on each skill domain (for example, a “+” on speed control indicates good speed control, and a “–” indicates poor speed control). Weight strengths are not indicated here but please see Supplemental Appendix A plots in Walsche et al for more detailed cluster loading information.
The risky, no control, jackrabbit skill cluster is distinguished from risky, poor control, jackrabbit, by higher negative loadings on throttle control and lane position.
Outcome Variables
The outcome is time to first crash from the first day of licensure, determined as the total number of days after a license was issued to the first involvement in a crash. Although measured to the day, total follow-up time was rounded to 15-day intervals to protect driver privacy.
Predictor Variables
The key predictor in this crash risk model is the 20 Skill Clusters and the 4 summary Driving Classes (see Table 2).
Adjustment Variables
The potentially confounding adjustment variables included in our model include the 5 licensing sites noted above, age at the time of licensure (16 to <17, 17 to <18, 18 to <19, and 19 to <25 years of age), sex (male/female), as well as sociodemographic summary measures of each county (as described previously32,33 ): income/education, transportation, and urbanicity.
Statistical Analysis
A Cox proportional hazard model was used to estimate the instantaneous risk of a police-reported crash during the follow-up period as a function of Skill Cluster and Driving Class membership. Effect coding was also used to estimate the change in risk associated with cluster assignment so that each category is compared with the overall average crash risk. All analyses are adjusted for licensing center location at which the test was taken; we treat these as “unadjusted” (that is, unadjusted for substantive confounders of interest). Analyses adjusting for individual-level age, sex and census-tract-level sociodemographic characteristics were also conducted. To account for potential correlations among drivers at each center, a generalized estimating equation approach was used.34 All analyses were completed by using SAS V9.4. The proportional hazard assumption was tested by using a bootstrap test of cumulative martingale residuals.35 For models that failed to follow the proportional hazard assumption, interactions with time were considered to assess how risk might have changed over follow-up within certain classes. Interaction with the age of licensure among the Driving Classes was also considered.
Results
Sample Characteristics
The majority (78%) of newly licensed drivers in this sample were tested in the Columbus area (Urban Sites 1, 2; Suburban Site 3), 14% were at a rural location (Site 5), and 8.4% were at another suburban location (Site 4; Table 3).
Variable . | Site 1 . | Site 2 . | Site 3 . | Site 4 . | Site 5 . |
---|---|---|---|---|---|
All . | n = 3929 (23%) . | n = 2786 (17%) . | n = 6412 (38%) . | n = 1422 (8%) . | n = 2365 (14%) . |
Age, n (%) | |||||
16 y | 1303 (33.2%) | 985 (35.4%) | 3163 (49.3%) | 842 (59.2%) | 1522 (64.4%) |
17 y | 595 (15.1%) | 392 (14.1%) | 923 (14.4%) | 201 (14.1%) | 265 (11.2%) |
18 y | 914 (23.3%) | 542 (19.5%) | 925 (14.4%) | 175 (12.3%) | 297 (12.6%) |
19–25 y | 1117 (28.4%) | 867 (31.1%) | 1401 (21.9%) | 204 (14.4%) | 281 (11.9%) |
Sex: male, n (%) | 1935 (49.3%) | 1492 (53.6%) | 3406 (53.1%) | 744 (52.3%) | 1227 (51.9%) |
Tract-level median household income, n (%) | |||||
Low: bottom 10th percentile | 566 (14.4%) | 442 (15.9%) | 577 (9.0%) | 144 (10.1%) | 339 (14.3%) |
10th–90th percentile: | 2975 (75.7%) | 1922 (69.0%) | 4709 (73.5%) | 1265 (89.0%) | 1986 (84.0%) |
High: top 10th percentile | 388 (9.9%) | 422 (15.1%) | 1125 (17.5%) | 13 (0.9%) | 40 (1.7%) |
Tract-level college degree education, n (%) | |||||
Low: bottom 10th percentile | 570 (14.5%) | 219 (7.9%) | 509 (7.9%) | 281 (19.7%) | 578 (24.5%) |
10th–90th percentile: | 2937 (74.8%) | 1855 (66.6%) | 4019 (62.7%) | 1126 (79.2%) | 1765 (74.6%) |
High: top 10th percentile | 422 (10.7%) | 712 (25.5%) | 1883 (29.4%) | 15 (1.1%) | 22 (0.9%) |
Variable . | Site 1 . | Site 2 . | Site 3 . | Site 4 . | Site 5 . |
---|---|---|---|---|---|
All . | n = 3929 (23%) . | n = 2786 (17%) . | n = 6412 (38%) . | n = 1422 (8%) . | n = 2365 (14%) . |
Age, n (%) | |||||
16 y | 1303 (33.2%) | 985 (35.4%) | 3163 (49.3%) | 842 (59.2%) | 1522 (64.4%) |
17 y | 595 (15.1%) | 392 (14.1%) | 923 (14.4%) | 201 (14.1%) | 265 (11.2%) |
18 y | 914 (23.3%) | 542 (19.5%) | 925 (14.4%) | 175 (12.3%) | 297 (12.6%) |
19–25 y | 1117 (28.4%) | 867 (31.1%) | 1401 (21.9%) | 204 (14.4%) | 281 (11.9%) |
Sex: male, n (%) | 1935 (49.3%) | 1492 (53.6%) | 3406 (53.1%) | 744 (52.3%) | 1227 (51.9%) |
Tract-level median household income, n (%) | |||||
Low: bottom 10th percentile | 566 (14.4%) | 442 (15.9%) | 577 (9.0%) | 144 (10.1%) | 339 (14.3%) |
10th–90th percentile: | 2975 (75.7%) | 1922 (69.0%) | 4709 (73.5%) | 1265 (89.0%) | 1986 (84.0%) |
High: top 10th percentile | 388 (9.9%) | 422 (15.1%) | 1125 (17.5%) | 13 (0.9%) | 40 (1.7%) |
Tract-level college degree education, n (%) | |||||
Low: bottom 10th percentile | 570 (14.5%) | 219 (7.9%) | 509 (7.9%) | 281 (19.7%) | 578 (24.5%) |
10th–90th percentile: | 2937 (74.8%) | 1855 (66.6%) | 4019 (62.7%) | 1126 (79.2%) | 1765 (74.6%) |
High: top 10th percentile | 422 (10.7%) | 712 (25.5%) | 1883 (29.4%) | 15 (1.1%) | 22 (0.9%) |
Tract-level percentiles are based on Ohio state-wide data (see previous paper35 ).
Table 4 presents sample characteristics. Overall, almost one-half of the newly licensed drivers were aged 16 years (46%), with fewer 19- to 25-year-old new drivers (22%). Relative to state-wide distributions of license applicants,18 this sample was from more urban tracts and was slightly more educated (and more so among drivers with no crash).
Variable . | All Drivers . | Drivers With No Crash . | Drivers With Crash . |
---|---|---|---|
All . | n = 16 914 (100%) . | n = 14 613 (100%) . | n = 2301 (100%) . |
Age, n (%) | |||
16 y | 7815 (46.2%) | 6876 (47.1%) | 939 (40.8%) |
17 y | 2376 (14.0%) | 1997 (13.7%) | 379 (16.5%) |
18 y | 2853 (16.9%) | 2388 (16.3%) | 465 (20.2%) |
19–24 y | 3870 (22.9%) | 3352 (22.9%) | 518 (22.5%) |
Sex: male, n (%) | 8804 (52.1%) | 7586 (51.9%) | 1218 (52.9%) |
Mean follow-up time, d (SD) | 465 (228) | 452 (227) | 548 (215) |
Time in learner permit | |||
Mean number, d (SD) | 190 (109) | 191 (108) | 183(111) |
TLP <14 d, n (%) | 1332 (7.9%) | 1121 (7.7%) | 211 (9.2%) |
TLP 14 d to <6 mo, n (%) | 4681 (27.7%) | 4004 (27.4%) | 677 (29.4%) |
TLP 6 to <12 mo, n (%) | 10901 (64.4%) | 9488 (64.9%) | 1413 (61.4%) |
Tract-level median household income, n (%) | |||
Low: bottom 10th percentile | 2068 (12.2%) | 1751 (12.0%) | 317(13.8%) |
Mid: 10th–90th percentile | 12857 (76.0%) | 11067 (75.7%) | 1790 (77.8%) |
High: top 10th percentile | 1988 (11.8%) | 1794 (12.3%) | 194 (8.4%) |
Tract-level college degree education, n (%) | |||
Low: bottom 10th percentile | 2157 (12.7%) | 1813 (12.4%) | 344 (14.9%) |
Mid: 10th–90th percentile | 11702 (69.2%) | 10048 (68.8%) | 1654 (71.9%) |
High: top 10th percentile | 3054 (18.1%) | 2751 (18.8%) | 303 (13.2%) |
Neighborhood urbanicity, n (%) | |||
Urban | 3901 (23.1%) | 3300 (22.6%) | 601 (26.1%) |
Suburban | 10561 (62.4%) | 9158 (62.7%) | 1403 (61.0%) |
Rural | 2452 (14.5%) | 2155 (14.7%) | 297 (12.9%) |
Variable . | All Drivers . | Drivers With No Crash . | Drivers With Crash . |
---|---|---|---|
All . | n = 16 914 (100%) . | n = 14 613 (100%) . | n = 2301 (100%) . |
Age, n (%) | |||
16 y | 7815 (46.2%) | 6876 (47.1%) | 939 (40.8%) |
17 y | 2376 (14.0%) | 1997 (13.7%) | 379 (16.5%) |
18 y | 2853 (16.9%) | 2388 (16.3%) | 465 (20.2%) |
19–24 y | 3870 (22.9%) | 3352 (22.9%) | 518 (22.5%) |
Sex: male, n (%) | 8804 (52.1%) | 7586 (51.9%) | 1218 (52.9%) |
Mean follow-up time, d (SD) | 465 (228) | 452 (227) | 548 (215) |
Time in learner permit | |||
Mean number, d (SD) | 190 (109) | 191 (108) | 183(111) |
TLP <14 d, n (%) | 1332 (7.9%) | 1121 (7.7%) | 211 (9.2%) |
TLP 14 d to <6 mo, n (%) | 4681 (27.7%) | 4004 (27.4%) | 677 (29.4%) |
TLP 6 to <12 mo, n (%) | 10901 (64.4%) | 9488 (64.9%) | 1413 (61.4%) |
Tract-level median household income, n (%) | |||
Low: bottom 10th percentile | 2068 (12.2%) | 1751 (12.0%) | 317(13.8%) |
Mid: 10th–90th percentile | 12857 (76.0%) | 11067 (75.7%) | 1790 (77.8%) |
High: top 10th percentile | 1988 (11.8%) | 1794 (12.3%) | 194 (8.4%) |
Tract-level college degree education, n (%) | |||
Low: bottom 10th percentile | 2157 (12.7%) | 1813 (12.4%) | 344 (14.9%) |
Mid: 10th–90th percentile | 11702 (69.2%) | 10048 (68.8%) | 1654 (71.9%) |
High: top 10th percentile | 3054 (18.1%) | 2751 (18.8%) | 303 (13.2%) |
Neighborhood urbanicity, n (%) | |||
Urban | 3901 (23.1%) | 3300 (22.6%) | 601 (26.1%) |
Suburban | 10561 (62.4%) | 9158 (62.7%) | 1403 (61.0%) |
Rural | 2452 (14.5%) | 2155 (14.7%) | 297 (12.9%) |
SD, standard deviation; TLP, time in learner permit.
Tract-level thresholds were derived by using state-wide data.33
Distribution of VDA Driving Classes and Skill Clusters
Table 2 reveals the distribution of the sample across the 4 Driving Classes and their subsets of Skill Clusters.
Crashes
In this newly licensed driver sample, 2 301 (13.6%) were observed to have at least 1 crash. Figure 2 reveals the Kaplan-Meier curve. In addition, 1.5% of these new drivers had 2 or more crashes (10.7% of those who had at least 1 crash) during follow-up. The maximum total number of crashes for an individual was 4. Table 5 reveals the mean number of crashes per year by VDA Driving Classes and the various Skill Clusters they represent. However, only time to first crash is used in our estimation of how the risk of crash varies by VDA classification.
. | . | . | Survival Model Hazard Ratios (CI 95%) . | ||
---|---|---|---|---|---|
Driving classes, skill clusters . | Mean length of follow-up (d) . | Mean number of crashes per year . | Adjusted for location . | Adjusted for location, age, sex . | Adjusted for location, age, sex, SES . |
No issues | 425 | 0.112 | 0.90 (0.87–0.94)* | 0.91 (0.87–0.95)* | 0.91 (0.87–0.95)* |
1. Good steering control | 419 | 0.112 | 0.88 (0.77–1.00) | 0.89 (0.77–1.02) | 0.89 (0.77–1.04) |
2. Cautious | 409 | 0.108 | 0.83 (0.72–0.96)* | 0.83 (0.72–0.96)* | 0.83 (0.73–0.95)* |
3. Good brake & steering | 446 | 0.115 | 0.91 (0.80–1.03) | 0.91 (0.81–1.03) | 0.91 (0.80–1.02) |
Minor issues | 518 | 0.119 | 1.04 (0.97–1.12) | 1.04 (0.97–1.12) | 1.04 (0.97–1.12) |
4. Skilled with hard throttle | 613 | 0.116 | 0.99 (0.94–1.04) | 0.99 (0.93–1.05) | 0.99 (0.94–1.05) |
5. Jerky braking | 615 | 0.116 | 1.03 (0.93–1.13) | 1.02 (0.92–1.13) | 1.03 (0.93–1.13) |
6. Speeder, tailgater, rule breaker | 506 | 0.121 | 1.07 1.02–1.14)* | 1.07 (1.00–1.13)* | 1.06 (1.01–1.12)* |
7. Skilled average | 410 | 0.118 | 0.93 (0.74–1.17) | 0.93 (0.75–1.16) | 0.93 (0.76–1.14) |
8. Below average control | 438 | 0.122 | 1.03 (0.73–1.46) | 1.03 (0.72–1.46) | 1.05 (0.74–1.48) |
9. Quick with controlled braking | 476 | 0.117 | 0.95 (0.75–1.19) | 0.94 (0.75–1.19) | 0.95 (0.76–1.19) |
10. Mild control issues | 453 | 0.132 | 1.11 (0.96–1.29) | 1.10 (0.96–1.26) | 1.07 (0.93–1.23) |
Major issues | 399 | 0.118 | 0.96 (0.88–1.04) | 0.96 (0.88–1.04) | 0.96 (0.88–1.04) |
11. Skilled rule breakers | 390 | 0.116 | 0.89 (0.71–1.11) | 0.88 (0.70–1.11) | 0.89 (0.70–1.13) |
12. Slow, poor control, rule breakers | 450 | 0.106 | 0.85 (0.61–1.19) | 0.86 (0.62–1.18) | 0.87 (0.64–1.19) |
13. Less control, rule breakers | 389 | 0.127 | 1.04 (0.95–1.15) | 1.03 (0.94–1.14) | 1.02 (0.92–1.14) |
14. Extremely slow, poor control, rule breaker | 398 | 0.131 | 1.01 (0.84–1.20) | 1.02 (0.85–1.22) | 1.00 (0.83–1.21) |
Major issues with dangerous behavior | 471 | 0.132 | 1.11 (1.01–1.22)* | 1.11 (1.00–1.23)* | 1.11 (1.00–1.23)* |
15. Aggressive tailgaters | 475 | 0.125 | 1.04 (0.79–1.38) | 1.04 (0.78–1.37) | 1.00 (0.83–1.21) |
16. Extremely aggressive, reckless | 484 | 0.139 | 1.08 (1.01–1.16)* | 1.07 (1.01–1.15)* | 1.10 (1.01–1.19)* |
17. Controlled jackrabbit | 448 | 0.130 | 1.05 (0.94–1.17) | 1.05 (0.95–1.15) | 1.05 (0.96–1.15) |
18. Less controlled jackrabbit | 495 | 0.125 | 1.08 (0.84–1.39) | 1.09 (0.84–1.41) | 1.08 (0.83–1.39) |
19. Risky, poor control, jackrabbit | 463 | 0.142 | 1.11 (0.95–1.29) | 1.11 (0.97–1.28) | 1.12 (0.94–1.33) |
20. Risky, no control, jackrabbit | 441 | 0.130 | 1.20 (1.02–1.42)* | 1.24 (1.06–1.44)* | 1.24 (1.00–1.46)* |
. | . | . | Survival Model Hazard Ratios (CI 95%) . | ||
---|---|---|---|---|---|
Driving classes, skill clusters . | Mean length of follow-up (d) . | Mean number of crashes per year . | Adjusted for location . | Adjusted for location, age, sex . | Adjusted for location, age, sex, SES . |
No issues | 425 | 0.112 | 0.90 (0.87–0.94)* | 0.91 (0.87–0.95)* | 0.91 (0.87–0.95)* |
1. Good steering control | 419 | 0.112 | 0.88 (0.77–1.00) | 0.89 (0.77–1.02) | 0.89 (0.77–1.04) |
2. Cautious | 409 | 0.108 | 0.83 (0.72–0.96)* | 0.83 (0.72–0.96)* | 0.83 (0.73–0.95)* |
3. Good brake & steering | 446 | 0.115 | 0.91 (0.80–1.03) | 0.91 (0.81–1.03) | 0.91 (0.80–1.02) |
Minor issues | 518 | 0.119 | 1.04 (0.97–1.12) | 1.04 (0.97–1.12) | 1.04 (0.97–1.12) |
4. Skilled with hard throttle | 613 | 0.116 | 0.99 (0.94–1.04) | 0.99 (0.93–1.05) | 0.99 (0.94–1.05) |
5. Jerky braking | 615 | 0.116 | 1.03 (0.93–1.13) | 1.02 (0.92–1.13) | 1.03 (0.93–1.13) |
6. Speeder, tailgater, rule breaker | 506 | 0.121 | 1.07 1.02–1.14)* | 1.07 (1.00–1.13)* | 1.06 (1.01–1.12)* |
7. Skilled average | 410 | 0.118 | 0.93 (0.74–1.17) | 0.93 (0.75–1.16) | 0.93 (0.76–1.14) |
8. Below average control | 438 | 0.122 | 1.03 (0.73–1.46) | 1.03 (0.72–1.46) | 1.05 (0.74–1.48) |
9. Quick with controlled braking | 476 | 0.117 | 0.95 (0.75–1.19) | 0.94 (0.75–1.19) | 0.95 (0.76–1.19) |
10. Mild control issues | 453 | 0.132 | 1.11 (0.96–1.29) | 1.10 (0.96–1.26) | 1.07 (0.93–1.23) |
Major issues | 399 | 0.118 | 0.96 (0.88–1.04) | 0.96 (0.88–1.04) | 0.96 (0.88–1.04) |
11. Skilled rule breakers | 390 | 0.116 | 0.89 (0.71–1.11) | 0.88 (0.70–1.11) | 0.89 (0.70–1.13) |
12. Slow, poor control, rule breakers | 450 | 0.106 | 0.85 (0.61–1.19) | 0.86 (0.62–1.18) | 0.87 (0.64–1.19) |
13. Less control, rule breakers | 389 | 0.127 | 1.04 (0.95–1.15) | 1.03 (0.94–1.14) | 1.02 (0.92–1.14) |
14. Extremely slow, poor control, rule breaker | 398 | 0.131 | 1.01 (0.84–1.20) | 1.02 (0.85–1.22) | 1.00 (0.83–1.21) |
Major issues with dangerous behavior | 471 | 0.132 | 1.11 (1.01–1.22)* | 1.11 (1.00–1.23)* | 1.11 (1.00–1.23)* |
15. Aggressive tailgaters | 475 | 0.125 | 1.04 (0.79–1.38) | 1.04 (0.78–1.37) | 1.00 (0.83–1.21) |
16. Extremely aggressive, reckless | 484 | 0.139 | 1.08 (1.01–1.16)* | 1.07 (1.01–1.15)* | 1.10 (1.01–1.19)* |
17. Controlled jackrabbit | 448 | 0.130 | 1.05 (0.94–1.17) | 1.05 (0.95–1.15) | 1.05 (0.96–1.15) |
18. Less controlled jackrabbit | 495 | 0.125 | 1.08 (0.84–1.39) | 1.09 (0.84–1.41) | 1.08 (0.83–1.39) |
19. Risky, poor control, jackrabbit | 463 | 0.142 | 1.11 (0.95–1.29) | 1.11 (0.97–1.28) | 1.12 (0.94–1.33) |
20. Risky, no control, jackrabbit | 441 | 0.130 | 1.20 (1.02–1.42)* | 1.24 (1.06–1.44)* | 1.24 (1.00–1.46)* |
SES, socioeconomic status.
Statistical significance at the α = 0.05 level.
Survival Model Results
Table 5 provides the hazard ratio estimates for each of the Skill Clusters and Driving Classes, relative to the overall average risk of crash. After adjusting only for licensing center location, the No Issues Driving Class had a significantly lower risk of crash, at 10% less than average (95% confidence interval [CI] 13% to 6%), whereas the Major Issues with Dangerous Behavior Class had a significantly higher risk of crash, 11% greater than average (95% CI 1% to 22%). In addition, all Skill Clusters in the No Issues Class had below-average risks of a crash, with the Cautious skill cluster drivers reaching statistical significance, with a crash risk 17% less than the average (95% CI 28% to 4%). On the other hand, all Skill Clusters in the Major Issues with Dangerous Behavior Class had above-average risk of crash, with 2 reaching statistical significance: the Risky, No Control, Jackrabbit drivers had a 20% increase in risk (95% CI 2% to 42%), and Extremely Aggressive, Reckless drivers had an 8% increase in risk (95% CI 1% to 16%). The Minor Issues and Major Issues Classes did not differ from the average crash risk. The Skill Clusters in the Minor Issues and Major Issues Classes were both above and below-average risk, although most were close to 1. The exception to this is the Speeder, Tailgater, Rulebreaker Skill Cluster in the Minor Issues Class, in which there was a 7% increase in risk relative to the overall risk (95% CI 2% to 14%). Adjusting for covariates (age, sex, and tract-level SES indicators) had little impact on crash hazard ratios.
Interactions With Licensing Age
There was evidence of an interaction between age and the Minor Issues Class (P = .005); in particular, drivers licensed at age 18 had an elevated crash risk 16% higher than average (95% CI 6% to 27%), but no other age group differed from average crash risk. Specifically, the Skilled Average Cluster within the Minor Issues Class revealed a strong interaction with the 18-year-old age group (who are not required to complete driver education and training before licensure): the risk of crash was 38% higher than average in this age group (95% CI 13% to 68%).
Nonconstant Hazard Ratios
The proportional hazards assumption, that the relative risk difference is constant over time, generally held for the Driving Classes, but this was not the case for 3 of the Skill Clusters: Below Average Control (P = .025), Quick with Controlled Braking (P = .045) and Risky, Poor Control, Jackrabbit (P = .005; See Fig 3). The Below Average Control Cluster had significantly elevated risk in 12 to 18 months after licensing (47% increase, 95% CI 7% to 101%). The Quick with Controlled Braking Cluster started out with near-average risk but, after 18 months of follow-up, had above-average risk (74% increase, 95% CI 35% to 126%). A similar pattern was seen for the Risky, Poor Control, Jackrabbit Cluster (112% increase at 12 to 18 months, 95% CI 64% to 174%; 53% increase after 18 months, 95% CI 11% to 162%).
Discussion
This population-level study confirmed that new driver skills at licensure (measured via a virtual driving assessment) were a predictor of time to crash early in licensure.
Specifically, compared with average risk, the best-performing No Issues Driving Class had a lower risk of crashing post-licensure, and the worst-performing Major Issues with Dangerous Behavior Driving Class was associated with higher crash risk. The more moderate Minor Issues Driving Class and the Major Issues Class did not have statistically significant differences from the overall risk of crash. However, within the Minor Issues Driving Class, the Speeder, Tailgater, Rule Breaker cluster had a significantly elevated risk of crash. Although this cluster had a small proportion of drivers (0.8%), it is a worrisome behavior pattern that could be targeted for interventions if identified at the time of licensure via the VDA. Given these findings, this Skill Cluster would likely be better placed in the Major Issues with Dangerous Behavior Driving Class and will be reclassified in subsequent studies and use of this tool.
Controlling for the covariates of age, sex, and tract-level sociodemographic indicators did not add any predictive power. However, there was an interaction between age and the Minor Issues Driving Class, whereby those aged 18 were at an elevated risk of a crash (especially prominent in the Skilled Average cluster). Of note, in Ohio, license applicants aged 18 years are only required to pass the licensing examination and are not subject to the rigorous licensing policies of those <18 years of age (a 6-month learner permit, mandated driver education and behind-the-wheel training). Future studies should explore the interaction between training, exposure and skills, and time-varying risk post-licensure.
Although the VDA does not capture an individual’s entire risk profile, it could be used in combination with other measures of crash-risk skills, attitudes, and behaviors to advance the development of better diagnostic prediction models of crash risk. With this information, programs can be developed that provide personalized and targeted interventions that mitigate skill deficits before licensure and in the early period of licensure when crash risk is highest for new drivers. For example, some online education interventions have been developed for improving hazard awareness and prevention that have been validated by observed reduction in crashes25,36 that could be provided at the time of licensure.
Limitations
One important limitation that we could not control for is miles driven, future studies should account for differences in driving exposure. In addition, the total follow-up time was rounded to 15-day intervals to protect data privacy, meaning the time to crash could have errors of up to 14 days. However, this error would be expected to be random and, thus, not likely to impact these results. Although the demographic profile of this analytical sample from 5 Ohio BMVs is rather similar to that of the state-wide license applicant population of Ohio in the same years (per our previous work33 ), this sample was skewed slightly older with slightly fewer 16-year-old newly licensed drivers (46% vs 54% state-wide), and more 19- to 25-year-old newly licensed drivers (23% vs 18% state-wide). Finally, the VDA database available to us from the Ohio BMV implementation did not contain visual-scanning metrics that could better indicate situational awareness, a key component for driver safety.
Conclusions
This population-level study reveals that new driver skills measured by the VDA at the time of licensure contribute to crash risk post-licensure. These findings take us a step closer to identifying the drivers at risk of crashing early in licensure and some of the skills that are protective of future crashes. These results can inform targeted interventions that could be delivered at the time of licensure, providing actionable feedback or training to newly licensed drivers on areas needing continued improvement post-licensure.
Acknowledgments
The authors would like to acknowledge Arcus and the Data Science and Biostatistics Unit at CHOP for their support in data management and geocoding the license address data.
Drs Walshe, Romer, and Winston conceptualized and designed the study, planned, oversaw and reviewed analysis, drafted the initial manuscript, and critically reviewed and revised the manuscript; Dr Elliot conceptualized the statistical analysis plan, executed the initial analyses, drafted the initial manuscript, and critically reviewed and revised the manuscript; Mr Cheng prepared the data set and assisted with conducting the analyses and drafting the initial manuscript; Dr Curry consulted on the study design and the analytical plan, reviewed the analyses conducted, and critically reviewed and revised the manuscript; Mr Grethlein prepared the data, reviewed the analysis conducted, assisted with drafting the manuscript, and reviewed the final manuscript; Mr Gonzalez prepared and honest brokered the research dataset and assisted with drafting the manuscript; and all authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.
FUNDING: Funded by the National Institutes of Health (NIH). Research reported in this publication was supported by the Eunice Kennedy Shriver Institute of Child Health and Human Development of the NIH under award number 1R21HD099635-01 (NICHD), by the US National Highway Traffic Safety Administration through the Ohio Traffic Safety Office (NHTSA), and by the State of Ohio’s Department of Administrative Services, managed by the Ohio Bureau of Motor Vehicles (Ohio), as well as a gift from New Jersey Manufacturers (NJM). Data collection was financed by NHTSA and Ohio. Data management and analyses were supported jointly by NICHD, Ohio, NHTSA, and the Annenberg Public Policy Center at the University of Pennsylvania. The preparation of this manuscript was supported jointly by NICHD, NHTSA, the Annenberg Public Policy Center, and NJM.
CONFLICT OF INTEREST DISCLOSURES: Flaura Winston and David Grethlein have an intellectual property and financial interest in Diagnostic Driving, Inc. The Children’s Hospital of Philadelphia (CHOP) has an institutional interest in Diagnostic Driving, Inc. Diagnostic Driving, Inc., created a virtual driving assessment system that is used in Ohio as an assessment at licensing centers and in driving schools to assess driver training programs. Flaura Winston serves as the chief scientific advisor of Diagnostic Driving, Inc. This potential conflict of interest is managed under a conflict-of-interest management plan from CHOP and the University of Pennsylvania whereby Flaura Winston has no interaction with participants (all field data collection procedures were conducted by Ohio Bureau of Motor Vehicles personnel) and all analyses were reviewed and approved by outside consultants with no intellectual or financial interest (John Bolte, a traffic injury researcher at the Ohio State University, and Nancy Kassam-Adams, a behavioral researcher at CHOP and the University of Pennsylvania). David Grethlein works as a Data Scientist at Diagnostic Driving Inc. His conflict is managed in the same manner as Dr Winston: David Grethlein has no interaction with participants and was not directly involved in the analysis of this paper, which has been reviewed and approved by outside consultants with no intellectual or financial interest (named above). The other authors have indicated they have no potential conflicts of interest relevant to this article to disclose.
Comments
Confounding Variables
I am writing in response to the thought-provoking journal article that proposes a correlation between driving skills measured at the time of licensure and subsequent crash rates during the initial period of licensure. While I commend the authors for their efforts in accounting for several variables such as sex, race, socioeconomic status, education, time in the learner permit phase, and driving classes skill clusters, it is crucial to note that there could be other confounders at play. Factors such as individual temperament, attitude towards risk-taking, and exposure to real-world driving conditions might impact crash rates during the early stages of licensure. Moreover, it would be prudent to consider the broader context of the overall driving score. The article seems to focus solely on the relationship between initial driving skills and crash rates, without incorporating a comprehensive evaluation of the overall driving competency. By neglecting the inclusion of the overall driving score, which considers performance across various skill clusters, the study may have overlooked significant insights into predicting early licensure crashes, rather than just individual skill clusters. Additionally, one must recognize that driving skills alone may not be the sole determinant of crash risk during the early stages of licensure. Other factors, such as previous experience with an adult, situational awareness, and decision-making abilities, develop and evolve over time. Therefore, it is essential to consider the dynamic nature of driving competence and its impact on crash rates. Furthermore, the article's emphasis solely on driving skills measured at the time of licensure may not accurately capture an individual's driving ability. Driving skills can be influenced by various factors, including anxiety, stress, or unfamiliarity with the testing environment. Also, whether they experienced inattention, neighborhood/areas, and time of accident should also be included to account for further modifications. It would be valuable to explore how these transient factors might impact the observed correlation between initial driving skills and early licensure crashes. In conclusion, while the journal article provides valuable insights into the relationship between initial driving skills and early licensure crashes, I urge caution in drawing definitive conclusions. The external validity of this study might not be as strong or be able to be applied widely if the sample only includes a few urban, rural, and suburban areas of Ohio. The presence of confounding variables, the exclusion of an overall driving score, and the dynamic nature of driving competence all warrant further investigation. Future research should consider a more comprehensive approach that incorporates these factors to develop a more accurate predictor of crash risk during the early stages of licensure.
Sincerely,