Background Traumatic cervical spine injuries are amongst the traffic injuries that can cause most harm to a person. Classifying subtypes of clinical presentations has been a method used in other pathologies to diagnose more efficiently and to address the appropriate treatment and the prognosis. The management of patients suffering from cervical injuries could be improved by classifying the severity of the impairment. This will allow clinicians to propose better treatment modalities according to the severity of the injury. Materials and methods The present study is a retrospective cohort study performed with the clinical data from 772 patients stored at Fisi-(ON) Health Group. All the patients treated for cervical spine injuries are evaluated using the EBI-5 ® system, which is based on inertial measurement unit (IMU) technology. The normalized range of motion of each patient was incorporated into a single index, the Neck Functional Holistic Analysis Score (NFHAS). Results Clustering analysis of the patients according to their NFHAS resulted in five groups. The Kruskal-Wallis H test showed that there were statistically relevant differences in the ROM values and NFHAS of the patients depending on the cluster they were assigned to: FE X 2 (4) = 551.59, p = 0.0005; LB ROM X 2 (4) = 484.58, p = 0.0005; RT ROM X 2 (4) = 557.14, p = 0.0005; NFHAS X 2 (4) = 737.41, p = 0.0005. Effect size with ηp 2 for the comparison of groups were: FE = 0.76, LB = 0.68, RT = 0.76 and NFHAS = 0.96.