Objectives:​ We aimed to develop and validate a fracture risk algorithm for the automatic identification of subjects at high risk of imminent and long-term fracture risk.

Research, design, and methods:​ A cohort of subjects aged 50–85, between 2007 and 2017, was extracted from the Catalan information system for the development of research in primary care database (SIDIAP). Participants were followed until the earliest of death, transfer out, fracture, or 12/31/2017. Potential risk factors were obtained based on the existing literature. Cox regression was used to model 1 and 5-year risk of hip and major fracture. The original cohort was randomly split in 80:​20 for development and internal validation purposes respectively. External validation was explored in a cohort extracted from the Spanish database for pharmaco-epidemiological research in primary care.

Results:​ A total of 1.76 million people were included from SIDIAP (50.7 % women with mean age of 65.4 years). Hip and major fracture incidence rates were 3.57 [95%CI 3.53 to 3.60] and 11.61 [95%CI 11.54 to 11.68] per 1000 person-years, respectively. The derived model included 19 risk factors. Internal validity showed good results on calibration and discrimination. The 1-year C-statistic for hip and major fracture were 0.851 (95%CI 0.853 to 0.864), and 0.717 (95%CI 0.742 to 0.749) respectively. The 5-year C-statistic for hip and major fracture were 0.849 (95%CI 0.847 to 0.852) and 0.724 (95%CI 0.721 to 0.727) respectively. External validation showed good performance for hip and major fracture risk prediction.

Conclusions:​ We have developed and validated a clinical prediction tool for 1- and 5-year hip and major osteoporotic fracture risks using electronic primary care data. The proposed algorithm can be automatically estimated at the population level using the available primary care records. Future work is needed on the cost-effectiveness of its use for population-based screening and targeted prevention of osteoporotic fractures.