Impact of adherence to biological agents on health care resource utilization for patients over the age of 65 years with rheumatoid arthritis.

OBJECTIVE: Poor adherence to therapy increases the patient and societal burden and complexity of chronic diseases such as rheumatoid arthritis (RA). In the past 15 years, biologic disease-modifying anti-rheumatic drugs (DMARDs) have revolutionized the treatment of RA. However, little data are available on the impact of adherence to biologics on health care resources. The objective of the study was to determine the long-term health care resource utilization patterns of RA patients who were adherent to biologic DMARD therapy compared to RA patients who were non-adherent to biologic DMARD therapy in an Ontario population and to determine factors influencing adherence.
METHODS: Patients were identified from the Ontario RA Database that contains all RA patients in Ontario, Canada, identified since 1991. The study population included RA patients, aged 65+ years, with a prescription for a biologic DMARD between 2003 and 2013. Exclusion criteria included diagnosis of inflammatory bowel disease, psoriatic arthritis or psoriasis in the 5 years prior to the index date and discontinuation of biologic DMARD, defined as no subsequent prescription during the 12 months after the index date. Adherence was defined as a medication possession ratio of ≥0.8 measured as the proportion of days for which a patient had biologic treatment(s) over a defined follow-up period. Adherent patients were matched to non-adherent patients by propensity score matching.
RESULTS: A total of 4,666 RA patients were identified, of whom 2,749 were deemed adherent and 1,917 non-adherent. The age (standard deviation) was 69.9 (5.46) years and 75% were female. Relative rates for resource use (physician visits, emergency visits, hospitalization, home care and rehabilitation) for the matched cohort were significantly lower (P⩽0.0001) in adherent patients. Non-adherent patients' use of oral prednisone (67%) was significantly higher (P⩽0.001) than that of the adherent cohort (56%).
CONCLUSION: RA patients adherent to biologic therapy have lower health care resource use and lower steroid use compared to non-adherent patients.

Introduction

Rheumatoid arthritis (RA) is a chronic autoimmune inflammatory disease characterized by persistent and progressive joint inflammation and destruction, resulting in pain, stiffness and functional disability. The primary presenting symptoms of RA include pain, weakness, stiffness, tenderness and swelling of the joints. The chronic inflammation leads to cartilage destruction, bone erosions and the formation of chronic granulation, scar tissue and bone fusions resulting in deformities of the hands, feet and other joints. Over time, and without appropriate treatment, this leads to loss of function and disability.

Conventional treatment of RA with traditional disease-modifying anti-rheumatic drugs (DMARDs) such as methotrexate (MTX) is recommended as first-line management of patients with persistent synovitis.1,2 During initial treatment with DMARD therapy, patients are often prescribed glucocorticoids or non-steroidal anti-inflammatory drugs (NSAIDs) as adjunctive therapy in addition to their DMARD(s) to help manage flares and symptoms if no other options are available. After failure with traditional DMARD therapy, treatment with a biologic DMARD is then often initiated.1,2

Adherence to prescribed drug therapies is necessary to prevent irreversible joint damage.3,4 However, studies indicate that adherence rate in RA patients varied between 49.5% and 98.5% depending on the definition and method used.4 Most of these studies evaluated DMARD therapy, and only a small number of them investigated adherence to anti-tumor necrosis factor (anti-TNF) therapies.5–7

The primary objective of this study was to determine the long-term health care resource utilization patterns of RA patients who were adherent to biologic DMARD therapy compared to RA patients who were non-adherent to biologic DMARD therapy in an Ontario population and to determine factors influencing adherence.

Methods

Administrative databases

This was a retrospective observational database study utilizing de-identified linked databases managed by the Institute for Clinical Evaluative Sciences (ICES). As part of a pilot project at the ICES, analyses were performed by ICES scientists on ICES administrative data to study the research question. The heath care resource outcomes measured included visits to general practitioner or family practitioner (GP/FP), rheumatologists, cardiologists and other medical doctors (MD), emergency department (ED) visits, hospitalizations and length of stay, same-day surgeries, home care services and rehabilitation services.

