Introduction

Musculoskeletal disorders accounted for 30 percent of the approximately 1.2 million workplace illnesses and injuries reported in 2006 in the United States [1]. A general population study in The Netherlands found that those with elbow or wrist/hand symptoms were over two times as likely to be disabled from work than those without such symptoms [2]. Similar results held for those with neck, shoulder or upper back complaints.

Although prevalent and potentially disabling, some disorders of the upper extremity are not readily classified. Both the Harrington et al. [3] and Sluiter et al. [4] criteria documents for clinical diagnoses or epidemiologic case definitions of upper extremity musculoskeletal disorders (UEMSDs) recognize a non-specific disorder. The diagnosis or characterization of a non-specific UEMSD, variously called “non-specific diffuse forearm pain”, “repetitive strain injury” (RSI), “cumulative trauma disorder” (CTD), etc. may be one of exclusion of more localized disorders (e.g., carpal tunnel syndrome, de Quervain’s disease). However, this non-specific disorder may in fact represent several separate disorders. This may have implications for prevention and treatment.

Subtypes in diseases including complex regional pain syndrome [5, 6], fibromyalgia [7], and Alzheimer’s disease [810] have been explored through factor analysis and cluster analysis. A similar analysis of patients with UEMSDs was possible through a questionnaire administered in 1999 through the Dutch RSI Patients Association (n = 806) from which detailed information on symptom location, quality, and intensity was obtained. The objective of this study was to determine if there were patterns of symptoms in these patients (who identified themselves as having RSI) that tended to occur together, with a long-term goal of discovering new upper extremity syndromes.

Specific UEMSDs include those that by definition are localized to particular anatomical regions. For instance, carpal tunnel syndrome as a compression of the median nerve in the carpal tunnel results in symptoms in the hand/wrist, symptoms in epicondylitis are confined to the elbow, and rotator cuff tendinitis is characterized by complaints in the shoulder region. Similarly, specific UEMSDs may be distinguished by quality of symptoms experienced. That is, the predominant symptom in the various localized tendinopathies is activity dependent pain, whereas paraesthesias and numbness are not likely in the absence of neuropathology. These distinguishing features of specific UEMSDs informed the study hypotheses, which were explored through factor analysis and subsequent cluster analysis of the patient association database. It was hypothesized that:

  1. 1.

    distinct UEMSD syndromes would be identified based on localization of symptoms affecting the following upper extremity regions:

    1. a.

      hand/wrist

    2. b.

      elbow

    3. c.

      shoulder

    4. d.

      neck

  2. 2.

    further discrimination of syndromes would be possible within the above distinct anatomical groups through classifying by present qualities of symptoms (pain, stiffness, paraesthesia, etc.) and by intensity of patient complaint.

Additionally, clusters were examined for their differences with respect to gender and age proportion, percentage of those working, and overall course of symptoms.

Methods

Description of Dataset

In the Netherlands, patients with the same disease are allowed and administratively helped by Governmental Services to join and start a society. Everyone who is diagnosed with or believes that they have RSI is free to become a member of the RSI-patient society. In February 1999, all registered 1,700 members of the Dutch National RSI-patient society were sent a questionnaire. All members received this survey with a response envelope at their home address as sent from the secretariat of the society. An accompanying letter was included from the researchers and the RSI-patient society board. Additionally, a general reminder was posted in the society’s official newsletter. A total of 53% responded. Five questionnaires were returned unopened, and 92 questionnaires were filled in insufficiently for analyses. Therefore, data from 806 questionnaires could be included in the analyses.

This study was conducted in accordance with the laws of the country in which it was performed. In The Netherlands, it is not necessary for study subjects to fill out a consent form in de-identified cross-sectional survey research inquiring about health status. If subjects do not want to participate, they just will not fill in the questionnaire.

