Introduction

Every year, between 15,000 and 20,000 individuals in Germany suffer a stroke caused by arterio-arterial embolisms due to atherosclerotic plaque of the carotid artery [1]. Standard revascularization procedures such as carotid endarterectomy (CEA) and carotid artery stenting (CAS) are available for both primary and secondary prevention. Evidence-based recommendations on the diagnosis, indications, choice of procedure, and follow-up have been published in national and international guidelines [1, 2].

All CEA and CAS procedures performed in German hospitals are subject to statutory external quality assurance (eQA), which includes indicators of outcome quality and details of inpatient treatment [3]. Secondary data analyses on specific questions, such as the link between age, gender, annual center volume, surgical technique, and the perioperative stroke and death rates have already been published [4,5,6,7,8]. Due to data protection regulations, it was not possible to carry out hospital-specific evaluations. In addition, eQA does not record information on patients’ place of residence, rendering it essentially impossible to make a differentiated analysis in terms of location.

An analysis of the frequency of surgical procedures in selected service areas revealed a moderate variation between the 402 German districts and cities in terms of procedure frequency for appendectomies and cesarean sections standardized for age and gender [9]. In contrast, coronary interventions and pacemaker implantation showed a high systematic variation [9]. So far, the small-area frequency of CEA and CAS procedures carried out in Germany was not analyzed. Although a certain regional variation should be statistically expected and conclusions should be drawn with caution [10], regional differences in the frequency of surgical procedures can be used to generate hypotheses in order to identify determinants of suboptimal care processes (e. g., access to the health system, supplier-induced increase in demand) [11,12,13,14]. Therefore, the aim of this study was: (a) to analyze the small-area frequency of CEA and CAS procedures for carotid stenosis in Germany in terms of place of residence, and (b) to identify regional characteristics (exploratory approach) associated with the regional frequency of CEA and CAS.

Methods

Data source and processing

The German diagnosis-related group (G-DRG) statistics compiled by the German Federal Statistical Office (Statistische Bundesamt, StBA) for the period 2012–2014 were evaluated using controlled remote data processing (CRDP; [15]). The methods used have already been described in detail elsewhere [16, 17]. The data are stored on StBA servers in accordance with the StBa data protection regulations. Staff at the Research Data Center ensured compliance with data protection regulations. The population consisted of all patients treated in a hospital on an inpatient basis. The study, which covered the whole of Germany, was approved by the Ethics Committee at the Technical University of Munich and was conducted in accordance with current guidelines [18, 19]. The assignment of regional variables (Tables 12 and 3) to the cases was carried out on a place of residence basis using the official municipality key. The case-specific DRG was used to create a link to the G‑DRG catalog. The sex and age group-specific number of inhabitants (www.genesis.destatis.de), the type of district and region in terms of residential areas (German Federal Institute for Research on Building, Urban Affairs and Spatial Development, Bundesinstitut für Bau, Stadt- und Raumforschung), data on health system infrastructure (www.inkar.de, German Association of Statutory Health Insurance Physicians, Kassenärztliche Bundesvereinigung), data on the prevalence of risk factors (www.versorgungsatlas.de, www.gbe-bund.de), the DRG casemix valuation ratio, and the average maximum length of stay were linked. The codes for identifying comorbidities, procedures, and secondary outcomes are listed in eTable 1.

Table 1 Characteristics of the patient cohort (aggregated 2012–2014)
Table 2 Diagnosis, management, and outcome (aggregated 2012–2014)
Table 3 Correlation between regional characteristics and procedure frequency, as well as a comparison of regions with high and low procedure frequencies

Case definition

All DRG cases in the reporting years 2012–2014 for which CEA or CAS was coded with the German Operations and Procedures Key (Operationen- und Prozedurenschlüssel, OPS), at the same time as being coded with the principal or secondary diagnosis carotid stenosis I65.2 (International Classification of Diseases ICD-10), were included. Cases with anatomically or etiologically unspecific codes (I65.3/8/9, I63.0/2) were excluded. eFigure 1 shows the patient flow chart, while eTables 2–3 list the characteristics of excluded patients. All analyses refer to one hospital case (analysis unit). With respect to the OPS codes, inclusion and exclusion criteria correspond to the official QA filter of the eQA (module 10/2, sqg.de, eTable 4–7), with the exception of minimal adjustments in the case of a single OPS code (5396.00–.03, blood vessel transposition).