The following administrative databases were used in this study: The Continuing Care Reporting System – Long-Term Care (CCRS LTC), which contains demographic, clinical, functional and resource utilization information on individuals receiving continuing care services in long-term care homes; The Discharge Abstract Database (DAD), which contains patient-level data for acute, rehabilitation, chronic and day surgery institutions in Ontario; The Home Care Database (HCD), which captured all services provided by or coordinated by Ontario’s Community Care Access Centers (CCACs); The National Ambulatory Care Reporting System (NACRS)/Same Day Surgery (SDS), which captures information on patient visits to hospital and community-based ambulatory care: day surgery, outpatient clinics and EDs; The National Rehabilitation Reporting System (NRS), which contains client data collected from participating adult inpatient rehabilitation facilities and programs; The Ontario Drug Benefit (ODB), which contains claims for prescription drugs received under the ODB program; The Ontario Health Insurance Plan (OHIP), which contains most claims paid by the OHIP (the data cover all health care providers who can claim under OHIP – this includes physicians, groups, laboratories and out-of-province providers); The Postal Code Conversion File (PCCF), which links six-character postal codes to standard geographic areas such as dissemination areas, census tracts and census subdivisions; The Registered Persons Data Base files (RPDB), which provided basic demographic information about anyone who has ever received an Ontario health card number and finally, The Ontario Rheumatoid Arthritis Database (ORAD), which is an ICES-derived validated population-based cohort that contains all Ontario RA patients identified since 1991. In the ORAD database, incident cases were identified from 1995/1996. Cases were defined by having at least one inpatient Canadian Institute for Health Information (CIHI) DAD diagnosis code (any type) for RA (International Classification of Diseases [ICD]-9 714; ICD-10 M05, M06) or three physician OHIP claims with a diagnosis code for RA (ICD-9 714) over a 2-year period, with at least one claim by a musculoskeletal specialist (according to ICES Physician Database MAIN SPECIALTY: Rheumatology, Orthopedic surgery, or Internal medicine).

The administrative databases employed in the study rely on coded data. Personal identifiers are not provided to researchers. Therefore, patient privacy was not a concern in this study.

Cohort definition

The study included all subjects previously diagnosed with RA and aged 65 years and older from the ORAD database, with a prescription for a biologic DMARD between January 1, 2003, and December 31, 2013. The index date was the first prescription for a biologic DMARD (Table 1) between January 1, 2003, and December 31, 2013. The exclusions were a diagnosis of Crohn’s disease, psoriatic arthritis, psoriasis or ulcerative colitis in the 5 years prior to the index date or the discontinuation of biologic drug defined as no subsequent prescription for any of the biologic drugs of interest during the 12 months after the index date.

Table 1

Biologic treatments of interest

Brand name	Generic name	DIN	Dose and dose frequency
Anti-TNF agents
Humira	Adalimumab	02258595	40 mg every 2 weeks
		09857294
		09857326
		09857327
Enbrel®	Etanercept	02242903	25 mg twice weekly or 50 mg once weekly
		02274728
		09857322
		09857394
Remicade®	Infliximab	02244016	3 mg/kg/dose at 0 weeks, 2 weeks and 6 weeks followed by maintenance therapy of
		09852956	3 mg/kg/dose every 8 weeks up to a maximum of six maintenance doses per year
Cimzia®	Certolizumab	02331675	400 mg at 0 weeks, 2 weeks and 4 weeks followed by maintenance therapy of 200 mg every 2 weeks. For maintenance dosing, 400 mg every 4 weeks may be considered
Simponi®	Golimumab	02324776	50 mg once a month
		02324784
Other mechanism of action
Kineret®	Anakinra	02245913	100 mg per day
Actemra®	Tocilizumab	023500920235010602350114	4 mg/kg/dose once every 4 weeks followed by an increase to 8 mg/kg/dose based on clinical response. For individuals whose body weight is ⩾100 kg, doses exceeding 800 mg per infusion are not recommended
Orencia®	Abatacept	02282097	IV: Year 1: 750 mg at weeks 0, 2 and 4 then every 4 weeks
		02402475	SC: 125 mg weekly
Rituxan®	Rituximab	0224192709857407	One course of treatment is 1,000 mg followed 2 weeks later by the second 1,000 mg dose. Two courses will be approved each year (courses should be at least 6 months)

Abbreviations: DIN, drug identification number; TNF, tumor necrosis factor; IV, intravenously; SC, subcutaneously.

Variable definition

Adherence was defined as the proportion of days for which a patient had biologic treatment(s) over their total overall follow-up period (Table 1 for the list of biologics included). Adherence was determined using the medication possession ratio (MPR): MPR = (Days supplied for the drug of interest/Follow-up time). Patients were considered adherent if they had an MPR of 0.8 or higher, a commonly used measure in studies on adherence.3,4,8 Days’ supply was captured in the ODB database. Additional variables reported at baseline or at index date were the year of index date, age, sex, income status, rural status and the resource utilization band (RUB). Further, heath care resource utilization in the year prior to the index date was reported, which included: GP/FP, rheumatologists, cardiologists and other MD visits, ED visits, hospitalizations and length of stay, same-day surgeries, home care services and rehabilitation services. In order to further characterize the study cohort, prior use (1 year prior to index date) of NSAIDs/COX-2 inhibitors (COX-2), corticosteroids, DMARDs and anti-hypertensives was also reported.

Statistical analyses

