Impact Factor 3.552

The 2nd most cited open-access journal in Clinical Neurology

This article is part of the Research Topic

Towards an Understanding of Tinnitus Heterogeneity

Original Research ARTICLE

Front. Neurol., 03 April 2017 | https://doi.org/10.3389/fneur.2017.00115

Cluster Analysis to Identify Possible Subgroups in Tinnitus Patients

  • 1Department of Otorhinolaryngology/Head and Neck Surgery, University of Groningen, University Medical Center Groningen, Groningen, Netherlands
  • 2Graduate School of Medical Sciences (Research School of Behavioural and Cognitive Neurosciences), University of Groningen, Groningen, Netherlands
  • 3Department of Neurosurgery, University of Groningen, University Medical Center Groningen, Groningen, Netherlands

Introduction: In tinnitus treatment, there is a tendency to shift from a “one size fits all” to a more individual, patient-tailored approach. Insight in the heterogeneity of the tinnitus spectrum might improve the management of tinnitus patients in terms of choice of treatment and identification of patients with severe mental distress. The goal of this study was to identify subgroups in a large group of tinnitus patients.

Methods: Data were collected from patients with severe tinnitus complaints visiting our tertiary referral tinnitus care group at the University Medical Center Groningen. Patient-reported and physician-reported variables were collected during their visit to our clinic. Cluster analyses were used to characterize subgroups. For the selection of the right variables to enter in the cluster analysis, two approaches were used: (1) variable reduction with principle component analysis and (2) variable selection based on expert opinion.

Results: Various variables of 1,783 tinnitus patients were included in the analyses. Cluster analysis (1) included 976 patients and resulted in a four-cluster solution. The effect of external influences was the most discriminative between the groups, or clusters, of patients. The “silhouette measure” of the cluster outcome was low (0.2), indicating a “no substantial” cluster structure. Cluster analysis (2) included 761 patients and resulted in a three-cluster solution, comparable to the first analysis. Again, a “no substantial” cluster structure was found (0.2).

Conclusion: Two cluster analyses on a large database of tinnitus patients revealed that clusters of patients are mostly formed by a different response of external influences on their disease. However, both cluster outcomes based on this dataset showed a poor stability, suggesting that our tinnitus population comprises a continuum rather than a number of clearly defined subgroups.

Introduction

Tinnitus is a prevalent condition, estimated to affect 5–18% of the adult population (1), which may lead to severe impairment in quality of life. Although many trials on tinnitus therapies have been conducted, hardly ever a treatment effect is demonstrated. A potential explanation for the lack of effectivity of these treatments might be the underlying heterogeneity of the disease. Therefore, consensus on the optimal treatment of tinnitus gradually shifts from a “one size fits all” approach to a more patient-tailored approach. Possibly, a particular group of patients would be more likely to respond to treatment, if a selection is made on etiology, tinnitus characteristics, or patient characteristics. It might be that in a specific subgroup of tinnitus patients a particular treatment is successful, while this treatment is not successful in another subgroup of tinnitus patients. Thus, insight in the heterogeneity of the tinnitus spectrum might improve the management of these patients.

Identification of tinnitus subgroups is also important with regard to concomitant mental distress. Hoekstra et al. demonstrated that patients who express certain characteristics (i.e., high percentage of experience of tinnitus during the day, self-reported depression or anxiety, and subjective experience of tinnitus loudness) are more at risk for a high tinnitus burden (2). This subgroup of patients with high tinnitus distress needs more extensive counseling and follow-up in order to prevent mental breakdown.

In an attempt to identify subgroups of tinnitus patients, cluster analysis was used in this study. Cluster analysis is a statistical technique that divides data into groups, or clusters, which are meaningful and/or useful. It is an explorative analysis that assigns patients to clusters based on certain characteristics, so that patients look very much alike within a cluster (high within-group homogeneity) and, at the same time, are very different from the other clusters (low between-group homogeneity) (3). In research, this cluster analysis method is not only used in medicine studies to identify groups of patients but also in marketing for finding customer segments for example.

