Skip to main content
  • Research article
  • Open access
  • Published:

Runners with patellofemoral pain demonstrate sub-groups of pelvic acceleration profiles using hierarchical cluster analysis: an exploratory cross-sectional study



Previous studies have suggested that distinct and homogenous sub-groups of gait patterns exist among runners with patellofemoral pain (PFP), based on gait analysis. However, acquisition of 3D kinematic data using optical systems is time consuming and prone to marker placement errors. In contrast, axial segment acceleration data can represent an overall running pattern, being easy to acquire and not influenced by marker placement error. Therefore, the purpose of this study was to determine if pelvic acceleration patterns during running could be used to classify PFP patients into homogeneous sub-groups. A secondary purpose was to analyze lower limb kinematic data to investigate the practical implications of clustering these subjects based on 3D pelvic acceleration data.


A hierarchical cluster analysis was used to determine sub-groups of similar running profiles among 110 PFP subjects, separately for males (n = 44) and females (n = 66), using pelvic acceleration data (reduced with principal component analysis) during treadmill running acquired with optical motion capture system. In a secondary analysis, peak joint angles were compared between clusters (α = 0.05) to provide clinical context and deeper understanding of variables that separated clusters.


The results reveal two distinct running gait sub-groups (C1 and C2) for female subjects and no sub-groups were identified for males. Two pelvic acceleration components were different between clusters (PC1 and PC5; p < 0.001). While females in C1 presented similar acceleration patterns to males, C2 presented greater vertical and anterior peak accelerations. All females presented higher and delayed mediolateral acceleration peaks than males. Males presented greater ankle eversion (p < 0.001), lower knee abduction (p = 0.007) and hip adduction (p = 0.002) than all females, and lower hip internal rotation than C1 (p = 0.007).


Two distinct and homogeneous kinematic PFP sub-groups were identified for female subjects, but not for males. The results suggest that differences in running gait patterns between clusters occur mainly due to sex-related factors, but there are subtle differences among female subjects. This study shows the potential use of pelvic acceleration patterns, which can be acquired with accessible wearable technology (i.e. accelerometers).

Peer Review reports


Patellofemoral pain (PFP) is the most common musculoskeletal overuse injury in runners, regardless of sex and age [1] and it has been suggested that atypical gait kinematics may play a role in its etiology [2,3,4]. However, a general consensus on the pathomechanics of this injury has yet to be reached [4] possibly due to the existence of more than a single atypical gait pattern [5,6,7,8].

Distinct running kinematic sub-groups have been identified in PFP patients, with a sub-group presenting lower peak hip adduction; another with greater peak knee abduction angles [6]; and a sub-group that presented an attempt to compensate for a greater initial hip internal rotation with an external rotation during midstance [7], suggesting the existence of multiple kinematic pathomechanical pathways or motor adaptations associated with PFP. It also has been shown that kinematic differences are influenced by sex-related factors, wherein males with PFP present lower angles of hip adduction and knee abduction during running [8]. These findings should be interpreted with caution, as they are from relatively small sample sizes (n = 16–22) and were based on visual inspection of the data, approaches which may not fully describe the etiology of PFP and related sub-groups.

The consensus statement from the 3rd International Patellofemoral Pain Research [9] concluded that “identification of sub-groups remains the ‘holy grail’ for PFP research”. Identification of sub-groups could provide insight into the pathomechanics associated with PFP as well as inform personalized treatment. One approach to identify homogenous sub-groups within a dataset is the use of cluster analyses. With the advance of technology and data science methods the use of machine learning techniques in gait analysis is growing and exploratory analysis of complex data such as gait kinematics is important to bring new insights in the field [10, 11]. Recent research from our laboratory [12] utilized a hierarchical cluster analysis (HCA) approach to successfully identify two distinct and homogeneous kinematic sub-groups among 121 healthy runners. However, because the acquisition of 3-dimensional (3D) kinematics data is time consuming, it usually relies on multiple assessors to collect data on larger sample sizes, introducing sources of imprecision into the data collection process, especially from marker placement errors [13,14,15]. Therefore, finding alternative methods for evaluating gait mechanics becomes important for clinical applications.

Recently, axial segment acceleration data has provided unique insight into running mechanics, discriminating between fatigue states [16] and training levels [17]. Therefore, the analysis of pelvic acceleration profiles could also be useful to identify sub-groups of runners with PFP, with the advantage of being less influenced by identification of anatomical landmarks when using optical motion capture systems. Furthermore, the study of segmental accelerations has the option to use wearable devices for data acquisition, which are becoming increasingly popular in both academia and industry, and there has been an effort to further investigate their potential applications in health systems [18,19,20]. Therefore, this approach may offer an accessible and objective method of assessment with clinical applicability.

