Article Text

Concise report
A matrix risk model for the prediction of rapid radiographic progression in patients with rheumatoid arthritis receiving different dynamic treatment strategies: post hoc analyses from the BeSt study
  1. K Visser1,
  2. Y P M Goekoop-Ruiterman2,
  3. J K de Vries-Bouwstra3,
  4. H K Ronday2,
  5. P E H Seys4,
  6. P J S M Kerstens5,
  7. T W J Huizinga1,
  8. B A C Dijkmans3,
  9. C F Allaart1
  1. 1Leiden University Medical Center, Leiden, The Netherlands
  2. 2HAGA Hospital, The Hague, The Netherlands
  3. 3VUMC, Amsterdam, The Netherlands
  4. 4Franciscus Hospital, Roosendaal, The Netherlands
  5. 5JBI, Amsterdam, The Netherlands
  1. Correspondence to Dr K Visser, Department of Rheumatology C1-R, Leiden University Medical Center, Albinusdreef 2, 2333 ZA Leiden, The Netherlands; k.visser{at}


Objectives To develop a matrix model for the prediction of rapid radiographic progression (RRP) in subpopulations of patients with recent-onset rheumatoid arthritis (RA) receiving different dynamic treatment strategies.

Methods Data from 465 patients with recent-onset RA randomised to receive initial monotherapy or combination therapy were used. Predictors for RRP (increase in Sharp-van der Heijde score ≥5 after 1 year) were identified by multivariate logistic regression analysis. For subpopulations, the estimated risk of RRP per treatment group and the number needed to treat (NNT) were visualised in a matrix.

Results The presence of autoantibodies, baseline C-reactive protein (CRP) level, erosion score and treatment group were significant independent predictors of RRP in the matrix. Combination therapy was associated with a markedly reduced risk of RRP. The positive and negative predictive values of the matrix were 62% and 91%, respectively. The NNT with initial combination therapy to prevent one patient from RRP with monotherapy was in the range 2–3, 3–7 and 7–25 for patients with a high, intermediate and low predicted risk, respectively.

Conclusion The matrix model visualises the risk of RRP for subpopulations of patients with recent-onset RA if treated dynamically with initial monotherapy or combination therapy. Rheumatologists might use the matrix for weighing their initial treatment choice.

Statistics from

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.


A major goal in the treatment of rheumatoid arthritis (RA) is to prevent joint damage progression and thus prevent long-term disability.1 At the group level, initial combination therapies including either prednisone or a tumour necrosis factor (TNF) inhibiting agent are more effective in achieving that goal than initial disease-modifying antirheumatic drug (DMARD) monotherapy.2,,5 Several predictive factors for radiographic progression have previously been reported to help identify the patients at risk of a more severe disease course who might benefit from such therapies.6,,9 However, in clinical practice it is still difficult to translate this evidence into treatment choices for individual patients with recently diagnosed RA. This is probably due to the fear of toxicity of DMARD plus prednisone combinations and the cost of TNF-blocking therapies.10 Furthermore, most existing algorithms assessing the risk for joint damage progression have not included treatment, yet treatment is one of the most important determinants of radiographic outcome.11,,14 A recent model which did include treatment was based on a static treatment study not reflecting the dynamics of clinical practice, where treatment is adjusted if the initial therapy fails to achieve low disease activity.15

We therefore set out to develop a matrix risk model for the prediction of rapid radiographic progression (RRP) based on a cohort of patients with recent-onset RA who were dynamically treated aimed at a low disease activity (disease activity score (DAS) ≤2.4).3 We identified baseline predictors for RRP in subsets of patients treated according to various treatment strategies which could help physicians in making their initial treatment choice.


Patients and treatment

Data from patients with recent-onset RA randomised to four different treatment strategies were used.3 All patients were DMARD-naïve and fulfilled the 1987 American College of Rheumatology criteria for RA.16 In groups 1 (n=126) and 2 (n=121), patients started with initial methotrexate monotherapy which could be switched to or extended with other DMARDs. In group 3 (n=133), patients started with a combination of methotrexate, sulfasalazine, hydroxychloroquine and a tapered high-dose prednisone. In group 4 (n=128), patients started with a combination of methotrexate and infliximab. Every 3 months the treatment was adjusted according to a fixed protocol, aiming at a DAS ≤2.4.17 Groups 1 and 2 were combined for the current analyses.

Statistical analysis

As outcome, RRP was defined as an increase in the Sharp-van der Heijde score (SHS) of ≥5 after 1 year, corresponding to the smallest detectable change and expert opinion on clinically relevant progression.18 19 Baseline characteristics were compared between patients with and without RRP per treatment group using univariate logistic regression analysis. In a multivariate analysis with backward selection including the variable treatment group, significant independent predictors for RRP were identified using a p value of >0.10 as removal criterion. Possible interactions with treatment were explored. More information on the baseline and outcome variables is available in the online supplement.

Matrix construction

To construct the matrix, the final regression model was fitted with all variables categorised, based on clinically relevant cut-offs and tertiles. The predicted risk of RRP was calculated for all possible risk factor combinations and presented in visual matrices for each treatment group. In addition, for the initial combination therapies compared with initial monotherapy, the risk differences were converted into numbers needed to treat (NNTs) and likewise presented in a matrix.