In 2008, Tyler et al. performed a preliminary cluster analysis on 153 patients with tinnitus (4). The cluster analysis of Tyler et al. identified distinct cluster characteristics, which were described as: (1) “constant distressing tinnitus,” (2) “varying tinnitus that is worse in noise,” (3) “tinnitus patients who are copers and whose tinnitus is not influenced by somatic modulation,” and (4) “tinnitus patients who are copers but whose tinnitus is worse in quiet environments.” Tyler et al. did not report a statistic value to identify the degree to which patients are clustered in these groups.

In this paper, we report on an exploratory cluster analysis of patients from the tinnitus database of the University Medical Center Groningen (n = 1,783 patients). We initially attempted to replicate the cluster analysis reported by Tyler et al (4); however, this was not possible as many of the variables used in their analysis were not identical or not available in our database. Instead, we report on two further cluster analyses. In the first analysis, the choice of variables that were entered in the cluster analyses was fully guided by the statistical techniques. In the second analysis, the selection of variables was based on the expert opinions in our tinnitus clinic. The aim of this study was to identify subgroups of tinnitus using cluster analysis, based on a very large dataset of tinnitus patients.

Materials and Methods

Tinnitus Population

This study was performed at the Otorhinolaryngology Department of the University Medical Center Groningen (The Netherlands), which has a specialized multidisciplinary care group for tinnitus patients since 2007. Patients with severe complaints of tinnitus can be referred to this care group for medical consultation and psychological support. Almost all patients who visit this care group have consulted an audiologist and/or otorhinolaryngologist earlier. However, these patients were referred to our specialized tertiary care group by these specialists, because of the severity and impact of the complaints. Consultation at our clinic consists of thorough evaluation by an otorhinolaryngologist, an audiologist, radiologist, a medical social worker, and/or a psychologist.

Variables

The variables that were available for this cluster analysis were demographic characteristics (e.g., sex and age), tinnitus characteristics (e.g., duration of tinnitus, onset, lateralization, pitch, variable loudness), factors of influence on patients tinnitus (e.g., influence of loud sounds, noisy environment, movement of head and neck), tinnitus and quality of life-related questionnaires [e.g., tinnitus handicap index (THI), visual analog scale (VAS), and hospital anxiety and depression scale (HADS)], and audiological characteristics [e.g., frequency matching, pure tone audiometry (PTAs), loudness matching of tinnitus]. Hearing loss was divided into categories based on the pure tone audiogram: (1) no or slight hearing loss (both ears thresholds <30 dB on PTA thresholds at 0.25–0.5–1–2–4–8 kHz), (2) asymmetrical hearing loss (≥30 dB difference between both ears on the mean PTA thresholds at 2–4–8 kHz), (3) bilateral high tone hearing loss (both ears thresholds ≥30 dB on PTA thresholds at 2–4–8 kHz), (4) bilateral severe hearing loss (PTA thresholds >30 dB on 0.25–0.5–1–2–4–8 kHz), and (5) others. The available variables are all listed in Table S1 in Supplementary Material. All patient-reported variables were completed by the patients in booklets during the visit at the tinnitus outpatient clinic. Physician reported data, such as audiological characteristics, were also reported in booklets by the physician. All these routinely collected data were anonymized and entered in a database. For the current analysis, these data were retrospectively analyzed. The collection of data was approved by the Institutional Reviewer Board of the UMCG. No full review was needed due to the retrospective nature of this study.

Selection of Variables for Cluster Analysis

All variables that were collected were entered in the database. However, not all of these variables could be entered in the cluster analysis. In cluster analysis, it is important to keep the sample size in mind when deciding how many variables to enter in the analysis. Formann recommends the number of variables (m) of 2m = sample size (5). In our study, the sample size is n = 1783, implying that the number of variables should be 10 or 11. There are two ways to select appropriate variables for cluster analysis: (1) a statistical approach with the use of principal component analysis (PCA) and (2) selection of variables based on “expert opinion,” i.e., variables that are presumed to be clinically relevant and thought to be discriminative in the total group. Both selection procedures were performed in this study, resulting in two different cluster analyses.

Variable Reduction by PCA

A PCA is a dimension reduction technique that condenses variables that are highly correlated into a set of factors, thereby removing overlap and redundancy. PCA with Varimax rotation was performed on all variables with missing values ≤20%. The PCA revealed several factors, and for each factor, the variable with the highest loading was selected for inclusion in the cluster analysis.