The questionnaire included demographic variables such as gender, age, height and weight. Work status was indicated as “full time”, “part time”, “working with adaptations”, “work pension” or “not working”. Other items included number of months since symptoms began, and course of complaints. Subjects were asked to assess the level of current complaints as compared with the level at symptom onset. Respondents could choose “no change”, “getting better”, “getting worse”, “getting better but more extensive”, and “getting worse and more extensive”. Symptom location, quality and intensity were also ascertained.

Description of Candidate Variables for Factor Analysis

Questionnaire respondents characterized their current or recent symptoms through 118 variables; these formed the candidates for factoring. These variables were of three types, distinguished by what attributes they were measuring and the scale by which the attributes were measured. Attributes included:

  1. 1.

    current symptom quality intensity (1–10 scale, 2 = “just sensible”, 10 = “maximal pain”) was ascertained separately for eleven body regions (neck, right(R)/left(L) shoulder, R/L upper arm, R/L elbow, R/L lower arm, and R/L wrist/hand). Possible qualities included: dull aching, stabbing/sharp, burning, toothache-like, or shooting pain, stiffness, paraesthesias, numbness, and coldness.

  2. 2.

    duration of symptoms in the last week by body region. Possible regions included: neck, shoulder, upper arm, elbow, lower arm, and hand/wrist. The left or right side of each region as well as dorsal or ventral (front or back) area could be specified. Possible duration categories included: 0 days, 1–2 days, 3–6 days and every day (4 point scale).

  3. 3.

    severity of symptoms in the last week. Questions included: severity of pain in the neck/shoulder/arm/hand, severity of pain in the neck/shoulder/arm/hand during a specific actvivity, severity of paraesthesias in the shoulder/arm/hand, severity of weakness in the neck/shoulder/arm/hand, and severity of stiffness in the neck/shoulder/arm/hand. Possible answers included: none, mild, average, severe, and extreme (5 point scale).

Factoring Procedure

The original 118 items were reduced using principal components analysis (PCA) as a first step in describing and identifying the clusters. Of the n = 806 individuals with usable questionnaires, twenty-nine subjects who reported diabetes, thyroid problems or rheumatoid arthritis were eliminated from the analysis to avoid confounding by these co-morbid conditions. This yielded a total of 777 subjects for factor analysis. Of these, 701 had complete data on the 118 items. One subject was later deleted at random during the clustering procedures so that the sample could be randomly split into five subsamples for replication analyses [11]. As a result, the final n = 700.

One of the more critical decisions in a PCA is to determine the correct number of factors to retain and rotate [1214]. The most common rule is to retain factors when eigenvalues are ≥1.0. This solitary criterion is the default procedure in most statistical packages. The shortcoming is that implementation of solitary criteria tends to under- or overestimate the number of true latent dimensions [1517]. Accordingly, each model was evaluated against the following four rules: (a) scree [18], (b) Glorfeld’s (1995) extension of parallel analysis (PA) [19, 20], (c) interpretability [12, 15], and (d) internal-consistency reliabilities required to be greater, or equal to, 0.70 [21, 22]. Although all four rules were applied, special emphasis was placed on PA and the scree test because results from several investigations demonstrated that PA is the best method for determining the correct number of factors to accept and that the scree test is a useful adjunct [16, 17, 19, 23].

Clustering Procedure

Upon identifying the underlying principal components, profiles (also called profile types or clusters) from the 700 participants were sorted according to level, shape, and dispersion among the factors. All factors were converted to z-scores (mean = 0.0, standard deviation = 1.0) prior to analysis in order to equalize their variances. Otherwise, scores with larger variances would contribute disproportionately to the formation of clusters [24]. The analyses required that individuals within each group be maximally similar to one another (maximum homogeneity) and maximally dissimilar to those in other groups (minimum overlap). The groups of similar profiles were also required to be reasonably replicable across groups rather than spurious mergers, as would occur by chance [25]. To this end, the sample was split into four random subsamples, each composed of 150 participants.