Multivariate logistic regression was used to identify factors that influence adherence. The equation is in the format: where patient adherence was the dependent variable and baseline characteristics at index and previous resource use in the year prior to index date were the independent variables. The multivariate logistic regression analysis was used to determine the adjusted odds ratios (ORs) for each variable, adjusted for all other patient characteristics. For example, log(b1) would capture the impact of age on determining adherence, adjusting for sex, income quintile and other variables. The variables evaluated in the model included age at index date, sex, income status at index date, rural status at index date, long-term care resident status at index date and the RUB at index date. Further, health care resource utilization in 1 year prior to the index date was included in the model as well, which included GP/FP, rheumatologists, cardiologists and other MD visits, ED visits, hospitalizations and length of stay, same-day surgeries, home care services and rehabilitation services.

To capture differences in the health care resource use of adherent and non-adherent patients, adherent and non-adherent patients were matched on a propensity score (PS) calculated using baseline characteristics including age, sex, income quintile, Local Health Integration Network of residence, rurality, long-term care residence, RUB status, index year of first biologic script and health services use 1 year prior to index date (visits to GP/FP, rheumatologists, cardiologists and other MDs, ED visits, hospitalizations, same-day surgeries, home care services and rehabilitation episodes). The PS model included all the 17 variables that were available in the data set since there were differences observed in variables between the adherent and the non-adherent cohort at baseline (Table 2). Specifically, it was important to balance the cohorts based on index year since the proportion of adherent to non-adherent patients varied depending on the index year. Further, the inclusion of all available variables would ensure that all potential confounders were corrected for. Once patients were matched, the unadjusted rate of each health care resource per person per 100 days was captured over each patient’s total follow-up time. The relative rates capturing the differences in the rate of each health care resource use was then modeled assuming that outcomes followed a negative binomial distribution.

Table 2

Baseline characteristics

Characteristics	Overall	Adherent (MPR ≥0.8)	Non-adherent (MPR ⩽0.8)	P-value
	N=4,667	n=2,196	n=2,471
Index year, n (%)				⩽0.001
2003	456 (9.8)	149 (6.8)	307 (12.4)
2004	237 (5.1)	78 (3.6)	159 (6.4)
2005	283 (6.1)	97 (4.4)	186 (7.5)
2006	279 (6)	97 (4.4)	182 (7.4)
2007	327 (7)	135 (6.1)	192 (7.8)
2008	406 (8.7)	166 (7.6)	240 (9.7)
2009	451 (9.7)	200 (9.1)	251 (10.2)
2010	349 (7.5)	157 (7.1)	192 (7.8)
2011	524 (11.2)	291 (13.3)	233 (9.4)
2012	658 (14.1)	385 (17.5)	273 (11)
Long-term care resident, n (%)	27 (0.6)	13 (0.7)	14 (0.5)	0.454
Demographics
Age, mean (SD)	69.94 (5.46)	70.60 (5.69)	69.47 (5.24)	⩽0.001
Female, n (%)	3,500 (75.0)	1,453 (75.8)	2,047 (74.5)	0.301
Income quintile, n (%)				0.211
Missing	16 (0.3)	8 (0.4)	8 (0.3)
1 – lowest	744 (15.9)	323 (16.8)	421 (15.3)
2	892 (19.1)	385 (20.1)	507 (18.4)
3	946 (20.3)	386 (20.1)	560 (20.4)
4	1,055 (22.6)	405 (21.1)	650 (23.6)
5 – highest	1,013 (21.7)	410 (21.4)	603 (21.9)
Geographics
Rural, n (%)	706 (15.1)	354 (16.1)	352 (14.2)	0.074
Health resource utilization
RUB, n (%)				⩽0.001
0–2	11 (0.2)	6–10*	1–5*
3	2,296 (49.2)	864 (45.1)	1,432 (52.1)
4	1,580 (33.9)	687 (35.8)	893 (32.5)
5	779 (16.7)	356–360*	419–423*
MD visits
GP/FP
Any visit, n (%)	4,562 (97.8)	1,873 (97.7)	2,689 (97.8)	0.798
Mean (SD)	12.7 (12.57)	13.5 (12.57)	12.1 (12.54)	⩽0.001
Median (IQR)	9 (5–16)	10 (6–17)	9 (5–15)	⩽0.001
Rheumatologist
Any visit, n (%)	2,682 (57.5)	1,045 (54.5)	1,637 (59.5)	⩽0.001
Mean (SD)	3.7 (4.89)	3.5 (4.89)	3.9 (4.89)	0.017
Median (IQR)	2 (0–6)	1 (0–6)	2 (0–6)	0.017
Cardiologist
Any visit, n (%)	1,144 (24.5)	512 (26.7)	632 (23)	0.004
Mean (SD)	0.6 (1.99)	0.7 (2.44)	0.5 (1.6)	0.001
Median (IQR)	0 (0–0)	0 (0–1)	0 (0–0)	0.001
Other
Any visit, n (%)	4,589 (98.3)	1,891 (98.6)	2,698 (98.1)	0.188
Mean (SD)	12.8 (10.25)	13.5 (11.52)	12.3 (9.22)	⩽0.001
Median (IQR)	11 (6–17)	11 (6–18)	10 (6–17)	⩽0.001
ED visits
Any visit, n (%)	1,706 (36.6)	777 (40.5)	929 (33.8)	⩽0.001
Mean (SD)	0.7 (1.42)	0.9 (1.6)	0.7 (1.27)	⩽0.001
Median (IQR)	0 (0–1)	0 (0–1)	0 (0–1)	⩽0.001
Same-day surgeries
Any visit, n (%)	1,180 (25.3)	502 (26.2)	678 (24.7)	0.239
Mean (SD)	0.4 (0.76)	0.4 (0.77)	0.4 (0.75)	0.521
Median (IQR)	0 (0–1)	0 (0–1)	0 (0–0)	0.521
Hospitalizations
Any visit, n (%)	845 (18.1)	410 (21.4)	435 (15.8)	⩽0.001
Mean (SD)	0.2 (0.59)	0.3 (0.64)	0.2 (0.55)	⩽0.001
Median (IQR)	0 (0–0)	0 (0–0)	0 (0–0)	⩽0.001
LOS
Mean (SD)	9.1 (12.18)	9.1 (11.6)	9.1 (12.72)	0.96
Median (IQR)	5 (3–10)	6 (3–10)	5 (2–10)	0.96
Home care
Any service, n (%)	338 (7.2)	170 (8.9)	168 (6.1)	⩽0.001
Mean (SD)	7.7 (47.32)	10.3 (57.08)	5.8 (39.01)	0.001
Median (IQR)	0 (0–0)	0 (0–0)	0 (0–0)	0.001
Rehabilitation
Any rehabilitation, n (%)	481 (10.3)	206 (10.7)	275 (10)	0.412
Mean (SD)	0.1 (0.47)	0.1 (0.49)	0.1 (0.45)	0.351
Median (IQR)	0 (0–0)	0 (0–0)	0 (0–0)	0.351

