Re-evaluating Safety and Effectiveness of Dabigatran Versus Warfarin in a Nationwide Data Environment: A Prevalent New-User Design Study.

INTRODUCTION: The new user cohort design is widely used to assess the effects of a new drug, such as dabigatran, but inherently excludes some users due to prior use of the comparator drug, for example warfarin. The prevalent new-user design offers a solution that includes all eligible users of the new drug.
OBJECTIVE: To evaluate the safety and effectiveness of dabigatran versus warfarin in non-valvular atrial fibrillation (NVAF) patients with prevalent new-user design.
METHODS: Taiwan National Health Insurance and mortality data from 2011 through 2015 were utilized. From an incident NVAF cohort, we identified dabigatran initiators as either incident or prevalent (switchers from warfarin) new users. Time- and prescription-based exposure sets were formed for dabigatran initiators to account for prior warfarin prescriptions. A comparable warfarin user was matched on the time-conditional propensity score to the dabigatran initiator in each set. The matched patients were followed for clinical outcomes, with Cox proportional hazards model used to estimate hazard ratios (HRs).
RESULTS: There were 10,811 dabigatran initiators, including 22% prevalent new users (switchers), who formed the exposure sets and were matched 1:1 to warfarin users. Dabigatran use was associated with lower risks of intracranial hemorrhage (HR 0.51; 95% confidence interval [CI] 0.39, 0.66) and gastrointestinal bleeding (HR 0.81; 95% CI 0.70, 0.92), compared with warfarin use. These effects were similar between the incident and prevalent new users.
CONCLUSION: Using a design that includes both incident and prevalent new users of dabigatran, the use of dabigatran is associated with lower major bleeding risk than warfarin use among patients with incident NVAF.

Key Points

Background

In early post-marketing safety assessment of newly marketed drugs, ideally all users of a new drug should be identified for evaluation of adverse outcomes. Some new drugs are novel therapies for disease conditions that do not have approved drug treatment and thus all users of the new drugs are genuine new users and would not have received any prior treatment. However, for disease conditions with existing drug treatment, newly approved drugs are either initiated in patients after switching from a prior drug (switchers) or initiated for patients without any prior medication for the disease condition (genuine new users). Safety assessment is further complicated by whether the new drug is approved as monotherapy, add-on therapy, or both. The new user design with active comparator has been the most common method in post-marketing safety and effectiveness research to study two drugs among comparable patients [1-4]. Studying exclusively new users of a new drug and new users of the comparator who have not received prior treatment would mimic a clinical trial with sufficient washout duration and to improve internal validity, but would exclude many users of the new drug who are switched from prior drugs in real-world practice. A pertinent example was a study evaluating the risk of angioedema associated with use of drugs targeting the renin-angiotensin-aldosterone system with beta-blocker as the active comparator; the study size was only 30% of all potential subjects when eligibility was restricted to new users who had no prior use of any study drugs [5]. Studying all eligible users of a new drug would provide more generalizable data and greater statistical power, but the heterogeneity of the new drug users would result in methodology challenges in identifying comparators and making valid inference on comparative safety and effectiveness.

The methodology challenges described above apply to post-marketing safety assessment of non-vitamin K oral anticoagulant (NOAC) for patients with non-valvular atrial fibrillation (NVAF) [6-12]. After a NOAC was introduced, some patients with prevalent NVAF might have been treated with warfarin suboptimally for years before switching to a NOAC; patients with newly diagnosed NVAF might be initially prescribed warfarin although the NOAC was available and were switched to the NOAC shortly afterwards, and other patients with newly diagnosed atrial fibrillation (AF) would initiate the NOAC without prior use of warfarin [11-15]. Despite the apparent heterogeneity of patients initiating NOAC, safety studies utilizing the active comparator new-user design of NOAC and warfarin have been implemented differently [6, 8, 16, 17]. In the early days after approval of NOACs when major bleeding risk was uncertain [18], regulatory agencies understandably would wish to include all NOAC initiators in safety studies rather than only studying the internally valid new users under the conventional research framework. Methodology concerns with including NOAC initiators who switched from warfarin use included under-ascertainment of outcomes of interest during prior warfarin exposure and change in confounder distribution as a result of prior warfarin use [4]. Another concern has been the “depletion of susceptibles” when the prevalent warfarin users who switched to NOAC had a lower bleeding risk [19, 20]. In prior studies of the newly introduced NOAC dabigatran, the proportion of dabigatran users with prior warfarin exposure ranged from 38 to 51% [11-15]. Including those who switched from warfarin in the dabigatran new users would result in the heterogeneity described above and excluding these dabigatran users would limit study size and generalizability of the study findings.