Cluster analysis was used to sort the 700 profiles [26, 27]. After evaluating numerous clustering algorithms, Ward’s [28] minimum-variance procedure was determined to best satisfy the research goals for first-stage clustering. The procedure began with Ward’s technique to produce first-stage clusters for each of the four independent subsamples. Three percent of the participants were trimmed from each of the four, first-stage samples, and thus were eliminated as outliers from the clustering process. Clustering during the second-stage analysis was completed again using Ward’s technique. Group centroids from the second-stage solution served as starting points for the third-stage, iterative-partitioning analysis conducted using K-means passes.

Results from the final solution were compared along three internal criteria. First, the overall solution was required to show a replicability rate ≥ 75% for the absorption of first-stage clusters into the final, third-stage taxonomy. Because cluster solutions may be unstable, replication can validate their integrity [29]. More specifically, a Monte Carlo simulation demonstrated that analyses are likely to identify the true number of clusters when replication rates are ≥75% [11]. Second, the typology was required to yield an average within-profile type homogeneity coefficient, \( \bar{H} \) ≥ 0.60 [30]. This property has been referred to as internal cohesion [31]. It requires small score dispersions within clusters and, thereby, for individuals within a cluster to be maximally similar. The third criterion required the typology to show external isolation [31]. This property dictates large score dispersions between clusters so that the clusters are maximally dissimilar from one another. Cattell’s r p coefficient [32] was calculated between profile types as our measure of external isolation and the typology needed to provide an average between-profile-types similarity coefficient, \( \bar{r}_{p} \) ≤ 0.40 [32]. \( \bar{H} \) and r p are each sensitive to similar profile levels and shapes, where a value of 1.00 indicates profiles identical in level and shape, 0.0 indicates chance similarity based on the full sample, and negative values indicate gross dissimilarity. Minimum criterion values have been established empirically for \( \bar{H} \) and r p (≥0.60 and ≤0.40, respectively) based on clustering and classification studies conducted with random samples larger than the current one, as well as with epidemiologically-representative data sets [3336].

Statistical Analysis

The clusters were compared with regard to gender, age, percentage working, and course of complaints. Subjects were regarded as working if they denoted their work status as full time, part time, or working with adaptations. Those on a work pension were regarded as not working. Symptom course was compared across three groups: (1) no change, (2) better, and (3) worse, or worse and more extensive, or better but more extensive. The Pearson’s chi-square test was used to separately compare the proportion of those working, percentage female, and distribution of course of symptoms within each cluster to that of the cohort overall. A one-way ANOVA followed by the Games-Howell post-hoc test (used because of unequal n in the clusters) determined differences among clusters in number of months since symptom onset and age. A P-value of ≤0.05 was regarded as statistically significant.

Clusters were regarded as mild if all factors were ≤0.3 standard deviations (SD) from the cohort mean, and moderate if all factors were <1.5 SD from the mean. If any factor in the cluster was ≥1.5 SD from the mean, the cluster was labeled severe. If any factor in the cluster was ≥5 SD from the mean, the cluster was labeled extreme. The interpretation of these levels was empirically based.

Results

Demographic Characteristics of Clustered Subjects

The analysis cohort consisted of n = 700 individuals of whom n = 462 (66%) were females. Mean age was 37.6 (SD = 8.6) years. Eighty-four percent had been diagnosed by a physician. Just under half (n = 321, 46%) listed their current occupation as “administrative activities/office work (secretary, accountant)” or “journalist/translator/text writer”. Information technology (IT) workers and those in industry/performing heavy physical work each comprised approximately 15% of respondents. Sixty percent (n = 417) were working at least part-time.

Morbidity Characteristics of Clustered Subjects

The 700 subjects who were clustered had mean values of symptom intensity just over “just sensible” of these types of symptoms: dull aching pain and stiffness in the neck and right shoulder, dull aching right lower arm pain, and right wrist symptoms including dull aching pain, stiffness, paraesthesias, and a cold feeling (Table 1). Additionally, the average participant experienced pain for at least 1–2 days during the previous week in each body region under question. Symptoms that were experienced for the greatest duration during the last week (between 1–2 and 3–6 days) included those on the right-side dorsal portion of the shoulder, those in the right dorsum of wrists, hands and fingers and those in the neck. Severity of pain during the last week was “average” in the neck/shoulder/arm/hand (mean = 3.0, SD = 1.0) and worse than average in the same region during specific activities (mean = 3.6, SD = 1.2) in the last week. Mean severity of paraesthesias, weakness, and stiffness during the last week was between “mild” and “average”.