The purpose of this exploratory study was to determine if running gait patterns in PFP runners could be clustered into homogeneous sub-groups using pelvic acceleration data, using a large dataset of males and females with PFP. Based on the results from previous studies, we hypothesized that more than one running gait pattern sub-group, or cluster, would be present in female PFP runners, since the studies suggesting the existence of sub-groups were mostly comprised of women [6, 7]. Furthermore, female runners with PFP would be different from their male counterparts, given that sex-related kinematic differences have been identified previously [8]. A secondary purpose was to analyze kinematic differences between the sub-groups, by comparing lower limb peak angles that are considered important in the pathomechanics of PFP, and thereby investigate the practical and clinical implications of clustering these subjects based on 3D pelvic acceleration data. Based on the kinematic sub-groups that has been described in the literature [6, 7], we expected female clusters to present differences in hip and knee frontal and transverse plane angles, and males to display lower peak angles of hip adduction and knee abduction [8].



Data from 110 physically active individuals with PFP with running as their primary exercise modality for at least 6 months, were analyzed in this cross-sectional study. The presence of PFP was confirmed by a licensed healthcare professional (i.e., athletic therapist, physical therapist or medical doctor) based on specific inclusion and exclusion criteria (Table 1). Subjects experiencing pain in other sites were also included in the study, however the primary complaint had to be PFP. Data was collected either at the University of Calgary or in clinical settings partnered with the Running Injury Clinic.

Table 1 Inclusion and exclusion criteriaa

Data collection

The data collection methods are described in detail elsewhere [21, 22]. Briefly, 8 high-speed digital video cameras (MX3/Nexus, Vicon, Oxford, UK) were used to film treadmill-running at 200 Hz. Spherical retro-reflective markers (9 mm diameter, Mocap Solutions, Huntington Beach, USA) were attached to the specific lower extremity anatomical landmarks bilaterally along with technical marker clusters on rigid shells placed to represent the pelvis and bilateral foot, shank, and thigh segments. Each participant wore the same shoes (Pegasus, Nike, Beaverton, USA) to standardize the footwear condition.

Following placement of all the anatomical and segment markers, each participant stood on a motorized treadmill (Bertec Corporation, Columbus, OH, USA) for a 1-s static trial. Upon completion of the static trial, the markers on the anatomical landmarks were removed while the technical marker clusters remained. The participants were instructed to warm-up on the treadmill for 2–3 min, and then ran on the treadmill at a comfortable self-selected pace (2.61 ± 0.20 m/s) for 20 s, in which approximately 60–80 consecutive running steps were collected for processing and analysis. All participants were experienced treadmill users and were permitted as much time as they required to familiarize themselves with treadmill running before beginning the data collection.

Data processing

Ankle, knee and hip joint sagittal plane angular accelerations were used for defining ground contact, using previously published event detection methods [23]. The position of the pelvis was measured using the centroid of the pelvic marker cluster [24] and pelvic acceleration was calculated by double differentiation of pelvis displacement using a modified Savitzky-Golay method [25]. Differentiation was performed at both stages using a time-window of 10 data points, and 4th order polynomial fitting. In order to emulate a wearable device, marker accelerations in the global coordinate frame were then converted to a local coordinate frame on the pelvis, using segment markers and rigid body transformations [26]. The local coordinate frame was aligned with the global frame during the static trial.

Each step cycle was normalized to 100 points, with 80 data points for stance and 20 data points for flight phase, since we are analyzing an axial segment. These normalized phases were then combined to represent 100% of the step cycle, averaged over all extracted steps, and standardized to zero mean and unit variance. The kinematic data (3 planes of motions × 100 time-normalized pelvic accelerations) were combined into one 300-dimensional row vector for each subject, creating a matrix of 110 subjects-by-300 data points.

Data analysis

The HCA method was used to identify homogeneous running gait patterns separately for males and females based on the pelvic acceleration time-series, by creating a cluster tree, or dendrogram for each sex-group. Agglomerative strategy or a “bottom up” approach was used, which consists of three steps: (1) a measure of dissimilarity between sets of subjects using the Euclidean distance, (2) subject linkage using the Ward’s minimum variance method [27], and (3) cluster determination using the variance ratio criterion [28].