Note:

Suppressed for privacy.

Abbreviations: MPR, medication possession ratio; SD, standard deviation; RUB, resource utilization band; MD, medical doctor; GP, general practitioner; FP, family practitioner; IQR, interquartile range; ED, emergency department; LOS, length of stay.

As a sensitivity analysis, three separate matched cohorts capturing patients with 3 years, 5 years and 7 years of follow-up were developed (1,095 days, 1,825 days and 2,555 days of follow-up, respectively, from the date of the first claim for any of the biologic treatments of interest) and reported the mean and median health care resource utilization per 100 days for adherent and non-adherent patients.

Results

A total of 4,666 RA patients were identified, of whom 2,749 (59%) were deemed adherent and 1,917 (41%) non-adherent. The age (standard deviation [SD]) was 69.9 (5.46) years and 75% were female (Table 2). There were a greater proportion of non-adherent patients from 2003 to 2010, while a larger number of adherent patients were identified from 2011 onward. Adherent patients tended to be significantly younger and from a higher income quintile than non-adherent patients. In the year prior to the index date, there were a lower number of adherent patients in the higher RUB (P⩽0.001; Table 2). Indeed, adherent patients had a slightly lower number of GP visits (mean [SD] of 12.1 [12.5] vs 13.5 [12.6]; P⩽0.001), ED visits (mean [SD] of 0.7 [1.3] vs 0.9 [1.6]; P⩽0.001), hospitalizations (mean [SD] of 0.2 [0.5] vs 0.3 [0.6]; P⩽0.001) and home care visits (mean [SD] of 5.8 [39.0] vs 10.3 [57.1]; P=0.001). In contrast, non-adherent patients had a lower number of rheumatologists visits (mean [SD] of 3.5 [5.0] vs 3.9 [4.9]; P⩽0.017).

Multivariate logistic regression was used to identify factors that predict the overall adherence to biologic drugs (Table 3). Increased age (OR 0.96 [0.95–0.97]; P⩽0.0001), previous use of NSAIDs/COXIBs (OR 0.84 [0.75–0.94]; P=0.0035), use of oral prednisone (OR 0.76 [0.67–0.85]; P⩽0.0001) and use of anti-hypertensives (OR 0.74 [0.66–0.83]; P⩽0.0001) were significantly associated with non-adherence. In contrast, being in a higher income quintile of 4 (OR 1.23 [1.02–1.49]; P=0.033), living in south east Ontario (OR 1.72 [1.19–2.50]; P=0.004) and visits to the rheumatologist (OR 1.02 [1.00–1.03]; P⩽0.0078) were associated with adherence.