Table 1 Mean intensity by present symptom quality in n = 700 members of the Dutch RSI Patients Association

The average duration of symptoms was 39 months (SD = 37.4). Somewhat over one-third of the respondents characterized the course of their symptoms when compared with symptom onset as “worse and more extensive”, and a slightly smaller number reported that their symptoms were improving. Approximately one-fifth stated that their complaints were “less but more extensive”.

Factor Results

Fourteen factors were interpreted according to the magnitude and meaning of their salient pattern coefficients. All coefficients greater than 0.45 were considered appreciable. Names for the 14 dimensions are provided below. Corresponding internal-consistency reliabilities (i.e., alpha (α) coefficients) are presented in parentheses: (1) left-sided shoulder through hand/wrist symptoms, chronic in lower arm (α = 0.92), (2) bilateral proximal (elbow, shoulder, and neck) cold and numbness (α = 0.91), (3) bilateral diffuse (defined as “neck-to-hand/wrist”) “like a toothache” pain (α = 0.93), (4) chronic bilateral stiff aching neck/shoulder (α = 0.84), (5) right-sided diffuse stabbing shooting pain (α = 0.89), (6) chronic (dorsal and ventral) wrist/hand/finger symptoms (α = 0.87), (7) bilateral diffuse burning pain (α = 0.89), (8) chronic bilateral distal (lower arm through hand/wrist) numbness and paraesthesias (α = 0.86), (9) chronic bilateral (dorsal and ventral) upper arm symptoms (α = 0.83), (10) bilateral arm (NOT neck and shoulder) stiffness (α = 0.85), (11) bilateral proximal (neck-to-elbow) paraesthesias (α = 0.88), (12) chronic bilateral shooting, stabbing, aching elbow pain (α = 0.86), (13) chronic bilateral lower arm aching pain (α = 0.73), (14) right-sided upper-to-lower arm aching pain (α = 0.79). Associations among the 14 retained components were low. Only one of the possible 91 non-redundant correlations was ≥0.40 (i.e., 14 (11 − 1)/2 = 91). This latter finding indicates that the 14 dimensions are essentially orthogonal and share less than 16% percent of their variance in common with one another (i.e., 0.402 = 0.16).

Cluster Results

First-stage clustering produced 41 clusters in total (an average of 8.2 profile groups per analysis). These were submitted to second-stage agglomerative clustering based on a 41 × 41 similarity matrix and the solution at all hierarchical steps was evaluated against the stated internal criteria. The second-stage, eight-cluster solution was the only one to satisfy all three of the internal criteria. Therefore, it was submitted to a third-stage, iterative-partitioning analysis which also resulted in a final typology of eight clusters.

Table 2 displays, for each of the cluster types, its estimated prevalence in the population, average coefficient for within-type homogeneity, between-types similarity, and replication rate. The average \( \bar{H} \) value of 0.99 is well above its a priori criterion of ≥0.60 and the average r p value of −0.095 satisfies the ≤0.40 criterion. The types replicated 90.6% of the time across the four independent experiments. Therefore, the obtained, eight-cluster solution also satisfied its a priori replication criterion of ≥75%.

Table 2 Cluster prevalences and statistical characteristics (n = 700)

In Table 3, each cluster is presented in terms of number of standard deviations (z-score) from the mean on a per factor basis. It is important to remember that the mean (where z-score = 0) symptomatic profile of the group under study consists of dull/aching neck and right arm (shoulder through wrist), stiff neck, shoulder, and wrist, and paraesthesias and coldness in the wrist (Table 1). Also, subjects were symptomatic in the last week, particularly in the right neck, shoulder and hand/fingers, with severity of pain reported as “average”.

Table 3 Mean score patterns for cluster types by factors, expressed as z-scores