Primary study variable and spatial resolution

The frequencies of CEA and CAS procedures (= hospital incidence, primary study variable) were calculated as indirect age and sex-standardized values [17]. All frequencies refer to patients’ place of residence (not the place of treatment) and to 100,000 inhabitants. Due to the expected low incidence at district and town level (NUTS3 Level, Nomenclature des Unités Territoriales Statistiques), the reporting years 2012–2014 were aggregated. In order to show CEA and CAS separately, it was necessary to reduce spatial resolution to spatial planning regions (so-called “Raumordnungsregionen”, http://www.bbsr.bund.de) in order to avoid blocks for data protection reasons and to strike a balance between spatial resolution, currentness of data, and statistical accuracy.

Statistical analysis

Categorical data are given as absolute frequencies with percentages. The median, 25% and 75% quantiles (Q1, Q3) were calculated uniformly for metric variables. The 95% confidence intervals (95% CI) for ranked forest plots were calculated using the method described by Sahai and Khursid [20]. Regional variation was calculated using the systematic component of variation (SCV) and the proportion of SCV in the total variation (SCV proportion; [21,22,23]). Global Moran’s I was calculated to analyze spatial autocorrelation and Getis-Ord Gi* was calculated to identify hot and cold spots. Standard ArcGIS algorithms (version 10.1, Environmental Systems Research Institute, Redlands, CA, USA) were used to this end. Data processing and analysis (in the form of CRDP, using the “NewVar Macro” version 1.2 provided by the StBA) was carried out with the SAS statistical program, version 9.2, for Microsoft Windows (Copyright © 2015 SAS Institute Inc., Cary, NC, USA). The graphic representation of data was prepared using Microsoft Excel and the statistics program R (version 3.2.1; The R Foundation, https://www.r-project.org). Statistical significance was set at α = 0.05.

Hypothesis-generating exploratory analyses

To investigate the relationship between regional characteristics (Table 3) and the frequency of CEA/CAS procedures, an exploratory (ecological study design) rank correlation analysis (Spearman) was performed on a NUTS3 level (n = 402 districts and towns). To enable the reader to estimate the correlations on the original scale, the raw values of NUTS3 regions with low procedure frequency (lower decile) and high procedure frequency (upper decile) were additionally compared. The choice of statistical method was guided by the empirical distribution evaluated with Q‑Q plots (t-test in the case of normal distribution, otherwise the Wilcoxon rank-sum test).

Results

Characteristics of patients included

A total of 88,182 procedures, 73,042 CEA and 15,367 CAS, were included in the study (Table 1) between 2012 and 2014. In 227 cases, both procedures were documented. Of the patients two thirds were male, and the median age of the total cohort was 72 years. The most common comorbidities included arterial hypertension (82%), diabetes mellitus (30%), and coronary heart disease (29%). The majority of patients lived in an urban district (39%) or a city (27%). Treatment most frequently took place in independent cities (48%), followed by urban districts (27%). eFigure 2 shows the cross tables on place of residence and place of treatment in the case of CEA and CAS. The median annual number of cases per hospital was 82 (Table 2).

Hospital mortality was 1.0%. Since the times at which a stroke, a transient ischemic attack, or other symptoms occur are not coded in DRG data, it was not possible to distinguish between index events and complications; therefore, this data was not evaluated. Patients spent a median of 6 days in hospital and were discharged home normally in 90% of cases or transferred to a rehabilitation facility in 5% of cases. The mean (arithmetic) case mix index was 1.51. Further details are provided in Tables 1 and 2.

Regional variation

The overall frequency (CEA + CAS at NUTS3 level) varied between 13.2 per 100,000 inhabitants (city of Augsburg) and 89.2 procedures (city of Wilhelmshaven). The CEA frequency varied at the planning region level from 13.5 (southwest Schleswig-Holstein, Germany) to 48.3 (east Upper Franconia, Germany) CEA per 100,000 inhabitants. The CAS procedures were most rarely carried out in the Schwarzwald-Baar-Heuberg region (1.55 per 100,000 inhabitants) and most frequently in the southwest region of Schleswig-Holstein (17.9 per 100,000 inhabitants). A cartographic representation is shown in Fig. 1 and a ranked forest plot in Fig. 2. Spatial statistical analysis revealed a clustered distribution pattern for both the total frequency (CEA + CAS) and the frequency of CEA and CAS (positive global autocorrelation, p-value < 0.001 each). Non-random variation was 95% (SCV 7.5) for all procedures, 95% (SCV 5.7) for CEA, and 97% (SCV 19.0) for CAS.

Fig. 1
figure 1

Age and sex standardized frequency of procedures per 100,000 inhabitants in the period 2012–2014. Spatial resolution refers to patients’ place of residence (districts and cities in a, spatial planning region in b and c). a Procedure frequency (CEA + CAS; global spatial autocorrelation: Moran’s I = 0.43; p < 0.001 clustered pattern). b Procedure frequency (CEA; global spatial autocorrelation: Moran’s I = 0.47; p < 0.001 clustered pattern). c Procedure frequency (CAS; global spatial autocorrelation: Moran’s I = 0.32; p < 0.001 clustered pattern)

Fig. 2
figure 2

Ranked “forest plot” of age and sex standardized procedure frequency per 100,000 inhabitants (inh.) in the period 2012–2014. Spatial resolution refers to patients’ place of residence (NUTS3 region in a, spatial planning region [SPR] in b and c)

The highest percentage of CEA procedures was found in Altmark (95.8%), while the lowest were found in southwest Schleswig-Holstein (43.1%), central Schleswig-Holstein (59.1%), and western Rhineland-Palatinate (59.4%). The distribution of the CEA percentage is shown in eFigure 3.

Exploratory analysis (to generate hypotheses)

The regional prevalence of risk factors was significantly and positively correlated with the overall frequency of CEA and CAS, as was the concentration of internal medicine beds (Table 3). Better income levels correlated with a lower frequency of procedures and debts with a higher frequency. Fewer CEAs and CASs were carried out in regions with high levels of inward and outward migration and high numbers of self-employed individuals, while more CEAs and CASs were performed in regions with high numbers of individuals in need of care. A comparison of regions with high (upper decile) and low (lower decile) procedure frequencies revealed significant differences for the same factors (with the exception of the concentration of internal medicine beds). Neither the type of district nor type of region in terms of residential areas was significantly associated with the frequency of procedures. There was also no significant correlation with characteristics of healthcare infrastructure, such as density of beds, resident vascular surgeons or angiologists. See Table 3 for further details.

Discussion

This study shows that the frequency of CEA and, in particular, CAS procedures for carotid stenosis in Germany is subject to significant regional variation. Spatial analysis revealed a clustered distribution pattern. In addition to the regional prevalence of diabetes mellitus, smoking, and obesity, socioeconomic factors, such as income and debts were associated with the overall frequency of procedures, but not indicators of health infrastructure.

A total of 73,000 CEA and 15,000 CAS procedures were included in the study period. When including the 7600 CEA and 4300 CAS performed in patients without specifically coded carotid stenosis (I65.2) (eTable 2, Cohort 2), these figures from the DRG statistics are consistent with eQA figures [3]. Although the two reporting channels are not independent of each other, they are subject to different control processes and incentives to document data. While it is not possible to assign the administrative DRG data to the indication groups A (asymptomatic), B (symptomatic, elective), and C (acute stroke, crescendo TIA, stroke-in-evolution, interventions under special conditions) clinically defined in eQA, it can be assumed that the majority of patients included in this study can be assigned to elective procedures. The higher mortality rate (4% vs. 1%) and higher Elixhauser comorbidity score (7 vs. 3) in the excluded patients (eTable 2, Cohort 2) points to this; however, due to different coding practices (e. g., crescendo TIA) or absent or non-specific administrative codes, it is not possible to rule out vagueness in the distinction between cohorts. Restricting patients to those that had both carotid stenosis and a matching procedure code results in a smaller but more clearly defined cohort and better reflects the elective care processes that do not overlap with emergency procedures.

Secondary data analysis of eQA data showed that nearly two thirds of CAS patients were classified in American Society of Anaesthesiology (ASA) categories I and II, while CEA patients were classified with ASA III in over two thirds of cases [3]. The latter is in line with the authors’ experience, whereas the former appears counterintuitive, particularly with respect to guideline recommendations [1]. In this study, however, the median Elixhauser comorbidity score was slightly higher for CAS compared with CEA (4 vs. 3). Although it is not possible to make a direct comparison between clinical evaluation with ASA categories and the administrative Elixhauser comorbidity score, a systematic underestimation of ASA category in CAS in the context of eQA seems probable. This should be the subject of critical discussion at the German Institute for Quality Assurance and Transparency in Healthcare (Institut für Qualitätssicherung und Transparenz im Gesundheitswesen, IQTiG) and the regional authorities for quality assurance, since ASA category is used for risk adjustment. Likewise, a consideration of relevant risk factors, such as hypertension, pre-existing cardiovascular diseases, or statin therapy, should be discussed [1].

The frequency of CEA varied between 14 and 48 per 100,000 residents (extremal coefficient 3.6). This value is comparatively low compared with the USA, for which a 9- to 10-fold variation between states has been described; it corresponds approximately to the variation in Canada [12,13,14]. The SCV for CEA was 5.7, which can be interpreted as a high variation [22] and lies (for comparison purposes) between that of cesarean sections (2.3), appendectomies (3.9), and percutaneous coronary interventions (7.8; [9]). In comparison, CAS has an extremely high variation with an SCV of 19. The fluctuation range (1.6–18 per 100,000 inhabitants) with an extremal coefficient of 12 is also higher than, e. g., in the USA (6–7; [13, 14]). The most recent study on regional distribution in the USA described an extremal coefficient of 8.6 for CEA and 13.9 for CAS, whereby these values relate to the Medicare population only [24].

Although standardization minimized sex and age effects in the present study, some variation still remained. Amongst others, this may be explained by differences in the prevalence of unobserved factors. The overall frequency of procedures was particularly high in, e. g., northern Bavaria, Mecklenburg-Western Pomerania, and North Rhine-Westphalia. A possible explanation for this might be the varying distribution of prevalence of risk factors, such as diabetes mellitus, smoking, or obesity, which could only be correlated on a regionally aggregated level (ecological study design), and the significant positive correlations identified (Table 3) must be interpreted with caution. As Anderson concluded, variation in the utilization of processes is not necessarily a bad thing per se [10] and can be the result of, e. g., unobserved factors.

The significant variation seen for CAS and the statistical outliers (3 out of 96 spatial planning regions) in the CAS share (eFigure 3) are, in the authors’ view, not only due to unobserved confounders, but presumably also due to differing indications; however, in the light of current guidelines, it is not possible to evaluate indication practices, particularly with administrative data. Different indication practices alone are probably not able to explain the spatially clustered distribution patterns, which, for instance, show a high frequency for CEA and a low frequency for CAS in north-eastern Bavaria and the opposite in Schleswig-Holstein, since the distance between the place of residence and treatment was less than 22 km linear distance for 75% of patients (Table 1).

An analysis of spatial variation could be used to generate hypotheses, e. g., using geocoded eQA data, to analyze where, how, and why carotid stenosis was treated—and with what outcome [10, 11]. It must be borne in mind that conservatively treated patients have not been covered as yet by eQA and only rarely by DRG data. Therefore, conclusions on the accuracy of the indication cannot be drawn with the current eQA procedure.

A positive correlation between smoking prevalence and CEA procedures was also observed in the USA [25]. Economic indicators such as debts were also positive, whereas household income correlated negatively with the frequency of procedures in Germany. The higher the affluence in a region, the fewer CEA and CAS procedures were carried out. A similar correlation has been described in the USA [25] and is consistent with other studies showing a positive link between high socioeconomic status and low cardiovascular morbidity [26].

In contrast, no correlation was seen between the frequency of CEA and CAS and parameters of health system infrastructure, such as the density of hospitals, general practitioners, vascular surgeons or angiologists. This is consistent with studies in the USA, which also found no correlation between the density of vascular surgeons and CEA frequency [25]. In that particular study, a link was found between a high concentration of beds and CEA frequency, which, in the German data was shown only for the density of internal medicine beds. Ultimately, however, there is no evidence of inadequate vascular health care services in Germany.

In summary, the study shows the extent of heterogeneity both in the indications and in the choice of method to treat carotid stenosis in Germany. Although natural variation and actual differences in incidence are likely, these cannot fully explain the observed variation. This raises the question of the extent to which national S3 guidelines on carotid stenosis have been adopted and the social question on nationwide, demand-oriented treatment for everyone.

Particularly with respect to the German Government’s “quality offensive” and the IQTiG’s efforts to improve eQA, this analysis, which in the authors’ opinion captures the most important, practice-relevant quintessence of the issue, raises the question of whether non-consideration of regional aspects and separating eQA analyses according to CEA and CAS represent a relevant weakness in the quality assurance process. Therefore, the inclusion of regional parameters in the evaluation of eQA data should be examined in order to assess outcome quality from the perspective of individual service providers, and not separately from the quality of indications, choice of procedure, and patient selection from the point of view of the regional patient population.

Limitations

The basic limitations of the data and methods used for this analysis have already been described [17, 27] and are discussed in more detail in the online supplement (eMethods). The most relevant limitations are the following:

  • The data are not findings made in the clinical setting, but administrative claims data for the purposes of hospital reimbursement.

  • Since clinical details, e. g., the degree of stenosis or initial neurological symptoms, are not coded in the DRG data, no conclusions can be drawn on the quality of the indications, choice of procedure or guideline conformity.

  • The StBA’s DRG statistics do not document which diagnoses were already present on admission, making it impossible to reliably differentiate between comorbidity and complication; similarly, it was not possible to measure neurological outcome.

  • The follow-up period covered only inpatient stays.

  • All analyses refer to patients’ place of residence; an analysis of the place of treatment on the level of NUTS3 or regional policy region was not possible for data protection reasons.

  • Exploratory analyses were performed to generate hypotheses and could only be carried out on an aggregated level. Therefore, it is not possible to rule out an ecological fallacy.

Practical conclusion

  • The total number of CEA and CAS procedures in the DRG statistics showed good agreement with eQA data. The fact that, although the two reporting channels are subject to different control processes yet are not independent of each other, points to complete data collection by eQA.

  • In relation to districts and towns, the overall age and gender-standardized incidence of CEA and CAS in carotid stenosis varied between 13 and 89 per 100,000 inhabitants.

  • The regional frequency of all CEA and CAS procedures demonstrated a positive spatial autocorrelation and, thus a clustered spatial pattern of distribution. The CEA and CAS were frequently performed in northern Bavaria, Mecklenburg-Western Pomerania, and North Rhine-Westphalia (44–89 per 100,000 inhabitants), whereas they were carried out less frequently in Baden-Württemberg, eastern Lower Saxony, and Schleswig-Holstein (13–27 per 100,000 inhabitants).

  • Only patients from the western Rhineland-Palatinate and central and south-west Schleswig-Holstein showed high percentages of CAS (41–57%), while other regions showed CAS percentages of between 4% and 33%.

  • This study shows the level of heterogeneity in both the indication and the choice of method to treat carotid stenosis in Germany. Although natural variation and differences in actual incidence are likely, they are unable to fully account for the variation observed. This raises the question of the extent to which the German S3 guidelines on carotid stenosis have been adopted, as well as the social question of nationwide, demand-oriented treatment for everyone.

  • Particularly in view of the German government’s quality offensive and the IQTiG’s efforts to improve eQA this article, which captures the quintessential, practice-relevant points, raises the question as to whether non-consideration of regional aspects and separating eQS analyses according to CEA and CAS represent a relevant weakness in the quality assurance procedure. Therefore, the inclusion of regional parameters and relevant risk factors in the evaluation of eQA data should be examined in order to assess outcome quality from the perspective of the individual service providers, and not separately from the quality of patient selection, indication, and choice of procedure from the point of view of the regional patient population.