Article Text

Download PDFPDF

Original article
Costs of early spondyloarthritis: estimates from the first 3 years of the DESIR cohort
  1. Stephanie Harvard1,2,3,
  2. Daphne Guh2,
  3. Nick Bansback2,3,
  4. Pascal Richette4,5,
  5. Maxime Dougados6,7,8,9,
  6. Aslam Anis1,2 and
  7. Bruno Fautrel3,10
  1. 1University of British Columbia, Vancouver, British Columbia, Canada
  2. 2Centre for Health Evaluation and Outcome Sciences, Vancouver, British Columbia, Canada
  3. 3Sorbonne Universités, UPMC Univ Paris 06, Paris, France
  4. 4Faculty of Medicine, University Paris Diderot, Paris, France
  5. 5AP-HP, Rheumatology Department, Lariboisiere University Hospital, Paris, France
  6. 6AP-HP, Rheumatology B Department, Cochin University Hospital, Paris, France
  7. 7Faculty of Medicine, University Paris Descartes, Paris, France
  8. 8EULAR Center of Excellence
  9. 9Clinical Epidemiology and Biostatistics, INSERM (U1153): PRES Sorbonne Paris-Cité, Paris, France
  10. 10AP-HP, Rheumatology Department, Pitié Salpétrière University Hospital, Paris, France
  1. Correspondence to Dr Stephanie Harvard; harvards{at}mail.ubc.ca

Abstract

Objectives To value health resource utilisation and productivity losses in DESIR, a longitudinal French cohort of 708 patients with early spondyloarthritis (SpA) enrolled between 2007 and 2010, and identify factors associated with costs in the first 3 years of follow-up.

Methods Self-reported clinical data from DESIR and French public data were used to value health resource utilisation and productivity losses in 2013 Euros. Factors associated with costs, including and excluding biological drugs, were identified in generalised linear models using the generalised estimating equations algorithm to account for repeated observations over participants.

Results The mean (±SD) annual cost per patient was €5004±6870 in year 1, decreasing to €4961±7457 in year 3. Patients who never received a biologic had mean 3-year total costs of €4789±6022 compared to €38 206±19 829 among those who received a biologic. Factors associated with increased total costs were peripheral arthritis (rate ratio (RR) 1.19; 95% CI 1.04 to 1.37; p<0.0001), time on biologics (RR 1.23 per month; 1.21, 1.24; p<0.0001), and average BASFI score (RR 1.18/10 point increase; 1.15, 1.25; p<0.0001). Factors associated with increased costs excluding biologics were baseline age (RR 1.10 per 5 year increase; 1.05, 1.16; p<0.0001), peripheral arthritis (RR 1.20; 1.02, 1.40; p<0.0133), time on biologics (RR 1.04 per month; 1.02, 1.05; p<0.0001), and average BASDAI score (RR 1.21 per 10 point increase; 1.16, 1.25; p<0.0001).

Conclusions In addition to biologics, factors like age, peripheral arthritis and disease activity independently increase SpA-related costs. This study may serve as a benchmark for cost of illness among patients with early SpA in the biologic era.

  • Spondyloarthritis
  • Ankylosing Spondylitis
  • Anti-TNF

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Key messages

What is already known about this subject?

  • Spondyloarthritis leads to healthcare costs as well as productivity loss, the latter having been shown to be the most significant cost driver. However, there is a lack of cost-of-illness studies among patients with early SpA in the biologics era.

What does this study add?

  • In the biologic era, medication is the greatest cost driver among patients with SpA, exceeding productivity loss. Other factors like age, peripheral arthritis and functional ability are independently associated with total costs. In a cohort of patients with early SpA using biologics for less than 3 years, biologic use was associated with increased costs outside of biologics.

How might this impact on clinical practice?

  • This study identifies factors associated with SpA-related costs in the biologic era, potentially pointing towards priority areas for biologics use and providing a new benchmark for cost of illness among patients with early SpA.

Introduction

Spondyloarthritis (SpA) is a family of chronic rheumatic diseases that includes ankylosing spondylitis (AS), psoriatic spondyloarthritis (PsA), reactive SpA, SpA with inflammatory bowel disease and undifferentiated SpA.1 In its early phases, SpA is often undifferentiated, and corresponds to the concept of inflammatory back pain (IBP), the key feature of the disease.2 To date, most economic studies of SpA have examined outcomes among patients with AS or PsA specifically, and have not included patients with other SpA subtypes or with early SpA symptoms that are still evolving. A severe form of SpA, AS has been shown to cause significant pain, functional disability and loss of mobility3 In turn, these physical outcomes are associated with substantial health resource utilisation and work productivity loss among patients with AS.4 ,5