Discriminative/predictive ability and internal validation

The area under the receiver operating characteristics (ROC) curve was calculated (area under the curve (AUC)) and the positive and negative predictive values (PPV, NPV) of the model were explored. The predicted risk for the monotherapy group was compared with the observed radiographic progression after 5 years of continued DAS-guided treatment. Additional information can be found in the online supplement.

SPSS Version 16.0 software (SPSS, Chicago, Illinois, USA) was used for all analyses.


The baseline characteristics of 465 of 508 patients (92%) for whom radiographs were available are shown in table S1 in the online supplement.

Univariate analysis

RRP was observed in 75/224 patients (33%) in the initial monotherapy group, in 16/120 patients (13%) in the initial combination with prednisone group and in 11/121 patients (9%) in the initial combination with infliximab group. Significant univariate predictors of RRP in all three treatment groups were C-reactive protein (CRP), erythrocyte sedimentation rate (ESR), SHS and the erosion score at baseline (see table S2 in the online supplement). The presence of rheumatoid factor (RF) and anti-citrullinated protein antibodies (ACPAs) were only significant univariate predictors in the monotherapy group.

Multivariate and fitted logistic regression analysis

In addition to the univariate predictors, age, gender, body mass index, symptom duration, swollen joint count, DAS, health assessment questionnaire and treatment group were also entered into the multivariate analysis. Because of colinearity, a merged variable including both RF and ACPA status was constructed and the erosion score, which is easier to perform, was included instead of the total SHS.

