- Open Access
Modeling rheumatoid arthritis using different techniques - a review of model construction and results
Health Economics Reviewvolume 4, Article number: 18 (2014)
Rheumatoid arthritis (RA) is a chronic, inflammatory disease with severe effects on the functional ability of patients. Due to the prevalence of 0.5 to 1.0 percent in western countries, new treatment options are a major concern for decision makers with regard to their budget impact. In this context, cost-effectiveness analyses are a helpful tool to evaluate new treatment options for reimbursement schemes.
To analyze and compare decision analytic modeling techniques and to explore their use in RA with regard to their advantages and shortcomings.
A systematic literature review was conducted in PubMED and 58 studies reporting health economics decision models were analyzed with regard to the modeling technique used.
From the 58 reviewed publications, we found 13 reporting decision tree-analysis, 25 (cohort) Markov models, 13 publications on individual sampling methods (ISM) and seven discrete event simulations (DES). Thereby 26 studies were identified as presenting independently developed models and 32 models as adoptions. The modeling techniques used were found to differ in their complexity and in the number of treatment options compared. Methodological features are presented in the article and a comprehensive overview of the cost-effectiveness estimates is given in Additional files 1 and 2.
When compared to the other modeling techniques, ISM and DES have advantages in the coverage of patient heterogeneity and, additionally, DES is capable to model more complex treatment sequences and competing risks in RA-patients. Nevertheless, the availability of sufficient data is necessary to avoid assumptions in ISM and DES exercises, thereby enabling biased results. Due to the different settings, time frames and interventions in the reviewed publications, no direct comparison of modeling techniques was applicable. The results from other indications suggest that incremental cost-effective ratios (ICERs) do not differ significantly between Markov and DES models, but DES is able to report more outcome parameters.
Given a sufficient data supply, DES is the modeling technique of choice when modeling cost-effectiveness in RA. Otherwise transparency on the data inputs is crucial for valid results and to inform decision makers about possible biases. With regard to ICERs, Markov models might provide similar estimates as more advanced modeling techniques.
In most industrialized countries the financial burden of financing health care systems has increased over the last years due to the rise of chronic diseases in consequence of demographic changes as well as the progress in the development of medical technologies. To counteract these developments, regulatory institutions and decision makers have a growing demand for transparent and valid information on the relative value of innovations in comparison to established technologies. Therefore data on relative efficacy as well as effectiveness with respect to specific patient populations today is expected by regulators and payers in the major markets to support the decision on whether to include a new treatment option into reimbursement schemes and treatment guidelines.
In this context different techniques for conducting cost-effectiveness analyses (CEAs) surfaced over the last 15 years, enabling a comparison of (incremental) costs and outcomes of new interventions vs. existing comparators or competitors. Thereby the field of modeling offers a variety of different methods allowing for an estimation of even long-term impacts of certain treatments on the course of a disease beyond the time of clinical studies. . These different approaches might have different pros and cons, e.g. with respect to the method itself and also in specific diseases, which needs to be considered by regulators and decision makers .
Over the past 15 years rheumatoid arthritis (RA) has been one of the most competitive disease areas with the introduction of new treatment strategies, e.g. various biologic therapies, leading to a broad body of literature on health economic research. Therefore, the aim of the present review is to explore the field of cost-effectiveness modeling in RA with respect to the different methods used. In this context strengths and weaknesses with respect to appropriate modeling approaches with respect to the course of RA itself as well as a comprehensive provision of treatment strategies are examined. Finally, some recommendations on the interpretation of the results from different modeling types or approaches and potential advantages of certain modeling techniques in future studies are given.
To identify developments for the use of modeling techniques in RA, a systematic literature search was conducted January 14, 2014 in PubMED via Medline using "Medical Subject Heading" (MeSH) terms. The search yielded 1,074 hits for ("Arthritis, Rheumatoid" [Mesh] AND "Economics" [Mesh]) with no filters or limitations applied. The review process is depicted in Figure 1. Publications were excluded if they were not using modeling techniques or not applied to specific treatment options. Publications were included if patients with RA were part of the study population.
Identified systematic reviews were used for cross validation of identified studies. The main focus of the review process was to identify analysis techniques applied in the modeling of the course of RA treatment. Firstly, the publications were categorized by modeling technique, country and year and screened for information on specific modeling parameters (e.g. cycle length). Secondly, clinical measures used to describe the course of the disease within the model were extracted and supplemented by the measurements used to reflect changes in those health states. Additionally, the sources of the efficacy data and the base-case results were extracted and are presented in Additional files 1 and 2.
Rheumatoid Arthritis (RA) is an inflammatory, auto-immune disease of chronic course with unknown etiology. The prevalence of RA is estimated to range between 0.5 to 1 percent in western populations whereby women are more often affected by the disease with a ratio of 3:1. Concerning the incidence, there is some evidence of a decreasing RA rate due to a later onset of the disease . Nevertheless, increased mortality of RA patients in comparison to the general population remains unchanged and the gap may be even increasing .
Etiology of the disease
The clinical course of the disease usually starts with unspecific symptoms like increasing exhaustibility, anorexia and weight loss, followed by specific symptoms of painfully swollen joints. Usually, the affection of joints occurs in a symmetrical manner and mostly affects joints of the hands and knees. As movement, stretching and pressure aggravate pain, patients tend to take a relieving posture often resulting in stiffness of the affected joints . The course of the disease is characterized by periodic phases of higher and lower disease activity, which can be measured via the Disease Activity Score (DAS) or the DAS with 28-joint count (DAS28) . The disease state can be classified using the scoring system published by the American College of Rheumatology (ACR)  (see also newest version from Alehata et al. (2010) ). In the long run, the chronic inflammation leads to damaged joints and destruction of cartilage . Thereby, the functional disability of RA patients increases on a diminishing scale and is often measured via the disability index of the Health Assessment Questionnaire (HAQ). The HAQ is one of the standard domain specific quality of life instruments in rheumatology that gives a detailed impression of specific functional impairments adding to the burden of disease ,. A premature mortality of RA is primarily cause-specific, but may also be influenced by adverse events of anti-rheumatic drugs .
Treatment options and treatment guidelines
In general, treatment options for RA include patient education, pharmacotherapy, physiotherapy and surgery . Within pharmacotherapy, drugs can be distinguished into Disease-Modifying-Anti-Rheumatic-Drugs (DMARDs) and symptom relief medications such as corticosteroids. DMARDs are recommended as monotherapy for patients with early RA in remission or with low disease activity or as combinational therapy in patients with established RA . Among the DMARDs, so-called biologics (e.g. anti-TNF compounds) play a major role today in the treatment algorithm of active RA patients. The American College of Rheumatology (ACR) treatment guidelines recommend different strategies for RA patients depending on the time since onset of disease (early vs. established RA), the grade of disease activity (low, moderate or high) and whether the patient is with or without poor prognosis . Disease activity should be reassessed at regular intervals and drug treatment should be modified or switched in the case of adverse events, non-response or loss of efficacy over time. More or less the same approaches to treatment are recommended by guidelines of The European League Against Rheumatism (EULAR).
Treatment of RA undoubtedly has become more efficacious after the introduction of biologic treatment alternatives but it has also become a major area of concern to health care payers when it comes to the budget impact. In such scenarios it is of paramount importance to explore relative cost effectiveness of different alternatives or approaches to treatment to sort out viable options not only for health care payers, but ultimately also for patients who usually have to finance the health care system one way or the other.
Health economic modeling techniques
Health economic models are common tools for decision analyses in health administrative bodies and are increasingly used in health economic research. Hereby, different model types might be distinguished from relative simple models such as decision trees to more advanced modeling exercises (e.g. agent-based models). Nowadays one of the most commonly used model type is that of a Markov-model.
The most simplistic model type used in health economics is the decision tree. Within a decision tree the course of a disease is displayed as a hierarchical, one-directional tree with either a chance node or a terminal node at the end of each branch. These models can either be constructed for single patients or cohorts, but without the possibility of interaction between individual subjects. Nevertheless, decision trees only allow for the analysis of fixed time frames, making the approach problematic in use for diseases with varying length in certain health states. Furthermore, to encounter heterogeneity between patient-groups, input parameters need to be changed and the model needs to be run separately for every subgroup.
In health economic Markov-chain modeling, the course of a disease is described by (various) discrete health states. Over time, cohorts move through or between these health states at the end of a cycle of a fixed time interval. Transition probabilities determine the number of patients remaining in a health state at the end of a cycle and also the number of subjects moving to other connected potential health states. Markov models are usually run until all patients of a cohort reached an absorbing state (e.g. death) or maximum number of cycles. The number of patients over all health states needs to stay constant in all cycles and patients are only allowed to be in one health state during a cycle . It is also important to note that there are no interaction or individual decisions modeled between the different subjects in a Markov model. As the disease is mimicked via the creation of a limited number of distinct health states, available data needs to be adjusted to represent the clinical course, costs and outcomes at those fixed health states. Similar to a decision tree Markov models can analyze subgroups only by adjustment of model parameters and repeated runs of the model to cover questions of heterogeneity. Although it is possible to cover small numbers of subgroups of to address questions of heterogeneity by additional health states, it is feasible to analyze high numbers of subgroups by repeated runs of the model with different parameter sets for technical boundaries.
Individual Sampling Methods (ISM)
Another method to address problems of heterogeneity is to use individual sampling models . In accordance to Markov models, individual sampling models or Monte Carlo simulations consist of different distinct states, but as only single patients are sent through the model, attributes of patients may be changed over time thus changing transition probabilities or other input parameters. This way individual sampling models overcome the need for a disproportional number of health states to track changes over time. For following references in this paper with respect to this method the main distinction of the individual sampling method and DES methods is defined as the use of fixed cycle lengths.
Discrete Event Simulations (DES)
In contrast to the above described modeling techniques, discrete event simulation (DES) focuses on the treatment pathway of a patient rather than modeling the course of the disease itself. In short, single entities (e.g. patients, health care professionals) with certain attributes (e.g. age, gender, disease severity) experience certain events (e.g. disease progression, hospital admission) based on their attributes, and, as a result of that, experience benefits or harms and consume resources. A strong emphasis of DES lies on the modeling of scarce resources and the behavior of entities competing or waiting for the availability of those resources . In this context, an interaction between entities, e.g. patients, is possible as they can for example form waiting queues and access resources following pre-defined rules. Also, a main advantage of DES is its capability of the inclusion of patient characteristics, thus accounting for heterogeneity .
Even though DES is one of the more advanced methods used in health economic evaluation it still has some disadvantages attached. Firstly, DES is rather process-oriented leaving out interactions between entities which are not connected to resources. Secondly, it leaves the study subjects in a rather passive role making them flow through the system instead of letting them decide independently . Additionally, the higher complexity of the model comes with a much higher demand on input data, which in many situation will not be easily accessible (e.g. patient-level data from clinical trials). Also, the sampling from distributions for each parameter for a large number of individuals involves model calibration, which is still under methodological development.
ISPOR Guidelines on good modeling practice (health transition modeling and DES-modeling)
Besides the decision tree method, the International Society for Pharmacoeconomics and Outcomes Research (ISPOR) provides guidelines for good modeling practice for each of the above mentioned modeling techniques. Concerning the choice of the model, the authors recommend a Markov cohort simulation, if the number of necessary health states which depict the health problem remains manageable. If some aspect of the health problem cannot be reproduced in an appropriate way, an ISM modeling exercise should be chosen, but sequential decisions on treatment options should not be considered in the same model . DES is recommended as first choice to model problems caused by constrained or limited resources. In addition, DES is seen as being also suitable to model time to one or multiple events stochastically and considering many characteristics of a patient influencing time. Thereby, the trade-off between the detailing of the model and the availability of data should be described and changes to the model structure or the influence of expert opinions filling data gaps should be explored and reported .
Cost-effectiveness modeling in RA
Summary over all identified analyses
Via the literature search n = 58 publications between 1996 and 2012 could be identified. 13 studies used a decision tree approach, 25 studies Markov-modeling, 13 studies an individual sampling method (ISM) and the DES-technique was applied seven times. As many of those studies are adaptations or enhancements of existing models, a family tree of all publications is displayed in Figure 2, where publications in bold letters represent the original model within a tree of references.
With 12 out of 25 studies, the highest number of independently developed models can be found for Markov-models with a similar share of six out of 13 for decision tree-studies. In contrast, four out of 13 ISMs and three out of seven DES-studies might be considered independently developed models.
On average the decision tree models have the shortest time horizons of all model types, besides Moore et al., Spiegel et al. and Zabinski et al. apply a life-time horizon -. The Markov models range from one year to life-time mostly applying a six month cycle length (three to twelve months), whereas the DES and ISM models, with one exception, take a life-time perspective. A summary of study characteristics can be found in Table 1 (sorted by model type).
As can be seen in Table 2 decision tree models tend to compare treatment strategies involving etanercept, abatacept, adalimumab and methotrexate. Markov models focus on methotrexate, an unspecified TNF-biologic or gold. Individual sampling models mostly involve methotrexate and rituximab/methotrexate. Compared to the overall number of models, DES models are evaluating a much higher number of treatment strategies as more often treatment patterns or pathways are modeled. Eleven of the DMARD strategies as well as 14 TNF-alpha-antagonist strategies compared are not further specified.
Description of Models by Modeling Technique
The earliest model included in this review is the decision tree analysis of Kavanaugh et al. , which has also been used by Choi et al.  and Choi et al. . The decision tree splits up after choice of treatment. For each treatment, side-effects (minor and major) may occur. If no side-effects emerge a patient can either (partially or fully) respond or not respond following ACR criteria. For patients with minor side-effects it is also possible to respond to the treatment. The time-horizon is six months. Choi et al.  and Choi et al.  performed only deterministic sensitivity analyses, whereas Kavanaugh et al.  used Monte-Carlo simulation for one-way sensitivity analysis.
The second family of decision tree models was built by Russell et al. , Beresniak et al. , Cimmino et al. , Puolakka et al. , Saraux et al.  and Beresniak et al. . The trees depict the treatment of patients within four six-month intervals starting with a TNF-alpha-antagonist. After each interval patients either show signs of remission, in which case they stick with the current treatment, or they change to the next intervention. In all publications the DAS28 score is used to define effectiveness using the low-disease activity state (DAS28 ≤ 3.2) as primary end-point describing cost-effectiveness. Besides Beresniak et al. , all models perform probabilistic sensitivity analysis and all six models report an ICER for costs per day in a low disease activity state with confidence intervals.
Chiou et al.  use an independent decision tree-model to evaluate the cost-effectiveness of several biologic response modifiers. Similar to the Russell-models, intervals of six months are used for a total time horizon of one year. After the first six months the effectiveness of an intervention is measured using ACR response criteria categorized into no ACR, ACR20, ACR50 and ACR70 responses. After the second interval it is possible that patients develop no, mild, moderate or severe complications. Only deterministic, one-way-sensitivity analysis is performed.
In the publication by Spiegel et al.  a decision tree is applied to the question whether rofecoxib and celecoxib are cost-effective compared to non-selective non-steroidal anti-inflammatory drugs (NSAIDs) in the treatment of chronic arthritis. Over a life-time time horizon, the number of gastrointestinal (GI) complications, i.e. dyspepsia, ulcer bleeding and ulcer perforation, is modeled for each of the alternatives. Hence, the model does not focus on the efficacy of the alternatives on RA, but on the reduction on adverse events. The cost-effectiveness measure is calculated as incremental costs per QALY. A meta-analysis is used to obtain parameters for clinical efficacy and PSA is performed.
In a similar study, Zabinski et al.  compared celecoxib alone vs. several NSAIDs in osteoarthritic (OA) and RA patients, also focusing on adverse events following these treatments. Thereby gastrointestinal events are supplemented by "anaemia with occult bleeding". The analysis is carried out for a time-horizon of 6 months and calculates costs from a third-party payer-perspective, but no ICER.
The decision tree constructed by Moore et al.  compares etoricoxib with several NSAID-strategies. Again, the main health states comprise gastrointestinal adverse events ("GI problems", "no GI problems"), but with a more complex depiction of the following treatment patterns. As in the two former studies, not only RA-patients, but also OA-patients were included in the population of the analysis. An ICER from the perspective of the NHS is reported including deterministic and probabilistic sensitivity analysis.
The most adopted Markov-model is from Kobelt et al. , which is based on two earlier works also from Kobelt et al. ,. Those studies are not included in this review as only hypothetical treatments are evaluated to justify the general use of Markov models in RA in those two studies. The authors built a model using seven states corresponding to HAQ-scores also including death with a time horizon of 10 years and a cycle-length of one year. In difference, the following studies by Kobelt et al.  and Kobelt et al.  reduced the number of disease states by one. Further, in the paper from 2011, Kobelt et al.  use a cycle length of six months and Lekander et al.  extended the time-horizon to 20 years. All following studies subdivided each health state by a high or low VAS-Score or high or low disease-activity measured by the DAS28-score. Whereas Kobelt et al.  analyze the cost-effectiveness of infliximab, the model from Kobelt et al.  does not evaluate the cost-effectiveness of a specific drug, but constitutes an advancement compared to Kobelt et al. . Kobelt et al. and Kobelt et al. investigate the cost-effectiveness of etanercept in combination with methotrexate and etanercept monotherapy. In addition, Lekander et al. extended this analysis by grouping TNF-alpha-inhibitors and Lekander et al.  analyzed the cost-effectiveness of infliximab. Another application of the model compared the results of the model for three different trial cohorts (Lekander et al. ). From Kobelt et al. on, probabilistic sensitivity analyses (PSA) has been used ,.
Welsing et al.  and Schipper et al.  are using the four Markov-states "remission", "low disease activity", "moderate disease activity" and "high disease activity" corresponding to DAS28-scores, although with different classifications. Patients may change the treatment after inadequate response. Both studies use a cycle length of three months and a time horizon of five years. PSA was performed using 1,000 iterations or 1,000 patients, respectively.
The original Markov-model from Wong et al.  has later been adapted by Barbieri et al.  for a longer time horizon. In both models, four different health states are defined using HAQ-scores for each of the five treatment options, thus resulting in 21 health states with the inclusion of death. The models extrapolate the results of the ATTRACT trial for infliximab from 54 weeks to a life-time horizon using a cycle length of six months. The first model has been validated in an earlier publication.  No PSA was performed.
Bae et al.  developed a first model for Korea to examine the cost-effectiveness of corticosteroids compared to NSAIDs and COX-2 inhibitors in RA patients. Over a life-time patients are moving through seven corticosteroid-specific or seven NSAID-specific complication-health states and the absorbing "death"-state in twelve month cycles. An ICER is given as costs per QALY and a deterministic sensitivity analysis was performed. The model has been adopted by Yun and Bae  to compare several NSAID and one COX-2 treatment strategy. Therefore, the health states concerning the NSAID-specific complications were kept in the model and supplemented by health states depicting treatment patterns.
In their first publication Maetzel et al.  focus on RA patients alone and the cost-effectiveness of adding leflunomide to a treatment sequence of DMARDs in Canada. The health states are divided into "ACR20-response" or "no ACR20-response" leading to continuation of current treatment and "severe AE" or "lack of efficacy" leading to the next treatment in the sequence. The result is given as an ICER (costs per QALY) and uncertainty is captured by probabilistic sensitivity analysis. In their second study, Maetzel et al.  analyze the cost-effectiveness of NSAIDs compared to COX-2 inhibitors. GI events are modeled for OA or RA patients over 12 health states, each supplemented by a decision tree reflecting the treatment pattern following the specific health state. The ICER is calculated as costs per QALY from a third-party payer perspective and only a deterministic sensitivity analysis is performed.
The model by Spalding et al.  includes three different Markov-states. Patients start with a specific first-line or initial drug therapy and move to a pooled treatment state or to death from either of the two other states. For each of the two treatment states certain costs and a specific reduction in the HAQ-score are assumed per 12 month cycle. Again, no PSA was conducted.
Tanno et al.  are presenting the only model for Japan. Similar to the model family around Kielhorn et al. , the Markov-states represent the treatment options and patients move in six months cycles to the next treatment in sequence, if they fail an ACR20 response. The authors include age- and sex-specific excess mortality due to RA as an exponential function of the HAQ-score. Patients are also assumed to have failed an initial therapy with bucillamine, whereas the main interest is in the cost-effectiveness of etanercept. Only one-way sensitivity analyses were performed and reported as tornado-plot.
The first intention of Marra et al.  lies in the cost-effectiveness of infliximab in combination with methotrexate compared to methotrexate monotherapy. Secondly, the authors also explore the influence of different quality of life questionnaires (HUI-2, HUI-3, SF-6D, EQ-5D) on health utility and therefore on cost-effectiveness. Within a time-horizon of 10 years patients move between 25 health states, as defined corresponding to HAQ-scores from zero to three in steps of 0.125 (similar to the BRAM approach described later) plus death in cycles of one week. PSA was performed and the differences between the different questionnaires are clearly demonstrated.
Nguyen et al.  present a similar but independently developed model to Spalding et al. . In the model, patients start with the treatment of one of five different TNF-α antagonists in the first health state and either move to the ACR50-response health state or a non-response health state. Patients in the non-responder health state are assumed to switch to tocilizumab as second-line treatment. Patients can die from all three health states. The cycle length of the model is three months and the time horizon is five years. One-way deterministic as well as probabilistic sensitivity analyses were applied.
The study of Bessette et al.  compares celecoxib as 1st, 2nd and 3rd line treatment option in addition to different NSAID-strategies. In similarity to other studies with a comparable research question, the health states consider mainly GI events supplemented by one health state for cardiovascular AEs. In the case of an AE, patients continue to the next treatment option in the sequence. Of all Markov-models Bessette et al. are using the shortest cycle length with one month (time-horizon: five years).
The German publication from Schӓdlich et al.  reports the results from an adaption of the "The Avara® interactive model", on which no further information or publication could be found. The cost-effectiveness of additional leflunomide in four different DMARD-treatment sequences is analyzed using ACR20-response to differentiate the health states. The model is run for six cycles of six months cycle length. The ICER for each sequence with leflunomide is calculated for the corresponding sequence without leflunomide.
Wu et al.  report the only cost-effectiveness model found for China comparing seven different treatment sequences. Health states are defined using ACR20/50/70 response criteria and the next treatment option in the particular sequence is selected after poor remission or adverse events. The results are presented as costs per QALY in a cost-effectiveness frontier. Deterministic as well as probabilistic sensitivity analyses are performed.
Individual sampling models
Brennan et al.  developed the first individual sampling model. Thereby, certain events for patients of a certain age, sex and HAQ-Score are calculated using fixed time intervals of six months for DMARDs and three months for etanercept, respectively. In detail, patients who failed two conventional DMARDs are assumed to receive a certain treatment which might lead to a response according to the ACR20 criteria, thus continuing treatment until loss of efficacy or an adverse event, no response or death. The responders may continue the treatment and non-responders move on to the next treatment in sequence. Costs and utilities for each patient are documented during the modeling process. Bansback et al.  describe a slightly modified model for Sweden, where adverse events and the loss of efficacy are explicitly included in the first cycle in which patients receive a new treatment. Brennan et al.  use their model from 2004 with data from the British Society for Rheumatology Biologics Registry (BSRBR) for 8,000 patients and Wailoo et al.  with data from the National Database of Rheumatic Diseases (NDB) for 10,000 US-patients. The latter model has been modified by Finckh et al.  to compare three different treatment strategies including DMARDs (conventional and biologic) and other pharmaceutical and non-pharmaceutical interventions.
In the individual sampling-model by Kielhorn et al.  patients start after inadequate response to two biological DMARDs and begin with states related to a certain treatment sequence. For each treatment patients may enter five different health states (including death) defined by the ACR-Score and move to the next treatment if they show no response. If they completed all treatment options patients move to a palliative care state. The time-horizon is 100 years with a cycle-length of 6 months. Merkesdal et al. , Hallinen et al.  and Soini et al.  are using the same model structure for country specific analyses. The analyses were conducted using 10,000 patients and 3,000 patients in both Finnish studies, respectively. Additionally, PSA was performed using 1,000 iterations in all models.
Vera-Llonch et al.  and Vera-Llonch et al.  lean their individual sampling models against Brennan et al. . Both studies are evaluating the cost-effectiveness of abatacept, using data from the clinical AIM or ATTAIN trials with methotrexate or an unspecified DMARD as comparator, respectively. Patients in the model who are assumed to start treatment with abatacept are considered to continue on treatment after 6 months of therapy if they show a HAQ-score improvement of at least 0.50, no side-effects (no cost and utility affects assumed), no co-morbid conditions and no surgery. The model analyses costs and QALYs of 1,000 female patients between 55 and 64 years over ten years and life-time with a cycle length of three months. Changes on HAQ-scores were considered as multiples of 0.125. The models were adopted by Yuan et al.  to analyze the cost-effectiveness of abatacept and rituximab.
The only non-adopted ISM was done by Diamantopoulus et al. . The model was developed to evaluate the cost-effectiveness of tocilizumab in addition to a sequence of DMARDs and TNF-inhibitors from an Italian third-party payer perspective. 10,000 patients are followed over a life-time. Treatment efficacy is measured via the ACR20/50/70 response criteria and patients continue to the next treatment option after no response. ICERs are calculated for all sequences and deterministic as well as probabilistic sensitivity analyzes are performed.
Discrete event simulations
Most DES models are developed for the United Kingdom and used within the Health Technology Assessment process. The first DES model has been developed by Jobanputra et al.  to assess the cost-effectiveness of infliximab and etanercept and is known as the Birmingham Preliminary Model (BPM). Within the model, the time of an individual patient on a certain DMARD in a sequence of DMARDs is calculated and all corresponding costs and health outcomes are documented. Afterwards, depending on the remaining lifetime of a patient, the patient moves to the next treatment option or dies. For some treatments effects of certain toxicities are included in the model which may lead to changes in the treatment cascade.
Barton et al.  extended the model by allowing patients to stop the current treatment not only by continuing with the next DMARD in sequence or death, but also by undergoing a joint replacement and/or an increase in the HAQ-Score. Thereby, in order to decide which event occurs and ends the treatment, the time of each event is calculated separately and the event with the shortest timeframe is taken into account. Changes in the HAQ-score are modeled with an interval of 0.125.
Clark et al.  are using the Birmingham Rheumatoid Arthritis Model (BRAM) from Barton et al. to evaluate the cost-effectiveness of anakinra. Chen et al.  also use the BRAM for adalimumab, etanercept and infliximab, although joint replacements are left out of both latter models. The later adaption by Chen also includes individual improvements on HAQ-scores and early withdrawal from a treatment by patients is considered. The last adaptation of the BRAM is used by Malottki et al.  which made further improvements to the model by allowing for a non-linear translation of HAQ-scores to utility values and enabling probabilistic sensitivity analysis.
With the work from Lindgren et al.  only one non-UK DES-model could be identified. The authors analyze the cost-effectiveness of rituximab as second biologic agent in different treatment sequences. As with the BPM and BRAM, patients may be on or off treatment or dead, although the evaluated treatment sequences only include TNF-alpha-antagonists thus being shorter than in other DES models. When patients are on treatment they are either in a high or low DAS28 state and change to other states by treatment discontinuation, re-initiation, changing disease activity or death. The model was filled with data from the REFLEX trial and the SSATG-database. The model has been adopted by Kobelt et al.  with no major changes to analyze costs and outcomes for one sequence.
The present review analyzed 58 studies reporting cost-effectiveness modeling exercises on RA. The number of models using decision trees, Markov or ISM methods does not vary significantly, with fewer models using DES, but no trends in time could be observed with respect to a preferential use of certain model-types. Nevertheless, a higher number of Markov models are reported as independently conducted CEAs, whereas decision trees, ISM and DES tended to be adopted more often. This might be due to the increasing complexity and the higher computational requirements of ISM and DES when compared to the relatively wide-spread use of Markov models in health economic research.
Besides Lindgren et al.  and Kobelt et al. , all DES studies primarily have been developed for the UK and were published as part of a health technology assessment (HTA), in which also models might be reviewed that are submitted by a pharmaceutical company to the National Institute for Clinical Excellence (NICE). The background given in the HTAs for the development of the Birmingham Preliminary Model (BPM)  and its updated versions – known as the Birmingham Rheumatoid Arthritis Model (BRAM) – is the ability of DES to more accurately depict comprehensive treatment sequences as they occur in RA patients (see Table 2) and a more flexible course of the disease. Barton et al.  also puts a strong emphasis on the timing on activities and how the BRAM deals with competing risks.  On the other hand an impending need for reliable data is stated, especially impact patterns on quality of life under certain treatments .
With regard to the differences in time and the different settings and perspectives, no valid direct comparisons between the results of the different modeling techniques could be performed and no publications could be identified performing a direct comparison between different methodologies with the same input data. However, two studies comparing DES and Markov models in HIV and depression could be found, respectively. Summarizing, Simpson et al.  report no considerable differences in the ICERs for a time-frame from five to ten years. However, as the DES model is capable to simulate a higher number of clinical parameters for HIV patients it seems to have a slightly better predictive ability than the Markov model, when the findings are compared to real life data. For depression, Le Lay et al.  compared a self-developed DES model with Markov models the authors identified in the literature. Though no ICERs are compared in the study, the authors also report the DES model to be able to incorporate more factors relevant to the treatment of depression than the Markov model, especially the patients' history and their attitude towards different treatment options.
These examples illustrate the advantages DES models might also provide in their further application for RA. Under the consideration of their weaknesses, an even higher need (and mostly lack thereof) for reliable data and the resulting increased incorporation of expert opinion or even assumptions , a wider spectrum of outcome parameters might be presented to regulators as well as other decision makers. This could be a benefit but also a disadvantage, depending on the willingness of individual decision makers to look at even more complex data sets. Nevertheless, in any case it will be possible to reduce the wider variety of results from a DES to the information depth that usually comes with a Markov model.
However, the selection of the most appropriate model type, and therefore which simplification of reality is chosen, is a lively discussion. Even as there are guidelines on model selection (see for example ,), following different guidelines may lead to different model types. Furthermore, questions on the relevance of certain parameters might be answered differently by modelers working on the same decision analytic problem, also leading to different approaches. Though simplicity is an advantage in models, if complex modeling techniques are available to tackle complex decision problems (e.g. life-time treatment sequences in RA) those model types should be applied . In this context, the consideration of patient heterogeneity in economic evaluations  further encourages the use of more complex modeling techniques like DES .
Besides the models following Russell et al. , which use costs per day of low-disease activity score (LDAS), all models reviewed use a costs per QALY approach. The use of QALYs as their main outcome is justified with recommendations of regulatory bodies or with better comparability with other studies. In this context the burden of the disease might be modeled more comprehensively by the inclusion of necessary non-pharmaceutical treatments and socioeconomic consequences of the disease. In future studies also structural limitations in the health care provision for RA patients may be implemented in DES models. For example, the implications of a decreasing number of rheumatologists might be included in a DES framework using its ability to simulate queues given a maximal capacity a rheumatologist is able to take care of. This might also be relevant for time consuming and resource blocking therapies, such as intravenous infusions including monitoring for several hours or physical examinations.
In conclusion, Markov models, ISM and DES are appropriate choices to model the cost-effectiveness of treatment options for RA patients, whereas decision trees are only able to produce results for short time horizons. Although ISM overcomes the lack of consideration of heterogeneity of RA patients in Markov (cohort) models, DES is able to produce even more detailed results. This might come at the cost of more assumptions needing to be incorporated in the model and decision makers should be aware of this. Therefore, extensive sensitivity analyses on input parameters as well as on structural components of the model and a clear declaration of data sources are necessary for a valid interpretation of the results.
The aim of this review was to provide an informative basis on the advantages and disadvantages of modeling techniques used to model RA and on the data sources for clinical parameters (see Additional files 1 and 2). If sufficient data are available, DES should be considered as the preferred modeling technique.
Barton P, Bryan S, Robinson S: Modelling in the economic evaluation of health care: selecting the appropriate approach. J Health Serv Res Policy 2004, 9(2):110–118. 10.1258/135581904322987535
Brennan A, Akehurst R: Modelling in health economic evaluation. What is its place? What is its value? Pharmacoeconomics 2000, 17(5):445–459. 10.2165/00019053-200017050-00004
Gabriel S, Michaud K: Epidemiological studies in incidence, prevalence, mortality, and comorbidity of the rheumatic diseases. Arthritis Res Ther 2009, 11(3):229. 10.1186/ar2669
Gonzalez A, Maradit Kremers H, Crowson CS, Nicola PJ, Davis JM, Therneau TM, Roger VL, Gabriel SE: The widening mortality gap between rheumatoid arthritis patients and the general population. Arthritis Rheum 2007, 56(11):3583–3587. 10.1002/art.22979
Hamilton K, Clair E: Tumour necrosis factor-alpha blockade_ a new era for effective management of rheumatoid arthritis. Expert Opin Pharmacother 2000, 1(5):1041–1052. 10.1517/14656518.104.22.1681
Anderson J, Caplan L, Yazdany J, Robbins ML, Neogi T, Michaud K, Saag KG, O'Dell JR, Kazi S: Rheumatoid arthritis disease activity measures: American College of Rheumatology recommendations for use in clinical practice. Arthritis Care Res 2012, 64(5):640–647. 10.1002/acr.21649
Arnett FC, Edworthy SM, Bloch DA, Mcshane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, Medsger TA, Mitchell DM, Neustadt DH, Pinals RS, Schaller JG, Sharp JT, Wilder RL, Hunder GG: The American rheumatism association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum 1988, 31(3):315–324. 10.1002/art.1780310302
Aletaha D, Neogi T, Silman AJ, Smolen J, Ward MM: 2010 Rheumatoid arthritis classification criteria: An American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum 2010, 62(9):2569–2581. 10.1002/art.27584
Welsing PMJ, van Gestel A, Swinkels H, Kiemeney L, van Riel P: The relationship between disease activity, joint destruction, and functional capacity over the course of rheumatoid arthritis. Arthritis Rheum 2001, 9: 2009–2017. 10.1002/1529-0131(200109)44:9<2009::AID-ART349>3.0.CO;2-L
Fries JF, Spitz P, Kraines RG, Holman HR: Measurement of patient outcome in arthritis. Arthritis Rheum 1980, 23(2):137–145. 10.1002/art.1780230202
Lubeck D: Health-related quality of life measurements and studies in rheumatoid arthritis. Am J Manag Care 2002, 8(9):811–820.
Riise T, Jacobsen BK, Gran JT, Haga H, Arnesen E: Total mortality is increased in rheumatoid arthritis. a 17-year prospective study. Clin Rheumatol 2001, 20(2):123–127. 10.1007/PL00011191
Akil M, Amos RS: ABC of Rheumatology: rheumatoid arthritis—II: treatment. BMJ 1995, 310(6980):652–654. 10.1136/bmj.310.6980.652
Singh JA, Furst DE, Bharat A, Curtis JR, Kavanaugh AF, Kremer JM, Moreland LW, O'Dell J, Winthrop KL, Beukelman T, Bridges SL, Chatham WW, Paulus HE, Suarez-almazor M, Bombardier C, Dougados M, Khanna D, King CM, Leong AL, Matteson EL, Schousboe JT, Moynihan E, Kolba KS, Jain A, Volkmann ER, Agrawal H, Bae S, Patkar NM, Saag KG: 2012 Update of the 2008 American College of Rheumatology recommendations for the use of disease-modifying antirheumatic drugs and biologic agents in the treatment of rheumatoid arthritis. Arthritis Care Res 2012, 64(5):625–639. 10.1002/acr.21641
Brennan A: Modelling the cost-effectiveness of etanercept in adults with rheumatoid arthritis in the UK. Rheumatology 2004, 43(1):62–72. 10.1093/rheumatology/keg451
Choi H, Seeger J, Kuntz K: A cost-effectiveness analysis of treatment options for methotrexate-naive rheumatoid arthritis. J Rheumatol 2002, 29(6):1156–1165.
Clark W, Jobanputra P, Barton P, Burls A: The clinical and cost-effectiveness of anakinra for the treatment of rheumatoid arthritis in adults: a systematic review and economic analysis. Health Technol Assess 2004, 8(18).
Siebert U, Alagoz O, Bayoumi AM, Jahn B, Owens DK, Cohen DJ, Kuntz KM: State-transition modeling: a report of the ISPOR-SMDM modeling good research practices task force-3. Value Health 2012, 15(6):812–820. 10.1016/j.jval.2012.06.014
Karnon J, Stahl J, Brennan A, Caro JJ, Mar J, Möller J: Modeling using discrete event simulation: a report of the ISPOR-SMDM modeling good research practices task force-4. Value Health 2012, 15(6):821–827. 10.1016/j.jval.2012.04.013
Moore A, Phillips C, Hunsche E, Pellissier J, Crespi S: Economic evaluation of etoricoxib versus non-selective NSAIDs in the treatment of osteoarthritis and rheumatoid arthritis patients in the UK. Pharmacoeconomics 2004, 22(10):643–660. 10.2165/00019053-200422100-00003
Spiegel BMR, Targownik L, Dulai GS, Gralnek IM: The cost-effectiveness of cyclooxygenase-2 selective inhibitors in the management of chronic arthritis. Ann Intern Med 2003, 138(10):795–806. 10.7326/0003-4819-138-10-200305200-00007
Zabinski RA, Burke TA, Johnson J, Lavoie F, Fitzsimon C, Tretiak R, Chancellor JV: An economic model for determining the costs and consequences of using various treatment alternatives for the management of arthritis in Canada. Pharmacoeconomics 2001, 19(Suppl 1):49–58. 10.2165/00019053-200119001-00004
Kavanaugh A, Heudebert G, Cush J, Jain R: Cost valuation of novel therapeutics in rheumatoid arthritis (CENTRA): a decision analysis model. Semin Arthritis Rheum 1996, 25: 297–307. 10.1016/S0049-0172(96)80016-3
Choi H, Seeger J, Kuntz K: A cost-effectiveness analysis of treatment options for patients with Methotrexate-resistant rheumatoid arthritis. Arthritis Rheum 2000, 43(10):2316–2327. 10.1002/1529-0131(200010)43:10<2316::AID-ANR20>3.0.CO;2-6
Russell A, Beresniak A, Bessette L, Haraoui B, Rahman P, Thorne C, Maclean R, Dupont D: Cost-effectiveness modeling of abatacept versus other biologic agents in DMARDs and anti-TNF inadequate responders for the management of moderate to severe rheumatoid arthritis. Clin Rheumatol 2009, 28(4):403–412. 10.1007/s10067-008-1060-4
Beresniak A, Ariza-Ariza R, Garcia-Llorente JF, Gossec L, Goupille P, Saraux A, Bamberger M, Bregman B, Dupont D: Modelling cost-effectiveness of biologic treatments based on disease activity scores for the management of rheumatoid arthritis in Spain. Int J Inflam 2011, 2011.
Cimmino MA, Leardini G, Salaffi F, Intorcia M, Bellatreccia A, Dupont D, Beresniak A: Assessing the cost-effectiveness of biologic agents for the management of moderate-to-severe rheumatoid arthritis in anti-TNF inadequate responders in Italy: a modelling approach. Clin Exp Rheumatol 2011, 29(4):633–641.
Puolakka K, Blafield H, Kauppi M, Luosujӓrvi R, Peltoma R, Leikola-Pelho T, Sennfalt K, Beresniak A: Cost-effectiveness modelling of sequential biologic strategies for the treatment of moderate to severe rheumatoid arthritis in Finland. Open Rheumatol J 2012, 6: 38–43. 10.2174/1874312901206010038
Saraux A, Gossec L, Goupille P, Bregman B, Boccard E, Dupont D, Beresniak A: Cost-effectiveness modelling of biological treatment sequences in moderate to severe rheumatoid arthritis in France. Rheumatology 2010, 49(4):733–740. 10.1093/rheumatology/kep434
Beresniak A, Baerwald C, Zeidler H, Kruger K, Neubauer AS, Dupont D, Merkesdal S: Cost-effectiveness simulation model of biologic strategies for treating to target rheumatoid arthritis in Germany. Clin Exp Rheumatol 2013, 31(3):400–408.
Chiou C, Choi J, Reyes CM: Cost-effectiveness analysis of biological treatments for rheumatoid arthritis. Expert Rev Pharmacoeconomics Outcomes Res 2004, 4(3):307–315. 10.1586/1473722.214.171.1247
Kobelt G: The cost-effectiveness of infliximab (Remicade(R)) in the treatment of rheumatoid arthritis in Sweden and the United Kingdom based on the ATTRACT study. Rheumatology 2003, 42(2):326–335. 10.1093/rheumatology/keg107
Kobelt G, Eberhardt K, Jönsson L, Jönsson B: Economic consequences of the progression of rheumatoid arthritis in Sweden. Arthritis Rheum 1999, 42(2):347–356. 10.1002/1529-0131(199902)42:2<347::AID-ANR18>3.0.CO;2-P
Kobelt G, Jönsson L, Lindgren P, Young A, Eberhardt K: Modeling the progression of rheumatoid arthritis: A two-country model to estimate costs and consequences of rheumatoid arthritis. Arthritis Rheum 2002, 46(9):2310–2319. 10.1002/art.10471
Kobelt G, Lindgren P, Lindroth Y, Jacobson L, Eberhardt K: Modelling the effect of function and disease activity on costs and quality of life in rheumatoid arthritis. Rheumatology 2005, 44: 1169–1175. 10.1093/rheumatology/keh703
Kobelt G, Lindgren P, Singh A, Klareskog L: Cost effectiveness of etanercept (Enbrel) in combination with methotrexate in the treatment of active rheumatoid arthritis based on the TEMPO trial. Ann Rheum Dis 2005, 64(8):1174–1179. 10.1136/ard.2004.032789
Kobelt G, Lekander I, Lang A, Raffeiner B, Botsios C, Geborek P: Cost-effectiveness of etanercept treatment in early active rheumatoid arthritis followed by dose adjustment. Int J Technol Assess Health Care 2011, 27(03):193–200. 10.1017/S0266462311000195
Lekander I, Borström F, Lysholm J, van Vollenhoven R, Lindblad S, Geborek P, Kobelt G: The cost-effectiveness of TNF-inhibitors for the treatment of rheumatoid arthritis in Swedish clinical practice. Eur J Health Econ 2012, 14.
Lekander I, Borgstrom F, Svarvar P, Ljung T, Carli C, van Vollenhoven RF: Cost-effectiveness of real-world infliximab use in patients with rheumatoid arthritis in Sweden. Int J Technol Assess Health Care 2010, 26(1):54–61. 10.1017/S0266462309990596
Lekander I, Kobelt G, Svarvar P, Ljung T, van Vollenhoven R, Borgstrom F: The comparison of trial data-based and registry data-based cost-effectiveness of infliximab treatment for rheumatoid arthritis in Sweden using a modeling approach. Value Health 2013, 16(2):251–258. 10.1016/j.jval.2012.11.002
Welsing PMJ, Severens JL, Hartman M, van Riel PLCM: Modeling the 5-year cost effectiveness of treatment strategies including tumor necrosis factor-blocking agents and leflunomide for treating rheumatoid arthritis in the Netherlands. Arthritis Rheum 2004, 51(6):964–973. 10.1002/art.20843
Schipper LG, Kievit W, den Broeder AA, van der Laar MA, Adang EM, Fransen J, van Riel PLCM: Treatment strategies aiming at remission in early rheumatoid arthritis patients: starting with methotrexate monotherapy is cost-effective. Rheumatology 2011, 50(7):1320–1330. 10.1093/rheumatology/ker084
Wong JB, Singh G, Kavanaugh A: Estimating the cost-effectiveness of 54 weeks of infliximab for rheumatoid arthritis. Am J Med 2002, 113: 400–408. 10.1016/S0002-9343(02)01243-3
Barbieri M, Wong J, Drummond M: The cost-effectiveness of infliximab for severe treatment-resistant rheumatoid arthritis in the UK. Pharmacoeconomics 2005, 23(6):607–618. 10.2165/00019053-200523060-00007
Wong JB, Ramey DR, Singh G: Long-term morbidity, mortality, and economics of rheumatoid arthritis. Arthritis Rheum 2001, 44(12):2746–2749. 10.1002/1529-0131(200112)44:12<2746::AID-ART461>3.0.CO;2-Z
Bae S, Corzillius M, Kuntz KM, Liang MH: Cost-effectiveness of low dose corticosteroids versus non-steroidal anti-inflammatory drugs and COX-2 specific inhibitors in the long-term treatment of rheumatoid arthritis. Rheumatology (Oxford) 2003, 42(1):46–53. 10.1093/rheumatology/keg029
Yun HR, Bae S: Cost-effectiveness analysis of NSAIDs, NSAIDs with concomitant therapy to prevent gastrointestinal toxicity, and COX-2 specific inhibitors in the treatment of rheumatoid arthritis. Rheumatol Int 2005, 25(1):9–14. 10.1007/s00296-003-0392-2
Maetzel A, Krahn M, Naglie G: The cost effectiveness of rofecoxib and celecoxib in patients with osteoarthritis or rheumatoid arthritis. Arthritis Rheum 2003, 49(3):283–292. 10.1002/art.11121
Maetzel A, Strand V, Tugwell P, Wells G, Bombardier C: Cost effectiveness of adding leflunomide to a 5-year strategy of conventional disease-modifying antirheumatic drugs in patients with rheumatoid arthritis. Arthritis Rheum 2002, 47(6):655–661. 10.1002/art.10793
Spalding J, Hay J: Cost-effectiveness of tumour necrosis factor-alpha inhibitors as first-line agents in rheumatoid arthritis. Pharmacoeconomics 2006, 24(12):1221–1232. 10.2165/00019053-200624120-00006
Tanno M, Nakamura I, Ito K, Tanaka H, Ohta H, Kobayashi M, Tachihara A, Nagashima M, Yoshino S, Nakajima A: Modeling and cost-effectiveness analysis of etanecerpt in adults with rheumatoid arthritis in Japan: a preliminary analysis. Mod Rheumatol 2006, 16: 77–84. 10.3109/s10165-006-0461-y
Marra CA, Marion SA, Guh DP, Najafzadeh M, Wolfe F, Esdaile JM, Clarke AE, Gignac MA, Anis AH: Not all "quality-adjusted life years" are equal. J Clin Epidemiol 2007, 60(6):616–624. 10.1016/j.jclinepi.2006.09.006
Nguyen C, Bounthavong M, Mendes M, Christopher ML, Tran JN, Kazerooni R, Morreale AP: Cost utility of tumour necrosis factor-alpha inhibitors for rheumatoid arthritis. an application of bayesian methods for evidence synthesis in a markov model. Pharmacoeconomics 2012, 30(7):575–593. 10.2165/11594990-000000000-00000
Bessette L, Risebrough N, Mittmann N, Roussy JP, Ho J, Zlateva G: Cost-utility of celecoxib use in different treatment strategies for osteoarthritis and rheumatoid arthritis from the Quebec healthcare system perspective. J Med Econ 2009, 12(3):246–258. 10.3111/13696990903288970
Schadlich PK, Zeidler H, Zink A, Gromnica-Ihle E, Schneider M, Straub C, Brecht JG, Huppertz E: Contribution of leflunomide to the cost effectiveness of sequential DMARD therapy of rheumatoid arthritis in Germany (Wirtschaftlichkeit von Leflunomid bei sequentieller Basistherapie der rheumatoiden Arthritis in Deutschland). Z Rheumatol 2004, 63(1):59–75. 10.1007/s00393-004-0570-y
Wu B, Wilson A, Wang F, Wang S, Wallace D, Weisman M, Lu L: Cost effectiveness of different treatment strategies in the treatment of patients with moderate to severe rheumatoid arthritis in china. PLoS One 2012, 7(10):e47373. 10.1371/journal.pone.0047373
Bansback N, Brennan A, Ghatnekar O: Cost effectiveness of adalimumab in the treatment of patients with moderate to severe rheumatoid arthritis in Sweden. Ann Rheum Dis 2005, 64(7):995–1002. 10.1136/ard.2004.027565
Brennan A, Bansback N, Nixon R, Madan J, Harrison M, Watson K, Symmons D: Modelling the cost effectivenessof TNF-alpha antagonists in the management of rheumatoid arthritis: results from the British Society of Rheumatology Biologics Registry. Rheumatology 2007, 46: 1345–1354. 10.1093/rheumatology/kem115
Wailoo AJ, Bansback N, Brennan A, Michaud K, Nixon RM, Wolfe F: Biologic drugs for rheumatoid arthritis in the medicare program: A cost-effectiveness analysis. Arthritis Rheum 2008, 58(4):939–946. 10.1002/art.23374
Finckh A, Bansback N, Marra CA, Anis AH, Michaud K, Lubin S, White M, Sizto S, Liang MH: Treatment of very early rheumatoid arthritis with symptomatic therapy, disease-modifying antirheumatic drugs, or biologic agentsa cost-effectiveness analysis. Ann Intern Med 2009, 151(9):612–621. 10.7326/0003-4819-151-9-200911030-00006
Kielhorn A, Porter D, Diamantopoulos A, Lewis G: UK cost-utility analysis of rituximab in patients with rheumatoid arthritis that failed to respond adequately to a biologic disease-modifying antirheumatic drug. Curr Med Res Opin 2008, 24(9):2639–2650. 10.1185/03007990802321683
Merkesdal S, Kirchhoff T, Wolka D, Ladinek G, Kielhorn A, Rubbert-Roth A: Cost-effectiveness analysis of rituximab treatment in patients in Germany with rheumatoid arthritis after etanercept-failure. Eur J Health Econ 2010, 11(1):95–104. 10.1007/s10198-009-0205-y
Hallinen TA, Soini EJO, Eklund K, Puolakka K: Cost-utility of different treatment strategies after the failure of tumour necrosis factor inhibitor in rheumatoid arthritis in the Finnish setting. Rheumatology 2010, 49(4):767–777. 10.1093/rheumatology/kep425
Soini E, Hallinen T, Puolakka K, Vihervaara V, Kauppi MJ: Cost-effectiveness of adalimumab, etanecerpt, and tocilizumab as first-line treatments for moderate-to-severe rheumatoid arthritis. J Med Econ 2012, 15(2):340–351. 10.3111/13696998.2011.649327
Vera-Llonch M, Massarotti E, Wolfe F, Shadick N, Westhovens R, Sofrygin O, Maclean R, Yuan Y, Oster G: Cost-effectiveness of abatacept in patients with moderately to severly active rheumatoid arthritis and inadequate response to methotrexate. Rheumatology 2008, 47: 535–541. 10.1093/rheumatology/ken007
Vera-Llonch M, Massarotti E, Wolfe F, Shadick N, Westhovens R, Sofrygin O, Maclean R, Yuan Y, Oster G: Cost-effectiveness of abatacept in patients with moderately to severly active rheumatoid arthritis and inadequate response to tumor necrosis factor-alpha antagonists. J Rheumatol 2008, 35(9):1–9.
Yuan Y, Trivedi D, Maclean R, Rosenblatt L: Indirect cost-effectiveness analyses of abatacept and rituximab in patients with moderate-to-severe rheumatoid arthritis in the United States. J Med Econ 2010, 13(1):33–41. 10.3111/13696990903508021
Diamantopoulos A, Benucci M, Capri S, Berger W, Wintfeld N, Giuliani G, Ricciardi W: Economic evaluation of tocilizumab combination in the treatment of moderate-to-severe rheumatoid arthritis in Italy. J Med Econ 2012, 15(3):576–585. 10.3111/13696998.2012.665110
Jobanputra P, Barton P, Bryan S, Burls A: The effectiveness of infliximab and etanercept for the treatment of rheumatoid arthritis: a systematic review and economic evaluation. Health Technol Assess 2002, 6(21).
Barton P, Jobanputra P, Wilson J, Bryan S, Burls A: The use of modelling to evaluate new drugs for patients with a chronic condition: the case of antibodies against tumour necrosis factor in rheumatoid arthritis. Health Technol Assess 2004, 8(11).
Chen Y-F, Jobanputra P, Barton P, Jowett S, Bryan S, Clark W, Fry-Smith A, Burls A: A systematic review of the effectiveness of adalimumab, etanercept and infliximab for the treatment of rheumatoid arthritis in adults and an economic evaluation of their cost-effectiveness. Health Technol Assess 2006, 10(42).
Malottki K, Barton P, Tsourapas A, Uthman AO, Liu Z, Routh K, Connock M, Jobanputra P, Moore D, Fry-Smith A, Chen YF: Adalimumab, etanercept, infliximab, rituximab and abatacept for the treatment of rheumatoid arthritis after the failure of a tumour necrosis factor inhibitor: a systematic review and economic evaluation. Health Technol Assess 2011, 15(14).
Lindgren P, Geborek P, Kobelt G: Modeling the cost-effectiveness of treatment of rheumatoid arthritis with rituximab using registry data from Southern Sweden. Int J Technol Assess Health Care 2009, 25(2):181–189. 10.1017/S0266462309090230
Kobelt G, Lindgren P, Geborek P: Costs and outcomes for patients with rheumatoid arthritis treated with biological drugs in Sweden: a model based on registry data. Scand J Rheumatol 2009, 38(6):409–418. 10.3109/03009740902865464
Simpson K, Strassburger A, Jones W, Dietz B, Rajagopalan R: Comparison of Markov model and discrete-event simulation techniques for HIV. Pharmacoeconomics 2009, 27(2):159–165. 10.2165/00019053-200927020-00006
Le Lay A, Despiegel N, Francois C, Duru G: Can discrete event simulation be of use in modelling major depression? Cost Eff Resour Alloc 2006, 4(1):19. 10.1186/1478-7547-4-19
Stahl J: Modelling methods for pharmacoeconomics and health technology assessment. an overview and guide. Pharmacoeconomics 2008, 26(2):131–148. 10.2165/00019053-200826020-00004
Grutters JC, Sculpher M, Briggs A, Severens J, Candel M, Stahl J, de Ruysscher D, Boer A, Ramaekers BT, Joore M: Acknowledging patient heterogeneity in economic evaluation. Pharmacoeconomics 2013, 31(2):111–123. 10.1007/s40273-012-0015-4
Caro JJ: Pharmacoeconomic analyses using discrete event simulation. Pharmacoeconomics 2005, 23(4):323–332. 10.2165/00019053-200523040-00003