key: cord-0702393-ebheb1t4 authors: Madrid-García, Alfredo; Pérez, Inés; Colomer, José Ignacio; León-Mateos, Leticia; Jover, Juan A; Fernández-Gutiérrez, Benjamín; Abásolo-Alcazar, Lydia; Rodríguez-Rodríguez, Luis title: Influence of colchicine prescription in COVID-19-related hospital admissions: a survival analysis date: 2021-03-26 journal: Ther Adv Musculoskelet Dis DOI: 10.1177/1759720x211002684 sha: ee38789bf4c39d3af9690b5987ebd7c36db68243 doc_id: 702393 cord_uid: ebheb1t4 AIMS: To analyze the association between colchicine prescription and COVID-19-related hospital admissions in patients with rheumatic and musculoskeletal diseases (RMDs). METHODS: Patients attending a rheumatology outpatient clinic from a tertiary care center in Madrid, Spain, from 1 September 2019 to 29 February 2020 were included. Patients were assigned as exposed or unexposed based on whether they were prescribed with colchicine in their last visit to the clinic during the 6 months before the start of the observation period. Treatment changes during the observation period were also considered. The primary outcome was COVID-19-related hospital admissions between 1 March and 20 May 2020. Secondary outcome included COVID-19-related mortality. Several weighting techniques for data balancing, based and non-based on the propensity score, followed by Cox regressions were performed to estimate the association of colchicine prescription on both outcomes. DISCUSSION: The number of patients entered in the study was 9379, with 406 and 9002 exposed and unexposed follow-up periods, respectively. Generalized Boosted Models (GBMs) and Empirical Balancing Calibration Weighting (EBCW) methods showed the best balance for COVID-19-related hospital admissions. Colchicine prescription did not show a statistically significant association after covariable balancing (p-value = 0.195 and 0.059 for GBM and EBCW, respectively). Regarding mortality, the low number of events prevented a success variable balancing and analysis. CONCLUSION: Colchicine prescription does not play a significant protective or risk role in RMD patients regarding COVID-19-related hospital admissions. Our observations could support the maintenance of colchicine prescription in those patients already being treated, as it is not associated with a worse prognosis. PLAIN LANGUAGE TITLE: Colchicine influence in COVID-19-related hospital admissions Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes a myriad of clinical signs and symptoms, together with typical laboratory abnormalities that manifest as the disease COVID-19. 1 The severity of COVID-19 symptoms can range from very mild to severe, with a variable disease progression among patients. Although the progression is usually characterized by symptoms such as fever, fatigue, dry cough or loss of taste or smell, in severe cases the course of the disease can be fatal, with pneumonia, ischemic events, or kidney complications. RNA viruses, such as SARS-CoV-2, trigger the innate immune system by destroying epithelial barriers, triggering the complement system response, and activating the toll-like receptors. In severe cases, an excessive immune response can lead to high levels of pro-inflammatory cytokines promoting a "cytokine storm" by overstimulating inflammatory-related genes. Pro-inflammatory circulating cytokines and chemokines are increased in SARS-CoV-2. 2 One of the immediate effects of this excessive response is the respiratory system's inflammation inducing respiratory failure. In order to prevent or modulate this exaggerated inflammatory response, several clinical trials with already approved and well known therapeutic agents 3 that inhibit the inflammatory response, such as hydroxicloroquine 4 or tocilizumab, 5 are being carried out. However, none of these drugs, often used in the management of patients with rheumatic and musculoskeletal diseases (RMDs), seems to achieve the desired results when treating COVID-19. Colchicine, a drug employed in rheumatology to treat diseases such as gout, Behçet's disease or other autoinflammatory syndromes, inhibits multiple pro-inflammatory mechanisms, mediating an anti-inflammatory response. 6 Colchicine effects have also been broadly discussed in cardiovascular disease prevention. 7 Together with its anti-inflammatory properties, colchicine's extensive clinical experience supports its potential theoretical application to help manage COVID-19. Several clinical trials to assess this drug's efficacy are under way, 8 although no result has been published at the time of this manuscript submission. The aim of our study was to analyze whether colchicine prescription plays a role in COVID-19 severity (defined as COVID-19-related hospital admission and mortality following admission) when affecting patients with RMDs. The study was carried out at the Hospital Clínico San Carlos (HCSC), a public tertiary care center of Madrid, Spain, with a catchment area of nearly 400,000 subjects. We performed a longitudinal retrospective observational study from 1 March 2020 (when the first COVID-19-related hospital admission was reported at the HCSC) to 20 May 2020. Eligible patients were selected from the HCSC-COVID-19 cohort, a cohort of patients with RMDs seen at least once at the HCSC rheumatology outpatient clinic between 1 March 2019 and 1 March 2020. 9, 10 In addition, patients had to fulfill the following criteria: HCSC Ethics Review Board approval was obtained as a retrospective study and waiver of informed consent was obtained for the use of deidentified clinical records. Furthermore, the study was conducted in accordance with the Declaration of Helsinki. Demographical and clinical data of relevance were retrieved from a departmental electronic health record (EHR; Penelope), from the HCSC Rheumatology Department. Outcomes data were extracted from the HIS. SARS-CoV-2 polymerase chain reaction (PCR) diagnostic assays results were obtained from the HCSC Microbiology Service. The primary outcome, COVID-19-related hospital admission, was defined as an admission to the HCSC during the observation period due to COVID-19 compatible symptoms, and a positive SARS-CoV2 PCR test, regardless of the duration of the admission or the department where the patient was admitted. The secondary outcome, COVID-19-related mortality, was defined as the decease of a subject previously journals.sagepub.com/home/tab 3 admitted to the HCSC due to COVID-19 compatible symptoms and a positive SARS-CoV2 PCR test. Decease had to occur during admission or in the 30 days after discharge. All deceases in these patients were considered related to COVID-19. The main covariable was "colchicine prescription": a patient was considered as having received a colchicine prescription if in the most recent visit registered in the EHR of the 6 months before the observation period, or during any visit of the observation period, they had: • Received by their attending rheumatologist a Spanish Drug and Medical Products Agency code for colchicine (703194 and 733204), and/or • Received a prescription with colchicine as a concomitant treatment, meaning that it had been prescribed by another physician (such as primary care physician). A further description of the identification of patients under colchicine therapy can be found in the Supplemental material file "Colchicine exposure assessment" online. These patients were allocated to the "Exposed group". Those that did not fulfill these criteria were allocated to the "Unexposed group". Treatment changes during the observation period were considered. Therefore, these patients were analyzed considering two (or more) distinct follow-up periods, as exposed and as unexposed. Demographic and clinical-related variables collected during the patient's visits using our EHR were also included. For both outcomes, covariable data were taken from the most recent visit before the beginning of the time at risk for each outcome: for the primary outcome, covariable data were taken from the patient's last visit of the 6 months before the observation period and from the visits during the observation period when colchicine treatment was modified. For the secondary outcome, covariable data were taken from the last patient's visit before the COVID-19-related hospital admission (regardless of whether it took place in the last 6 months before the observation period or during the observation period). The covariables considered are described in the Supplemental file "Covariables". Continuous variables were summarized using median and interquartile range (IQR) and categorical variables using proportions. COVID-19-related admissions and mortality crude incident rates and 95% confidence intervals (CI) were presented as the number of events per 100 patient-years and were estimated by dividing the number of events that occurred during follow-up by the number of person-years of exposure. For COVID-19-related admissions, the total follow-up time comprised the elapsed time from 1 March 2020 until 20 May 2020. Patients contributed risk time until any of the following censoring points: (a) a change in the colchicine treatment, including the prescription of colchicine in a patient not treated with this drug, regardless of the cause, or colchicine discontinuation, also regardless of the cause; or (b) hospital admission at the HCSC. This means that for those who were admitted, the follow-up time would be defined as the elapsed time between the beginning of the study (1 March 2020) and the hospital admission. For those who were not admitted, the follow-up time would be defined as the elapsed time between the beginning of the study and the end of the study (20 May 2020). Those changing treatment had the same follow-up and they contributed to the group of exposed or unexposed depending on whether they were prescribed or not with colchicine. For COVID-19-related mortality, the total follow-up time comprised the elapsed time from the beginning of COVID-19-related hospital admission (which took place between 1 March 2020 and 20 May 2020) until the end of the study (20 June 2020). Patients contributed risk time until any of the following censoring points: (a) death; or (b) 30 days after discharge. Kaplan-Meier cumulative incidence curves were estimated to account for COVID-19-related admissions and mortality. Cox regression models were fitted to estimate the influence of colchicine prescription in the hazard of COVID-19-related admissions and mortality. We considered the influence of other covariables: those with a prevalence ⩾5% in the "exposed group" and a percentage difference between both groups ⩾30%; those covariables already associated with prognosis and/or mortality risk [such as hypertension, diabetes, cancer, dyslipidemia, cardiovascular, lung (interstitial lung disease/chronic obstructive lung disease) or thyroid diseases], regardless of their prevalence or percentage difference; and those covariables related to the use of disease-modifying antirheumatic drugs. For diagnoses and treatments covariables with a prevalence <5%, we combined those that were related to similar pathologies or treatments, in order to achieve this criterion. Supplemental Excel file "Initial variables" summarizes all of the initially considered variables and the merging process into groups. Details on the methods used to balance the distribution of variables can be found in the Supplemental file "Statistical analyses". Finally, Cox regression models including the variable "Colchicine prescription" were calculated and adjusted by the variables that remained unbalanced. We limited the number of variables in the multivariable model following the rule of 10 events per variable. [11] [12] [13] Our null hypothesis was that there was not association, in our cohort, between colchicine prescription and our primary endpoints. Results from the Cox models were expressed as hazard ratios (HRs) with 95% CIs. In order to verify if the proportional hazards assumption held for our variable of interest, the Schoenfeld residuals and the scaled Schoenfeld residuals were used. Statistical analyses were performed with R version 4.0.2 and WeightIt 14 and Cobalt 15 R packages. Figure 1 shows the inclusion criteria and patient flow. By 1 March 2020, 9379 patients (385 exposed and 8994 unexposed) were identified. During the observation period, 28 patients underwent a medication change (20 were prescribed with colchicine, seven discontinued the treatment, and one suspended and then resumed the medication). We considered that these patients had two/three follow-up periods, as exposed and as unexposed. None of the patients were readmitted to the HCSC during the following-up period. For the primary outcome COVID-19-related hospital admissions, 9408 treatment periods were analyzed. The median (IQR) number of days that exposed patients were under colchicine prescription by 1 March 2020 was 413 (7-1379) days. Discontinuations in the treatment of less than 3 months, taking place before the observation period, were disregarded. More than half of the patients had been taking colchicine for more than a year (204, 53%); one-quarter had been recently prescribed (⩽7 days; 99, 26%); the other quarter had been using colchicine for less than a year (82, 21%). Demographic and clinical-related variables at the last rheumatology outpatient visit in the 6 months before the observation period or at the visit during the observation period when colchicine treatment was modified can be found in Table 1 . We identified 132 COVID-19-related hospital admissions (at the HCSC): 12 (2.96%) in patients from the exposed group and 120 (1.34%) in patients from the unexposed group. The median (IQR) time to event (TTE) for COVID-19-related hospital admission (time elapsed from 1 March until outcome or censoring) was 28.5 (20-38) days; 28.5 (20-38) days in the unexposed and 28 (19.5-36.5) days in the exposed groups. The median (IQR) length-of-stay of admitted patients was 9 (5-16) days; 8.5 (5-15.25) and 15 (7.75-20.25 ) days in the unexposed and exposed groups, respectively. Demographic and clinical-related variables at the last visit before COVID-19-related hospital admission can be found in Table 2 . Regarding COVID-19-related deaths, 34 were identified; four (0.99%) and 30 (0.33%) in the exposed and unexposed groups, respectively. The median (IQR) TTE for COVID-19-related mortality following hospital admission (time elapsed from admission until outcome or censoring) was 6 (2-9.75) days; 6 (2.25-9.75) days in the unexposed group and 4.5 (0.75-10.5) days in the exposed group. There was only one patient with a COVID-19-related decease 19 days after discharge. Supplemental Tables S1 and S2 show the mean differences for the selected covariables, unadjusted and balanced (adjusted, trimmed, and stabilized) using different methods based and not based on the Propensity Score (PS), respectively. For the methods based on the PS, Generalized Boosted Models (GBMs) showed the lowest number of unbalanced variables (n = 3). In addition, we selected the stabilized weights as they showed journals.sagepub.com/home/tab 5 the lowest coefficient of variation 0.296. Regarding non-PS methods, the three tested were able to balance all variables. We selected the Empirical Balancing Calibration Weighting, as it showed the lowest maximum mean difference compared with the other two methods. In addition, we selected the stabilized weighting as it showed the lowest coefficient of variation 5.706. In the unbalanced analysis, colchicine prescription was associated with a statistically significant higher hazard of COVID-19-related hospital admission [HR (95% CI) = 2.237 (1.236-4.049), p-value = 0.008]. Kaplan-Meier curves for COVID-19-related hospital admission can be found in Figure 2 . Table 3 shows the association of colchicine prescription and COVID-19 after balancing using PS and non-PS based methods with stabilized weights. Proportional hazard assumption test results are shown through the Schoenfeld residuals p-values. As can be seen, all models fulfilled the proportional hazard assumption (p-value > 0.05). In the selected balancing methods, we observed that colchicine prescription was still associated with a higher HR although statistical significance was lost. Supplemental Table S3 shows the complete models, including the unbalanced covariables. In the unbalanced analysis, colchicine prescription was not associated with a statistically significant higher hazard of COVID-19-related decease [HR (95% CI) = 1.371 (0.482-3.899), p-value = 0.554]. Kaplan-Meier curves for COVID-19-related death can be found in Figure 3 . Supplemental Table S4 shows the mean differences for the selected covariables, unadjusted and balanced (adjusted, trimmed, and stabilized) using different methods based on PS. Singularity and convergence issues appeared when testing non-PS based methods, likely due to the low number of events (n = 34 deaths). After balancing, the number of unbalanced variables was much higher, compared with the primary outcome analysis. The method able to balance the highest number of variables was the trimmed GBM, leaving 13 covariables unbalanced. The overall number of events and the number of covariables that would have to be included in the multivariate analysis [11] [12] [13] prevented us from carrying out a proper analysis of the secondary outcome. Our study aimed to assess the role that colchicine could play in a cohort of RMD patients in terms of COVID-19-related hospital admission and mortality. In this setting, colchicine prescription was not associated with a statistically significant higher hazard of admission after covariable balancing. Regarding COVID-19-related mortality, the number of unbalanced covariables remained too high with regard to the number of events, and, therefore, the colchicine effect was not analyzed for this outcome. Colchicine is one of the several drugs that are being currently tested for efficacy in COVID-19 16, 17 due to its anti-inflammatory effects: colchicine is able to inhibit the microtubule polymerization, interfering with inflammatory pathways related to the innate immune system, such as neutrophil recruitment and adhesion, superoxide production, and inflammasome, rho-associated protein kinase or NF-κB protein complex activation. 6 Several observational or uncontrolled studies have shown a beneficial effect of colchicine when used with the intention to treat COVID-19. Sandhu et al. 18 showed a lower mortality, a lower rate of intubations, and a higher number of discharges among hospitalized patients who received colchicine (n = 53) compared with standard of care (n = 144). In addition, Scarsi et al. 19 also observed a better survival rate in 122 patients admitted with severe COVID-19 (pneumonia and acute respiratory distress syndrome) and treated with colchicine, compared with 140 patients treated with standard of care. Brunetti et al. 20 performed a propensity score matched comparison of admitted patients receiving colchicine (n = 33) versus standard care (n = 33). Those receiving colchicine showed a lower 28-days all-cause mortality compared with those not receiving this medication. Although in the three studies the allocation of colchicine was not random, both the intervention and control groups were comparable. Colchicine has also been used in a community setting: Della-Torre et al. 21 reported a favorable outcome in an uncontrolled case series of nine domiciliary consecutive COVID-19 patients treated with this drug. In addition, data from two randomized controlled trials has been released: Lopes et al. 22 presented an interim analysis of a single-center randomized, double-blinded, placebo-controlled clinical trial including 35 patients (18 placebo and 17 colchicine). Subjects in the colchicine group showed a shorter time of need for supplemental oxygen, shorter hospitalization time, and a higher discharge rate. Finally, the results from the COLCORONA study have been released, 17 apparently showing a significant reduction in mortality. Comparing these results with those in our study, we did not observe a significant difference in the crude mortality rate among those exposed and unexposed to colchicine. This finding could be related to the fact that we were not able to properly balance the baseline characteristics of those exposed and unexposed to colchicine for this outcome, and, therefore, the propensity score analysis could not be carried out. On the other hand, we did not use colchicine as a specific treatment for COVID-19, and therefore it may have been discontinued after admission. Unfortunately, we do not have data regarding who continued treatment or who did not. We want to point out that the use of this drug was not included in the COVID-19 treatment protocols used at our center, and therefore its prescription during admission was not promoted. Finally, it is possible that colchicine's positive effect is related to the use of higher dosages, similar to those used to treat gout flares (⩾1 mg/day). In our case, most patients were receiving a lower dosage (0.5-0.6 mg/day) when outpatients, and among those who were admitted and continued their treatment, it is very likely that colchicine was not prescribed at higher dosages, and patients just continued with their previous dosages. Other studies have shown that the chronic use of colchicine was not associated with a higher risk of severe COVID-19: in a cohort of patients with familial Mediterranean fever, 23 a similar proportion of severe SARS-CoV-2 infection compared with the general population was reported. Furthermore, a group of five patients with iatrogenic allogenosis 24 seemed to present with a milder COVID-19 course compared with their close contacts who also developed the disease. In our study, analyzing the association between colchicine and COVID-19 in terms of chronic exposure (at least for most of the patients that had been receiving this medication by 1 March 2020), we observed similar results, in the sense that colchicine use was not associated with a significantly different COVID-19 prognosis (in our case defined as COVD19-related admission) compared with those patients not receiving this drug. As previously mentioned for COVID-19-related mortality, we could hypothesize that the dosage used by those patients taking colchicine was not high enough to improve COVID-19 prognosis. Although the published American College of Rheumatology guidelines 25 do not mention nor provide colchicine-related recommendations, based on our results and previous evidence we suggest not discontinuing colchicine therapy in patients who were prescribed with this drug, as the previous use of colchicine does not seem to constitute a risk factor for severe COVID-19. This recommendation may be of particular relevance for patients with vascular disease, since COVID-19 affectation and severity in this group of patients has been demonstrated and discussed since the beginning of the pandemic. [26] [27] [28] The unjustified colchicine discontinuation could negatively affect the patient's health and render them prone to potential rheumatic flares or ischemic complications. Regarding the former, we have to take into account the possible ramifications of a flare in a quiescent or controlled disease considering the current situation due to the pandemic. Not only were our healthcare systems almost at full capacity (emergency departments and hospitalizations) and with reduced or even at times suspended outpatient care, but also citizens were subjected to sometimes stringent limitations of movement (including home confinement). 29 As an example, if we consider autoinflammatory syndromes, flares can be associated with very severe and life-threatening manifestations. In addition, typically they are associated with fever, one of the more common symptoms of COVID-19. In other diseases, such as gout, an otherwise mild flare could be associated with severe limitations, as patients are unable to access healthcare and/or family caregivers are not able to provide support. Our study has a series of limitations, the main one being its retrospective nature. Treatment allocation was not carried out at random, and although we carried out different methods to increase the comparability of both groups, prescription bias cannot be completely compensated. However, we observed similar results regardless of the method used to balance the covariables between patients prescribed or not with colchicine, which increases our confidence regarding the validity of our results. On the other hand, the weighting strategies employed do not prevent bias from unmeasured confounders. Some studies have highlighted different ways to reduce this bias in observational studies, such as doing sensitivity analyses with external validation data. 30, 31 Although the considered variables in our study were only those accessible in our outpatient EHR, and we did not carry out an external validation, our results were in line with previously published observations in populations chronically exposed to colchicine. Another limitation is that the outcome data were extracted from a single center. On the one hand, there is the issue of the generalizability of the results, being obtained from a homogeneous group of patients with RMDs from a tertiary care center. It is important to consider that the effect of colchicine may be different in patients with other diseases that are also chronically exposed to this drug. On the other hand, patients may have been admitted to other hospitals due to COVID-19, or they may have died of COVID-19 also in other centers or before attending the HCSC. Because we do not have access to those data, we may have wrongfully classified patients as not having developed the outcomes. This may also have introduced a selection bias in our study, as we cannot assess whether the characteristic of those admitted to the HCSC and those that may have been admitted to other centers were similar or not. Furthermore, due to the same reasons, readmissions of the studied population in other hospitals during the follow-up period could not be assessed. Another limitation of our study is that patient adherence to colchicine treatment was not assessed. We were able to analyze prescription, but not use, and we assumed that those prescribed with colchicine were taking the medication. Furthermore, we have prescription data at the time of the visits from the EHR. Colchicine treatment may have been discontinued or initiated, but not registered at the EHR. Nevertheless, our study also presents several strengths: we have included patients with a board spectrum of RMDs, reflecting the clinical practice from our health area. In addition, more than 400 episodes of colchicine treatment were analyzed, and several weighting techniques applied, with consistent results. Our results suggest that colchicine does not play a significant protective either a risk role in COVID-19 severity, in a population of RMD patients. Whereas in the unbalanced analysis colchicine prescription was associated with a statistically significant higher hazard of COVID-19-related hospital admission, in the balanced analyses, the statistical significance was lost. Therefore, we conclude that there is no significant association between colchicine exposure and COVID-19 severity. These results should be interpreted as an initial approach to the possible role that colchicine could play in a COVID-19 exposed population. In light of our observations, colchicine treatment seems safe in these circumstances, and it likely should be continued in RMD patients. LRR and LAA conceived and designed the study. IP, JIC, LLM, JAJ, and BFG collected data. AMG, LAA, and LRR performed the data analysis and interpreted the data. All of the authors were involved in the drafting and/or revising of the manuscript. The authors declare that there is no conflict of interest. The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Instituto de Salud Carlos III (ISCIII), Ministry of Health, Spain (CPII17/00014; PI18/01188; CP16/00916; and RD16/0012/0014) and cofounded by el Fondo COVID-19 with pulmonary involvement. An autoimmune disease of known cause The cytokine storm in COVID-19: an overview of the involvement of the chemokine/chemokinereceptor system 15 drugs being tested to treat COVID-19 and how they would work Observational study of hydroxychloroquine in hospitalized patients with Covid-19 Effective treatment of severe COVID-19 patients with tocilizumab Colchicine update Colchicine for prevention of cardiovascular events Studies covid19 and colchicine Hospital admissions in inflammatory rheumatic diseases during the COVID-19 pandemic: incidence and role of disease modifying agents Risk factors for hospital admissions related to COVID-19 in patients with autoimmune inflammatory rheumatic diseases Importance of events per independent variable in proportional hazards analysis. I. Background, goals, and general strategy Importance of events per independent variable in proportional hazards regression analysis. II. Accuracy and precision of regression estimates A simulation study of the number of events per variable in logistic regression analysis WeightIt: weighting for covariate balance in observational studies covariate balance tables and plots Characteristics of clinical trials evaluating cardiovascular therapies for coronavirus disease 2019 registered on ClinicalTrials.gov: a cross sectional analysis. Am Heart J. Epub ahead of print 26 Colchicine Coronavirus SARS-CoV2 Trial (COLCORONA) (COVID-19) A case control study to evaluate the impact of colchicine on patients admitted to the hospital with moderate to severe COVID-19 infection Association between treatment with colchicine and improved survival in a single-centre cohort of adult hospitalised patients with COVID-19 pneumonia and acute respiratory distress syndrome Colchicine to weather the cytokine storm in hospitalized patients with COVID-19 Treating COVID-19 with colchicine in community healthcare setting Beneficial effects of colchicine for moderate to severe COVID-19: an interim analysis of a randomized, double-blinded, placebo controlled clinical trial Clinical course of COVID-19 in a cohort of 342 familial Mediterranean fever patients with a long-term treatment by colchicine in a French endemic area Epub ahead of print 16 American College of Rheumatology guidance for the management of rheumatic disease in adult patients during the COVID-19 pandemic: Version 2. Arthritis Rheumatol. Epub ahead of print COVID-19 and the cardiovascular system COVID-19 and cardiovascular disease Asociación entre el daño miocárdico y el pronóstico de pacientes hospitalizados por COVID-19 con y sin cardiopatía. Registro CARDIOVID. Rev Española Cardiol. Epub ahead of print 16 Treatment adherence during the COVID-19 pandemic and the impact of confinement on disease activity and emotional status: a survey in 644 rheumatic patients Adjusting effect estimates for unmeasured confounding with validation data using propensity score calibration Sensitivity analysis and external adjustment for unmeasured confounders in epidemiologic database studies of therapeutics