Abstract Background Preference-based health-state utility values (HSUVs), such as the EuroQol five-dimensional questionnaire (EQ-5D-5L), are needed to calculate quality-adjusted life-years (QALYs) for cost-effectiveness analyses. However, these are rarely used in clinical trials of interventions in axial spondyloarthritis (axSpA). In these cases, mapping can be used to predict HSUVs. Objective To develop mapping algorithms to estimate EQ-5D-5L HSUVs from the Bath Ankylosing Disease Activity Index (BASDAI) and the Bath Ankylosing Spondylitis Functional Index (BASFI). Methods Data from the British Society for Rheumatology Biologics Register in Ankylosing Spondylitis (BSRBR-AS) provided 5122 observations with complete BASDAI, BASFI, and EQ-5D-5L responses covering the full range of disease severity. We compared direct mapping using adjusted limited dependent variable mixture models (ALDVMMs) and optional inclusion of the gap between full health and the next feasible value with indirect response mapping using ordered probit (OPROBIT) and generalised ordered probit (GOPROBIT) models. Explanatory variables included BASDAI, BASFI, and age. Metrics to assess model goodness-of-fit and performance/accuracy included Akaike and Bayesian information criteria (AIC/BIC), mean absolute error (MAE) and root mean square error (RMSE), plotting predictive vs. observed estimates across the range of BASDAI/BASFI and comparing simulated data with the original data set for the preferred/best model. Results Overall, the ALDVMM models that did not formally include the gap between full health and the next feasible value outperformed those that did. The four-component mixture models (with squared terms included) performed better than the three-component models. Response mapping using GOPROBIT (no squared terms included) or OPROBIT (with squared terms included) offered the next best performing models after the three-component ALDVMM (with squared terms). Simulated data of the preferred model (ALDVMM with four-components) did not significantly underestimate uncertainty across most of the range of EQ-5D-5L values, however the proportion of data at full health was underrepresented, likely due in part to model fitting on a small number of observations at this point in the actual data (4%). Conclusions The mapping algorithms developed in this study enabled the generation of EQ-5D-5L utilities from BASDAI/BASFI. The indirect mapping equations reported for the EQ-5D-5L facilitate the calculation of the EQ-5D-5L utility scores using other UK and country-specific value sets.