In the paper, crash prediction models for estimating the safety of the rural motorways are presented. Separate models were developed for total crashes and severe (fatal plus all injury) crashes. Generalized linear modeling techniques were used to fit the models, and a negative binomial distribution error structure was assumed. The study made use of a sample of 2,245 crashes (728 severe crashes) that occurred in the period 2001–2005 in the Motorway A16 Naples-Canosa in Italy. Many characteristics of the motorway are sub-standard. It allowed to investigate a wide spectrum of geometric configurations. The models were developed by the stepwise forward procedure using explanatory variables related to: traffic volume and composition, horizontal a...