Following identification of homogeneous clusters (sub-groups) of PFP runners, differences in demographics, injury characteristics, vertical displacement of the pelvic centroid and peak joint angles were examined using one-way ANOVA (Tukey test for post-hoc analyses) and chi-squared test (α = 0.05), and effect sizes were calculated based on η2 and Cramer’s V indices, respectively. In case the data did not present a normal distribution (Shapiro-Wilk test) or a homogeneous variance between sub-groups (Levene test), the Kruskal-Wallis test was performed (Dunn’s test for post-hoc analyses). Differences in pelvic acceleration patterns were examined after applying a principal component analysis (PCA) to the standardized data matrix, and they were identified based on the interpretation of principal components (PCs) that presented a large effect size (η2 > 0.14) [12, 29], which were used to reconstruct the acceleration waveforms for a better mechanical interpretation [30]. The squared coefficients of correlations between the PC scores and the raw acceleration data (squared loading) [31] were used to calculate the relative loading of the PCs in the vertical (VT), antero-posterior (AP) and medio-lateral (ML) directions to aid in the interpretation of the PCs [32].

We also selected joint angles that are considered important in PFP pathomechanics and that have been suggested to differ between PFP sub-groups [6, 33], to compare between sub-groups. The analyzed peak joint angles were: ankle eversion; knee flexion, knee abduction and knee external rotation; and hip adduction and internal rotation.

A Pearson’s correlation coefficient was calculated between the significant PCs and demographic, injury characterization and kinematic variables that presented differences between sub-groups to determine whether these latter factors were significantly correlated with the acceleration patterns. All data processing and statistical analysis were performed on MATLAB 9.1 (The MathWorks Inc., Natick, MA,USA).


Identification of PFP sub-groups

For the female subjects, the variance ratio criterion determined the optimal number of clusters to be two sub-groups (C1 and C2) (Fig. 1a), whereas for the male subjects, no sub-groups could be identified (Fig. 1b) in the HCA.

Fig. 1
figure 1

Dendrogram of the hierarchical cluster analysis. Clustering of PFP patients produced by the hierarchical cluster analysis. a Male subjects; (b) Female subjects

Subject clinical and demographic characteristics for each cluster are presented in Table 2. There was a significant difference in height (H = 49.9; df = 2; p < 0.001; η2 = 0.32), with male subjects being taller than both female clusters (p < 0.001 for males vs C1, and males vs C2); and also in mass (H = 50.2; df = 2; p < 0.001; η2 = 0.32), with male subjects being heavier than both female clusters (p < 0.001 for males vs C1, and males vs C2). There were no significant differences between clusters with respect to age (F = 2.4; df = 2; p = 0.095; η2 = 0.04), running speed (H = 1.3; df = 2; p = 0.521; η2 = 0.01), years running (H = 0.9; df = 2; p = 0.629; η2 = 0.01); ratio of unilateral to bilateral involvement (χ2 = 1.7; p = 0.434; V = 0.12) or subjects with multiple injury sites (χ2 = 0.5; p = 0.781; V = 0.07).

Table 2 Number of PFP participants and subject specifications (Mean and (SD)) for the determined clusters

Table 3 presents descriptive statistics for each cluster regarding pelvic acceleration components, vertical displacement and lower limb kinematics. There were significant differences with large effect sizes between sub-groups only for the following principal components of pelvic acceleration: PC1 (F = 39.7; df = 2; p < 0.001; η2 = 0.43), with the distinct group being the females in C2 (p < 0.001 for males vs C2, and C1 vs C2); and PC5 (F = 19.8; df = 2; p < 0.001; η2 = 0.27), wherein C1 was the sub-group with significant difference (p < 0.001 for males vs C1, and C2 vs C1). However, these PCs explained less than 25% of the variance in the dataset (17.0% and 7.7%, respectively). Height presented a significant correlation with PC1 (p = 0.016), but not with PC5 (p = 0.064), and the correlation coefficients were weak (r < 0.30). Although, body mass was different between subgroups, it had no significant correlation (p > 0.05) with either of the selected PCs.

Table 3 Mean and standard deviation of PC scores, vertical displacement and peak joint angles

Differences in running kinematics between sub-groups

There was a significant difference in peak ankle eversion (H = 15.1; df = 2; p < 0.001; η2 = 0.12), wherein male PFP subjects presented greater angles than C1 (p = 0.003) and C2 (p = 0.004), and this joint angle had a low but significant correlation with PC1 (p < 0.036, r = − 0.20), but not with PC5. Peak knee abduction was also significantly different between clusters (H = 12.3; df = 2; p = 0.002; η2 = 0.09), with males exhibiting lower angles when compared to C1 (p = 0.019) and C2 (p = 0.005), and this joint angle was also only correlated with PC1 (p < 0.001, r = 0.38). There were also differences in hip adduction (F = 6.5; df = 2; p = 0.002; η2 = 0.11), with lower angles for male subjects in comparison to C1 (p = 0.043) and C2 (p = 0.003). The same tendency occurred for hip internal rotation (F = 5.2; df = 2; p = 0.007; η2 = 0.09), however the difference was only significant for males compared to women in C1 (p = 0.006). While peak hip adduction was correlated with PC1 (p = 0.005, r = − 0.27), hip internal rotation demonstrated a correlation with PC 5 (p = 0.014, r = − 0.23). There were no significant differences for knee flexion (H = 2.2; df = 2; p = 0.331; η2 < 0.01) and external rotation (F = 1.4; df = 2; p = 0.251; η2 = 0.03). Vertical displacement of the pelvis presented differences between sub-groups (F = 11.1; df = 2; p < 0.001; η2 = 0.17), wherein females from C1 displayed lower magnitudes of displacement when compared to males (p < 0.001) and females in C2 (p = 0.011); and this variable was only correlated with PC5 (p = 0.026, r = 0.21).