Prevalent new user design [21, 22] has been proposed recently to account for prior comparator drug use experience and to retain as many users of the newly introduced drug as possible in post-marketing observational studies. In this study design, users of a new drug (some had not used any comparator drug of interest before and others had used a comparator drug for the same indication) are compared against users of a comparator with similar prior exposure experience in the same analysis, accurately reflecting the mix of patients using the new drug in real-world settings. The objective of this study is to implement this novel design to evaluate safety and effectiveness of dabigatran among NVAF patients after its introduction to Taiwan, using warfarin as the comparator and taking into account the methodologic challenges when a substantial portion of dabigatran initiators were switched from warfarin.

Methods

Data Source

Taiwan National Health Insurance (NHI) is a publicly-funded single-payer compulsory program started in 1995 with coverage of more than 99% of the population in Taiwan [23]. A majority of healthcare providers contract with the Taiwan NHI Administration to provide a wide range of medical services. Clinical diagnoses are coded according to the International Classification of Disease Ninth Revision Clinical Modification (ICD-9-CM) through the end of 2015 and transitioned to the Tenth Revision starting in 2016. Diagnoses and drug prescriptions in NHI have been validated particularly in the field of cardiovascular diseases, demonstrating high positive predictive values and consistency for related co-morbidities as well as concomitant medications [24, 25]. For research purposes the NHI claims can be linked to Taiwan National Death Registry to provide information on date and cause of death [26]. For this study NHI and Death Registry data from 2011 to 2015 were accessed through the Health and Welfare Data Science Center of the Department of Statistics, Ministry of Health and Welfare in Taiwan.

Base Cohort

We identified patients with incident diagnosis of NVAF (ICD-9-CM 427.3 in ambulatory care or hospital admission claims at any claim position) from 2012 through 2015 with a look-back period of at least 1 year to ascertain prior diagnoses. Patients with mitral stenosis, valvular heart diseases, valve replacement, commissurotomy, heart transplantation, extracorporeal circulation or venous thromboembolism diagnoses 12 months before incident AF or age under 20 years at the time of first AF diagnosis were excluded. A base cohort was formed with patients having first warfarin or dabigatran prescribed on or after the initial AF diagnosis.

The first warfarin or dabigatran prescribed on or after AF was the cohort entry drug and its prescription date was the cohort entry date. Patients with more than one oral anticoagulant prescribed on cohort entry date were excluded. Patients who received rivaroxaban or apixaban (alternate NOACs) prior to cohort entry date were excluded (Fig. 1). Edoxaban was approved in February 2016 and not reimbursed by NHI until September 2016, and thus is not considered in this study (Table 1).

Fig. 1

Development of base cohort

Table 1

Date of approval for non-vitamin K oral anticoagulants in Europe and reimbursement in Taiwan

Non-vitamin K oral anticoagulants	Market approval			Reimbursement
	US FDA	EMA	TFDA	NHIA
Dabigatran	2010/10/19	2008/03/18	2011/07/13	2012/06/01
Rivaroxaban	2011/11/04	2008/09/30	2012/03/08	2013/02/01
Apixaban	2012/12/28	2011/05/18	2013/08/14	2014/06/01
Edoxaban	2015/08/01	2017/04/20	2016/02/02	2016/09/01

US FDA US Food and Drug Administration, EMA European Medicine Agency, TFDA Taiwan Food and Drug Administration, NHIA National Health Insurance Administration