Variable Reduction by “Expert Panel”

After excluding variables with missing values >20%, variables were selected by a group of tinnitus care professionals and investigators [Minke J. C. van den Berge (otolaryngology resident, Ph.D. candidate in tinnitus research), Pim van Dijk (medical physicist, audiologist, involved in the tinnitus care group), and Emile de Kleine (medical physicist, audiologist, involved in the tinnitus care group)]. Based on clinical experience and knowledge, those variables were selected that were deemed important in discriminating subgroups of tinnitus.

Cluster Analysis

The “two-step” cluster analysis method was used as the analyses that contained both categorical and continuous variables (6). Continuous variables were standardized by default. For distance measures, the log-likelihood method was used, as both continuous and categorical variables were entered in the analysis. The number of clusters to be formed was not specified in advance. The “silhouette measure of cohesion and separation” is a measure for the overall goodness-of-fit of the cluster structure that was found. It ranges from −1 to 1 (<0.25: no substantial structure; 0.26–0.50: weak structure and could be artificial; 0.51–0.70: reasonable structure; 0.71–1.0: strong structure) (7).

Differences in characteristics between clusters were compared according to the cluster membership variable, using one-way ANOVA for continuous variables and Pearson Chi-square tests for categorical variables. SPSS version 23.0 (Chicago, IL, USA) was used for all tests. The significance level was set at α = 0.05, and all tests were two tailed.

Results

Subject Characteristics

For this study, data from 1,783 consecutive patients who visited the UMCG tinnitus clinic between July 2007 and June 2016 were collected. The baseline characteristics of this study population are shown in Table 1. Variables that had >20% missing values are not shown in this table. In this population, 39.3% were females and the mean age was 53.6 ± 13.5 years. Tinnitus was unilateral in 50.7% of the cases and bilateral or central in 48.2%. The mean THI in the total patient group was 42.5 ± 23.2.

TABLE 1
www.frontiersin.org

Table 1. Demographic and tinnitus-related characteristic of included tinnitus patients (n = 1,783).

Outcome of Cluster Analysis with Variables Selected by PCA

The PCA was performed to obtain eigen values for each factor. The Kaiser–Meyer–Olkin measure of sampling adequacy was 0.681 and the Bartlett’s test of sphericity was significant [χ2 (6,127) = 325, p < 0.001], both indicating an appropriate factor model. A total of eight factors were extracted (based on the eigenvalue >1 rule), which together explained 55% of the total variance. Variables with the highest loading on each factor were selected. Subsequently, these variables (n = 8) were entered in the cluster analysis. The clustering revealed a four-cluster solution. As the analysis excludes every case when there is any variable with a missing value (listwise exclusion), the analysis was based on n = 976 patients. The cluster outcome showed a “silhouette measure of cohesion and separation” of 0.20, indicating that it is a “no substantial” cluster solution (7). Characteristics of these four identified clusters are shown in Table 2. The variables in the table are ranked from most discriminative between groups (top of the table) to less discriminative (bottom of the table). All variables differed statistically significant between the four clusters, except for the variables “VAS tinnitus annoyance” and “frequency of the tinnitus” (p = 0.925 and p = 0.478, respectively).

TABLE 2
www.frontiersin.org

Table 2. Characteristics of the four clusters identified by clustering with variable selection based on principal component analysis.

Cluster 1 (n = 293) is characterized by the fact that tinnitus is not easily influenced: loud sounds, sleep deprivation, and nap in the afternoon have no effect on their tinnitus. These patients have a relatively high difference between hearing loss in the right and left ear. These patients have relatively low HADS-depression scores.

Cluster 2 (n = 259) is distinguished by a gradual onset of the tinnitus. Also in this group, tinnitus is easily negatively influenced, especially by loud sounds and sleep deprivation. Both make their tinnitus louder.

Cluster 3 (n = 197) is a group of patients who report that their tinnitus is less loud when they hear loud sounds. Sleep deprivation and a nap in the afternoon mostly have no effect on their tinnitus.

Cluster 4 (n = 227) is typically a group with tinnitus of acute onset. They report that their tinnitus is easily negatively influenced by loud sounds or sleep deprivation. They show relatively high HADS-depression scores.

Outcome of Cluster Analysis with Variables Selected by Expert Panel