PC1 presented a high relative loading in the VT direction (47.9%), representing variations in the peak acceleration and the magnitude at early stance phase (Fig. 2a). There was a lower relative loading of PC1 in the AP direction (28.9%), wherein it represented a phase shift of the posterior acceleration peak in early stance (Fig. 2b). In the ML direction, PC1 also represented phase shifts in ML peak accelerations towards the stance and swing limbs during the first half of stance phase (Fig. 2c), but it was the lowest relative loading (23.2%).

Fig. 2
figure 2

Time-normalized pelvic accelerations. a Vertical acceleration, (b) Anteroposterior acceleration, and (c) Mediolateral acceleration for males and female sub-groups C1 and C2 during stance phase (1%–80%) and flight phase (81%–100%; gray area) of running. Regions represented by the significant principal components are indicated in the graphs

PC5 also had relatively high loadings in the VT axis (43.0%), denoting a difference in the rate of magnitude decrease after the peak acceleration (Fig. 2a), although these differences are subtle. The AP relative loading was the lowest for PC5 (25.0%) and it indicated a magnitude difference in the forward acceleration after weight acceptance (Fig. 2b). In the ML direction, there was a low relative loading (32.1%), representing a difference in magnitude variation during the first half of stance phase (Fig. 2c).

Overall, when comparing the pelvic acceleration patterns, males had similar acceleration patterns to females in C1 in the VT and AP directions, but the latter presented higher and delayed peaks ML accelerations. Females in C2 displayed lower acceleration magnitudes in early stance and a higher peak acceleration in the VT direction; a greater forward peak in early stance; and delayed peak accelerations in the AP and ML directions.


PFP sub-groups based on pelvic acceleration

The first purpose of the present study was to determine if running gait patterns in individuals experiencing PFP at the time of testing could be clustered into homogeneous sub-groups based on combinations of pelvic acceleration components. In support of our hypothesis, two distinct and homogenous sub-groups (clusters) were present in females with PFP, and these clusters were different when compared to PFP males. These results are similar to previous studies that also reported two to three different running patterns based on visual inspection of 3D kinematic data [6, 7] and mechanical differences between males and females with PFP [8].

There were no significant differences in running speed between sub-groups, which is a factor that has been shown to affect axial segment acceleration [34], especially in the ML axis [35]. Male subjects were significantly taller and heavier than females and these anthropometric differences are known to influence 3D kinematics during running [36]. However, there was a very weak correlation for height, and no correlation for body mass with the acceleration PCs that presented differences between sub-groups suggesting that the relationship with those factors was minimal.

The advantage of investigating pelvic acceleration as a measure of running mechanics is that it is less influenced by marker placement errors and is a much simpler method than a full 3D gait assessment, as it depends only on the trajectory of a single pelvic marker cluster. Additionally, these factors allow for the use of data from multiple research centres, allowing for the application of ‘big data’ analytics and a better understanding of the interaction between biomechanical factors and musculoskeletal injuries [10, 11]. Furthermore, the results of the present study opens the possibility for the use of wearable devices for data acquisition, such as a single triaxial accelerometer on the pelvis, an approach which is becoming increasingly popular in industry and health care [18, 20]. Therefore, the current work identifying sub-groups of PFP patients is a novel finding that can guide future studies in providing better context that can hopefully improve clinical practice.

Identification of differences in running gait patterns between sub-groups

A secondary purpose was to analyze peak joint angles between clusters to better understand the practical and clinical implications of clustering subjects with PFP based on 3D pelvic acceleration data. In general, differences in joint kinematics were sex-related, since there were no significant differences between female clusters, except for peak hip internal rotation. Moreover, the magnitude of mean differences were within the threshold for detectable kinematic changes reported by Osis et al. [15] for knee abduction (3.4o) and hip internal rotation (5.6o). However, the differences in ankle eversion and hip adduction between males and females are greater than the error margins caused by marker placement errors, confirming the findings of Willy et al. [8] who reported males with PFP to have less hip adduction than their female counterparts.