Cluster 1 consists of most factors between 0.5 and 1.0 standard deviations below the mean value for the cohort, while Cluster 8 is characterized by all but one factor being one or more standard deviations above the mean (Table 3). Again, a z-score of zero does not denote no symptoms, but rather the mean value of symptoms in the cohort.

Cluster Membership (Tables 4, 5)

Cluster 1: Mild Diffuse MSD (Prevalence = 37.9%)

This group was characterized by mild symptoms in the upper extremity with no one anatomical area predominating. A significantly greater proportion of subjects in this group were working (75.5%) than in the cohort overall. Their symptom course was significantly different from that of the entire population, with almost half (47%) reporting fewer symptoms from symptom onset through the date of filling out the questionnaire. They had symptoms for a mean of 2.5 years; this was a shorter duration than most of the other clusters (see below). Also, significantly more of these cluster members were male (44.9%) than in the other clusters.

Table 4 Cluster description, symptom intensity and localization
Table 5 Demographics, percentage working, mean symptom duration and symptom course, by cluster (n = 700)

Cluster 2: Mild Bilateral Aching Stiff Neck and Shoulders (Prevalence = 25.4%)

These cluster members had localized bilateral aching and stiffness in the neck and shoulder. The mean duration of symptoms was significantly greater than that of cluster 1 (40.6 vs. 30.3 months). Other characteristics were similar to that of the overall population.

Cluster 3: Moderate Bilateral Arm and Wrist/Hand/Finger MSD (Prevalence = 15.5%)

These individuals had chronic bilateral upper arm and hand/wrist/finger symptoms. Significantly fewer were working (46%) in comparison to the cohort as a whole. Seventy-four percent characterized their symptoms as getting worse or more widespread.

Cluster 4: Severe Diffuse MSD (Prevalence = 4.6%)

These cluster members had bilateral widespread symptoms of burning pain combined with distal numbness. Bilateral elbow symptoms and aching stiff arms were also prominent. Duration of symptoms (60.7 months) was longer than in any other cluster, and significantly different in duration than cluster 1. Almost sixty percent of these subjects reported their symptom course as “worse and more extensive.” Three-quarters of the subjects were female.

Cluster 5: Moderate MSD with Highlighted Bilateral Elbow Shooting Stabbing Pain (Prevalence = 11.0%)

This group experienced chronic bilateral elbow symptoms, with diffuse right-sided arm pain. Fewer subjects were working in comparison to the cohort as a whole (P ≤ 0.05), although 44% of this group was working. Most (83%) reported worsening or spreading symptoms.

Cluster 6: Severe Diffuse MSD with Highlighted Coldness, Numbness, and Paraesthesias (CNP) (Prevalence = 2.9%)

Chronic overall coldness, numbness, and paraesthesias characterized this cluster with lesser intensities of pain and stiffness. Although not statistically significant due to low power, almost all (85%) of the members were female. Slightly fewer than half of the cluster members were not working.

Cluster 7: Extreme Diffuse MSD with Highlighted Pain (Prevalence = 1.9%)

This group was characterized by high intensity chronic bilateral symptoms throughout the upper extremity. Symptom qualities included aching, burning, “like a toothache”, and stabbing/sharp pain as well as stiffness and numbness. Although there were only 13 members in this cluster, several of the examined characteristic proportions were significantly different from that of the overall patient population. Only one of the 13 cluster members was working (P ≤ 0.05) and all were females (P ≤ 0.05). Their symptom course was significantly different from that of the entire cohort with almost one-fourth reporting no change in their symptom course. Non-significantly different from the mean symptom duration of the cohort as a whole, their symptom duration was almost 5 years.

Cluster 8: Extreme Diffuse MSD with Highlighted CNP (Prevalence 0.9%)

These individuals experienced a great intensity of chronic bilateral symptoms. Coldness, numbness, and paraesthesias were predominant. Stiffness, which was worse distally than proximally, and overall pain were also experienced. Although the sample size was small, interestingly, 50% of these individuals were working and they had the shortest mean symptom duration (27.8 months) of any cluster. Five of the six individuals were female and for most, symptoms were getting worse, with none reporting symptom improvement.