Independent significant predictors for RRP were the baseline erosion score, RF/ACPA, CRP and treatment group. A comparable model was found if ESR was included instead of CRP. No significant interaction was present between each independent predictor and the variable treatment group. The final model with all variables categorised is shown in table 1. Categorising the variables decreased the explained variance (Nagelkerke's R2) from 0.33 to 0.31 in the CRP model and from 0.30 to 0.29 in the ESR model. Treatment alone accounted for one-third of the prediction of RRP (R2=0.11 in both models). The R2 of the model with RF or ACPA as dichotomous variables instead of combined was 0.28.

Table 1

Independent predictive variables for rapid radiographic progression (RRP) in the final fitted multivariate logistic regression model with all variables categorised

Matrix models

The matrices show the percentage predicted risk of RRP after 1 year of dynamic treatment with initial monotherapy or initial combination therapy for patients with specific risk profiles, indicated by the coloured boxes in figure 1. Matrices including ESR instead of CRP had comparable results (data not shown). The matrices of the initial combination therapy groups showed an overall reduced risk compared with the initial monotherapy matrix, especially in patients with an unfavourable risk profile. For example, the upper right box in figure 1A–C shows that, for a patient who is RF and ACPA positive and has a high CRP level and already has erosions at baseline, the predicted risk for RRP is 78% with initial monotherapy which is reduced to 42% or 34% with initial combination therapy including prednisone or infliximab, respectively.

Figure 1

Matrix model including C-reactive protein (CRP) with predicted risk (in percentages) of rapid radiographic progression (RRP) in subpopulations of patients with recent-onset rheumatoid arthritis (RA) per treatment strategy. (A) Initial methotrexate monotherapy. (B) Initial disease-modifying antirheumatic drug (DMARD) combination therapy with prednisone. (C) Initial combination therapy with infliximab (IFX). Red boxes represent a high risk (≥50%), orange boxes an intermediate risk (20–50%), yellow boxes a low risk (10–20%) and green boxes a very low risk (<10%) of RRP. ACPA, anti-citrullinated protein antibodies; RF, rheumatoid factor.

NNT matrix

The risk reduction associated with treatment choice can be demonstrated more clearly by calculating the corresponding NNT with initial combination therapy to prevent one patient from RRP. For the example given above, the NNT is 1/(78–42)×100=2.8 (95% CI 2 to 5) patients for the combination with prednisone and 1/(78–34)×100=2.3 (95% CI 2 to 4) patients for the combination with infliximab. Figure 2 shows that the NNTs are low for patients with bad prognostic factors, indicating that the risk of RRP for these patients can be largely reduced with initial combination therapy. The NNT is 2–3 for patients with a high risk (≥50%) for RRP, 3–7 for patients with an intermediate risk (20–50%) and 7–25 for patients with a low predicted risk (<20%).

Figure 2

Matrix model with numbers needed to treat (NNTs) with initial combination therapy to prevent one patient from progression if treated with initial monotherapy. (A) Initial combination therapy with prednisone. (B) Initial combination therapy with infliximab (IFX). The NNT is <5 in the uncoloured boxes, 5–10 in the light grey boxes and >10 in the dark grey boxes. ACPA, anti-citrullinated protein antibodies; CRP, C-reactive protein; RF, rheumatoid factor.

Discriminative/predictive ability and internal validation

An extended description of the following results can be found in the online supplement. For the model including CRP and ESR, the AUCs of the ROC curves were 0.81 (95% CI 0.77 to 0.86) and 0.80 (95% CI 0.75 to 0.85), respectively, indicating a reasonable discriminative ability (see figure S1 in online supplement). The model reliably classified 52% of the patients with a PPV of 62% and NPV of 91% (see table S3 in online supplement). The validity of the prediction after 5 years of continued DAS-steered treatment is shown in figure S2 in the online supplement.


To help rheumatologists decide how to start treatment in patients with recently diagnosed RA, we developed a matrix model from which the risk for RRP in patients with specific combinations of risk factors, if treated with different treatment strategies, can be assessed with a limited number of easily accessible clinical variables. In addition, the matrix gives insight into the risks of initial overtreatment by showing NNTs.

The current matrix is based on a dynamic treatment cohort where 3-monthly treatment adjustments were aimed at achieving and maintaining a DAS ≤2.4, closely resembling clinical practice. Our results showed that, in a tight control setting with frequent treatment adjustments, previously reported predictive variables such as the baseline erosion score, CRP, ESR, RF and ACPA still predict RRP after 1 year.7 8 11 14 20 It is also clear that the initial treatment choice is a main determinant of RRP and therefore should be included in a useful prediction model.

Our matrix confirms that, with initial combination therapy with infliximab or prednisone, the risk of RRP is markedly reduced in comparison with initial monotherapy, especially in patients with an unfavourable risk profile.21 This observation carries even more weight considering that, in the dynamic treatment design, patients on initial monotherapy could receive more intensive treatment within a relatively short time period.3 Moreover, initial combination therapy could be tapered and stopped in approximately 50% of the patients once a low disease activity was achieved.22 In this respect, the NNT allows consideration of whether to prevent RRP in one patient; initial temporary overtreatment in others is acceptable.

External validation in other dynamic treatment cohorts will be needed to increase the generalisability of the results. The general RA population probably represents a group of patients with a better prognosis and a lower prevalence of RRP than those included in the BeSt study.23 Consequently, and because the matrix is better at predicting who will not have RRP, the NPV will be even higher in the general population of patients with recent-onset RA.

In conclusion, our matrix model offers the possibility of assessing the risk of RRP for each patient with a specific combination of risk factors. It requires only a few easily accessible variables, includes relevant initial treatment options, it has been developed in a cohort resembling the dynamics of daily practice and is visually attractive and easy to use. It could therefore be a useful tool to help rheumatologists choose the optimal initial treatment for their patients with recently diagnosed RA.


The authors would like to thank all patients as well as the following rheumatologists (other than the authors) who participated in the Foundation for Applied Rheumatology Research (all locations are in The Netherlands): W M de Beus (Medical Center Haaglanden, Leidschendam); C Bijkerk (Reinier de Graaf Gasthuis, Delft); M H W de Bois (Medical Center Haaglanden, The Hague); F C Breedveld (Leiden University Medical Center, Leiden); G Collée (Medical Center Haaglanden, The Hague); J A P M Ewals (Haga Hospital, The Hague); A H Gerards (Vlietland Hospital, Schiedam); B A M Grillet (De Honte Hospital, Terneuzen); J H L M van Groenendael (Franciscus Hospital, Roosendaal); K H Han (Medical Center Rijnmond-Zuid, Rotterdam); J M W Hazes (Erasmus University Medical Center, Rotterdam); H M J Hulsmans (Haga Hospital, The Hague); M H de Jager (Albert Schweitzer Hospital, Dordrecht); J M de Jonge-Bok (retired); M V van Krugten (Walcheren Hospital, Vlissingen); H van der Leeden (retired); W F Lems (Slotervaart Hospital, Amsterdam); M F van Lieshout-Zuidema (Spaarne Hospital, Hoofddorp); A Linssen (retired); P A H M van der Lubbe (Vlietland Hospital, Schiedam); C Mallée (Kennemer Gasthuis, Haarlem); E T H Molenaar (Groene Hart Hospital, Gouda); H C van Paassen (Sint Franciscus Gasthuis, Rotterdam); A J Peeters (Reinier de Graaf Gasthuis, Delft); H K Markusse (deceased); R M van Soesbergen (retired); P B J de Sonnaville (Oosterschelde Hospital, Goes); I Speyer (Bronovo Hospital, The Hague); K S S Steen (Kennemer Gasthuis, Haarlem); J Ph Terwiel (Spaarne Hospital, Hoofddorp); A E Voskuyl (VU Medical Center, Amsterdam); M L Westedt (Bronovo Hospital, The Hague); S ten Wolde (Kennemer Gasthuis, Haarlem); J M G W Wouters (Sint Franciscus Gasthuis, Rotterdam); D van Zeben (Sint Franciscus Gasthuis, Rotterdam). The authors would also like to thank all other rheumatologists and trainee rheumatologists who enrolled patients in this study and all research nurses for their contributions.


View Abstract

Supplementary materials


  • Funding The BeSt study was supported by the Dutch College of Health Insurances. Schering-Plough and Centocor provided additional funding.

  • Competing interests None.

  • Ethics approval This study was conducted with the approval of the Leiden University Medical Center ethics committee.

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