Phinyomark et al. [12] reported the existence of two different sub-groups of asymptomatic runners based on a HCA of lower limb joint kinematics, and when they compared the peak knee abduction angles of those clusters with a sample of subjects with PFP, group differences were dependent on the cluster of healthy individuals that was used as reference. Interestingly, all PFP sub-groups from the current study presented greater values of knee abduction when compared to the ones reported for healthy runners (healthy C1: 8.0o; healthy C2: 4.4o). However, there is a tendency for a progressively greater alteration in knee frontal plane angles when comparing males to females in C1 and C2, although there was no significant difference between the female clusters. This could be related with distinct pathomechanical pathways or differences in response to treatment. For example, in a previous work [37] we found that non-responders to exercise treatment protocol presented greater knee abduction angles during late stance and swing phases of running gait, and the current findings suggest that this could be identified by pelvic acceleration data.

To our knowledge, this is the first study to investigate pelvic acceleration profiles in runners with PFP, and the identification of sub-groups could generate insights about differences in pathomechanics or adaptations to pain. Additionally, the analysis of segmental acceleration profiles minimizes measurement imprecisions originating from marker placement errors that propagate into the calculation of joint angles in 3D kinematics [14, 15]. Furthermore, the results of the current study suggest that accelerations acquired using wearable devices [24] may utilise this method in a clinical setting as an evidence-informed method to improve patient care and rehabilitation decisions.

The pelvic acceleration data can provide some clinical insight that can help clinicians make decisions regarding treatment options. For example, peak resultant pelvic acceleration is related to center of mass acceleration during 10 to 75% of stance phase [38]. Therefore, pelvic accelerations can provide some insights on shock absorption and lower limb stiffness. Nevertheless, this connection must be made with caution, since accelerations based on segmental measures overestimate the behavior of center of mass [38]. Women in C2 presented a higher VT peak acceleration, suggesting a diminished capacity for shock absorption. Since no differences in peak knee flexion angles were detected, this could be an indication of greater leg stiffness in these subjects, which is partially supported by the findings that women present higher leg stiffness during running [39] and drop jump landing tasks [40] when compared to males. In contrast, females in C1 were similar to males regarding VT acceleration patterns, which could be explained by the lower VT displacement.

Women also presented higher and delayed peak accelerations in the ML direction, suggesting differences in the control of side-to-side body movement during the first half of the stance phase, when these oscillations occur. This pattern could be related to the larger hip adduction angles exhibited during running, which led to increases in ML accelerations. In addition, females in C2 displayed a delay in peak AP accelerations in early stance, causing a prolonged period of deceleration. It is possible that this finding is related to strength differences between males and females [41, 42], as stronger individuals may be able to exert shorter impulses to achieve the same net change in momentum, however, strength differences were not quantified in the current study.

Although the identification of sub-groups among the female subjects with PFP did not coincide with significant differences in peak lower limb joint angles, there seems to be a progression of values in knee abduction and hip internal rotation depending on the cluster of female subjects. Specifically, there is a tendency for C1 to have lower knee abduction and higher hip internal rotation than C2. These factors could be related to symptom severity or differences in response to treatment, but would need further investigation.


In addition to the differences in height and weight between males and females that were already discussed, other limitations to the current research study are acknowledged. First, this study included both subjects with uni- or bilateral involvement and with secondary pain symptoms besides PFP, which could have also modified running mechanics. However, there was no significant difference in the distribution of those variables between the two subgroups, leading us to believe that it was not an important factor for this clustering. Additionally, these types of patients are frequently seen in clinical practice, therefore these PFP patients are important to include in research studies.

Second, we did not have access to other clinical variables that could influence running mechanics and explain the differences that were found between sub-groups. For example, Selfe et al. [43] has described 3 clusters of PFP patients that were grouped based on clinical measures of strength, flexibility and joint alignment and mobility. Additionally, experimental pain induction in the knee joint has been shown to cause reductions in peak torque in maximal voluntary contraction of knee flexors and extensors [44] and increased sway displacement during quiet stance [45], indicating that pain level could be a driver of changes in motor control. Therefore, future studies should include the aforementioned clinical variables to investigate whether they are related to the differences in running pattern found between sub-groups to have a better understanding in a clinical context.