In 2012, a systematic review identified 15 economic studies of SpA conducted worldwide since 2002, nine of which evaluated health resource utilisation and work productivity loss. The median cost of health resource utilisation in 2012 US dollars was $3764 per patient per year, compared to a median cost of productivity loss of $4999.6 More recently, biological TNF blockers have been made available for the treatment of moderate to severe SpA, increasing the cost of illness compared to previous years in which non-steroidal anti-inflammatory drugs (NSAIDs) were the only available drug therapy.7

The objective of this study was to value health resource utilisation and productivity losses among patients in the DESIR cohort,8 a longitudinal, multicentre study of early SpA in France. We further aimed to determine factors associated with total costs and with costs excluding biologics (‘non-biologic costs’) among DESIR patients.

Methods

Participants

The DESIR cohort9 includes 708 patients aged 18–50 with early IBP lasting more than 3 months but less than 3 years, and suggestive of SpA according to the rheumatologists' assessment (score ≥5 on a 0–10 Numerical Rating Scale where 0=not suggestive and 10=very suggestive). Patients were required to fulfil the Calin10 or Berlin11 IBP criteria. Patients with a definitive diagnosis of non-SpA back pain, a history of previous biological TNF-blocker use, or conditions that could affect informed consent and/or compliance (eg, alcoholism, psychiatric disorders) were excluded. All patients attended their baseline study visit between October 2007 and May 2010. Follow-up visits collecting self-report questionnaire data occurred every 6 months in the first 2 years and every year thereafter, and covered the health resource use since the previous visit. Our analyses included the first 3 years of follow-up, that is, baseline plus follow-up visits at months 6, 12, 18, 24, 36. Data up to month 12 were from the DESIR database locked in February 2014, and data from months 18 to 36 from that locked in April 2015.

Costing methods

Our analysis estimated the total cost from the societal perspective (all payers combined) of all-cause health resource utilisation, specifically health practitioner visits, hospitalisations (including emergency room visits and surgeries), medical workups, and medications), and work productivity losses. Although we aimed to include as many societal costs as possible, other patient out-of-pocket costs (eg, transportation, devices, caregiver expenses and other costs sometimes termed ‘non-medical’) were not included.

Valuing health resource use

Detailed costing methods are described in online supplementary tables S1–S4. Unit costs for health practitioner visits were estimated by adding base tariffs from the French National Health Insurance12 to the average cost of extra billings based on data from the IRDES (Institut de Recherche et de Documentation en Economie de la Santé)13 Eco-Santé database.14 Data on medical workups included laboratory analyses on blood and urine and functional/diagnostic tests. On the basis of clinical expertise, the cost of blood tests was estimated assuming the following standard laboratory analyses among patients with SpA: erythrocyte sedimentation rate (ESR); C-reactive protein (CRP); haemogram including platelets; thrombocytes; transaminases; creatinine; and creatinine clearance. Blood collection was valued at 15 min of nurse care; the unit cost of blood tests was the cost of analysis plus the cost of collection as per the French public insurance fee schedule. Urine tests were assumed to include urine protein only and no collection fee was applied. The following tests were attributed a single standard cost: mammography, MRI, bone densitometry, scintigraphy, respirometry, colonoscopy and fibroscopy. For X-rays, ultrasounds, scans and ‘other exams’, free-text data were reviewed individually and assigned specific examination codes which were linked to unit costs.

Data on medications included NSAIDs, DMARDs (conventional synthetic—csDMARDs—or biologic—bDMARD), corticosteroids (oral, intramuscular, intra-articular and intravenous) and analgesics. For abatacept, infliximab and certolizumab, prices per milligram were derived from list prices in France using the Vidal dictionary.15 For all other drugs, cost data from the French public drug program16 were used to derive a price per milligram for each drug. Drugs were valued by multiplying price per milligram by number of milligrams per day by number of days of use (except injectable corticoids, which were valued by multiplying price per milligram by number of milligrams in the reported number of injections). Where dose data were missing, standard dosages were imputed taking into account the patient's weight where applicable. For infliximab only, a standard cost for administration in day hospital was applied every 6 weeks for the duration of the therapy.