Table 3

Logistic regression results predicting overall adherence to biologic drugs among patients with RA in Ontario

	Odds ratio (95% CI)
Patient characteristics
Age	0.97 (0.96–0.98)
Female	0.93 (0.81–1.07)
Income quintile (REF =1 – lowest)
2	1.03 (0.84–1.26)
3	1.12 (0.91–1.36)
4	1.22 (1.0–1.49)
5 – highest	1.14 (0.94–1.39)
Rural	1.18 (0.98–1.43)
Long-term care resident	1.23 (0.55–2.75)
Drug use prior to index date
NSAIDs/COXIBs	0.83 (0.73–0.94)
Corticosteroids
Oral prednisone	0.89 (0.78–1.01)
Intra-articular steroid injections	0.96 (0.84–1.10)
Other	0.48 (0.19–1.25)
DMARDs
MTX	1.19 (1.05–1.36)
Others	1.16 (1.01–1.33)
Anti-hypertensives	0.84 (0.79–0.96)
Health services use prior to index date
RUB (REF =0–2)
3	2.81 (0.8–9.85)
4	2.42 (0.69–8.53)
5	2.63 (0.74–9.36)
Number of MD visits
GP/FP	1.00 (0.99–1.00)
Rheumatologist	1.02 (1.00–1.03)
Cardiologist	0.97 (0.94–1.01)
Others	1.00 (0.99–1.00)
Number of ED visits	0.96 (0.91–1.01)
Number of same-day surgeries	1.04 (0.96–1.13)
Number of hospitalizations	0.95 (0.83–1.08)
Number of home care services	1.00 (1.00–1.00)
Number of rehabilitation services	1.05 (0.91–1.20)

Note: Bold represent statistically significant values.

Abbreviations: RA, rheumatoid arthritis; REF, reference; NSAIDs, non-steroidal anti-inflammatory drugs; COXIBs, COX-2 inhibitors; DMARDs, disease-modifying anti-rheumatic drugs; MTX, methotrexate; RUB, resource utilization band; MD, medical doctor; GP, general practitioner; FP, family practitioner; ED, emergency department.

The impact of adherence on concomitant medication use is presented in Table 4. Since ODB only covers patients 65 years and older, the cohorts only included patients aged 66 years and older at the index date to allow the evaluation of drug use prior to the index date. In this cohort of 3,181 patients, 1,418 (45%) were labeled as non-adherent and 1,763 (55%) were labeled as adherent. Prior to the index date, a lower proportion of adherent patients were taking oral prednisone (52.4% vs 56.1%, P=0.033) and anti-hypertensives (61.5% vs 65.8%, P=0.012). In contrast, adherent patients were more likely to be taking MTX (66.8% vs 59.7%, P⩽0.001) or other DMARDs (79.4% vs 73.5%, P⩽0.001). After the index date, when patients were taking biologic therapy, adherent patients were less likely to be taking concomitant NSAIDs/COXIBs (53.5% vs 58.5%, P=0.006), corticosteroids (72.3% vs 79%, P⩽0.001) including oral prednisone (56.3% vs 67.3%, P⩽0.001) and intra-articular steroid injections (51.3% vs 55.4%, P=0.024) and anti-hypertensives (74.8% vs 80.7%, P⩽0.001). Concomitant non-biologic DMARD use was similar between the two groups.

Table 4

Medication use, among those aged 66 years and older, before and after the index date

Drug group, n (%)	Overall	Adherent (MPR ≥0.8)	Non-Adherent (MPR ⩽0.8)	P-value
	N=3,181	n=1,418	n=1,763
Prior to index date
NSAIDs/COXIBs	1,607 (50.5)	724 (51.1)	883 (50.1)	0.585
Corticosteroids
Any use	2,136 (67.1)	960 (67.7)	1,176 (66.7)	0.552
Oral prednisone	1,719 (54)	796 (56.1)	923 (52.4)	0.033
Intra-articular steroid injections	1,433 (45)	623 (43.9)	810 (45.9)	0.258
Others	16 (0.5)	11–15	1–5	0.051
DMARDs
Any use	2,869 (90.2)	1,258 (88.7)	1,611 (91.4)	0.012
MTX	2,024 (63.6)	846 (59.7)	1,178 (66.8)	⩽0.001
Others	2,442 (76.8)	1,042 (73.5)	1,400 (79.4)	⩽0.001
Anti-hypertensives	2,017 (63.4)	933 (65.8)	1,084 (61.5)	0.012
Index date to the end of follow-up
NSAIDs/COXIBs	1,773 (55.7)	829 (58.5)	944 (53.5)	0.006
Corticosteroids
Any use	2,394 (75.3)	1,120 (79)	1,274 (72.3)	⩽0.001
Oral prednisone	1,947 (61.2)	955 (67.3)	992 (56.3)	⩽0.001
Intra-articular steroid injections	1,690 (53.1)	785 (55.4)	905 (51.3)	0.024
Others	109 (3.4)	65 (4.6)	44 (2.5)	0.001
DMARDs
Any use	2,656 (83.5)	1,174 (82.8)	1,482 (84.1)	0.338
MTX	1,913 (60.1)	832 (58.7)	1,081 (61.3)	0.13
Others	1,693 (53.2)	797 (56.2)	896 (50.8)	0.002
Anti-hypertensives	2,463 (77.4)	1,145 (80.7)	1,318 (74.8)	⩽0.001