Finally, this investigation used an HCA approach, which is an unsupervised machine learning technique suitable for exploratory analyses, to determine whether this type of data could be useful in the identification of subgroups within a cohort of runners with PFP. Overall, our hypothesis was supported by the findings and suggest that a supervised analysis could also be applied to identify specific subgroups with specific clinical relevance. For example, recent work from our laboratory used a supervised machine learning method to classify runners with PFP into responders or non-responders to exercise treatment based on running kinematic data, achieving 78% of classification accuracy [37]. Thus, a similar approach could be applied in this context, using pelvic acceleration data to develop an objective method for the identification of such subgroups with greater accessibility in a clinical setting. Regardless, the present study is an important first step to verify the utility of simple measures, like pelvic accelerations, for the objective assessment of gait biomechanics.


In conclusion, using a hierarchical cluster analysis, the present study is the first to identify distinct pelvic acceleration patterns during running gait in a large group of PFP runners. Two homogenous female sub-groups were identified based on pelvic accelerations with one sub-group demonstrating a delay in the posterior and mediolateral acceleration peaks compared to the other. However, both female sub-groups presented greater acceleration peaks than males in all directions. Further analysis of peak kinematic angles provided clinical context to these sub-groups and revealed that gender-differences hip internal rotation, an important factor related to PFP, is distinct among the female sub-group. These results suggest that the variability observed in running gait patterns for PFP runners occur mainly due to sex-related factors, but there are subtle differences among females that could influence the interpretation of kinematic data. The findings also highlight potential for the use of data acquired with accessible wearable technology in the identification of sub-groups in PFP patients. Future research can use this approach in order to classify PFP patients and develop targeted intervention and injury prevention strategies.







Hierarchical cluster analysis




Principal component analysis


Principal components


