Stata release 15

To assess the impact of differences in methods with previous studies and between the two Spanish databases, four sensitivity analyses were performed. In BIFAP, an additional and broader case definition was considered, including all recorded skin cancer cases identified by codes and without any validation process. Additionally, cumulative duration was calculated similarly to the SIDIAP database, this is, regardless of the duration of the prescription. In SIDIAP socioeconomic status and alcohol abuse were further considered, as in BIFAP the former was not available and the information on the latter is scarce. Finally, the main analysis was conducted excluding ever users of amiloride in order to explore the specific contribution of this drug as it is the principal combination with HCTZ.
Secondary analyses were performed in BIFAP for other diuretics with suggested photosensitizing properties such as furosemide, indapamide, and amiloride 1 , but also for other antihypertensives (calcium-channel blockers (CCB), angiotensin II receptor blockers (ARB) and angiotensin-converting enzyme inhibitors (ACEi) 2 . These analyses were adjusted for HCTZ use.
All analyses were performed using STATA Release 15.1 (StataCorp, College Station, TX, USA) and R 35 . The study was approved by the scientific committees of BIFAP and SIDIAP, as well as by the Ethics Committee of IDIAPJGOL.

The analysis was designed to examine the effects of chronic medical conditions at baseline on the subsequent occurrence of pneumococcal disease and determine if these effects differed between the age groups. We calculated the number of pneumococcal disease cases for each chronic medical condition during the observational period and estimated the incidence rates per 10 0,0 0 0 patient-years. Cox proportional hazards models were constructed to estimate the hazard ratio (HR) of each chronic medical condition for pneumococcal disease occurrence. Separate models were used to analyze patients in the following age groups: 0-18 years, 19-49 years, 50-64 years, and ≥65 years. Kaplan-Meier curves were generated to examine the pneumococcal disease occurrence according to the number of chronic medical conditions in each patient (0 vs 1 vs 2 vs ≥3 conditions). Statistical analyses were performed using Stata Release 15.1 (StataCorp LLC, College Station, Texas). Twotailed P values below 0.05 were considered significant.

We used Microsoft SQL Server 2016 for data linkage and extraction, Stata release 15 for statistical analysis, and Microsoft Excel 2016 for producing tables and figures. At first, we performed the univariate analysis and calculated no-show rates in each variable. Then, we used the multivariate logistic regression to calculate the AORs of missed appointments in each variable. Within each variable, the item with the minimal value served as the reference. For nonfirst appointments, we further built a logistic regression model with the personal attendance pattern of the previous three appointments and those variables having item(s) with an AOR > 1.5 in the previous regression model. In the current study, we set the heuristic threshold of effect size of AOR at 1.5 to select both statistical significance and clinical meaningful predictors in our model. The 95% CIs were calculated based on normal approximation to represent statistical significances. The adjusted no-show rates were calculated for appointments of each characteristic as they were average on all other characteristics. Those rates were interpreted as the predicted probabilities of missed appointments if they were average in all other respects. A two-tailed level of 0.05 was considered statistically significant.