Note: Bold represent statistically significant values.

Abbreviations: MPR, medication possession ratio; NSAIDs, non-steroidal anti-inflammatory drugs; COXIBs, COX-2 inhibitors; DMARDs, disease-modifying anti-rheumatic drugs; MTX, methotrexate.

Finally, to evaluate health care resource utilization and to limit potential biases, adherent patients and non-adherent patients were matched on a PS calculated using 17 confounding variables. All available variables were utilized in the calculation of the PS since the adherent and non-adherent patient cohorts differed with respect to year of index date, age, income quintile and health care resource use in the year prior to index date. Further, it would be particularly important to match the cohorts on the prior health care resource use as this would ensure that any divergence observed in the follow-up period in health care resource use was an impact of adherence to medication.

The matched cohort included a total 3,686 patients, 1,843 per group (Table 5). In contrast to the whole cohort, there was no significant difference in follow-up time between the two matched groups. The relative rate of visits to GPs, cardiologists and other physicians were, respectively, 13%, 22% and 18% lower in adherent patients. Interestingly, the mean number of visits to a rheumatologist was 6% higher in adherent patients. Also, the relative rate of hospitalization, ED visits, same-day surgeries and home care visits were 32%, 25%, 9% and 28% lower, respectively, in adherent patients. Moreover, the median length of stay was shorter by 8 days, in adherent patients. Overall, the relative rates for all health care resource use end points were significantly lower in adherent patients except for rheumatologist visits (Table 6).

Table 5

Resource use for the matched* cohort

Resource	Resource use			P-value
	Total	Adherent (MPR ≥0.8)	Non-adherent (MPR <<0.8)
	N=3<686	n=1<843	n=1,843
Follow-up time (days)
Mean (SD)	1<943 (1<112)	1<919 (1<085)	1,967 (1,138)	0.187
Median (IQR)	1<831 (1<003–2<736)	1<828 (974–2<673)	1,836 (1,018–2,786)	0.187
GP/FP
Any visit, n (%)	3<670 (99.6)	1<834 (99.5)	1,836 (99.6)	0.616
Rate per 100 follow-up days
Mean (SD)	4.34 (3.80315)	4.7 (4.06444)	3.97 (3.48451)	⩽0.001
Median (IQR)	3.25 (1.9753–5.3486)	3.58 (2.1277–5.8315)	2.93 (1.8328–4.798)	⩽0.001
Rheumatologist
Any visit, n (%)	3<161 (85.8)	1<562 (84.8)	1,599 (86.8)	0.081
Rate per 100 follow-up days
Mean (SD)	1.37 (1.44)	1.25 (1.38)	1.49 (1.49)	⩽0.001
Median (IQR)	0.96 (0.34–1.93)	0.84 (0.26–1.75)	1.09 (0.46–2.08)	⩽0.001
Cardiologist
Any visit, n (%)	2<588 (70.2)	1<352 (73.4)	1,236 (67.1)	⩽0.001
Rate per 100 follow-up days
Mean (SD)	0.34 (0.68)	0.39 (0.75)	0.29 (0.59)	⩽0.001
Median (IQR)	0.12 (0–0.37)	0.15 (0–0.44)	0.097 (0–0.31)	⩽0.001
Other
Any visit, n (%)	3<682 (99.9)	1<841 (99.9)	1,841 (99.9)	1
Rate per 100 follow-up days
Mean (SD)	4.59 (3.88)	5.00 (4.31)	4.2 (3.35)	⩽0.001
Median (IQR)	3.7 (2.27–5.6)	3.9 (2.38–6.17)	3.47 (2.13–5.26)	⩽0.001
ED visits
Any visit, n (%)	3<059 (83)	1<601 (86.9)	1,458 (79.1)	⩽0.001
Rate per 365 follow-up days
Mean (SD)	1.71 (1.75)	1.35 (2.00)	1.00 (1.50)	⩽0.001
Median (IQR)	0.69 (0.23–1.51)	0.84 (0.33–1.64)	0.55 (0.16–1.20)	⩽0.001
Same-day surgeries
Any visit, n (%)	2<386 (64.7)	1<214 (65.9)	1,172 (63.6)	0.148
Rate per 365 follow-up days
Mean (SD)	0.37 (0.58)	0.39 (0.58)	0.35 (0.57)	0.037
Median (IQR)	0.23 (0–0.52)	0.25 (0–0.54)	0.22 (0–0.50)	0.037
Hospitalizations
Any visit, n (%)	2<390 (64.8)	1<293 (70.2)	1,097 (59.5)	⩽0.001
Rate per 365 follow-up days
Mean (SD)	0.48 (0.78)	0.58 (0.80)	0.39 (0.75)	⩽0.001
Median (IQR)	0.24 (0–0.62)	0.33 (0–0.78)	0.18 (0–0.48)	⩽0.001
Length of stay (in days)
Mean (SD)	29.93 (40.5)	37.69 (48.62)	20.79 (25.18)	⩽0.001
Median (IQR)	15 (6–38)	19 (6–51)	11 (5–28)	⩽0.001
Home care
Any service, n (%)	1<220 (33.1)	723 (39.2)	497 (27)	⩽0.001
Rate per 365 follow-up days
Mean (SD)	23.86 (75.92)	28.10 (76.76)	19.60 (75.19)	⩽0.001
Median (IQR)	0 (0–1.2195)	0 (0–3.2164)	0 (0–0.1628)	⩽0.001
Rehabilitation
Any rehabilitation, n (%)	533 (14.5)	327 (17.7)	206 (11.2)	⩽0.001
Rate per 365 follow-up days
Mean (SD)	0.04 (0.135)	0.05 (0.16)	0.02 (0.01)	⩽0.001
Median (IQR)	0 (0–0)	0 (0–0)	0 (0–0)	⩽0.001