Patellofemoral pain




  1. Taunton JE. A retrospective case-control analysis of 2002 running injuries. Br J Sports Med. 2002;36:95–101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Ferber R, Kendall KD, Farr L. Changes in knee biomechanics after a hip-abductor strengthening protocol for runners with patellofemoral pain syndrome. J Athl Train. 2011;46:142–9.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Ferber R, Bolgla LA, Earl-Boehm JE, Emery C, Hamstra-Wright K. Strengthening of the hip and core versus knee muscles for the treatment of patellofemoral pain: a multicenter randomized controlled trial. J Athl Train. 2015;50:366–77.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Barton CJ, Levinger P, Menz HB, Webster KE. Kinematic gait characteristics associated with patellofemoral pain syndrome: a systematic review. Gait Posture. 2009;30:405–16.

    Article  PubMed  Google Scholar 

  5. Dierks TA, Manal KT, Hamill J, Davis IS. Proximal and distal influences on hip and knee kinematics in runners with Patellofemoral pain during a prolonged run. J Orthop Sport Phys Ther. 2008;38:448–56.

    Article  Google Scholar 

  6. Dierks TA, Manal KT, Hamill J, Davis I. Lower extremity kinematics in runners with patellofemoral pain during a prolonged run. Med Sci Sports Exerc. 2011;43:693–700.

    Article  PubMed  Google Scholar 

  7. Noehren B, Pohl MB, Sanchez Z, Cunningham T, Lattermann C. Proximal and distal kinematics in female runners with patellofemoral pain. Clin Biomech. 2012;27:366–71.

    Article  Google Scholar 

  8. Willy RW, Manal KT, Witvrouw EE, Davis IS. Are mechanics different between male and female runners with patellofemoral pain? Med Sci Sports Exerc. 2012;44:2165–71.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Witvrouw E, Callaghan MJ, Stefanik JJ, Noehren B, Bazett-Jones DM, Willson JD, et al. Patellofemoral pain: consensus statement from the 3rd international patellofemoral pain research retreat held in Vancouver, September 2013. Br J Sports Med. 2014;48:411–4.

    Article  PubMed  Google Scholar 

  10. Ferber R, Osis ST, Hicks JL, Delp SL. Gait biomechanics in the era of data science. J Biomech. 2016;49:3759–61.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Phinyomark, A, Petri, G, Ibáñez-Marcelo, E, Osis ST, Ferber R. Analysis of Big Data in Running Biomechanics: Current Trends and Future Directions. Journal of Medical and Biological Engineering – Special Issue: Recent Advances in Biomedical Engineering. 2018;38:244–260.

  12. Phinyomark A, Osis S, Hettinga BA, Ferber R. Kinematic gait patterns in healthy runners: a hierarchical cluster analysis. J Biomech. 2015;48:3897–904.

    Article  PubMed  Google Scholar 

  13. Della Croce U, Cappozzo A, Kerrigan DC. Pelvis and lower limb anatomical landmark calibration precision and its propagation to bone geometry and joint angles. Med Biol Eng Comput. 1999;37:155–61.

    Article  CAS  PubMed  Google Scholar 

  14. Gorton GE, Hebert DA, Gannotti ME. Assessment of the kinematic variability among 12 motion analysis laboratories. Gait Posture. 2009;29:398–402.

    Article  PubMed  Google Scholar 

  15. Osis ST, Hettinga BA, Macdonald S, Ferber R. Effects of simulated marker placement deviations on running kinematics and evaluation of a morphometric-based placement feedback method. PLoS One. 2016;11:e0147111.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Schütte KH, Maas EA, Exadaktylos V, Berckmans D, Venter RE, Vanwanseele B. Wireless tri-axial trunk Accelerometry detects deviations in dynamic Center of Mass Motion due to running-induced fatigue. PLoS One. 2015;10:e0141957.

    Article  PubMed  PubMed Central  Google Scholar 

  17. McGregor SJ, Busa MA, Yaggie JA, Bollt EM. High resolution MEMS accelerometers to estimate VO2 and compare running mechanics between highly trained inter-collegiate and untrained runners. PLoS One. 2009;4:e7355.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Chan M, Estève D, Fourniols J-Y, Escriba C, Campo E. Smart wearable systems: current status and future challenges. Artif Intell Med. 2012;56:137–56.

    Article  PubMed  Google Scholar 

  19. Li RT, Kling SR, Salata MJ, Cupp SA, Sheehan J, Voos JE. Wearable performance devices in sports medicine. Sports Health. 2016;8:74–8.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Haghi M, Thurow K, Stoll R. Wearable devices in medical internet of things: scientific research and commercially available devices. Healthc Inform Res. 2017;23:4.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Pohl MB, Lloyd C, Ferber R. Can the reliability of three-dimensional running kinematics be improved using functional joint methodology? Gait Posture. 2010;32:559–63.

    Article  PubMed  Google Scholar 

  22. Osis ST, Hettinga BA, Macdonald SL, Ferber R. A novel method to evaluate error in anatomical marker placement using a modified generalized Procrustes analysis. Comput Methods Biomech Biomed Engin. 2015;18:1108–16.

    Article  PubMed  Google Scholar 

  23. Osis ST, Hettinga BA, Ferber R. Predicting ground contact events for a continuum of gait types: an application of targeted machine learning using principal component analysis. Gait Posture. 2016;46:86–90.

    Article  PubMed  Google Scholar 

  24. Gullstrand L, Halvorsen K, Tinmark F, Eriksson M, Nilsson J. Measurements of vertical displacement in running, a methodological comparison. Gait Posture. 2009;30:71–5.

    Article  PubMed  Google Scholar 

  25. D’Errico J. Movingslope [Internet]. MATLAB Cent. File Exch. 2007 [cited 2016 Jun 15]. Available from:

  26. Söderkvist I, Wedin PA. Determining the movements of the skeleton using well-configured markers. J Biomech. 1993;26:1473–7.

    Article  PubMed  Google Scholar 

  27. Ward JH. Hierarchical grouping to optimize an objective function. J Am Stat Assoc. 1963;58:236–44.

    Article  Google Scholar 

  28. Calinski T, Harabasz J. A dendrite method for cluster analysis. Commun Stat Theory Methods. 1974;3:1–27.

    Article  Google Scholar 

  29. Cohen J. Statistical power analysis for the behavioral sciences. 2nd ed. Hillsdale, NJ: Lawrence Erlbaum Associates; 1988.

    Google Scholar 

  30. Brandon SCE, Graham RB, Almosnino S, Sadler EM, Stevenson JM, Deluzio KJ. Interpreting principal components in biomechanics: representative extremes and single component reconstruction. J Electromyogr Kinesiol. 2013;23:1304–10.

    Article  PubMed  Google Scholar 

  31. Abdi H, Williams LJ. Principal component analysis. Wiley Interdiscip Rev Comput Stat. 2010;2:433–59.

    Article  Google Scholar 

  32. Phinyomark A, Hettinga BA, Osis S, Ferber R. Do intermediate- and higher-order principal components contain useful information to detect subtle changes in lower extremity biomechanics during running? Hum Mov Sci. 2015;44:91–101.

    Article  PubMed  Google Scholar 

  33. Noehren B, Hamill J, Davis I. Prospective evidence for a hip etiology in patellofemoral pain. Med Sci Sports Exerc. 2013;45:1120–4.

    Article  PubMed  Google Scholar 

  34. Kawabata M, Goto K, Fukusaki C, Sasaki K, Hihara E, Mizushina T, et al. Acceleration patterns in the lower and upper trunk during running. J Sports Sci Routledge. 2013;31:1841–53.

    Article  Google Scholar 

  35. Lin S-P, Sung W-H, Kuo F-C, Kuo TBJ, Chen J-J. Impact of center-of-mass acceleration on the performance of ultramarathon runners. J Hum Kinet De Gruyter Open. 2014;44:41–52.

    Google Scholar 

  36. Grau S, Maiwald C, Krauss I, Axmann D, Horstmann T. The influence of matching populations on kinematic and kinetic variables in runners with Iliotibial band syndrome. Res Q Exerc Sport. 2008;79:450–7.

    Article  PubMed  Google Scholar 

  37. Watari R, Kobsar D, Phinyomark A, Osis S, Ferber R. Determination of patellofemoral pain sub-groups and development of a method for predicting treatment outcome using running gait kinematics. Clin Biomech. 2016;38:13–21.

    Article  Google Scholar 

  38. Nedergaard NJ, Robinson MA, Eusterwiemann E, Drust B, Lisboa PJ, Vanrenterghem J. The relationship between whole-body external loading and body-worn Accelerometry during team-sport movements. Int J Sports Physiol Perform. 2017;12:18–26.

    Article  PubMed  Google Scholar 

  39. Sinclair J, Shore HF, Taylor PJ, Atkins S. Sex differences in limb and joint stiffness in recreational runners. Hum Mov De Gruyter Open. 2015;16:137–41.

    Google Scholar 

  40. Lyle MA, Valero-Cuevas FJ, Gregor RJ, Powers CM. Control of dynamic foot-ground interactions in male and female soccer athletes: females exhibit reduced dexterity and higher limb stiffness during landing. J Biomech. 2014;47:512–7.

    Article  PubMed  Google Scholar 

  41. Lephart SM, Ferris CM, Riemann BL, Myers JB, Fu FH. Gender differences in strength and lower extremity kinematics during landing. Clin Orthop Relat Res. 2002;(401):162-9.

  42. Nakagawa TH, Moriya ÉTU, Maciel CD, Serrão FV. Trunk, pelvis, hip, and knee kinematics, hip strength, and gluteal muscle activation during a single-leg squat in males and females with and without Patellofemoral pain syndrome. J. Orthop. Sport. Phys. Ther. 2012;42:491–501.

    Article  Google Scholar 

  43. Selfe J, Janssen J, Callaghan M, Witvrouw E, Sutton C, Richards J, et al. Are there three main subgroups within the patellofemoral pain population? A detailed characterisation study of 127 patients to help develop targeted intervention (TIPPs). Br. J. Sports Med. 2016;bjsports-2015-094792.

  44. Henriksen M, Rosager S, Aaboe J, Graven-Nielsen T, Bliddal H. Experimental knee pain reduces muscle strength. J Pain Churchill Livingstone. 2011;12:460–7.

    Article  Google Scholar 

  45. Hirata RP, Arendt-Nielsen L, Shiozawa S, Graven-Nielsen T. Experimental knee pain impairs postural stability during quiet stance but not after perturbations. Eur J Appl Physiol. 2012;112:2511–21.

    Article  PubMed  Google Scholar 