Discussion

Eight distinct clusters comprised of UEMSD symptoms were formulated in a dataset of RSI patients. Most clusters consisted of diffuse rather than localized symptoms, in contrast to our initial hypothesis. The exception was mild bilateral aching stiff neck and shoulders (cluster 2). By and large the clusters were distinguished by severity and quality, rather than locality, of symptoms.

On average, the study patients experienced a symptom intensity of slightly more than “just sensible” dull aching pain in the neck, and right shoulder, lower arm and wrist, stiffness in the neck and right shoulder and wrist, and wrist paraesthesias and coldness. Cluster severity ranged from mild through extreme as compared with these overall mean symptom intensities. Most of the RSI patients experienced mild symptoms; cluster sizes were smaller as symptom severity increased.

In the largest cluster, consisting of almost 40% of subjects (cluster 1), members had mild symptoms of relatively short duration (a mean of 2.5 years). At the time of the survey, three-fourths in this cluster were working and almost half were improving. Another 25% of the patient association members had localized mild bilateral neck and shoulder symptoms consistent with tension neck syndrome [37] and possibly trapezius myalgia [38] (cluster 2). Clusters 3 and 5 were characterized by symptoms of moderate intensity (total prevalence = 26.5%).

Based on symptom qualities, there appeared to be two types of non-specific more severe MSDs (total prevalence = 10.1%). The first was a type which was dominated by pain (cluster 7—extreme diffuse MSD with highlighted pain). The second was a type dominated by coldness, numbness and paraesthesias (cluster 6—severe diffuse MSD with highlighted CNP, cluster 8—extreme diffuse MSD with highlighted CNP). Interestingly, among those with the most extreme symptoms, clusters 7 and 8, the latter group is characterized by relatively short (2.3 years) symptom duration with half of the cluster members working. It is possible that the predominance of coldness, numbness, and paraesthesias rather than pain has allowed many of these patients to continue their employment.

The diffuse nature of the clusters we identified indicates the possible incompleteness of conventional case definitions for work-related upper extremity morbidity, which tend to be localized to particular anatomical areas. Mackinnon and Novak [39] describe a non-specific upper extremity pain syndrome upon which symptoms specific to a particular cumulative trauma disorder (such as CTS) are superimposed, of which cluster 5 appears to be an example. Other physicians have noted diffuse symptoms in office workers, an occupational group which accounts for nearly half of the present sample [40, 41]. Bilateral pain, widespread pain, and pain in contiguous anatomical sites greater than to be expected given an assumption of independence between sites have been found in the general population [2, 42] and in workers [43]. Moreover, practicing clinicians have observed that UEMSD symptoms seem to become more diffuse with time [44, 45]. Various theories have been proposed to explain such widespread symptoms, including a positive feedback mechanism of spreading muscle overuse/underuse [39] and peripheral or central sensitization [4648]. Supporting evidence has been found in an animal model [49].

This point prevalence study lends support to prior clinical reports describing anatomically diffuse symptoms among patients with UEMSDs. In our study, diffuse symptoms are found in a majority of patients at a mean duration of 3.25 years from symptom onset. In light of these findings, development of patient assessment tools in both clinical and research settings should routinely incorporate information on symptom pattern and duration. Such information can help to reformulate clinical case definitions, diagnostic categories, prognostic indicators and, ultimately, help to improve clinical management, including occupational rehabilitation.

To our knowledge, this study is unique in examining symptom qualities in UEMSD patients in such detail. We are aware of only one other author who has conducted cluster analyses in UEMSD patients. Reading et al. [50, 51] performed three separate cluster analyses based predominantly on physical examination signs in the wrist/hand, shoulder and elbow in working age adults with and without upper extremity pain. In the hand/wrist as well as in the elbow region, cluster analyses were based primarily on severity and location, whereas the shoulder clusters were formed on the basis of severity only. It is unknown how many of the patients in their sample had involvement in more than one anatomical region. As our data relied on symptoms only, it is difficult to make a comparison with these findings.