Dabigatran initiators in the base cohort included both incident and prevalent new users. An incident new user, in some studies referred to as incident user or new user [2], was operationally defined in this study as a new user of dabigatran who did not use oral anticoagulants before (oral anticoagulant naïve), whereas a prevalent new user (switcher) was a dabigatran initiator who had prior use of another oral anticoagulant (oral anticoagulant non-naïve or oral anticoagulant experienced) [12, 13, 27]. In some reports the term initiators is used, which may represent new users of a certain drug regardless of prior use of drugs of the same class [28]. Different terms used in the literature to describe different drug-use scenarios are summarized in Table 2. In this report we use the terms “incident new users” and “new users” interchangeably, as well as “prevalent new users” and “switchers”.

Table 2

Terminology used to describe incident and prevalent use of a drug

Used in this report	Described in published literature [reference]	Definition
Initiator	New user [4, 41, 42] Naïve user [4] Incident user [2, 42] Initiator [28]	New users of the drug regardless of prior use of other drugs for the same disease condition
Incident new user or new user	New user [4, 21] Naïve user [4] Incident user [21] Incident new user [21]	New user of a drug without prior use of any drug for the same disease condition
Prevalent new user or switcher	New user [4] Switcher [21] Prevalent new user [21]	New user of a drug with prior use of other drug(s) for the same disease condition
Prevalent user (not considered in this report)	Prevalent user [4, 41, 42] Non-naïve user [4] Continuing user [2] Current user [2]	Already using the drug at the start of study period

Co-morbidities of interest and their corresponding ICD-9-CM codes and NHI specific procedure codes are described in Supplementary Table 1, and concomitant drugs of interest with corresponding Anatomical Therapeutic Chemical Classification System (ATC) codes are described in Supplementary Table 2.

Time- and Prescription-based Exposure Sets

Warfarin and dabigatran prescriptions were individual units to be evaluated to construct exposure sets from the base cohort, while in traditional cohort studies, the unit of interest would be individual patients. An exposure set was formed for every new dabigatran prescription; number of prior warfarin and time since cohort entry to new dabigatran prescription were enumerated among patients initiating dabigatran. Patients who had alternate NOACs prescribed before cohort entry were excluded from the exposure sets. When constructing the exposure sets, we also excluded patients with two oral anticoagulants prescribed on the same day after cohort entry. For each exposure set, there was one dabigatran initial prescription (dabigatran index prescription) and multiple corresponding warfarin prescriptions. Eligible warfarin prescriptions in the exposure set were those who fulfilled the following criteria: (1) time since cohort entry to the warfarin prescription was within 15 days of the interval between cohort entry and dabigatran index prescription, (2) there were same number of prior warfarin prescriptions before the dabigatran index prescription, and (3) they occurred during the same calendar year as the dabigatran index prescription (Fig. 2). Initiating dabigatran prescriptions were then categorized as either an incident new-user exposure set—no warfarin prescription before the initial dabigatran prescription, or a prevalent new user (switcher) exposure set—having at least one warfarin prescription prior to the initial dabigatran prescription. A patient could contribute warfarin prescriptions to multiple exposure sets.

Fig. 2

Time- and prescription-based exposure set (matched on calendar year). For a comparator drug prescription to be included in an exposure set of a study drug, the comparator drug must have: (1) within 15-day difference of time from cohort entry, (2) the same number of prior warfarin prescriptions, and (3) the same calendar year as the new dabigatran prescription is included in the exposure set of the new dabigatran prescription

Time-conditional Propensity Score (PS)

Propensity score (PS) estimation with anticoagulant prescription (dabigatran or warfarin) as the dependent binary variable was carried out with conditional logistic regression within each exposure set. Covariates in the regression model included age, gender, co-morbidities, and concomitant medications among the initial dabigatran and warfarin prescriptions at the time of forming or being included in the exposure sets. We ascertained co-morbidities and concomitant medications during the 90 days before the dabigatran or warfarin prescription of interest, as most patients with chronic conditions like NVAF in Taiwan visit physicians every 3 months for chronic medication refills, so diagnoses of both chronic and acute conditions would be captured in the health insurance claims. We estimated PS for every dabigatran and warfarin prescription included in the exposure sets by using the coefficients obtained from the conditional logistic regression model. For new users, PS would be the probability of a patient initiating a dabigatran prescription and among prevalent new users (switchers), PS would be the probability of a patient switching from warfarin prescription to dabigatran.

PS Matching and Study Cohort