Notes:

Adherent patients were matched to non-adherent patients on a PS calculated using age, sex, income quintile, LHIN of residence, rurality, RUB, index year of first biologic script and health services use prior to index date (visits to GP/FP, rheumatologists, cardiologists and other MDs, ED visits, hospitalizations, same-day surgeries, home care services and rehabilitation episodes).

Abbreviations: MPR, medication possession ratio; SD, standard deviation; IQR, interquartile range; GP, general practitioner; FP, family practitioner; ED, emergency department; PS, propensity score; RUB, resource utilization band; MD, medical doctor; LHIN, Local Health Integration Network.

Table 6

Relative rates for resource use for the matched* cohort, using negative binomial distribution

Resource use	Adherence over entire follow-up period
	Relative rate	Lower CL	Upper CL
MD visits
GP/FP	0.88	0.83	0.93
Rheumatologist	1.06	1.00	1.12
Cardiologist	0.78	0.69	0.88
Others	0.82	0.78	0.86
ED visits	0.75	0.69	0.82
Same-day surgeries	0.91	0.84	0.98
Hospitalizations	0.68	0.63	0.74
Home care	0.72	0.55	0.93
Rehabilitation	0.60	0.49	0.74

Notes:

Adherent patients were matched to non-adherent patients on a PS calculated using age, sex, income quintile, LHIN of residence, rurality, RUB, index year of first biologic script and health services use prior to index date (visits to GP/FP, rheumatologists, cardiologists, and other MDs, ED visits, hospitalizations, same-day surgeries, home care services and rehabilitation episodes).

Abbreviations: CL, confidence limit; MD, medical doctor; GP, general practitioner; FP, family practitioner; ED, emergency department; PS, propensity score; RUB, resource utilization band; LHIN, Local Health Integration Network.

Discussion

The results of this study clearly show that, in a population cohort of RA patients, adherence to biologic DMARD therapy might be associated with lower health care resource utilization and clearly demonstrates the benefit of adherence to therapy when treating chronic inflammatory diseases such as RA.

There was a larger proportion of non-adherent patients in earlier index years (2003–2010) over later index years that could be explained either by better education on adherence or, more likely, by a longer duration of follow-up in the use of the MPR, which can be associated with a greater chance of being qualified as non-adherent. Additional differences at baseline were noted between adherent and non-adherent patients, most notably in age and income quintile that, although statistically significant, were very small numerically and are unlikely to be clinically relevant. Small numerical, yet statistically significant differences in health resource use were noted prior to the index date. Although the effect size is much smaller than what we observed after the index date in a matched cohort of patients, the similar pattern could reflect non-adherence to conventional therapy, although that was not examined in this study.