For the alternative method of choosing variables for clustering, 11 variables were selected by a panel of experts in the field. The selected variables (see Table 3) were entered in the cluster analysis. The outcome was a three-cluster solution, with a “silhouette measure of cohesion and separation” of 0.20, again indicating a poor solution. Because of listwise exclusion as described earlier, this analysis was based on n = 761 patients. About 527 of these patients were also included in the first cluster analysis. Also in this table, variables are ranked according to their degree of discriminative value. All variables differed significantly between the clusters (all p-Values <0.001).

TABLE 3
www.frontiersin.org

Table 3. Characteristics of the four clusters identified by clustering with variables selected by expert opinion.

Cluster 1 (n = 287) is a group of patients whose tinnitus is not easily influenced: loud sounds, stress, or movement of head and neck have no effect on their tinnitus loudness. Patients prefer a noisy environment. Sounds are never to seldom experienced as uncomfortably loud. The tinnitus is mostly unilateral. Although most patients in this group have no or slight hearing loss, other types of hearing loss are present in this group as well. They are not very much bothered or depressed by their tinnitus, as the THI and HADS-depression scores are low.

Cluster 2 (n = 247) is a predominantly male group, whose tinnitus gets worse by stress, loud sounds, and movement of head and neck. These patients prefer to be in a noisy environment. Sometimes, sounds are experienced as uncomfortably loud. Most of the patients have no or slight hearing loss. Tinnitus is bilateral, and the loudness of the tinnitus is variable.

Cluster 3 (n = 227) is characterized by the fact that their tinnitus is easily negatively influenced: loud sounds and stress clearly make their tinnitus louder. These patients prefer a silent environment. Often, patients find sounds uncomfortably loud. Tinnitus is often bilateral with most patients having no or slight hearing loss or asymmetrical hearing loss. The loudness of the tinnitus is variable.

Discussion

Key Results

In this study, we performed cluster analysis with the aim to identify subgroups in a population of tinnitus patients. Variable selection for cluster analysis was performed in two ways: by a strict methodological approach based on PCA, and by expert opinion. These analyses identified four- and three-patient clusters, where the clusters showed clearly different characteristics. However, the clustering solution in both analyses was not substantial, as indicated by a poor cluster solution quality.

Although both cluster analyses gave different outcomes, there were also interesting similarities. In both cluster solutions, the effects of “stress” and “loud sounds” on tinnitus have a relatively high discriminative value between groups. In each analysis, a group was revealed in which patients report that their tinnitus gets louder from loud sounds, and there was a group that reported that their tinnitus got less loud. In an earlier cluster analysis by Tyler et al., it is described that their found clusters differed by the effect of external factors on patients’ tinnitus: some patients are easily and negatively influenced by external factors and, in others, this has no effect (4). On the contrary, Tyler et al. describe a group that is characterized by high scores on tinnitus questionnaires and the HADS depression and anxiety scale. However, this was not reflected in our cluster solutions.

In the cluster analysis based on variables selected by experts, there was a clear distinction between a group that preferred a silent environment for their tinnitus and another group that had a preference for a noisy environment. The fact that some patients with tinnitus prefer noise and others prefer silence has been described earlier (8). This is interesting, as one might speculate that the latter group may have a higher change of responding well to sound therapy than the other group.

When interpreting these results, it must be kept in mind that the “silhouette measure” of both analyses was only 0.2. This is lower than the critical boundary of 0.25, which implies that there was no substantial clustering in this patient cohort. A lack of clustering indicates that the transition from one cluster to another is relatively smooth, without clear-cut boundaries. As a comparison, consider a group of cities, where the coordinates of the cities would go into a cluster analysis. If one group of cities is clearly separated from another group by a stretch of open land, the silhouette value will be large (when viewed from a distance, the cities will have a distinct silhouette of their skyline). However, if there is no such open land between the clusters, the silhouette value is low, consistent with the absence of substantial clustering. In our patient cohort, there were clearly no distinct “open stretches of land” between the clusters, suggesting that patient form a continuum rather than a clear clustering. As discussed above, the cluster analysis of Tyler et al. identified clusters with characteristics that show some resemblance to the clusters reported here. Unfortunately, Tyler et al. do not report a silhouette value or other measure of clustering. Hence, at present, it is not possible to discuss the clustering strength in their cohort.