The formulated clusters were not identical with respect to gender composition. As UEMSD severity increased, so generally did the proportion of females in the associated cluster. Musculoskeletal symptom severity (due to more prevalent exposure to work and non-work risk factors) was found to be greater in females than in males in a study of pubic service employees [52]. However, since females may tend to report symptoms in general more easily than males [53], further research would be necessary to confirm an association between MSD symptom severity and gender.

Although small numbers in several clusters prevented adequate comparisons, several clusters differed both by percentage working and by symptom course as compared with the surveyed patient population as a whole. Of these, clusters 5 (moderate MSD with highlighted bilateral elbow shooting stabbing pain) and 7 (extreme diffuse MSD with highlighted pain) were characterized by worsening or stable symptoms and by fewer individuals working. As mentioned previously, cluster 1 (mild diffuse MSD) was composed of more members working with improving symptoms. It should be noted that each of these clusters were characterized by diffuse symptoms.

Study Limitations

There are several limitations to the study. Brauer et al. [54] found that subjects were able to recall their pain severity over a period of 3 months. Hence, recall of symptoms within the last week would presumably be reliable. Nonetheless, recall bias may have affected respondents’ recollection of symptom intensity at the time of onset, depending on whether they were better or worse at the time of the administration of the survey. Current symptom status affects individuals’ recall of prior pain; those currently symptomatic may have the tendency to over-recall previous pain at a particular point in time whereas those without symptoms may “forget” or underestimate previous pain at a given occasion [55]. However, it is unknown how this would affect subjects’ reporting of their overall symptom course.

The low response rate (53%) may indicate a selection bias affecting the results. It is likely that the non-respondents were more impaired than those answering the survey. Although the whole range in intensity or severity of complaints was represented in the final sample, it is likely that this survey may under-represent severe MSDs or that it may underestimate the prevalence of MSDs that impair writing ability (such as MSDs particularly affecting the fingers/hand/wrist), as these subjects may have been prevented by their symptoms from filling out the questionnaire.

The findings may be specific to the occupational exposures of this particular cohort, which was composed primarily of office workers, with lesser percentages of IT specialists and industrial employees. Clusters based on symptoms alone may be more sensitive to disease states than those based on both signs and symptoms and useful for secondary preventive purposes, yet they may not be specific enough to determine suitable primary prevention in the workplace or treatment options in the clinic.

A strength of the research presented in this paper (that was undertaken to define an empirical classification system for upper extremity MSD symptoms) was the employment of cluster analysis, a classification technique for forming homogenous groups within complex data sets [56]. Cluster analysis aims at the classification of individuals or objects based solely on an internal analysis of similarities and differences in multivariate data patterns [57]. The methodology has found numerous applications in psychology [26, 33, 58]. This approach has also been applied in characterizing Alzheimer’s disease [8, 9], systematic lupus erythematosus [59], complex regional pain syndrome [5, 6], Parkinson’s disease [60] and other physical disorders.

Notwithstanding such wide use, methodological ambiguities exist in the cluster analysis literature with respect to the choice of similarity measure, recovery capabilities, and replicability of structures. [57]. Our clustering strategy employed the multistage Euclidean grouping (MEG) method described by McDermott [25] that employed Ward’s similarity measure and explicitly evaluates its recovery capability and structure replication rate.

Conclusions

Eight clusters based primarily on symptom quality and severity were formed from a survey of upper extremity MSD patients. The clusters, most of which consisted of diffuse symptomatology, differed with respect to percentage working, gender proportion and symptom course. Future studies should formulate clusters based on detailed symptoms as well as signs in comparable populations. Since the clusters were determined through reports of present symptoms and of symptoms within 1 week of filling out the questionnaire, it is not possible to predict a probable symptom course. Future studies could collect data from health care providers prospectively to determine the possible prognostic value of these clusters with respect to natural history, chronicity, and return to work.