Data on hospitalisations were assigned a diagnosis or act code and linked to Diagnosis-Related Group (DRG) codes using 2012 data from the French national agency for hospital information (ATIH, ‘Agence Technique de l'Information sur l'Hospitalisation’).17 Selected DRG were then linked to costs using 2012 ATIH data.18 The base-case analysis used the 6-digit DRG selected as the best fit; sensitivity analyses were performed using DRG sharing the same first five digits as the best-fit DRG (see online supplementary table S3). All unit costs were adjusted to 2013 Euros based on the Consumer Price Index19 and multiplied by frequency of use.

Valuing work productivity loss

Data on patients' profession were collected in eight categories (table 1). Average daily wage data were obtained from the French Ministry of Agriculture20 and from the French National Institute of Statistics and Economic Studies (‘INSEE’).20 Productivity loss was valued by multiplying the number of self-reported workdays lost per period by the average daily wage by professional category over the entire population of French workers as estimated by these public data sources. The age and sex distribution of the DESIR cohort was compared to that of the population of French workers from which average daily wages were obtained and wages were not further adjusted for age and sex.

Table 1

Baseline characteristics of the DESIR cohort

Handling missing clinical and cost data

Missing cost and clinical data were imputed using the Monte Carlo Markov Chain (MCMC) multiple imputation procedure, the last observation carried forward (LOCF) method, probabilistic imputation, or with negative values based on clinical expertise. Specifically, for clinical variables in which fluctuations are normal and commonly observed clinically (ie, ASDAS-CRP, BASDAI, BASFI and HAQ), missing data were handled with MCMC imputation, which was considered the best strategy for representing this variability. For clinical variables observed to have greater stability, either the LOCF method or probabilistic imputation was used in lieu of MCMC imputation, as these methods provide better control over the variability of values imputed; LOCF was used if all patients had baseline data, while probabilistic imputation was used if baseline data were missing. The imputation model is described in detail in the online supplementary appendix.

Statistical analyses

Patient characteristics were described by the mean±SD or number (%) at baseline. For descriptive purposes, patients were divided into four subsets based on imaging status: X-ray damage only; MRI inflammation only; X-Ray damage and MRI inflammation; neither. Costs of health resource utilisation and productivity losses were expressed as yearly costs. For health resource utilisation, year 1 costs were the sum of costs reported at months 6 and 12; year 2 costs were the sum of costs reported at months 18 and 24; year 3 costs were those reported at month 36. For productivity losses, yearly costs were those reported at the 12-month, 24-month and 36-month visits. Yearly costs were described by the median and mean±SD.

To describe factors associated with annual total costs and non-biologic costs over the 3 years, we used generalised linear models with a γ distribution of log link and a generalised estimating equations algorithm to account for repeated measures within participants. The independent variables year, baseline age and sex were chosen to be included in all models throughout development. Other independent variables of interest included baseline education; profession; presence of peripheral arthritis at baseline; DESIR study centre; number of months on biologics; baseline imaging status; baseline BASDAI and BASFI values; and mean BASDAI and BASFI values (where year 1 mean=mean of months 0, 6, 12; year 2 mean=mean of months 12, 18, 24; and year 3 mean=mean of months 24, 36). At the outset, baseline and mean BASDAI and BASFI values were tested for collinearity and examined alone and together in models. Goodness of fit assessed by the quasi-likelihood statistic (QIC) was used to compare the model fit among models with various BASDAI and BASFI values (baseline, mean, or baseline and mean together).21–23 The model with the best fit was chosen for inclusion in subsequent models. All other independent variables of interest were examined separately in models adjusting for year, age and sex, and variables associated with the outcome at p<0.20 were selected for inclusion in subsequent models. Model selection was then done in a backward stepwise manner, beginning with all variables and removing those that were not associated with the outcome at p<0.05 to increase goodness of fit based on the QIC. All analyses were performed using SAS 9.4.

Results

Patient characteristics