Cluster analysis has been upcoming in medical research. Recently, an interesting cluster analysis on bilateral Meniere disease was published to define clinical subgroups with potential similar etiologies. In this study, five clinical variants of bilateral Meniere disease were found based on six clinical variables and with a high silhouette measure of 0.8 (9). This study is not only beneficial to improve the selection of patients but also can explain the negative treatment effects of several treatment trials, as results can be biased by a heterogeneous patient group based on etiology (9). The difficulty in cluster analysis is that it is a type of analysis that is very sensitive to change of variables. The selection of variables is critical for the outcome of the cluster analysis (6). Generally, highly correlated variables should be avoided and it is important to select variables that can make a clear-cut differentiation between clusters (6). The systematic statistical approach of selecting variables using the highest factor loading on extracted factors by PCA is often used and has the advantage of choosing variables in a reproducible, transparent way. A downside of this technique is that the factor solution only explains a certain amount of variance and, therefore, much information is discarded. Eliminating factors with low loadings on the extracted factors has the same effect (10). This may lead to a reduced success of a subsequent cluster analysis. On the other hand, a disadvantage of selecting variables based on clinical knowledge or “gut feeling” is that it is less transparent. Also, unrecognized highly discriminating variables may remain undiscovered.

Strengths of the Study

For this study, a very large database of tinnitus patients was used with almost 1,800 patients. Even after exclusion of patients with missing values, n = 976 and n = 761 could be included in the cluster analyses. We expect that, if clear clustering would have existed with these variables, we would have been able to find it in these groups. There was an overlap of 527 patients who were included in both the first cluster outcome and the second cluster outcome. This is a substantial overlap, pointing out that it does not seem likely that the differences between both cluster analyses are caused by the differences in included patients.

Limitations

This study explored the patient cohort of a tertiary tinnitus referral center. Thus, the population described here consists of a group of tinnitus patients who were persistent in their search for treatments for their tinnitus. Our patient cohort may, therefore, be biased with a certain type of tinnitus patients. Potentially, a study including also less-persistent help seeking and non-help-seeking subjects would have identified a clearer clustering.

Although we had access to a large database of tinnitus patients (n = 1,783), over the years, there were changes in the variables that were collected because of changes made to the diagnostic protocol. Since the cluster analyses required a complete set of data for each patient, not all 1,783 patients could be included in the analysis, but 976 and 761.

Furthermore, it is debated internationally whether tinnitus is a disease or a symptom. One can look at it in both ways: when tinnitus is a result of an acoustic neuroma, then tinnitus can be a symptom. However, if we look at tinnitus as the result of defect on a cellular level of the auditory cortex, then tinnitus can be regarded as a disease. In most patients visiting our clinic, the etiology of the tinnitus is unclear. The fact is that these patients included in our dataset experience bothersome tinnitus. Within this group, we aimed to find subgroups such as patients with continuous central, loud tinnitus tend to have a high score on THI and VAS and find that their tinnitus gets worse in noisy environments. If we are able to find such patterns, may be we can adjust our treatment strategy to that (in this example, hearing aids might not be successful). Although the raised issue about tinnitus being a symptom or disease is important, we believe that this analysis looking for clusters of patients based on tinnitus characteristics transcends this issue.

Finally, the low silhouette value indicates that this patient cohort represents a heterogeneous group without clear clustering. Obviously, any cluster analysis outcome highly depends on the variables that we entered into the clustering algorithm. Our patient data consisted of mainly audiometry and questionnaire metrics. In these cluster analyses, tinnitus patients appear to represent a continuum rather than clearly defined subgroups, based on a low silhouette measure. However, it is possibly that other metrics (e.g., fMRI/EEG, genetic evaluation) are able to identify tinnitus subgroups. In other words, the lack of clustering in our analyses does not imply that clusters do not exist. However, if clusters exist, they cannot be identified with the variables that were considered here.

Conclusion