One warfarin prescription was matched to the initial dabigatran prescription within the same exposure set according to the estimated PS. A caliper of 0.06 standard deviation of logit PS within each exposure set was used to define the range of estimated PS within which to select the warfarin prescription. Warfarin prescription within the caliper of the corresponding dabigatran prescription and having the nearest PS to the dabigatran prescription was selected. The matching process was carried out in chronological order of when the exposure sets were formed. Once the patient with a warfarin prescription was matched to an initial dabigatran prescription, other warfarin prescriptions of the same patient in other exposure sets would no longer be eligible for matching to the yet to be matched new dabigatran patients. A patient may be matched as a warfarin patient to a dabigatran patient and later initiated with dabigatran, which was described as a prevalent new user (switcher) who was being “reused” [21]. The selected warfarin prescription and the corresponding new dabigatran prescription were index prescriptions and the prescription date became the index date for individual patients. While prescription was the unit utilized in formation of exposure sets and PS modeling within the base cohort, the corresponding individuals became the unit of analysis after 1-to-1 prescription matching. The study cohort comprised matched dabigatran and warfarin users, including both dabigatran incident and prevalent new users.

Clinical Outcomes

Outcomes of interest included intracranial hemorrhage, gastrointestinal bleeding, and ischemic stroke that resulted in hospital admission and all-cause mortality. Diagnosis codes that have been evaluated with Taiwan NHI data were used [12]. Diagnosis code for ischemic stroke was validated with data from a single medical center in southern Taiwan, where 97.9% of NHI ischemic stroke cases were confirmed through review of medical records [24].

Statistical Analysis

One warfarin prescription representing a warfarin user of an exposure set was randomly sampled from each exposure set before PS matching to evaluate comparability between the dabigatran users and warfarin users. Clinical diagnoses and drug use 1 year prior to the dabigatran index prescriptions and the randomly selected warfarin index prescriptions were evaluated in order to assess the effect of balancing potential confounders between dabigatran and warfarin users with PS matching. Standardized mean differences for prevalence of covariates were calculated with a threshold of less than 0.1, implying well balanced characteristics [29].

Patients in the study cohort were followed from index date to the earliest of first event of interest, switching to another oral anti-coagulant, death, or 31 December 2015. Cox proportional hazard analysis was undertaken using the maximum partial likelihood method for estimation with a robust sandwich covariance matrix to account for the same patient contributing dabigatran-exposed person time and warfarin-exposed person time [30]. Hazard ratios (HRs) and the corresponding 95% confidence intervals (CIs) were reported. Stratified analysis for incident and prevalent new users was performed to explore whether bleeding, ischemic stroke and mortality risks would be modified by prior warfarin use. All analyses were performed using SAS software, version 9.4 (SAS Institute, Inc, Cary, NC, USA).

Results

We identified 143,326 incident NVAF patients in the Taiwan NHI from 2012 through 2015, with 30,458 patients having at least one warfarin or dabigatran prescription on or after the date of incident AF diagnosis. After excluding patients with incomplete NHI enrollment information, patients who received rivaroxaban or apixaban before receiving initiating dabigatran or warfarin, and those who were prescribed dabigatran and warfarin on the cohort entry date, the base cohort comprised 28,562 patients.

Exposure Sets

Among the base cohort AF patients, 528 were ineligible for evaluation in the prevalent new-user design due to having more than one oral anticoagulant prescribed on the same day after cohort entry. A total of 10,830 base cohort patients received dabigatran; 77% were incident new users and 23% were prevalent new users. As a result, 10,830 exposure sets could potentially be formed by an initial dabigatran prescription. After considering time and number of warfarin prescriptions from cohort entry up to the initial dabigatran prescriptions, 10,811 (99.8%) initiating dabigatran prescriptions were able to find at least one eligible warfarin prescription to form time- and prescription-based exposure sets.

Study Cohort after PS Matching

After estimating PS for all dabigatran or warfarin prescriptions in each exposure set, a warfarin prescription having estimated PS being nearest within the caliper of that for the dabigatran prescription in the same exposure set was identified for 10,781 of dabigatran prevalent new users; 8385 (78%) were dabigatran incident new users and 2396 (22%) were prevalent new users.

Using one randomly selected warfarin user from each exposure set to evaluate the warfarin users’ and dabigatran users’ baseline characteristics, we found that the process of formulating the exposure sets according to prior warfarin use and time since first diagnosis of NVAF resulted in moderate comparability between the two groups (Supplementary Table 3). Dabigatran incident new users were older (mean age 75 years) than dabigatran prevalent new users (mean age 73 years) and had a lower prevalence of congestive heart failure, prior history of stroke, valvular heart disease, hemiplegia, and less concomitant use of statins (Supplementary Table 4). For the study cohorts after PS matching, mean age was 74 years and about 60% were male. Approximately 70% of patients had hypertension, 30% had diabetes, and 10% had renal impairment. Prevalence of all measured baseline covariates was balanced between the two study cohorts (Table 3).

Table 3

Characteristics of study subjects, propensity score-matched dabigatran and warfarin users in time- and prescription-based exposure sets (1 year before index prescription)

	After propensity score match
	Dabigatran		Warfarin		Standardized difference
n, (%)	10,781		10,781
Mean age, year (SD)	74.52	10.16	73.79	10.98	0.07
Male	6329	59%	6276	58%	0.01
CHA2DS2-VASc
Mean (SD)	3.90	1.64	3.77	1.85	0.07
Median (Q1–Q3)	4.00	(3–5)	4.00	(2–5)
Congestive heart failure	3209	30%	3301	31%	− 0.02
Hypertension	7768	72%	7674	71%	0.02
Chronic pulmonary disease	2308	21%	2248	21%	0.01
Diabetes	3237	30%	3255	30%	0.00
Renal failure	1170	11%	1287	12%	− 0.03
Liver disease	920	9%	924	9%	0.00
Peptic ulcer disease excluding bleeding	1260	12%	1258	12%	0.00
Ischemic stroke	3229	30%	3010	28%	0.04
Intracranial hemorrhage	210	2%	209	2%	0.00
Gastrointestinal bleeding	584	5%	612	6%	− 0.01
Myocardial infarction	350	3%	383	4%	− 0.02
Vascular disease	706	7%	735	7%	− 0.01
Hemiplegia or paraplegia	579	5%	605	6%	− 0.01
Dementia	605	6%	605	6%	0.00
Gout	1013	9%	1066	10%	− 0.02
Percutaneous coronary intervention	215	2%	228	2%	− 0.01
Coronary artery bypass surgery	7	0%	25	0%	− 0.04
Antiplatelet	3271	30%	3211	30%	0.01
Low-dose aspirin	7144	66%	6923	64%	0.04
NSAID	7101	66%	7188	67%	− 0.02
Antiepileptics	292	3%	290	3%	0.00
Antifungal	130	1%	126	1%	0.00
Diuretics	2677	25%	2658	25%	0.00
Beta-blockers	6938	64%	6979	65%	− 0.01
CCB	7341	68%	7306	68%	0.01
ACEI	2369	22%	2404	22%	− 0.01
ARB	6284	58%	6164	57%	0.02
Statin	3812	35%	3694	34%	0.02
Other lipid-lowering agents	680	6%	693	6%	0.00
Insulin	1200	11%	1271	12%	− 0.02
Other oral antidiabetics	2967	28%	2970	28%	0.00
Dabigatran new users	8385	78%
Dabigatran switchers	2396	22%

SD standard deviation, Q1–Q3 quartile 1-quartile 3, NSAID nonsteroidal anti-inflammatory drug, CCB calcium channel blocker, ACEI angiotensin-converting enzyme inhibitor, ARB angiotensin receptor blocker

Median duration of follow-up among the study cohort was 0.9 year. Incidence rates among dabigatran prevalent new users for intracranial hemorrhage, gastrointestinal bleeding, and ischemic stroke were 0.70, 3.17, and 4.91 per 100 person-year, respectively, and were 1.38, 3.96, and 5.44 per 100 person-year, respectively, among warfarin users. There were 1472 prevalent new users whose prior warfarin-exposed person time before initial dabigatran use was matched to another dabigatran user. The intracranial hemorrhage, gastrointestinal bleeding, and ischemic stroke incidence during the prior warfarin-exposed person time were 1.50, 3.67, and 11.48 per 100 person-year.

Decreased risk of intracranial hemorrhage (HR = 0.51; 95% CI 0.39, 0.66) and gastrointestinal bleeding (HR = 0.81; 95% CI 0.70, 0.92) were observed among dabigatran initiators compared with the warfarin users (Table 4). There was no strong evidence to suggest a lower risk of ischemic stroke among dabigatran and warfarin users (HR = 0.91; 95% CI: 0.81, 1.02). Dabigatran initiators had a reduced risk for all-cause mortality (HR = 0.59; 95% CI 0.54, 0.64) (Table 5).

Table 4

Bleeding outcome analyses

Outcomes	Total	Number of events	Total person-year	Incidence rate (per 100-person year)	Matched hazard ratio (95% CI)
Intracranial hemorrhage
Dabigatran	10,781	85	12,159.64	0.70	0.51 (0.39, 0.66)
Warfarin	10,781	161	11,664.71	1.38	1
Incident new users
Dabigatran	8385	65	9304.59	0.70	0.50 (0.37, 0.68)
Warfarin	8385	122	8726.07	1.40	1
Prevalent new users
Dabigatran	2396	20	2855.05	0.70	0.53 (0.31, 0.91)
Warfarin	2396	39	2938.63	1.33	1
Gastrointestinal bleeding
Dabigatran	10,781	379	11,959.77	3.17	0.81 (0.70, 0.92)
Warfarin	10,781	452	11,416.87	3.96	1
Incident new users
Dabigatran	8385	304	9136.81	3.33	0.84 (0.72, 0.98)
Warfarin	8385	341	8551.01	3.99	1
Prevalent new users
Dabigatran	2396	75	2822.96	2.66	0.69 (0.51, 0.92)
Warfarin	2396	111	2865.86	3.87	1

CI confidence interval

Table 5

Effectiveness outcome analyses

Outcomes	Total	Number of events	Total person-year	Incidence rate (per 100-person year)	Matched hazard ratio (95% CI)
Ischemic stroke
Dabigatran	10,781	580	11,820.41	4.91	0.91 (0.81, 1.02)
Warfarin	10,781	614	11,297.11	5.44	1
Incident new users
Dabigatran	8385	438	9055.96	4.84	0.84 (0.74, 0.96)
Warfarin	8385	490	8437.16	5.81	1
Prevalent new users
Dabigatran	2396	142	2764.45	5.14	1.19 (0.94, 1.50)
Warfarin	2396	124	2859.95	4.34	1
All-cause mortality
Dabigatran	10,781	872	12,214.31	7.14	0.59 (0.54, 0.64)
Warfarin	10,781	1437	11,748.9	12.23	1
Incident new users
Dabigatran	8385	707	9344.38	7.57	0.59 (0.53, 0.65)
Warfarin	8385	1138	8793.46	12.94	1
Prevalent new users
Dabigatran	2396	165	2869.93	5.75	0.57 (0.47, 0.69)
Warfarin	2396	299	2955.44	10.12	1

CI confidence interval

In the subgroup analyses stratified by incident new users and prevalent new users, similar levels of lower risk of major bleeding with dabigatran use were observed compared with warfarin. For gastrointestinal bleeding, HR was 0.84 for the dabigatran incident new users and it was 0.69 among the prevalent new users (Table 4). On the other hand, ischemic stroke risk was lower among incident new users of dabigatran (HR = 0.84; 95% CI 0.74, 0.96), but not observed among prevalent new users (HR = 1.19; 95% CI 0.94, 1.50). The reduction in all-cause mortality rate with dabigatran was similar in incident (HR = 0.59) and prevalent new users (HR = 0.57) when compared to warfarin users (Table 5).

Discussion

In this population-based study conducted within an incident NVAF base cohort, we identified 8385 dabigatran incident new users (new users) and 2396 prevalent new users (switchers) in a nationwide data environment and implemented the prevalent new-user design to evaluate comparative safety and effectiveness of dabigatran and warfarin. The novel design allows for studying incident new users and prevalent new users in the same analysis through time- and prescription-based exposure sets. Comparable warfarin users were identified for both dabigatran incident and prevalent new users, and the dabigatran users had a significant reduction in intracranial hemorrhage and gastrointestinal bleeding risks. Dabigatran was more effective in the prevention of ischemic stroke among the incident new users but the same effect was not observed among the prevalent new users.