At baseline, the mean age of patients was 33.7±8.6 years and mean disease duration was 1.5±0.8 years. Just under half of the patients were male (46.2%) and the majority were Caucasian (89.6%) (table 1). Most patients were HLA-B27-positive (57.9%). The mean BASDAI score (on 100) at baseline was 44.7±20.0 and the mean BASFI score was 30.5±22.8. The most common extrarheumatological manifestations were psoriasis (15.8%) and uveitis (8.5%). At baseline, 286 patients satisfied the Assessment of SpondyloArthritis international Society (ASAS imaging criteria (40.4%) and 404 the ASAS clinical criteria (57.1%); 486 patients satisfied one or the other (68.6%). Table 1 shows the characteristics of patients by imaging status, that is, X-ray damage alone (7.2%), MRI inflammation (13.8%); X-ray damage and MRI inflammation (18.6%), neither (56.9%); a total of 24 patients had no data available to describe imaging status. Attendance at follow-up visits declined over time: of the 708 patients enrolled at baseline, 704 attended the 6-month visit, 698 the 12-month visit, 691 the 18-month visit, 692 the 24-month visit, and 631 the 36-month visit.

Costs

Health resource utilisation

Table 2 shows the estimated costs of health resource utilisation and productivity losses among all DESIR patients over 3 years. Medication was the largest cost component in all years, representing over 50% of all costs in all years. The absolute mean cost of medication increased from €2680.2±5339.6 in year 1 to €3339.2±6224.1 in year 2 to €3396.3±6476.5 in year 3. The proportion of patients incurring costs from biologics was approximately a quarter each year. The cost of biologics accounted for approximately 95% of medication costs each year.

Table 2

Cost components as proportion of total costs by year of follow-up and biologics use

Health practitioner visits, the second largest cost component of health resource utilisation in all years, declined from €733.2±788.9 in year 1 to €615.0±783.4 in year 2 and €481.0±676.2 in year 3. Hospitalisation costs were incurred by a quarter to a fifth of patients, representing under 10% of costs in each year. The mean hospitalisation cost per patient was €498.4±1318.9 in year 1, €526.2±1325.7 in year 2 and €423.3±1100.8 in year 3. Medical workups were consistently the smallest component of health resource utilisation, never exceeding 5% of yearly costs.

Productivity losses

Approximately a third of patients (30.1%) incurred productivity losses in year 1; this decreased to a quarter of patients in year 2 (23.4%) and year 3 (24.4%). The mean cost of productivity losses was €843.4±2899.1 in year 1, €768.9±3368.0 and €497.0±1951.0 in year 3. In year 1, productivity losses represented 16.8% of total costs, compared to 14.1% of total costs in year 2 and 10.0% of total costs in year 3.

Total costs

The mean annual total cost per patient was €5004.1±6870.2 in year 1, €5444.1±7936.9 in year 2 and €4960.6±7457.4 in year 3. Mean annual total costs rose with biologic use: patients who did not receive a biologic had mean total costs of €1867.9±3056.6 in year 1, €1628.2±3090.3 in year 2 and €1292.9±2059.6 in year 3, compared to €11 736.4±7882.6 in year 1, €13 635.4±8915.1 in year 2 and €12 833.9±8667.4 in year 3 among patients who received a biologic. In all years, biologics represented over 50% of total costs, up to 66.3% in year 3. Overall, the mean 3-year total cost among patients was €15 408.7±19 793.5. Patients who never received a biologic had mean 3-year total costs of €4789.04±6021.8 compared to €38 205.74±19 829.3 among those who received a biologic. The estimated total costs in the cohort over 3 years amounted to €10 909 399.9±19 574.0. Of this, patients who received a biologic accounted for an estimated €8 596 291.6±4860.9 or 78.8% of all costs in the cohort. Figure 1 compares the trends in costs over the study period overall and for patients who received and did not receive a biologic.

Figure 1

Yearly costs among biologics users.

Factors associated with costs

Total costs

In the final model (table 3), each month of biologics use was associated with a 23% increase (rate ratio (RR): 1.23; 95% CI 1.21 to 1.24) in total costs. The presence of peripheral arthritis at baseline was associated with a 19% increase (RR: 1.19; 95% CI 1.04 to 1.37) in total costs, while every 10-point increase in the average BASFI score was associated with an 18% increase in total costs (RR: 1.18; 95% CI 1.15 to 1.22). Although absolute costs were higher in year 2 than year 1 (table 4), when adjusting for time on biologics, costs in years 2 and 3 were significantly lower than in year 1 (table 3). Imaging status was not significantly associated with total costs and was excluded from the final model (figure 2).

Table 3

Models of total cost and non-biologic costs

Table 4

Health resource use and productivity loss costs in 2013 Euros by year of follow-up and biologics use

Figure 2

Yearly costs among non-biologics users.

Cost excluding biologics

In the final model (table 3), each month of biologics use was associated with a 4% increase (RR: 1.04; 95% CI 1.02 to 1.05) in non-biologic costs. The presence of peripheral arthritis at baseline was associated with a 20% increase (RR: 1.20; 95% CI 1.02, 1.40) in non-biologics costs; every 10-point increase in the average BASDAI score was associated with a 21% increase in non-biologic costs (RR: 1.21; 95% CI 1.16 to 1.25). Every 5-year increase in age was associated with a 10% increase (RR: 1.10; 95% CI 1.05 to 1.16) in non-biologic costs, while males had 31% less costs relative to females (RR: 0.69; 95% 0.58 to 0.81). Imaging status was not significantly associated with costs excluding biologics and was excluded from the final model.

Discussion

This study described costs related to health resource utilisation and work productivity losses among patients with early SpA over 3 years. Our findings reflect the cost impact of biologics, which represented by far the largest cost component among DESIR patients. In the prebiologic era among patients with AS, work productivity loss was observed to account for the majority of total costs, for example, from 62% of total annual costs among patients with AS in Hong Kong24 to 74% of total annual costs among patients with AS in the UK.25 In our study, an analysis of patients with early SpA in the biologic era, productivity loss may be expected to represent a lower proportion of costs given patients' better functional status and higher drug costs; that said, productivity loss costs among DESIR patients still represented between 10% and 17% of annual costs over the study period. In terms of absolute costs, work productivity loss costs among DESIR patients (which ranged from €843 in year one to €497 in year 3) are consistent with those in a Dutch cohort of patients with early SpA (approximately €422 over 1 year).26 We note that our estimates of productivity loss costs are conservative: we did not calculate the costs of presenteeism, that is, the cost of reduced performance at work, which is another potentially important source of costs.27

Currently, there is increased interest in non-radiographic axial SpA, defined as SpA in the absence of definite sacroiliac changes on X-ray, which includes both patients who satisfy only the clinical criteria for SpA and those with active inflammation in the sacroiliac joints visible only on MRI.28 To the best of our knowledge, this is the first longitudinal study of patients with SpA that has explored the impact of radiographic damage on cost, distinguishing between damage on X-Ray with or without MRI inflammation. Interestingly, while preventing progression to radiographic SpA remains a high priority from the clinical perspective, we found no significant independent effect of imaging status on total costs or costs excluding biologics.

In this study, the most significant cost component was biologic drugs, and patients who ever received a biologic accounted for over 75% of total costs in the cohort. While this was not surprising, the finding that approximately a quarter of DESIR patients had received a biologic by the 6-month visit (biologics use at baseline was not permitted9) was higher than anticipated, suggesting that biologics use was already relatively common even among patients with early SpA in 2010. Whether biologics provide good value for money is being explored in a separate study, which considers the impact of biologics on health-related quality of life. Our study also found peripheral arthritis, age and functional ability to be related to total costs, indicating subgroups where biologic use may be prioritised. In this study, we found that, in the first 3 years, treatment with biologics leads to increased costs in other health resource use domains, not a reduction as is sometimes expected; follow-up is needed to determine whether this trend holds over a longer period.

This study has some limitations. The DESIR cohort includes only patients with early SpA in France between 2007 and 2014; the results cannot necessarily be generalised to patients with SpA with more advanced disease or in other settings at other times. Indeed, it is unknown whether similarly high rates of biologic use exist among patients with early SpA in other countries; another recently established early SpA cohort (‘SPACE’) is ongoing in the Netherlands,29 but to the best of our knowledge the prevalence of biologics use in this cohort has not as yet been reported. A recent Canadian study of patients with SpA seen clinically between 2003 and 2014 found that 48% had been treated with a biologic, yet the mean disease duration (8 years) was much longer in this sample than in DESIR.30 Generalisable to patients with early SpA in France, our findings will allow for future comparisons with other cohorts of patients with early SpA, as well as the observation of trends within DESIR over time.

In terms of costing methods employed, we did not include non-SpA-related drugs, transportation or presenteeism, or include other patient out-of-pocket costs. Attendance declined over time and data imputation was required. Resource utilisation was self-reported and the potential exists for either underreported or overreported utilisation. In order to exclude the excess resource utilisation associated with study participation, we excluded the baseline visit from cost estimates. However, study participation could have also resulted in underestimation of resource utilisation following the baseline visit, if study visits replaced normal resource utilisation but were not reported as such. Finally, work productivity losses were estimated using self-reported workdays lost (rather than hours lost) and average daily wages by professional category, estimated from the entire population of French workers. We used profession-specific daily wages and did not further adjust wages for age and sex, although there are small differences in the age and sex distribution in DESIR compared to the French working population from which income data were derived (48% female and 77% aged 25–54 among French workers vs 54% female and 81% aged 25–54 in DESIR31). By not adjusting for the 6% higher proportion of females in DESIR, it could be argued that we have slightly overvalued productivity losses, as daily wages among females in France are lower than among males; however, we note that wage inequality functions more generally to minimise the value of productivity loss. Although our valuation method provides a good estimate of productivity losses, we acknowledge that a greater level of accuracy could be achieved with direct self-report data on wage and productivity losses expressed in hourly terms. Despite these limitations, our findings are consistent with another recent description of patients with early SpA26 and reflect patterns of health resource utilisation and work productivity loss within one of the largest early SpA cohorts to date.

Acknowledgments

The authors acknowledge Sabrina Dadoun, Violaine Foltz, Fanny Roure, Guillaume Bertrand and Huiying Sun. The DESIR cohort is conducted under the control of Assistance publique–Hopitaux de Paris via the Clinical Research Unit Paris-Centre and under the umbrella of the French Society of Rheumatology and Inserm (Institut national de la santé et de la recherche médicale). The database management is performed within the department of epidemiology and biostatistics (Professor Jean-Pierre Daurès, D.I.M., Nîmes, France). An unrestricted grant from Wyeth Pharmaceuticals was allocated for the first 5 years of the follow-up of the recruited patients. We also wish to thank the different regional participating centres: Pr Maxime Dougados (Paris–Cochin B), Pr André Kahan (Paris–Cochin A), Pr Olivier Meyer (Paris–Bichat), Pr Pierre Bourgeois (Paris–La Pitié-Salpetrière), Pr Francis Berenbaum (Paris–Saint Antoine), Pr Pascal Claudepierre (Créteil), Pr Maxime Breban (Boulogne Billancourt), Dr Bernadette Saint-Marcoux (Aulnay-sous-Bois), Pr Philippe Goupille (Tours), Pr Jean-Francis Maillefert (Dijon), Dr Xavier Puéchal (Le Mans), Pr Daniel Wendling (Besancon), Pr Bernard Combe (Montpellier), Pr Liana Euller-Ziegler (Nice), Pr Philippe Orcel (Paris–Lariboisière), Pr Pierre Lafforgue (Marseille), Dr Patrick Boumier (Amiens), Pr Jean-Michel Ristori (Clermont- Ferrand), Dr Nadia Mehsen (Bordeaux), Pr Damien Loeuille (Nancy), Pr René-Marc Flipo (Lille), Pr Alain Saraux (Brest), Pr Corinne Miceli (Le Kremlin Bicêtre), Pr Alain Cantagrel (Toulouse), Pr Olivier Vittecoq (Rouen).

References

Footnotes

  • Contributors SH collected the economic data, cleaned and merged the clinical and economic data sets, undertook the costing methods described, developed the economic models, and drafted the manuscript. DG led the data management effort, provided expert statistical input (both theoretical and in computer programming), and designed the economic models. NB made substantial contributions to the development of methods, economic aspects of model development, and the interpretation of results. PR made substantial contributions to the conception of the work and acquisition of data. MD made substantial contributions to the conception of the work, the acquisition of data, clinical aspects of model development, and the interpretation of results. AA made substantial contributions to the design of the work, economic aspects of model development, and the interpretation of results. BF led the development of costing methods, the conception and design of the work, and contributed substantially to the acquisition of data, model development and the interpretation of results. All authors reviewed the work and revised it critically for important intellectual content. All authors give final approval of the version published and agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

  • Funding The DESIR cohort is financially supported by unrestricted grants from both the French Society of Rheumatology and Pfizer Ltd, France.

  • Competing interests None declared.

  • Ethics approval University of British Columbia Research Ethics Board.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement Additional unpublished data from this study are available to peer reviewers from the researchers on request.