This study identified a number of factors associated with adherence including increased age.6,7,9 With respect to medications, another study has reported a strong positive association between a prescription for DMARDs 6 months prior to anti-TNF treatment and adherence.5 We did not find that association although it is likely caused by the reimbursement criteria for biologic DMARD therapy in Ontario that requires prior failure of one to three DMARDs, and that 90.2% of the cohort was already on DMARDs. However, we did identify an association between use of NSAIDs/COXIBs and oral prednisone and non-adherence. A possible explanation is that these patients have more active disease (which could be due to non-adherence) or that the higher number of concomitant medications increases the risk for non-adherence. Living in a rural environment and more frequent visits to a rheumatologist were also associated with better adherence.

Being adherent to biologic DMARD therapy was associated with a significant reduction in all health care resource utilization outcomes measured, except for visits to a rheumatologist, which was slightly greater in adherent patients. A previous study investigating adherence to anti-TNF therapy also found that adherent patients had 21.5% reduction in mean annual non-pharmacological expenditure for each patient,7 a rate almost identical to what we found in this study.

There was a significantly higher proportion of non-adherent RA patients who were using NSAIDs/COXIBs and corticosteroids such as oral prednisone and intra-articular steroid injections. This difference was noted in the pre-index date but increased by over two-fold once those patients were under biologic DMARD therapy, while no change was noted in the use of DMARDs and anti-hypertensives. It is possible that non-adherence to biologic DMARDs lead to a greater incidence of disease flare, which is then managed with NSAIDs/COXIBs and corticosteroids.1

Although utilizing health administrative databases has many benefits, a limitation is that data are collected for administrative or billing purposes and are therefore subject to limitations. Health administrative databases can be subject to inaccuracies such as in the identification of health conditions (eg, RA) reflecting a documented diagnosis by a physician or an institution. The data elements from CIHI DAD used in this study for the purpose of defining the diagnosis of RA, however, have been found to be very reliable (CIHI Data Quality Study of the 2009–2010 Discharge Abstract Database, Ottawa: CIHI, 2012). The study excluded other commonly occurring immunological conditions; however, certain conditions including axial spondyloarthropathy, systemic lupus erythematosus and Sjogren’s disease were not excluded and biologics could be prescribed for these indications even though not specifically indicated. The OHIP data set contains most physician billing claims paid by OHIP, including records by physicians that are not under a fee for service remuneration model. These records, referred to as shadow billing, are known to be incomplete since the physicians are salaried and the incentive to bill is reduced. This has an implication to the definition of RA, since every encounter with a physician may not be recorded in the OHIP data set. Since the ODB routinely captures claims for prescription drugs for people 65 years and older only, the study cohort is limited to that age group, which represents 44% of RA patients in Ontario as of 2010.10 Furthermore, ICES does not have data on prescriptions paid out of pocket or by private insurance companies, which could underestimate the use of biologic DMARDs and also makes distinguishing between biologic naïve and biologic experienced patients impossible. The RPDB provides basic demographic information about anyone who has ever received an Ontario health card number. Thus, the cohort of people diagnosed with RA is limited to those with health insurance coverage. Also, the ICES databases utilized did not directly capture marital status and educational level that can be considered a limitation of the study, since these aspects may reflect on adherence. However, patients’ income levels were considered, and this could reflect the level of education of the patients included in the study. There was a larger number of non-adherent than adherent patient group at the earlier index date and the former had a longer duration of follow-up. This could be explained either by the possibility of better adherence over time or that patients with longer follow-up were more likely to become non-adherent as time goes by. To limit this potential for bias, sensitivity analyses were carried out using three separate cohorts capturing patients with 3 years, 5 years and 7 years of follow-up. There was no difference in the data between any of those cohorts and the whole cohort used in this study (data not shown). Finally, retrospective studies using a database can only explore the relationship between factors in that database and adherence. Potentially relevant factors that are not in this database will remain undiscovered and it is not possible to find a causal relationship.4

Finally, like most studies evaluating adherence, the definition of adherence is based on MPR, which is calculated based on filled claims and does not guarantee that patients have actually taken the medication. In addition, this study has inherent limitations of administrative database studies that include the lack of clinical outcome data for RA.

A central strength of this study is that it was conducted in a large diversified population cohort that included all RA patients aged 65 years and older treated with biologic DMARDs within the province of Ontario. Another strength is that the study is linked to many administrative databases, including health care resource utilization and drug use, allowing for a comprehensive analysis with the power to identify small effects with confidence along with minimization of potential biases. This study also included all biologic DMARDs (not limited to anti-TNFs), and future research can examine potential differences.

Conclusion

We found that adherence to biological DMARD therapy is associated with significantly lower health care resource utilization and the need for concomitant therapies such as NSAIDs/COXIBs and corticosteroids. These findings have implications on patients, physicians and health care payers and highlight the importance of adherence to therapy in the management of a chronic disease such as RA.