Download references


Funding for this research was provided by the Ministry of Education of Brazil, Coordination for the Improvement of Higher Education Personnel (CAPES) [grant #9408–13-4]. Funding was also provided by the National Athletic Trainers Association – Research and Education Foundation through the Outcomes Grant program [Grant #808OUT003R], Alberta Innovates: Health Solutions [Grant #200700478], and a Discovery Grant (Grant No. 1028495) and Accelerator Award (Award No. 1030390) through the Natural Sciences and Engineering Research Council of Canada (NSERC). The funders had no role in the study design, collection, analysis and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.

Availability of data and materials

Data used for the present study will not be made available, as they come from a database managed by the Running Injury Clinic. Further inquiries about the database can be sent to RF.

Author information

Authors and Affiliations



RW was responsible for most of the data processing and analysis for the PCA and HCA, and interpretation of the results. SO manages the Running Injury Clinic database and devised methodologies for processing of the 3D gait kinematic data. AP contributed to the interpretation of the results. RF, along with the first author, was the idealizer of the study, participated in the interpretation of the results and was a major contributor in writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Reed Ferber.

Ethics declarations

Ethics approval and consent to participate

The University of Calgary’s Conjoint Health Research Ethics Board (CHREB) approved the collection of the data, storage within a research database, and extraction of anonymous data for analysis (Ethics ID E-24519). All participants provided written informed consent to participate.

Competing interests

There are no financial or personal relationships with other people or organizations that could potentially and inappropriately influence (bias) the submitted work and conclusions.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Watari, R., Osis, S.T., Phinyomark, A. et al. Runners with patellofemoral pain demonstrate sub-groups of pelvic acceleration profiles using hierarchical cluster analysis: an exploratory cross-sectional study. BMC Musculoskelet Disord 19, 120 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: