- Research article
Patterns of health care utilization related to initiation of amitriptyline, duloxetine, gabapentin, or pregabalin in fibromyalgia
Arthritis Research & Therapyvolume 17, Article number: 18 (2015)
Several pharmacologic treatments are available for fibromyalgia, but little is known about the comparative effectiveness of these treatments on health care utilization.
Using US commercial insurance claims data (covering 2007 to 2009), we conducted a cohort study to examine the comparative effectiveness of amitriptyline, duloxetine, gabapentin, and pregabalin on health care utilization in patients with fibromyalgia. We measured patients’ medication adherence using the proportion of days covered (PDC) and estimated multivariable rate ratios (RRs) for outpatient visits, prescriptions, hospitalization, and emergency department (ED) visits in propensity score (PS)–matched cohorts.
Cohorts of 8,269 amitriptyline, 9,941 duloxetine, and 18,613 gabapentin initiators were compared with their PS-matched pregabalin initiators. During the baseline 180-day period, patients had, on average, seven to nine physician visits, including six to eight specialist visits, and received eight prescription drugs. The mean PDC up to 180 days varied from 38.6% to 67.7%. The number of outpatient visits, prescriptions, and hospitalizations decreased slightly after initiating one of the study drugs, but the number of ED visits increased after treatment initiation. Compared to pregabalin, duloxetine was associated with decreased outpatient visits (RR, 0.94; 95% confidence interval (CI), 0.88 to 1.00), prescriptions (RR, 0.94; 95% CI, 0.90 to 0.98), hospitalizations (RR, 0.75; 95% CI, 0.68 to 0.83), and ED visits (RR, 0.85; 95% CI, 0.79 to 0.91). Little difference in health care utilization rates was noted among amitriptyline and gabapentin initiators compared to those who were started on pregabalin.
Fibromyalgia patients had high health care utilization before and after initiation of amitriptyline, duloxetine, gabapentin, or pregabalin. Medication adherence was suboptimal. Overall, fibromyalgia treatment had little impact on reducing health care utilization, but duloxetine initiators had less health care utilization than those started on pregabalin.
Fibromyalgia is a common, highly disabling syndrome characterized by chronic widespread pain, fatigue, and problems with sleep, memory, and mood [1-3]. Treatment of fibromyalgia is challenging and often requires both pharmacologic and nonpharmacologic treatments [4,5]. For pharmacologic treatment, tricyclic antidepressants (TCAs) such as amitriptyline, selective serotonin reuptake inhibitors (SSRIs), serotonin and noradrenaline reuptake inhibitors (SNRIs) such as duloxetine and milnacipran, and monoamine oxidase inhibitors are used to reduce pain, depression, fatigue, and health-related quality of life [6,7]. Gabapentin and pregabalin are also commonly prescribed to ameliorate pain and sleep problems associated with fibromyalgia [6,7].
Patients with fibromyalgia are known to be heavy users of health care systems, including outpatient visits, hospitalizations, and prescription drug use [8-10]. Despite the number of available drugs for fibromyalgia, few studies have directly compared the effects of different pharmacologic agents for fibromyalgia on health care utilization. In studies comparing pre- and posttreatment, pregabalin initiators had increased inpatient and outpatient visits after initiating pregabalin, whereas duloxetine initiators had decreased inpatient and outpatient visits after initiating duloxetine [11-13]. In another study, researchers reported increases in total health care costs after initiating treatment with pregabalin or duloxetine for fibromyalgia . In a study comparing pregabalin to TCAs, pregabalin initiators had increased total health care costs from pretreatment to follow-up . Although all of these studies were performed using US claims data, the effects of fibromyalgia treatment on health care utilization may vary across different health care systems, particularly across different countries.
The objective of this study was to examine the comparative effectiveness of four commonly used drugs for fibromyalgia—amitriptyline, duloxetine, gabapentin, and pregabalin—on the rate of health care utilization in a US population-based cohort of commercially insured patients with fibromyalgia. Patients’ adherence to these drug regimens (that is, continuation or discontinuation) was also assessed. We hypothesized that the effect of amitriptyline, duloxetine, or gabapentin on health care utilization would differ from the effect of pregabalin.
We conducted a cohort study using claims data obtained from UnitedHealthcare, a US commercial health insurance plan, which insures primarily working adults and their family members, for the period from 1 January 2007 to 31 December 2009. The details of the data source are presented elsewhere . The study protocol was approved by the Institutional Review Board of the Brigham and Women’s Hospital. Patients’ informed consent was not required by the Institutional Review Board, as the dataset was deidentified to protect subject confidentiality.
Among patients ages 18 years and older who had at least one visit coded with the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD 9-CM), code 729.1x, for fibromyalgia, new users of a study drug—amitriptyline, duloxetine, gabapentin, or pregabalin—were identified. Milnacipran users were initially considered, but they were subsequently excluded because of the small number of patients taking milnacipran (less than 0.5% of all subjects). To be included in the study cohort, we required patients to have at least 180 days of continuous health plan eligibility before receiving the first prescription of a study drug. Patients were required to be naïve to all four drugs during the 180-day baseline period. We excluded patients who were prescribed more than one study drug at the same index date.
Follow-up time started on the index date, defined as the dispensing date of the first study drug, and ended at the first of any of the following censoring events: discontinuation or switching of study drugs, loss of health plan eligibility, end of study database, or death.
Health care utilization
Changes in health care utilization factors (for example, the number of outpatient visits to any physicians, primary care physicians and any specialists; number of prescription drugs; acute care hospitalizations; emergency department (ED) visits; physical therapy) were assessed during follow-up.
Variables potentially related to fibromyalgia symptoms, initiation of the study drugs, or health care utilization were assessed using data from the 180-day baseline period. These variables included age, sex, fibromyalgia-related comorbidities (back pain, chronic headache, depression, anxiety, chronic pain, fatigue, sleep disorder, and abdominal pain condition), other comorbidities (hypertension, diabetes, stroke, cardiovascular disease, malignancy, and inflammatory arthritis), use of other medications (opioids, benzodiazepines, anticonvulsants excluding gabapentin and pregabalin, TCAs excluding amitriptyline, SSRIs, SNRIs excluding duloxetine, other antidepressants, muscle relaxants, nonbenzodiazepine sleep disorder drugs, oral glucocorticoids, nonsteroidal anti-inflammatory drugs (NSAIDs), topical analgesics, migraine drugs, and gastrointestinal protective drugs), and health care utilization factors (number of outpatient visits to any physicians, primary care physicians, and specialists; number of prescription drugs; physical therapy; ED visits; and acute hospitalizations) [15,16]. Additional files 1 and 2 contain the lists of diagnoses and procedure codes and medications. To quantify patients’ comorbidities, we also calculated the Deyo-adapted Charlson Comorbidity Index based on the ICD-9-CM [17,18].
We compared the baseline characteristics among the initiators of each study drug. To control for potential confounding by indication, we used propensity score (PS) estimation . In multivariable logistic regression models using pregabalin initiators as the common reference, we estimated a PS. Pregabalin was chosen as the common reference because (1) we wanted to directly compare one treatment approved by the US Food and Drug Administration (FDA) to another FDA-approved treatment (duloxetine) and other non-FDA-approved treatments (amitriptyline and gabapentin), and (2) the number of pregabalin initiators was greater than that of duloxetine initiators. All the aforementioned baseline covariates except health care utilization factors were included in the logistic models. Each pair—amitriptyline versus pregabalin, duloxetine versus pregabalin, and gabapentin versus pregabalin—was then matched on the PS using the greedy matching algorithm . To measure patients’ adherence to their fibromyalgia treatments, we calculated the proportion of days covered (PDC) as the total number of days’ supply provided by prescriptions for a given study drug up to 180 days from the index date divided by 180 days.
Adjusting for baseline health care utilization factors, we used multivariable Cox models to estimate rate ratios (RR) for ED visits, acute care hospitalizations, and physical therapy during the follow-up period in all three PS-matched pairs . Multivariable Poisson regression, adjusted for baseline health care utilization factors, was used to estimate RR for outpatient visits to any physician, primary care physicians, and any specialists, as well as the number of prescription drugs, during the follow-up period in all three PS-matched pairs. To deal with the potential overdispersion in Poisson models, the standard errors were adjusted with the scale option in SAS 9.3 software (SAS Institute, Cary, NC, USA). To assess a long-term effect of treatment, a subgroup analysis using multivariable Poisson and Cox regression models was performed in patients with treatment duration of at least 180 days. All analyses were done using SAS 9.3 software.
Among 116,183 patients who had at least one diagnosis of fibromyalgia between 2007 and 2009, we identified 13,404 amitriptyline, 18,420 duloxetine, 23,268 gabapentin, and 19,286 pregabalin initiators . Of those, 8,269 amitriptyline, 9,941 duloxetine, and 18,613 gabapentin initiators were matched to pregabalin initiators on their PS with a ratio of 1:1. Thus, different subsets of pregabalin users were compared for each pairwise analysis.
Demographic characteristics, comorbidities, and medications at baseline were well balanced and statistically not different between the PS-matched cohorts (Table 1). The mean age ranged from 47.7 to 50.0 years, depending on the treatment pairs. The majority of patients were women (77% to 82%). The median starting daily doses were 25 mg for amitriptyline, 60 mg for duloxetine, 300 mg for gabapentin, and 75 mg for pregabalin. Back pain was the most common comorbidity (52% to 63%), followed by hypertension (30% to 36%). During the 180-day baseline, more than one-half of patients received at least one prescription for opioids and one-third used benzodiazepines. Use of SSRIs, non-SSRIs/SNRIs antidepressants, muscle relaxants, oral steroids, and NSAIDs was also common.
On average, patients had 7.4 to 8.8 outpatient visits and 5.9 to 7.7 specialist visits during the 180-day baseline period. The mean number of prescription drugs at baseline was also high, between 7.9 and 8.4. Approximately 10% to 12% of patients had at least one ED visit and 10% to 14% had at least one acute hospitalizations during the 180-day baseline period.
Continuation of fibromyalgia medication
The mean (standard deviation (SD)) duration of treatment was 77 (92) days for the amitriptyline-pregabalin pairs, 94 (110) days for the duloxetine-pregabalin pairs, and 82 (97) days for the gabapentin-pregabalin pairs. The mean (SD) PDCs up to 180 days of follow-up for the amitriptyline and pregabalin pair were 38.6% (30.1) and 50.0% (31.2), respectively. Among the amitriptyline and pregabalin initiators, 7.9% remained on the treatment for at least 180 days. For the duloxetine and pregabalin pair, the mean (SD) PDC values were 55.0% (34.0) and 67.8% (33.4), respectively. Among duloxetine and pregabalin initiators, 16.5% and 8.3%, respectively, remained on the treatment for at least 180 days. The mean PDCs (SD) up to 180 days of follow-up were 45.0% (32.5) for gabapentin initiators and 55.4% (32.4) for pregabalin. Among gabapentin and pregabalin initiators, 9.0% and 9.9%, respectively, remained on the treatment for at least 180 days (Table 2). Of the patients who remained on the treatment for 180 days, the dose was increased in only 5.7% of amitriptyline, 8.7% of duloxetine and gabapentin, and 13.9% of pregabalin initiators.
Changes in health care utilization rates
The mean number of outpatient visits and the proportion of patients with at least one hospitalization numerically decreased slightly after treatment initiation compared to before treatment initiation. The number of prescription drugs did not change, and the proportion of patients with at least one ED visit or at least one physical therapy visit increased after treatment initiation in all groups. In the amitriptyline and pregabalin pair, 4.2% of patients taking amitriptyline and 4.5% of patients taking pregabalin were started on another study drug during follow-up. Similarly, 5.6% of duloxetine and 7.1% of pregabalin initiators were prescribed another study drug during follow-up. In the gabapentin-pregabalin pair, 11.5% of gabapentin initiators and 13.8% of pregabalin initiators received another study drug during follow-up. Table 3 shows patients’ unadjusted health care utilization patterns before and after treatment initiation.
In multivariable Poisson and Cox regression analyses (Table 4), duloxetine initiation was associated with significantly decreased outpatient visits (RR, 0.94; 95% CI, 0.88 to 1.00), prescription drug use (RR, 0.94; 95% CI, 0.90 to 0.86), hospitalizations (RR, 0.75; 95% CI, 0.68 to 0.83), ED visits (RR, 0.85; 95% CI, 0.79 to 0.91), and physical therapy visits (RR, 0.81; 95% CI, 0.71 to 0.92), compared to pregabalin. Amitriptyline initiators had similar rates of health care utilization, except for the significantly lower rate of hospitalization (RR, 0.70; 95% CI, 0.62 to 0.80) and outpatient visits (RR, 0.91; 95% CI, 0.84 to 0.99), compared to pregabalin. Little difference in the rate of health care utilization existed between gabapentin and pregabalin initiators.
In a subgroup of patients who received at least 180 days of treatment, amitriptyline, duloxetine, and gabapentin initiators had lower rates of outpatient visits and specialist visits compared to pregabalin initiators (Table 5).
In this large, US population-based cohort study, patients newly prescribed amitriptyline, duloxetine, gabapentin, or pregabalin for fibromyalgia generally had high health care utilization before and after treatment initiation. None of the four fibromyalgia treatments seemed to have much impact on health care utilization; however, duloxetine initiators had significantly fewer outpatient and ED visits and hospitalizations, less prescription drug use, and fewer physical therapy appointments compared to pregabalin initiators. As reported in prior research, depression is a predictor of high health care use . The beneficial effect of duloxetine versus pregabalin seen in this study may be attributable to the fact that duloxetine is also an antidepressant and thus manages depressive mood in patients with fibromyalgia. Amitriptyline initiators also had fewer hospitalizations and physical therapy visits during follow-up than pregabalin initiators did.
Our study highlights several important issues in the management of fibromyalgia. First, fibromyalgia drugs appear to have little effect on reducing health care utilization by patients with fibromyalgia, as seen in previous studies of pregabalin, duloxetine, and TCAs [10-13]. Whether the lack of effect on health care utilization is due to inadequate benefits or concurrent benefit with new side effects of treatment is unclear. It is also difficult to determine whether the use of health care was inappropriately high or clinically necessary. Second, discontinuation rates are high across all study drugs. In previous studies, researchers have found that the effectiveness of these medications is limited , and many patients discontinue these medications due to side effects . Third, pharmacologic treatment alone may not be effective in fibromyalgia . Prior meta-analyses of various nonpharmacologic treatments such as massage therapy, aquatic physical therapy, balneotherapy, and hydrotherapy showed mixed results on fibromyalgia symptoms [25-28], although beneficial effects of aerobic exercise and cognitive behavioral therapy on pain and mood were consistently noted [29,30]. Nonetheless, in future longitudinal study, investigators should explore the utility of multifaceted approaches to the treatment of fibromyalgia, including patient education, exercise, and mental health services, with the goal of reducing health care utilization.
There are limitations to our study. First, selection of a drug was not random in this observational cohort. We used rigorous pharmacoepidemiologic methods in the study design and data analysis, such as the new user design and PS-matched analysis, to minimize confounding by indication [19,31]. However, confounding by unmeasured variables, such as severity or activity of fibromyalgia, functional status, other comorbidities, use of nonpharmacologic treatments, and socioeconomic status, can still be an issue. For example, differences existed in comorbidities prior to PS matching, such that duloxetine initiators had a higher proportion of diagnosis and treatment for depression (24%) and anxiety (21%) than the pregabalin group (12% and 13%, respectively) . If the presence or the severity of such conditions, which are related to use of health care systems, were incompletely controlled for, this could lead to residual confounding in our analysis. Second, there was also potential for misclassification, as we relied on diagnosis codes and drug-dispensing records to select patients with fibromyalgia and identify their comorbidities. Third, in this study, the reasons for use of health care resources (for example, outpatient and ED visits, hospitalization) were not available. Fourth, some of the visits to physicians or EDs may not have been related to fibromyalgia. Furthermore, the number of physician visits may not be an appropriate indicator of the effectiveness of fibromyalgia treatment. Fifth, we did not directly compare health care utilization between amitriptyline, duloxetine, and gabapentin initiators. Sixth, although it is possible that pharmacologic treatment has a beneficial effect on health care utilization in patients with fibromyalgia compared to no treatment, our study was designed to examine the comparative effectiveness of the four study drugs. Seventh, because access to health care is different in each health insurance plan in the United States, our results may not be generalizable to populations covered by different health insurance plans, particularly in different countries. Eighth, we did not examine the effect of combination therapy, but less than 10% of patients were put on combination therapy during the study period .
Fibromyalgia patients had substantial health care utilization before and after initiation of amitriptyline, duloxetine, gabapentin, or pregabalin and generally had high discontinuation rates of these medications. Fibromyalgia treatment did not substantially decrease health care utilization. Compared to pregabalin initiators, duloxetine initiators had less health care utilization during follow-up. In future studies, researchers should examine the effectiveness of pharmacologic treatments combined with nonpharmacologic interventions in treating chronic pain, as well as the use of health care resources in fibromyalgia.
US Food and Drug Administration
International Classification of Diseases, Ninth Revision, Clinical Modification
Nonsteroidal anti-inflammatory drug
Primary care physician
Proportion of days covered
Serotonin and noradrenaline reuptake inhibitor
Selective serotonin reuptake inhibitor
Wolfe F. The epidemiology of fibromyalgia. J Musculoskelet Pain. 1993;1:137–48.
Wolfe F, Ross K, Anderson J, Russell I. Aspects of fibromyalgia in the general population: sex, pain threshold, and fibromyalgia symptoms. J Rheumatol. 1995;22:151–6.
Wolfe F, Clauw DJ, Fitzcharles MA, Goldenberg DL, Katz RS, Mease P, et al. The American College of Rheumatology preliminary diagnostic criteria for fibromyalgia and measurement of symptom severity. Arthritis Care Res (Hoboken). 2010;62:600–10.
Clauw D, Crofford L. Chronic widespread pain and fibromyalgia: what we know, and what we need to know. Best Pract Res Clin Rheumatol. 2003;17:685–701.
Hassett AL, Gevirtz RN. Nonpharmacologic treatment for fibromyalgia: patient education, cognitive-behavioral therapy, relaxation techniques, and complementary and alternative medicine. Rheum Dis Clin North Am. 2009;35:393–407.
Häuser W, Bernardy K, Uçeyler N, Sommer C. Treatment of fibromyalgia syndrome with gabapentin and pregabalin: a meta-analysis of randomized controlled trials. Pain. 2009;145:69–81.
Häuser W, Walitt B, Fitzcharles MA, Sommer C. Review of pharmacological therapies in fibromyalgia syndrome. Arthritis Res Ther. 2014;16:201.
White L, Robinson R, Yu A, Kaltenboeck A, Samuels S, Mallett D, et al. Comparison of health care use and costs in newly diagnosed and established patients with fibromyalgia. J Pain. 2009;10:976–83.
Robinson RL, Kroenke K, Williams DA, Mease P, Chen Y, Faries D, et al. Longitudinal observation of treatment patterns and outcomes for patients with fibromyalgia: 12-month findings from the reflections study. Pain Med. 2013;14:1400–15.
Gore M, Tai KS, Chandran A, Zlateva G, Leslie D. Clinical comorbidities, treatment patterns, and healthcare costs among patients with fibromyalgia newly prescribed pregabalin or duloxetine in usual care. J Med Econ. 2012;15:19–31.
Peng X, Sun P, Novick D, Andrews J, Sun S. Real-world comparison of health care utilization between duloxetine and pregabalin initiators with fibromyalgia. J Pain Res. 2014;7:37–46.
Sun P, Peng X, Sun S, Novick D, Faries DE, Andrews JS, et al. Direct medical costs and medication compliance among fibromyalgia patients: duloxetine initiators vs. pregabalin initiators. Pain Pract. 2014;14:22–31.
Zhao Y, Sun P, Watson P, Mitchell B, Swindle R. Comparison of medication adherence and healthcare costs between duloxetine and pregabalin initiators among patients with fibromyalgia. Pain Pract. 2011;11:204–16.
Gore M, Tai KS, Chandran A, Zlateva G, Leslie D. Clinical characteristics, pharmacotherapy, and healthcare resource use among patients with fibromyalgia newly prescribed pregabalin or tricyclic antidepressants. J Med Econ. 2012;15:32–44.
Kim SC, Landon JE, Solomon DH. Clinical characteristics and medication uses among fibromyalgia patients newly prescribed amitriptyline, duloxetine, gabapentin, or pregabalin. Arthritis Care Res (Hoboken). 2013;65:1813–9. doi:10.1002/acr.22071.
Gore M, Sadosky A, Zlateva G, Clauw D. Clinical characteristics, pharmacotherapy and healthcare resource use among patients with fibromyalgia newly prescribed gabapentin or pregabalin. Pain Pract. 2009;9:363–74.
Deyo R, Cherkin D, Ciol M. Adapting a clinical comorbidity index for use with ICD-9-CM administrative databases. J Clin Epidemiol. 1992;45:613–9.
Schneeweiss S, Seeger J, Maclure M, Wang P, Avorn J, Glynn R. Performance of comorbidity scores to control for confounding in epidemiologic studies using claims data. Am J Epidemiol. 2001;154:854–64.
Rubin D. Estimating causal effects from large data sets using propensity scores. Ann Intern Med. 1997;127:757–63.
Parsons LS. Reducing bias in a propensity score matched-pair sample using greedy matching techniques. Sugi 26 Proceedings. Long Beach, CA; April 22-25, 2001 Cary, NC: SAS Institute; 2001. p. 214–26. http://www2.sas.com/proceedings/sugi26/p214-26.pdf. Accessed 30 January 2015.
Cameron AC, Trivedi PK. Regression analysis of count data. Cambridge, UK: Cambridge University Press; 1998.
McCracken C, Dalgard OS, Ayuso-Mateos JL, Casey P, Wilkinson G, Lehtinen V, et al. Health service use by adults with depression: community survey in five European countries. Evidence from the ODIN study. Br J Psychiatry. 2006;189:161–7.
Wolfe F, Walitt BT, Katz RS, Lee YC, Michaud KD, Häuser W. Longitudinal patterns of analgesic and central acting drug use and associated effectiveness in fibromyalgia. Eur J Pain. 2013;17:581–6.
Häuser W, Jung E, Erbslöh-Möller B, Gesmann M, Kühn-Becker H, Petermann F, et al. The German fibromyalgia consumer reports–a cross-sectional survey. BMC Musculoskelet Disord. 2012;13:74. doi:10.1186/1471-2474-13-74.
Nüesch E, Häuser W, Bernardy K, Barth J, Jüni P. Comparative efficacy of pharmacological and non-pharmacological interventions in fibromyalgia syndrome: network meta-analysis. Ann Rheum Dis. 2013;72:955–62.
Li YH, Wang FY, Feng CQ, Yang XF, Sun YH. Massage therapy for fibromyalgia: a systematic review and meta-analysis of randomized controlled trials. PLoS One. 2014;9:e89304.
Lima TB, Dias JM, Mazuquin BF, da Silva CT, Nogueira RM, Marques AP, et al. The effectiveness of aquatic physical therapy in the treatment of fibromyalgia: a systematic review with meta-analysis. Clin Rehabil. 2013;27:892–908.
Naumann J, Sadaghiani C. Therapeutic benefit of balneotherapy and hydrotherapy in the management of fibromyalgia syndrome: a qualitative systematic review and meta-analysis of randomized controlled trials. Arthritis Res Ther. 2014;16:R141.
Häuser W, Klose P, Langhorst J, Moradi B, Steinbach M, Schiltenwolf M, et al. Efficacy of different types of aerobic exercise in fibromyalgia syndrome: a systematic review and meta-analysis of randomised controlled trials. Arthritis Res Ther. 2010;12:R79.
Bernardy K, Klose P, Busch AJ, Choy EH, Häuser W. Cognitive behavioural therapies for fibromyalgia. Cochrane Database Syst Rev. 2013;9:CD009796.
Walker A. Confounding by indication. Epidemiology. 1996;7:335–6.
This study was supported by an investigator-initiated grant from Pfizer. SCK is supported by National Institutes of Health (NIH) grant K23 AR059677 and a research grant from Pfizer. YCL is supported by NIH grants K23 AR57578 and R01 AR064850 and a research grant from Forest.
This study was supported by an investigator-initiated grant from Pfizer. Pfizer had no role in the study design, data collection, data analysis, or writing of the manuscript, nor did it have a role in approval of the content of the submitted manuscript. Publication of this article was not contingent on the approval of Pfizer. The study was conducted by the authors independently of the sponsor. The sponsor was given the opportunity to make nonbinding comments on a draft of the manuscript, but the authors retained the right of publication and to determine the final wording of the article. YCL received research grants from Forest. She holds stock in Cubist, Merck, Perrigo, and Express Scripts. SCK received a research grant from Forest. JEL has no competing interests. None of the authors had any nonfinancial competing interests related to this study.
SCK and YCL conceived of the study. SCK, JEL, and YCL participated in the study design and data interpretation. SCK and JEL performed data analysis. SCK drafted the manuscript. All authors read and approved the final manuscript.