Hernan and Robins proposed the “target trial” concept, and within that framework [31], the prevalent new-user design would emulate a trial that allows for patient enrollment regardless of prior use of the comparator drug. Within a pool of eligible subjects with the same usage experience of the comparator drug with regard to the number of prescriptions within a defined baseline timeframe, one subject would be randomly assigned to the new drug and another would be randomly assigned to the comparator.

If dabigatran is evaluated in a conventional new-user design without considering prior warfarin use [1, 4], the nearly 10,000 dabigatran initiators in Taiwan would form a heterogeneous group, introducing potential confounding. One the other hand, if the study population only comprised dabigatran initiators who did not use warfarin previously, study size would be reduced by 22% and there would be no information on the prevalent new users (switchers). With the use of prevalent new user design in the post-marketing assessment of dabigatran in a healthcare system, virtually all incident NVAF patients prescribed dabigatran regardless of prior warfarin use were evaluated in the same analysis, providing comprehensive safety experience of dabigatran use as well as the capacity to evaluate drug effects separately among incident and prevalent new users. The design is most relevant in a therapeutic area where the safety of a novel medication is assessed in comparison with an existing regimen.

Our findings on intracranial hemorrhage, gastrointestinal bleeding, ischemic stroke, and all-cause mortality in the dabigatran versus warfarin comparison were within the range of those reported from prior studies from Denmark and the USA [8, 10, 11, 27, 32]. The longitudinal nature of the Taiwan NHI allows for extended look-back period to evaluate warfarin use before dabigatran initiation, which may not be feasible in health insurance claims data in the USA. In the US studies, look-back periods ranged from 3 months to 12 months [8, 16, 33], and may result in mis-classification of prior warfarin use and limited ascertainment of prior warfarin use duration. The longitudinal nature of the Taiwan NHI is more similar to the Danish data environment, where ascertainment of prior warfarin use may be comprehensively addressed either by the design of exposure sets or other analytical methods. In a Danish study that evaluated bleeding risks among NOACs and warfarin users, prior warfarin use experience was evaluated as a covariate in multiple regression models [13]. In applying the prevalent new user design, the use of exposure sets for index prescriptions in the same calendar year could simultaneously control for the duration and frequency of prior warfarin use in a prospective manner, which offered a transparent alternative to conventional cohort design with time-dependent analyses in a Cox model. A similar approach could be applied for the evaluation of rivaroxaban, another NOAC [17]. However, the real-world drug use scenario is more complex for subsequently approved NOACs, such as apixaban and edoxaban, as there would be switching from warfarin and between NOACs. The prevalent new user design is also applicable in head-to head comparisons amongst different NOACS [33-36]. However, the design of exposure sets and PS estimation would be more complicated, as the design needs to take into account for various patterns of prior warfarin and NOACs use. Furthermore, the implementation of prevalent new user design would be more complicated in disease conditions for which a drug of interest may be used in monotherapy and add-on therapy. For NVAF only monotherapy with an oral anticoagulant is clinically indicated, but that is not the case for other diseases. For instance, long-acting inhaled corticosteroids may be added to a beta-adrenergic agonist for combination treatment for chronic obstructive pulmonary disease [37] or a sulfonylurea as a second-line oral antidiabetic agent would be added to first-line metformin in diabetic patients [22]. When a new drug of interest is added on to the existing regimen, defining additional exposure group is suggested to complement the incident and prevalent new users for cohort studies [22].