Two cluster analyses of a large patient cohort identified three and four groups of tinnitus patients, respectively. The clustering was not substantial, as a low silhouette measure of the cluster solutions was found. This indicates that in this particular cohort, tinnitus patients appear to represent rather a continuum than clearly defined subgroups. This finding may have consequences for future treatments: if clear subgroups would have been present, clearly distinct treatment might be developed in the future. However, for a continuum of patients, it may be necessary to use a number of treatments to find the optimum for each individual patient. Obviously, our conclusion is based on the set of variables that were at our disposal. Possibly, new future ways to characterize tinnitus patients may be able to find distinct subgroup in tinnitus patients.

Author Contributions

Conception of study and design; writing of manuscript: MB and PD. Intellectual contribution to analyses: MB, PD, and EK. Critical revision of the manuscript: RF, RA, PD, and JD. Data acquisition: MB, PD, RA, RH, and EK.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer CM and handling Editor declared their shared affiliation, and the handling Editor states that the process nevertheless met the standards of a fair and objective review.

Funding

There were no funding sources for this study.

Supplementary Material

The Supplementary Material for this article can be found online at https://www.frontiersin.org/article/10.3389/fneur.2017.00115/full#supplementary-material.

References

1. Savage J, Waddell A. Tinnitus. BMJ Clin Evid (2014) 2014:0506.

PubMed Abstract | Google Scholar

2. Hoekstra CE, Wesdorp FM, van Zanten GA. Socio-demographic, health, and tinnitus related variables affecting tinnitus severity. Ear Hear (2014) 35(5):544–54. doi: 10.1097/AUD.0000000000000045

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Clatworthy J, Buick D, Hankins M, Weinman J, Horne R. The use and reporting of cluster analysis in health psychology: a review. Br J Health Psychol (2005) 10(Pt 3):329–58. doi:10.1348/135910705X25697

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Tyler R, Coelho C, Tao P, Ji H, Noble W, Gehringer A, et al. Identifying tinnitus subgroups with cluster analysis. Am J Audiol (2008) 17(2):S176–84. doi:10.1044/1059-0889(2008/07-0044)

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Formann AK. Die Latent-Class-Analyse: Einführung in Theorie und Anwendung. Weinheim, Germany: Beltz (1984). Available from: http://catalog.hathitrust.org/api/volumes/oclc/16851238.html

Google Scholar

6. Mooi E, Sarstedt M. Cluster analysis. A Concise Guide to Market Research. Berlin, Heidelberg: Springer-Verlag (2011). p. 237–84.

Google Scholar

7. Kaufman L, Rousseeuw PJ, editors. Finding Groups in Data: An Introduction to Cluster Analysis. Hoboken, NJ: John Wiley & Sons, Inc. (1990).

Google Scholar

8. Stouffer JL, Tyler RS. Characterization of tinnitus by tinnitus patients. J Speech Hear Disord (1990) 55(3):439–53. doi:10.1044/jshd.5503.439

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Frejo L, Soto-Varela A, Santos-Perez S, Aran I, Batuecas-Caletrio A, Perez-Guillen V, et al. Clinical subgroups in bilateral Meniere disease. Front Neurol (2016) 7:182. doi:10.3389/fneur.2016.00182

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Dolnicar S, Grunn B. Challenging “factor-cluster-segmentation”. J Travel Res (2008) 47:63–71. doi:10.1177/0047287508318910

CrossRef Full Text | Google Scholar

Keywords: tinnitus, cluster analysis, subgroup identification, heterogeneity of tinnitus, principle component analysis

Citation: van den Berge MJC, Free RH, Arnold R, de Kleine E, Hofman R, van Dijk JMC and van Dijk P (2017) Cluster Analysis to Identify Possible Subgroups in Tinnitus Patients. Front. Neurol. 8:115. doi: 10.3389/fneur.2017.00115

Received: 29 November 2016; Accepted: 14 March 2017;
Published: 03 April 2017

Edited by:

Jose Antonio Lopez-Escamez, Granada University Hospital, Spain

Reviewed by:

Eduardo Martin-Sanz, Hospital de Getafe, Spain
Andrés Soto-Varela, Complejo Hospitalario Universitario de Santiago, Spain
Coral Del Val Muñoz, University of Granada, Spain

Copyright: © 2017 van den Berge, Free, Arnold, de Kleine, Hofman, van Dijk and van Dijk. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Minke J. C. van den Berge, m.j.c.van.den.berge@umcg.nl