Findings from this study were mostly consistent with those from prior reports from Taiwan using the same nationwide data over different time periods [10, 11, 38]. Dabigatran use was associated with approximately 50% reduced risk of intracranial hemorrhage and all-cause mortality [10, 11]. The 20% risk reduction in gastrointestinal bleeding observed in our study is similar to that of another report (HR = 0.77; 95% CI 0.59, 1.02) [38], although two earlier studies reported a similar gastrointestinal bleeding risk among dabigatran and warfarin users [10, 11]. As for ischemic stroke, we did not find a favorable ischemic stroke prevention effect for dabigatran, and previous studies from Taiwan reported 20–40% ischemic stroke risk reduction associated with dabigatran [10, 11, 38]. A methodology concern for prior studies was that for recurrent events such as bleeding or ischemic stroke, event rates during warfarin-exposed person-time before dabigatran initiation among prevalent new users were not ascertained in conventional retrospective cohort studies [11]. Studying prevalent new users without addressing comparable exposed person-time and not ascertaining outcome events during such person-time may introduce bias in prevalent user designs [4]. In this study, rates of outcomes of interest during the pre-dabigatran warfarin-exposed person-time were ascertained and major bleeding rate was similar to that during overall warfarin-exposed person time and ischemic stroke incidence rate was higher during the prior warfarin use person-time. These findings suggested under-ascertainment of early events in the prior warfarin use among dabigatran prevalent new users in previous study [4].

With the use of time- and prescription-based exposure sets in prevalent new user design, we found that distributions of potential confounders, except for renal failure and age, between dabigatran and warfarin users were empirically balanced before PS matching, illustrating an efficient means to control for potential confounders without using complex statistical analysis. After matching on time-conditional PS, distribution of all measured confounders was balanced between dabigatran and warfarin users. In this study with prevalent new-user design, given that comparable warfarin users were selected within the exposure sets and the subgroup of interest was prior warfarin use status, subgroup analysis of incident and prevalent new user could be performed without breaking the matched sets, while retaining balanced characteristics from PS matching. A methodology review suggested that most propensity matched studies (33% of the 83 studies reviewed) appeared to perform subgroup analysis according to any variable of interest directly from the matched cohorts, which would break the PS-matched sets and patient characteristics may no longer be balanced within the subgroups [39]. We suspect that similar issues may also be present in subgroup analysis along with inverse probability of treatment weighting where the cohort was subdivided directly from the weighted cohort. Although a recent simulation study suggested that limited bias was observed for breaking matched sets in subgroup analysis, the authors concluded that it was due to the low prevalence of broken matches. If more covariates were added to the PS model or high dimensional PS models were used, prevalence of broken matches may not be low in subgroup analysis [40]. In the context of prevalent new-user design studies, incident or prevalent new-user status were of primary interest and thus evaluation of these subgroups would be very straightforward. However, subgroup analyses on variables other than prior active comparator use would face similar concerns on imbalanced confounder distribution after stratification.

While the public health insurance program in Taiwan generated extensive longitudinal data with minimal loss to follow-up, our limited study budget only allowed for access to NHI data starting from 2011. As such, dabigatran users for whom AF onset was before 2011 were not included in the study. It did not affect the validity of the inference, but the number of dabigatran prevalent new users could be increased if an earlier data access time was available. Patients who had AF onset before 2011, received warfarin, and subsequently switched to dabigatran would be identified, resulting in a larger number of dabigatran prevalent new users in the study. As in any study based on health insurance claims, lack of lifestyle factors and laboratory examination results, especially international normalized ratio (INR) for monitoring warfarin use, are the limitations. Although INR tests are required for warfarin dose adjustment, when both frequency and duration of prior warfarin treatment are accounted in the exposure sets, the warfarin utilization information captured in the exposure sets may indirectly control for anticoagulant treatment effects and mitigate the limitations of no INR data in the study.

Conclusion

In conclusion, we utilized a novel design that took into account prior drug-use history in post-marketing safety evaluation of a new drug (dabigatran) after it was introduced to Taiwan. The method could be applied to the evaluation of all new drugs in similar post-marketing settings.

Below is the link to the electronic supplementary material.

Supplementary material 1 (PDF 335 kb)

The prevalent new user design allows for inclusion of both dabigatran incident new users and prevalent new users in the same analysis where prior use of warfarin was accounted for through the time- and prescription-based exposure sets.

Dabigatran use was associated with lower risk of intracranial hemorrhage, gastrointestinal bleeding and all-cause mortality compared to warfarin use.

Outcomes among dabigatran incident new users and prevalent new users could be evaluated directly from subgroup analysis, where lower ischemic stroke risk compared to warfarin use was observed among dabigatran incident new users but not prevalent new users.