Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Unsupervised Event Characterization and Detection in Multichannel Signals: An EEG application
Next Article in Special Issue
Reliability of Sleep Measures from Four Personal Health Monitoring Devices Compared to Research-Based Actigraphy and Polysomnography
Previous Article in Journal
Color-Space-Based Visual-MIMO for V2X Communication
Previous Article in Special Issue
Measurements of Generated Energy/Electrical Quantities from Locomotion Activities Using Piezoelectric Wearable Sensors for Body Motion Energy Harvesting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electromyographic Patterns during Golf Swing: Activation Sequence Profiling and Prediction of Shot Effectiveness

1
Intelligent Systems Laboratory, Centre for Applied Intelligent Systems Research, Halmstad University, Kristian IV:s väg 3, PO Box 823, Halmstad S-30118, Sweden
2
Department of Electrical Power Systems, Kaunas University of Technology, Studentu 50, Kaunas LT-51368, Lithuania
3
Department of Information Systems, Kaunas University of Technology, Studentu 50, Kaunas LT-51368, Lithuania
4
School of Business, Engineering and Science, Halmstad University, Kristian IV:s väg 3, PO Box 823, Halmstad S-30118, Sweden
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(4), 592; https://doi.org/10.3390/s16040592
Submission received: 3 March 2016 / Revised: 11 April 2016 / Accepted: 17 April 2016 / Published: 23 April 2016
(This article belongs to the Special Issue Wearable Biomedical Sensors)

Abstract

:
This study analyzes muscle activity, recorded in an eight-channel electromyographic (EMG) signal stream, during the golf swing using a 7-iron club and exploits information extracted from EMG dynamics to predict the success of the resulting shot. Muscles of the arm and shoulder on both the left and right sides, namely flexor carpi radialis, extensor digitorum communis, rhomboideus and trapezius, are considered for 15 golf players (∼5 shots each). The method using Gaussian filtering is outlined for EMG onset time estimation in each channel and activation sequence profiling. Shots of each player revealed a persistent pattern of muscle activation. Profiles were plotted and insights with respect to player effectiveness were provided. Inspection of EMG dynamics revealed a pair of highest peaks in each channel as the hallmark of golf swing, and a custom application of peak detection for automatic extraction of swing segment was introduced. Various EMG features, encompassing 22 feature sets, were constructed. Feature sets were used individually and also in decision-level fusion for the prediction of shot effectiveness. The prediction of the target attribute, such as club head speed or ball carry distance, was investigated using random forest as the learner in detection and regression tasks. Detection evaluates the personal effectiveness of a shot with respect to the player-specific average, whereas regression estimates the value of target attribute, using EMG features as predictors. Fusion after decision optimization provided the best results: the equal error rate in detection was 24.3% for the speed and 31.7% for the distance; the mean absolute percentage error in regression was 3.2% for the speed and 6.4% for the distance. Proposed EMG feature sets were found to be useful, especially when used in combination. Rankings of feature sets indicated statistics for muscle activity in both the left and right body sides, correlation-based analysis of EMG dynamics and features derived from the properties of two highest peaks as important predictors of personal shot effectiveness. Activation sequence profiles helped in analyzing muscle orchestration during golf shot, exposing a specific avalanche pattern, but data from more players are needed for stronger conclusions. Results demonstrate that information arising from an EMG signal stream is useful for predicting golf shot success, in terms of club head speed and ball carry distance, with acceptable accuracy. Surface EMG data, collected with a goal to automatically evaluate golf player’s performance, enables wearable computing in the field of ambient intelligence and has potential to enhance exercising of a long carry distance drive.

1. Introduction

The capability of hitting long distances in golf always had a lot of attention in both the research and coaching literature. Recent studies highlight the importance of driving distance in relation to golf performance [1], and the strategic advantage gained by hitting irons further has been theorized to give advantageous ball trajectory conditions, which increase the possibility of higher accuracy [2]. Identifying factors that influence the ability to hit further with the irons is an important area for developing golf performance. Initial ball velocity is dependent on the centeredness of impact, club head velocity (i.e., magnitude and direction) and club face orientation [3,4]. Understanding the muscle activation patterns that characterize an effective golf swing has potential to facilitate golf coaching [5], as well as injury prevention [6].
Surface electromyography (EMG) data-based predictive models have been applied in a variety of fields. Prostheses control by mapping EMG signals into several classes of forearm motions [7,8,9,10,11,12,13,14,15,16,17,18,19], evaluation of local muscle fatigue in static [20,21,22,23,24], as well as dynamic exercises [25,26,27,28,29], screening for neuromuscular disorders [30,31,32,33] and analysis of the gait of arthritic patients [34,35] are probably the main application fields. Prediction of externally-applied forces to human hands by relating the EMG signals and experimentally known forces [36], prediction of forces on the lumbar spine [37], estimation of the iso-kinetic knee torque [38], recognition of simple and complex finger flexions [39], anticipation of head motion [40] and finger joint angle [41] for virtual-environment applications, distinguishing between normal subjects and subjects that experienced stroke of several degrees of severity (severe, moderate and mild) [42], measuring age-related reduction in number of motor units [43], a computer interface for text typing based on EMG signals induced by wrist movements [44], hand gesture sensing [45], assessing the subject’s vigilance level [46], stride detection and stride length estimation during walking [42] and classification of different tremor types [47] are other recent examples of using EMG data-based predictive models.
While predictive modeling using EMG features as predictors is a popular topic in various EMG signal analysis-based applications, prediction of shot success from EMG activity during golf swing is rarely researched. All EMG studies relevant to the golf swing were synthesized in [48], enabling the explanation of why certain muscles are active during specific phases of the swing in relation to the movement of the body during the swing. A recent review in [49] found 19 articles regarding EMG variables of the golf swing, but none of the reviewed articles is dedicated to studying relations between EMG signals and data characterizing biomechanical swing kinematics and/or the effectiveness parameters of the resulting shot. A comparison of two different threshold methods for EMG activity onset timing during the golf swing [50] found that club type has no influence on onset detection, and there is a similar activation time between muscles on the right and left body sides.
Based on the swing kinematics, some studies [51,52] find that neither upper torso nor pelvis rotation alone, but rather the X-factor (upper torso and pelvis separation) and delayed initiation of movement of the arms and wrists in the downswing are significantly and strongly correlated with ball velocity. A biomechanical analysis of professional swings in [53] also concludes that the effectiveness of the swing for a long drive appears to be related to a sequential activation of muscle groups that optimizes transfer of momentum through the kinetic chain, and high club velocity can be achieved. It is widely agreed that a biomechanically-correct transition, as far as the golf swing progression is concerned, works up from the hips (pelvis) to the chest (thorax), then the lead arm and, finally, the club.
Noting the lack of research on EMG temporal patterns generated during a golf swing and the potential of using information on EMG dynamics for the prediction of shot success, we try to address these aspects in our study. Predictive models are constructed here to predict club head speed or ball carry distance from EMG features. Data from golf players are described in Section 2. Proposed and employed methods are outlined in Section 3. The performed experiments are summarized in Section 4. Results are reviewed in Section 5, and final remarks are made in Section 6.

2. Equipment and Data

Fifteen elite golf players, 10 men and five women (mean age = 22 ± 2 years), with a handicap of −1.5 or better, were considered for the study. All participants were right-handed. Each subject was instructed to take five “normal” swings with their own 7-iron and to hit into a net, eight meters from the tee area. The total swing time ranged between 1.5 and 2.2 s.
Surface EMG signals were collected with bipolar surface electrodes collected by an eight-channel EMG recorder (ME Muscle Tester 6000 Megawin Electronics, Kuopio, Finland), recorded at a 1000-Hz sampling frequency. Electrodes were placed with a 20-mm inter-electrode distance along the muscle belly of each muscle, with a reference electrode placed perpendicular to the measurement electrodes to minimize EMG cross-talk. The muscles investigated were the left and right hands’ extensor digitorum communis (EDC), flexor carpi radialis (FCR) and the left and right side of rhomboideus and trapezius. Exact locations and mapping of sensors on the body are illustrated in Figure 1. More details on the placement of sensors are in Table A1.
The Flightscope X2 Doppler-radar launch monitor system (Flightscope, Version 6.0.9, Stellenbosch, South Africa) measured ball carry distance, lateral distance (left/right), club and ball speed, swing height, flight time, horizontal and vertical launch angles and vertical impact angle. Club head speed and ball carry distance parameters were selected as measures of shot effectiveness. Figure 2 shows the effectiveness of each shot with respect to these parameters. Pearson’s correlation coefficient between speed and distance was statistically significant at 99% confidence (p-value < 0.01) and positive, which conforms to previous research [54]. Shots of men had moderate (0.538) and shots of women had strong (0.647) positive correlation. Men had a noticeably higher sex-specific average for both parameters. The most effective player was G8 for men and G4 for women, whereas the least effective player was G11 for men and G3 for women.

3. Methodology

The most popular techniques used to extract features characterizing EMG signals are autoregressive modeling, higher order frequency moments, multi-resolution wavelet analysis, the maximum eigenvalues of the time-varying covariance patterns between EMG channels, signal variance values computed in overlapping windows, techniques based on fractal analysis, principle component analysis, independent component analysis, non-negative matrix factorization, short time Fourier transform and the orthogonal fuzzy neighborhood discriminant analysis.
We use time domain features in this work, assuming that such features are better at describing the underlying muscle orchestration than frequency domain features. By muscle activation pattern, we refer to the orchestration of EMG activity onsets in eight muscles of the arm and shoulder, namely right and left sides for extensor digitorum communis (EDC), flexor carpi radialis (FCR), rhomboideus and trapezius, during a golf swing. A comparison of onset sequences through EMG channels between swings of the same player and between different players is done visually.
Channel-specific statistics of filtered EMG signal, pairwise correlations between filtered EMG dynamics in various muscles during swing, the details of EMG signal peaks in each muscle and the concurrency of peak positions are computed as predictive features for the random forest model. Using computed features, the regression task seeks to estimate the target attribute value (club head speed or ball carry distance), and the classification task seeks to detect an effective shot, where effectiveness corresponds to being above the personal average of the player in question.

3.1. Muscle Activation Sequence Profiling

We expect to obtain some insights into golf swing patterns by studying muscle activation profiles, characterized by EMG signals recorded from different muscles. Figure 3 present two examples of EMG signals recorded from two golf players. On the X axis is shown a sample number (1000 samples/s is the sampling frequency), while the Y axis reflects the absolute value of an EMG signal amplitude. Since noise is always present in EMG signals, some preprocessing of EMG signals is needed to obtain robust characterization of the muscle activation sequence. Therefore, in addition to the raw EMG signal, information resulting from signal preprocessing is also shown in Figure 3.
We define a muscle activation profile as a vector, x , of time differences between the start of swing-related activation in the reference muscle and the other seven muscles. Thus, x consists of seven components x = x 1 , . . . , x 7 . Channel 8 (left trapezius) was used as the reference for activation times in this work. The start of swing-related activations, marked by a dashed vertical line ending with a circle in Figure 3, are to be determined in all eight channels, to obtain components of x . This is not an easy task, since rather substantial activation can be observed in some muscles, even before the golf swing start (see Figure 3). To detect the beginning of swing-related activations in different channels, the following processing steps are performed.
  • Filter the signal | s ( t ) | in each channel i by four sliding Gaussian filters G 1 , G 2 , G 3 and G 4 of the following width: 256, 512, 1024, 2048 samples. The resulting signals v 1 i ( t ) , . . . , v 4 i ( t ) are shown in Figure 3 by the red, green, blue and black curves.
  • For each channel i, determine a threshold value T i :
    T i = α c a r d { } t = t v 4 m a x i - δ t = t v 4 m a x i + 2 δ | s i ( t ) | , i = 1 , . . . , 8
    where c a r d { } means the number of elements in the sum, α and δ are parameters selected experimentally and t v 4 m a x i stands for the position of the maximum value of the i-th channel signal filtered by the G 4 filter. Values of α = 1 . 2 and δ = 500 worked well in our study. The determined thresholds are shown on the Y axis in Figure 3.
  • The start of swing-related activation in the i-th channel t s i is then given by:
    t s i = arg max t = m i n t ( | t - t T i v 2 , 3 , 4 | ) v 1 i ( t ) t
    where t T i v 2 , 3 , 4 stands for three time moments, where v 2 i ( t ) = T i , v 3 i ( t ) = T i and v 4 i ( t ) = T i .
  • Define λ as the maximum difference in activation time between the reference channel and all remaining channels:
    λ = m a x i = 1 , . . . , 7 ( | t s 8 - t s i | )
    Inspect and post-process the results manually if λ < λ m i n or λ > λ m a x , where global λ m i n and λ m a x are determined experimentally from all EMG channels and swings of all participants.
The determined t s i values constitute components x i and x . In Figure 3, the determined t s i values are shown by dashed vertical lines with a circle on the top.

3.2. EMG Features for Predictive Modeling

Pre-processing of an EMG signal usually corresponds to full-wave rectification, followed by filtering of the rectified signal to obtain a smooth linear envelope of the EMG. Most software packages for EMG analysis in sports biomechanics use a fourth-order Butterworth filter [55]. We apply a second-order Butterworth filter for both high- and low-pass filtering. To eliminate phase lag, resulting from a second-order filter, the signal is filtered again in the reverse direction, which effectively converts the filter to a fourth-order filter. Low-pass filtering follows after demeaning the high-pass filtering result, and the cut-off frequency is 20 and 400 Hz for high- and low-pass, respectively. An example of the filtered signal is provided in the bottom part of Figure 4. The Butterworth high-pass filter with a corner frequency of 20 Hz is also recommended by [56] for general use.
Peak detection, which is at the core of feature extraction here, is done using a MATLAB function ipeak from the package ipeak7.zip (Version 7.1), created and maintained by Prof. Thomas C. O’Haver [57]. A program window with the ipeak function result for the filtered EMG signal from a single channel is displayed in Figure 5. To find peaks and extract a swing segment from the surface EMG signal, the following steps are performed.
  • Normalize the rectified signal | s ( t ) | in each channel i to the peak activation level (see Figure 4):
    n o r m a l i z e ( s ) = s ( t ) max ( s ( t ) )
  • Filter the normalized signal of each channel i using the high- and low-pass Butterworth filters (the result can be seen as the resulting black curve in the bottom part of Figure 4).
  • In each channel i, detect all peaks in the filtered signal using the ipeak function (with maximum peak density parameter PeakD = 2), and retain the two highest peaks, denoting them in the order of appearance as P 1 i and P 2 i . In Figure 5, Peak # 10 would correspond to P 1 i , and Peak # 11 would correspond to P 2 i .
  • Pool absolute time-based positions of retained peaks from all channels into a single vector of positions p . Enumerate all possible i values into two identical lists idx 1 and idx 2 of channel indexes.
  • Check for outliers among the elements of the vector p using the box-plot rule, finding values exceeding the 1.5 inter-quartile range above/below the upper/lower quartile. If outliers are found, remove the numbers of the corresponding channels from idx 1 and/or idx 2 .
  • Compute swing start ( t 1 ) and end ( t 2 ) times by using the position and width of the marginal P 1 and P 2 peaks:
    k = arg min i idx 1 ( p o s ( P 1 i ) ) and t 1 = p o s ( P 1 k ) - w i d t h ( P 1 k ) 2
    k = arg max i idx 2 ( p o s ( P 2 i ) ) and t 2 = p o s ( P 2 k ) + w i d t h ( P 2 k ) 2
    where k is the channel number of the marginal peak and functions p o s and w i d t h return the position and width of the peak.
  • Retain only the segment from t 1 to t 2 .
  • If outliers were found (in Step # 5), proceed from Step # 3.
  • Convert an absolute position of each peak within the signal into a relative position within the swing by taking into account the swing duration ( t 2 - t 1 ).
After extracting the swing segment, the feature vector x is determined by the feature set, where the very first feature x 0 is always sex (0, woman; 1, man), and the remaining features are related to: (1) channel statistics and correlations; (2) properties of the first peak; (3) properties of the second peak; (4) properties derived from both peaks.
Table 1 presents a list of feature sets used in this work. Statistics, used to enrich Feature Sets # 1, 2, 5, 6, 8, 9, 11, 12, 15–22 with various aspects of the distribution in question, encompass the following 12 characteristics: minimum, maximum, mean, median, lower quartile ( Q l o ), upper quartile ( Q u p ), standard deviation, inter-quartile range, lower range (median, Q l o ), upper range ( Q u p , median), skewness and kurtosis. Since there are 4 channels on one side and the EMG dynamics of each channel is compressed into 12 characteristics, Feature Sets # 1, 2 (ChanStats) contain 48 (4 × 12) features plus an extra feature for the duration in seconds; therefore, 49 features in total.
The main properties of a peak, estimated using the ipeak function (see the top part of Figure 5), are the following: position, height, width, area and fitting error. All peak properties are used in Feature Sets # 7, 10, 13, 14, while Feature Sets # 8, 9, 11, 12, 15–22 are calculated using peak position only. Feature Sets # 13, 14 are based on the mathematical operation between properties of the first and second peaks, where Feature Set # 13 uses subtraction and Feature Set # 14 uses the division of the corresponding properties. Features Sets # 7, 10, 13, 14 are constructed considering 5 properties for each of the 8 channels; therefore, they contain 40 (5 × 8) features in total.
Feature Sets # 3–6, 8, 9, 11, 12, 15–22 are derived from all possible channel pairs. Having 8 channels, we obtain 28 arrangements. Correlation-related feature sets are based on the arithmetical operation between Spearman and Pearson correlations for each pair of channels. Feature Sets # 3, 5 (containing CorrDeltain their name) consider subtraction, and Feature Sets # 4, 6 (containing CorrRatioin their name) consider the division of correlations for the corresponding channel pair.
Peak position-based Feature Sets # 8, 9, 11, 12, 15–22 (containing Sync in their name) evaluate peak position concurrency for each pair of channels. The concurrency is measured here as a signed difference for Feature Sets # 8, 11, 15–18 or an absolute difference for Feature Sets # 9, 12, 19–22 (containing Absin their name). Feature Sets # 8, 9, 11, 12 are constructed using the difference between positions of the 1st and 2nd peaks separately. These separate concurrences are further combined by the basic arithmetical operations: subtraction in Feature Sets # 15, 19; addition in Feature Sets # 16, 20; multiplication in Feature Sets # 16, 21; division in Feature Sets # 17, 22.

3.3. Prediction Using Random Forest

Regression models seek to predict club head speed or the ball carry distance of a shot. Binary classification models try to detect the personal effectiveness of speed or distance, that a resulting shot for a player in question is effective or ineffective in terms of speed or distance. Regression has a continuous, while detection has a discreet target attribute. Discretization into labels for the detection task is done with respect to the player-specific average: a shot resulting in being above the personal average club head speed or ball carry distance is considered as effective, and a shot resulting in being below the personal average club head speed or ball carry distance is considered as ineffective.
We used random forest (RF) [58] for predictive modeling. The MATLAB port of the original RF algorithm was obtained from [59]. RF can be used for both classification and regression. In classification, RF is a committee of decision trees, where the final output is derived from voting (see Figure 6). In regression, voting by a committee is replaced with averaging. Random forests have been successfully applied in a variety of fields [60]. The core idea of RF is to combine many (B in total) decision trees, built using different bootstrap samples of the original dataset and a random subset (of predetermined size q) of features x 1 , , x p (see Figure 6). RF is known to be robust against over-fitting, and as the number of trees increases, the generalization error converges to a limit [58]. Low bias and low correlation are essential for the robust generalization performance of the ensemble. To get low bias, trees are unpruned (grown to the maximum depth). To achieve the low correlation of trees, randomization is applied. RF is constructed in the following steps:
  • Choose the forest size B as a number of trees to grow and the subspace size q p as a number of features to provide for each tree node.
  • Draw a bootstrap sample (random sample with replacement) of the learning dataset, which generally results in 2 3 · n unique observations for training, thus leaving 1 3 · n for testing as the out-of-bag (OOB) dataset for that particular tree.
  • Grow an unpruned tree using the bootstrap sample. When growing a tree, at each node, q variables are randomly selected out of the p available.
  • Repeat Steps 2 and 3 until the size of the forest reaches B.
Because of the “repeated measures” aspect in golf data, where each player is represented by several swings, sampling part of the RF had to be modified to ensure that all swings of each player are either included in a bootstrap sample or left aside as OOB. The list of passed parameters to the external interface function (C++ source code to build MATLAB mex files) was expanded by the additional indexList parameter, containing subject’s unique ID for each swing. The implemented modification of bootstrap sampling converts leave-one-out validation into leave-one-subject-out validation. Furthermore, the RF setting, allowing to perform stratified sampling, was configured in the classification task to preserve the sex ratio of the full dataset in each drawn bootstrap sample.

3.4. Decision-Level Fusion

Individual RFs were built independently on bootstrap sets, and the decisions of these individual experts were combined in a meta-learner fashion. RF was used in this work both as a base learner and as a meta learner. This implies that outputs from RF models from the first stage are treated as inputs (meta-features) for another RF in the second stage. The set of meta-features was slightly expanded by always including sex and intentionally avoiding removal of sex in the decision optimization stage.
For the classification task, an input to the meta-learner is the difference between class a posteriori probabilities obtained from the base-learner. Given a trained RF, this difference is estimated as:
d ( { t 1 , . . . , t L } , x ) = i = 1 L f ( t i , x , q = 2 ) L - i = 1 L f ( t i , x , q = 1 ) L
where x is the object being classified, L is the number of trees t 1 , . . . , t L in the random forest for which observation x is OOB, q is a class label and f ( t i , x , q ) stands for the q-th class frequency in the leaf node, into which x falls in the i-th tree t i of the forest:
f ( t i , x , q ) = n ( t i , x , q ) j = 1 Q n ( t i , x , q j )
where Q is the number of classes and n ( t i , x , q ) is the number of training data from class q and falling into the same leaf node of t i as x .
For the regression task, input to the meta-learner is the averaged OOB prediction of the base-learner. Given a trained RF, this estimate is provided by:
d ( { t 1 , . . . , t L } , x ) = i = 1 L f ( t i , x ) L
where x is the data instance represented by the feature set in question, L is the number of trees t 1 , . . . , t L in the random forest for which the observation x is OOB and f ( t i , x ) corresponds to the predicted value of target attribute.
Meta-features were also investigated by performing permutation-based variable importance analysis. The mean decrease in accuracy (for the detection task) or the mean decrease in root mean squared error (for the regression task) were considered here as variable importance measures. Each meta-feature is processed by permuting its values several times, and the mean difference in RF performance on OOB data is estimated.

3.5. Decision Optimization

In ensemble learning, non-generative ensemble or decision optimization [61,62] refers to model selection through tuning of the meta-learner and/or selection of meta-features. RF parameter B was fixed throughout experiments, but parameter q was selected in a grid search (for individual RFs) or brute force (for combiner RF) fashion. All features were used for individual RFs. Decision optimization to search out a subset of meta-features was done only for the RF combiner through a specific variant of the evolutionary algorithm. Each generated meta-feature subset was evaluated with all possible q values, and the best evaluation, summarized using the cost of log-likelihood ratio (Cllr, see Section 3.6.1 for details), is returned as a result of the fitness function.
We have chosen angle-modulated differential evolution (AMDE) as the decision optimization technique for wrapper-based meta-feature selection in the meta-learner model. AMDE was initially proposed in [63] as a homomorphous mapping technique enabling the differential evolution algorithm to operate within a binary space. In [64], AMDE was used to select features for support vector machine, while in [65], AMDE was applied for decision optimization of meta-learner RF.
The angle modulation (AM) technique uses a generating function to transform a 4-dimensional continuous-valued space into a larger binary-valued n-dimensional Hamming space (n is the total number of individual feature sets), which corresponds to convenient tuning of a feature inclusion mask by optimizing 4 continuous parameters. The transformation is achieved using the following trigonometric function:
g ( x ) = s i n ( 2 π ( x - a ) × b × c o s ( 2 π × c ( x - a ) ) ) + d
where the coefficients a, b, c and d are four optimized parameters determining the shape of the generating function (a, horizontal shift; b, maximal s i n frequency; c, maximal c o s frequency; d, vertical shift) and x is an element from a set (size n) of evenly-spaced values, for example, starting from 0.1 with a step of 0.1 until 0.1*n is reached. We generate a binary selection vector g by:
g k = 1 if g ( x k ) > 0 0 if g ( x k ) 0
where k = 1 p and the g ( x k ) value is calculated using Equation (10). By proceeding with k, a vector g of length n is generated. The vector g represents the potential solution within the binary problem space; in our case, a binary mask of the meta-features to select.

3.6. Evaluating Predictive Performance

Generalization performance of RF was evaluated by OOB validation. It is well known that for large RF, each observation of a learning set is used as an OOB observation, and the OOB performance estimate corresponds to an unbiased estimate of a test set error.

3.6.1. Assessing Detection

To evaluate the goodness of detection, detector’s scores for OOB data were obtained. Votes of RF were converted to a proper score vector by normalizing votes for a specific class through division by the total number of times the case was OOB. Using soft decisions (scores), instead of hard decisions, enables more precise evaluation of the detection. We evaluated scores from RF using several measures: (1) the cost of the log-likelihood ratio (Cllr); (2) the detection error trade-off (DET) curve and the equal error rate (EER); (3) the receiver operating characteristic (ROC) curve and the area under the curve (AUC).
Various thresholds are possible to convert a soft decision into hard, resulting in a multitude of confusion matrices, and the overall performance of the detector can be visually summarized by the ROC or DET curve. Due to the logarithmic scale, the DET curves allow a comparison of several systems at a glance more easily than the ROC curves [66]. A quick way to compare the accuracy of detectors with different DET (ROC) curves is the equilibrium point, widely known as the equal error rate (EER). EER is the point where the DET (ROC) curve intersects the diagonal and: (1) false positive rate (miss rate) = false negative rate (false alarm rate), for DET; or (2) true positive rate (sensitivity) = true negative rate (specificity), for ROC.
The cost of the log-likelihood-ratio (Cllr) is the most comprehensive goodness-of-detection metric, used here as the main criterion for model selection and ranking of predictive models. The log-likelihood-ratio is the logarithm of the ratio between the likelihood that the shot was effective and the likelihood that the shot was ineffective. More details on Cllr are provided in [67]. The Cllr, DET, EER, ROC and AUC measures were estimated using the ROC convex hull method, available in the BOSARIS toolkit [67]. In a detection task, a well-calibrated and useful detector should have Cllr < 1, EER < 50% and AUC > 0.5.

3.6.2. Assessing Regression

Two types of errors were used: the root mean squared error (RMSE) and the mean absolute percentage error (MAPE). Pearson’s correlation coefficient between actual and predicted values was also computed. RMSE is used here as the main criterion for model selection and ranking of predictive models.
The reference values for RMSE, MAPE and correlation were estimated by considering a variant of naive prediction, which simply outputs a sex-specific average instead of a real attempt at prediction. Satisfactory performance, which is better than using a sex-specific average as prediction, can be confirmed when: (1) RMSE < 6.91, MAPE < 4.12 and correlation > 0.742 for club head speed; (2) RMSE < 13.37, MAPE < 8.15 and correlation > 0.738 for ball carry distance.

4. Experimental Investigations

4.1. Muscle Activation Profiles

To get deeper insights into the variations of the muscle activation profiles, the profile representations x = x 1 , . . . , x 7 were visualized. Figure 7 presents the resulting map of the profiles. Activation profiles indicate that each individual player exhibits a rather stable activation pattern. Stable refers to most players having both stable activation sequences and time differences, except for players G3 and G10 (see Figure 7). Several players also exhibit almost identical activation patterns. Therefore, it is interesting to investigate the relation of the activation profile to parametrically- or biomechanically-effective transition. Nonetheless, profiles themselves can hardly be suitable as EMG features for prediction of player-specific effectiveness and exact club head speed or ball carry distance.
Interestingly, for golf player G15, correct biomechanically (see Appendix A), but not among the most effective players of the same sex, EMG activity in muscles related to the Channels 2 (right EDC), 3 (left FCR) and 5 (right rhomboideus) starts before the reference muscle (left trapezius) is activated, but EMG activity in muscle related to Channel 7 (left rhomboideus) starts after activation of the reference muscle (left trapezius).
Regarding the relation of muscle activation sequence to sex-specific effectiveness, the following insights were derived. Effective men (G5, G8 and G9) and women (G4 and G14) players had the reference muscle as lagging, and no channels were leading the reference; all channels were activated noticeably before the reference or simultaneously with it (hence, no activation times with high positive values). Profiles for effective men (G5, G8 and G9) had almost an identical shape, where Channels 1, 4, 8 were activated in a synchronized fashion and remaining channels in a specific order earlier. Profiles for effective women (G4 and G14) were distinct with regard to Channels 3–4, where player G14 with more scattered shots had Channel 4 synchronized to Channels 1 and 8 (as in effective men), and Channel 3 was activated a bit earlier; but similarly, effective player G4 had a rather unique profile shape with less synchronization to the reference channel. Most profiles (with the exception of G3, G6, G15 and G17 players), and especially effective men and women, revealed a consistent sequential order: avalanche effect, where Channels 5–7 formed a chain of activations with respect to the reference (Channel 8): first 5, then 6, then 7 and, finally, 8.

4.2. Predictive Performance

The total number of trees in RF was 5000. Several specific values of q ( p , 2 · p , p) in prediction using individual feature sets and all possible values of q (2–22) in decision-level fusion were tested and the best performing q setting retained.

4.2.1. Detection Performance

The detection task was more successful at predicting club head speed than ball carry distance, both when using individual feature sets (see Table 2) and when combining them through decision-level fusion (see Table 3 and Figure 8). The lowest EER using a single feature set was 32% for the speed and 39% for the distance (see Table 2).
The best individual feature sets were quite different with respect to the target attribute: the difference between correlations (CorrDelta and StatsCorrDelta) for speed and the difference between absolute concurrencies of peaks (PeakAbsSyncDelta), besides the absolute concurrency of the first peak (PeakP1AbsSync) for distance prediction.
Decision-level fusion together with decision optimization achieved EER of 24.32% for the club head speed and 31.71% for the ball carry distance target attributes (see Table 3 and Figure 8). The decision optimization step reduced the need for various feature sets by more than half (see left side of Figure 9): ∼41% (nine out of 22) of meta-features were retained for speed and ∼36% (eight out of 22) of meta-features were retained for distance prediction.
Ranking of meta-features revealed similar tendencies as individual rankings of feature sets for speed, but not for distance. For example, the top three rankings for speed are identical (compare Rank of the club head speed in Table 2 with the left side of Figure 9 or see the detection part of Table B1 for convenience), but the top three rankings for distance are distinct (compare Rank of ball carry distance in Table 2 with the left side of Figure 9 or see the detection part of Table B2 for convenience). It appears that when combining feature sets through decision-level fusion, additional feature sets can become more important than individually best ones. For the distance case, such complimentary feature sets are the statistics of left channels (LeftChanStats) and the absolute concurrency of the second peak (PeakP2AbsSync).

4.2.2. Regression Performance

The regression models were more successful in predicting club head speed than ball carry distance, both when using individual feature sets (see Table 4) and when combining them (see Table 5 and Figure 10). The lowest MAPE using a single feature set was 4.7% for the speed and 8.7% for the distance (see Table 4). Such poor performance of individual feature sets for the regression task is unsatisfactory, since all performance metrics (RMSE, MAPE and correlation) appear to be worse than a naive guess that the effectiveness parameters of any swing in question simply are equal to the sex-specific average.
Best individual feature sets were quite different with respect to the target attribute: properties of the first peak (PeakP1), besides the difference between correlations (StatsCorrDelta) for the speed, and the ratio between absolute concurrencies of first and second peaks (PeakSyncDiv), besides statistics of left channels (LeftChanStats) for the distance prediction.
Decision-level fusion together with decision optimization proved to be successful in improving regression up to the satisfactory performance level and achieved MAPE of 3.22% for club head speed and 6.38% for ball carry distance target attributes (see Table 5). Nonetheless, the club head speed prediction (see the left side of Figure 10) defaults to the sex-specific average for women players, which could imply that exact club head speed prediction is not an easy task, either due to the task itself or due to the small sample of women participants. The decision optimization step reduced the need of various feature sets, but the reduction was more prominent for the speed than for the distance (see the right side of Figure 9): ∼41% (nine out of 22) meta-features retained for speed versus∼64% (14 out of 22) meta-features retained for distance prediction.
Ranking of the meta-features revealed no similar tendencies with individual rankings of feature sets, irrespective of the target attribute. For convenient comparison, see Appendix B (the regression part in Table B1 and Table B2). Here, we reiterate that when combining feature sets through decision-level fusion, additional feature sets can reveal themselves as more important than individually best ones. For the speed case, such complimentary feature sets are the statistics of the channels corresponding to the muscles on the right side of the body (RightChanStats) and multiplication of concurrencies for the first and the second peaks (PeakSyncProd). For the distance case, such complimentary feature sets are the difference between correlations (StatsCorrDelta) and the absolute concurrency of the first peak (PeakP1AbsSync).

4.2.3. Summary of Performance

The top feature sets for the club head speed target attribute were statistics of activity in the muscles on the right side of the body (RightChanStats), difference between correlations (StatsCorrDelta and CorrDelta) and the ratio between absolute concurrencies of first and second peaks (PeakAbsSyncDiv). The top feature sets for the ball carry distance target attribute were the statistics of activity in the muscles on the left side of the body (LeftChanStats), absolute concurrency of the first peak (PeakP1AbsSync), properties of the second peak (PeakP2) and difference between absolute concurrencies of the first and second peaks (PeakAbsSyncDelta). Interestingly, even though muscle activity on both sides of the body is important, the importance of body side is target-related: the right side muscles are useful in predicting club head speed; meanwhile, the left side muscles are useful in predicting ball carry distance. For more details, see Appendix B.
In summary, channel statistics from both sides (RightChanStats and LeftChanStats) could be considered as the most informative features. The difference between Pearson and Spearman correlations for each channel pair with the statistical summary of these differences (StatsCorrDelta) is the next important feature set. Both peaks proved to be useful separately, but in rather different aspects: absolute concurrency of the first peak (PeakP1AbsSync) versus basic properties of the second peak (PeakP2). Feature sets combining position concurrency of both peaks, such as the ratio of absolute concurrencies (PeakAbsSyncDiv) or the difference between signed concurrencies (PeakSyncDelta), also appear to be valuable.

5. Discussion

The introduced muscle activation profiles helped to summarize an EMG activity onset sequence in the monitored channels during the golf swing and revealed that the profile pattern of each player is very stable. Moreover, effective players, especially among men, had almost identical profiles. For effective profiles, EMG activity onsets in right FCR, left EDC and left trapezius muscles happened almost simultaneously, whereas activity in right EDC, left FCR, right rhomboideus and right trapezius muscles was detected earlier and constituted a specific avalanche effect. The effect is defined as the following sequence of EMG activity onsets: right rhomboideus, then right trapezius and then left rhomboideus. Such a sequence is supported by previous studies [48,49], where right rhomboideus and right trapezius are reported as most active in the back swing. Furthermore, the delayed activation of the left trapezius, right FCR and left ECR muscles in the present study also conforms to previous findings, where high levels of activity in left trapezius [48] and forearm [49] muscles during forward swing were detected. Interestingly, we observed similar activation patterns of the right FCR and left EDC muscles among effective men and for one of the more effective women; however, the most effective woman had a rather unique activation pattern. Due to the small sample size, 10 men and five women, further research into sex-specific differences in forearm muscle activity is required. Increased knowledge on these differing and unique neuromuscular patterns may aid in the characterization of an effective movement pattern for the optimization of golf swing performance.
Developing and maintaining a golf swing characterized by biomechanically-correct transition is not an easy task. As can be seen in Appendix A, only one participant managed to sustain the correct order in all of his shots. This experimental data may indeed partially explain what many coaches experience and why they consider sequential movement of body segments a key technical parameter [5]. Nonetheless, even if correct biomechanically, the player was not among the most effective, possibly due to the lack of synchronization after the avalanche effect, which was revealed by activation sequence profiles.
Ranking of the feature sets found a difference in body side relevance with respect to the club head speed and ball carry distance target attributes. Results show that club head speed is better predicted from the statistics of the dynamics in the right side muscles, whereas ball carry distance is better predicted from the statistics of dynamics in the left side muscles.
Some EMG feature sets revealed their predictive power only when decisions from these feature sets were used in combination with decisions arising from other feature sets. Therefore, in Appendix B, rankings from various modeling aspects were estimated and averaged into several overall rankings, namely target-and-task-specific and target-specific ranks. Ultimate ranking, averaged over targets, revealed channel statistics from both body sides and correlation-related feature sets as the most informative. The peak-related feature sets were also valuable, especially for the ball carry distance target attribute, where of notable importance are the absolute concurrency of the first peak and the basic properties of the second peak. The ratio of absolute concurrencies and the difference between signed concurrencies of the peaks were useful for both target attributes.

6. Conclusions

The two highest peaks in the filtered EMG signal of each channel were used to extract swing segment, the part of EMG stream with the most pronounced muscle activity, potentially containing all swing phases, and calculate various feature sets. Peak-based modeling of EMG dynamics proved to be useful for predicting shot effectiveness. The detection of personal effectiveness was an easier task than regression, where an exact club head speed or ball carry distance had to be predicted. Despite the task, the individual feature sets were clearly outperformed by the decision-level fusion. Club head speed was more predictable than ball carry distance: the lowest EER in detection was 24.3% for speed and 31.7% for distance; the lowest MAPE in regression was 3.2% for speed and 6.4% for distance.
From the ranking of feature sets, we speculate that for right-handed golfers right side upper body muscle activity is directly related to club head speed, whereas left side upper body muscle activity could be related to adjustments of club face orientation directly before the moment of impact to achieve a successful long shot. EMG channel statistics from both body sides, correlations between channels and peak-related features (concurrency of the first peak, properties of the second peak, ratio of absolute concurrencies of the peaks and difference between signed concurrencies of the peaks) were among the most important predictors of golf shot effectiveness.
Muscle activation profiles helped to reveal individual sequences, which appeared to be stable from shot to shot, especially for the most effective players. EMG activity for effective profiles started with a common avalanche effect from a triad of muscles (right rhomboideus, right trapezius, left rhomboideus) followed by the simultaneous activation of another triad of muscles (right FCR, left EDC and left trapezius). The level of synchronicity in EMG activity onsets in right FCR, left EDC and left trapezius muscles at the moment of impact during golf swing appears to influence shot effectiveness more than the biomechanical correctness of swing kinematics.
Analysis of EMG data from more participants would be beneficial to strengthen the results and conclusions. Decision-level fusion of EMG features with biomechanical features, automatically extracted from video kinematics, could be an interesting future direction. For example, combining the prediction of club head speed with an estimated angle of ball trajectory could potentially improve the accuracy of ball carry distance prediction. Such contributions would expand the scope of wearable computing applications in golf, provide insights for golfers and instructors alike and help to improve carry distance drives.

Acknowledgments

We would like to thank Thomas C. O’Haver for sharing ipeak software publicly. This work has been partially supported by research Grant No. 2012/0319 from the Knowledge Foundation (KK-stiftelsen) of Sweden. The study was approved by the regional Swedish ethics committee (Lund, Dnr 2016/12). All of the participants were informed of the study design and were made aware of the fact that they could terminate their participation at any time without any questions being asked. Thereafter, subjects signed a written consent to participate in the study. The authors are grateful to all the participants.

Author Contributions

Antanas Verikas performed the literature review, supervised the research, contributed to discussion and conclusions, and helped in revising the manuscript. Evaldas Vaiciukynas designed EMG features, experimented with prediction, interpreted experimental results and prepared the manuscript. Adas Gelzinis devised muscle activation sequence profiling, interpreted experimental results and contributed to the manuscript. James Parker and M. Charlotte Olson were responsible for EMG recording setup, Doppler-radar launch monitor system, electromagnetic motion capture system, collected golf swing data from study participants and contributed to the manuscript and its revision.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A.

Data on swing kinematics were collected using a four-sensor Polhemus Liberty electromagnetic motion capture system (Polhemus Inc., Colchester, VT, USA), operating at 120 Hz. Fifteen landmarks were digitized to define segment lengths, orientations and joint axes (see last column of Table A1). The swing events were determined from sensors on the club; the topmost position of back swing was determined as the moment the sensor on the club started to rotate with a positive direction around the global Y axis; and impact was determined as the moment at which a rapid deceleration of the club sensor in the down swing was noticed. Angular velocities and displacements of the pelvis, thorax and lead arm were calculated using standard biomechanical principles in conjunction with advanced motion measurement software (AMM 3D, Phoenix, AZ, USA).
In Figure A1, negative speeds reflect movement away from the ball, the back swing, while positive speeds indicate movement towards the ball, the down swing. The transition phase is defined as the interval between two points of rotation direction change, namely the direction of the pelvis and the club. The sequence of colors, violet, green, red and cyan, corresponds to the correct order of transitions, hips (pelvis), chest (thorax), arm and club.
Table A1. Sensor and landmark positions.
Table A1. Sensor and landmark positions.
Body PartSensor PlacementLandmarks Used for Digitization
ClubBelow gripTop of grip
Hozel
Club head, bottom groove at heel
Club head, bottom grove at toe
Club head, top groove at toe
Left armBehind the elbowLeft acromion process
ThoraxOn T5Lateral epicondyle
Medial epicondyle
Thorax/upper-bodyOn T5Left acromion process
Right acromion process
Right side mid thorax, high
Right side mid thorax, low
PelvisSacrumLeft greater trochanter
Right greater trochanter
Point above left greater trochanter
Figure A1. Illustration of a biomechanically-correct transition in a golf swing. The golf swing progresses in time on the X (horizontal) axis from left to right, and the speed of player’s body parts is on the Y (vertical) axis.
Figure A1. Illustration of a biomechanically-correct transition in a golf swing. The golf swing progresses in time on the X (horizontal) axis from left to right, and the speed of player’s body parts is on the Y (vertical) axis.
Sensors 16 00592 g011
Figure A2 presents experimental data illustrating transition sequences for fourteen experienced golf players. In Figure A2, color change means the change of rotation direction, and ’0’ on the X axis stands for the moment of impact. The swings of two players (G8 and G10) were not recorded, and the swings of one extra player (G13) were included instead; also, data for some shots (of G2 and G7) were missing. Only one player (G15) maintains the correct order in all five shots.
Figure A2. Illustration of transition sequences for 14 golf players (∼5 shots each), where color change means the change of rotation direction. The sequence of colors (violet, green, red and blue) for player G15 corresponds to the correct transition sequence (pelvis, thorax, arm and club). Note: position in samples on the X (horizontal) axis is shown, where 0 corresponds to the moment of impact.
Figure A2. Illustration of transition sequences for 14 golf players (∼5 shots each), where color change means the change of rotation direction. The sequence of colors (violet, green, red and blue) for player G15 corresponds to the correct transition sequence (pelvis, thorax, arm and club). Note: position in samples on the X (horizontal) axis is shown, where 0 corresponds to the moment of impact.
Sensors 16 00592 g012

Appendix B.

Table B1. Ranking of feature sets for club head speed prediction. Ind., individual; Imp., importance; Avg., average.
Table B1. Ranking of feature sets for club head speed prediction. Ind., individual; Imp., importance; Avg., average.
#Feature SetRanking for Detection TaskRanking for Regression TaskOverall Ranking
Ind.Imp.Avg.RankInd.Imp.Avg.RankAvg.Rank
1RightChanStats475.54714.014.751
2LeftChanStats15912.081078.5310.256
3CorrDelta111.0121915.0158.003
4CorrRatio333.0322614.0138.504
5StatsCorrDelta222.0222212.097.002
6StatsCorrRatio62214.01142213.01113.5011
7PeakP1162219.01712211.5815.2514
8PeakP1Sync222222.0221639.5515.7515
9PeakP1AbsSync112216.5151449.0412.759
10PeakP252213.51013810.5712.008
11PeakP2Sync72214.512192220.52117.5019
12PeakP2AbsSync82215.013202221.02218.0020
13PeakDelta172219.51852213.51216.5017
14PeakRatio212221.52132212.51017.0018
15PeakSyncDelta988.5662214.01311.257
16PeakSyncSum192220.520152218.51919.5022
17PeakSyncProd142218.0161729.5513.7512
18PeakSyncDiv102216.01492215.51615.7515
19PeakAbsSyncDelta182220.019112216.51718.2521
20PeakAbsSyncSum1369.57122217.01813.2510
21PeakAbsSyncProd1248.05182220.02014.0013
22PeakAbsSyncDiv20512.59856.529.505
Ranking based on the performance of the individual feature sets (denoted as Ind. in Table B1 and Table B2) was averaged with ranking based on variable importance of the meta-learner RF (denoted as Imp. in Table B1 and Table B2) to reveal the final ranking (denoted as Rank in Table B1 and Table B2). Obtained task-specific rank averages (denoted as Avg. in Table B1 and Table B2) were further averaged to reveal target-specific overall ranking for club head speed (Table B1) and ball carry distance (Table B2).
Table B2. Ranking of feature sets for ball carry distance prediction.
Table B2. Ranking of feature sets for ball carry distance prediction.
#Feature SetRanking for Detection TaskRanking for Regression TaskOverall Ranking
Ind.Imp.Avg.RankInd.Imp.Avg.RankAvg.Rank
1RightChanStats192220.52115610.5815.5016
2LeftChanStats22111.55295.548.502
3CorrDelta20814.012221418.02016.0018
4CorrRatio122217.017211015.51616.2519
5StatsCorrDelta212221.5221116.0513.7514
6StatsCorrRatio17712.0616410.0711.007
7PeakP1162219.02091210.5814.7515
8PeakP1Sync52012.5918712.51012.5011
9PeakP1AbsSync22212.06624.028.001
10PeakP262113.511433.518.502
11PeakP2Sync152218.5191258.5613.5012
12PeakP2AbsSync1126.53201115.51611.007
13PeakDelta32212.5972214.51413.5012
14PeakRatio454.5282215.0159.756
15PeakSyncDelta1438.5452213.51211.007
16PeakSyncSum92215.515172219.52117.5021
17PeakSyncProd132217.518141313.51215.5016
18PeakSyncDiv72214.513184.539.505
19PeakAbsSyncDelta142.51102216.0189.254
20PeakAbsSyncSum82215.014132217.51916.2519
21PeakAbsSyncProd102216.016192220.52218.2522
22PeakAbsSyncDiv18612.0632212.51012.2510
Table B3. Ultimate ranking of feature sets for the prediction of golf shot effectiveness.
Table B3. Ultimate ranking of feature sets for the prediction of golf shot effectiveness.
#Feature SetClubDistAvg.Rank
1RightChanStats4.7515.5010.132
2LeftChanStats10.258.509.381
3CorrDelta8.0016.0012.008
4CorrRatio8.5016.2512.3810
5StatsCorrDelta7.0013.7510.384
6StatsCorrRatio13.5011.0012.259
7PeakP115.2514.7515.0018
8PeakP1Sync15.7512.5014.1314
9PeakP1AbsSync12.758.0010.384
10PeakP212.008.5010.253
11PeakP2Sync17.5013.5015.5020
12PeakP2AbsSync18.0011.0014.5015
13PeakDelta16.5013.5015.0018
14PeakRatio17.009.7513.3812
15PeakSyncDelta11.2511.0011.137
16PeakSyncSum19.5017.5018.5022
17PeakSyncProd13.7515.5014.6316
18PeakSyncDiv15.759.5012.6311
19PeakAbsSyncDelta18.259.2513.7513
20PeakAbsSyncSum13.2516.2514.7517
21PeakAbsSyncProd14.0018.2516.1321
22PeakAbsSyncDiv9.5012.2510.886
Ultimate ranking (denoted as Rank in Table B3) of feature sets is estimated after combining the target-specific overall rank averages (replicated in the Club and Dist columns of Table B3) into the ultimate average (denoted as Avg. in Table B3).

References

  1. Hellström, J.; Nilsson, J.; Isberg, L. Drive for dough. PGA Tour golfers’ tee shot functional accuracy, distance and hole score. J. Sports Sci. 2014, 32, 462–469. [Google Scholar] [CrossRef] [PubMed]
  2. Joyce, C.; Burnett, A.; Cochrane, J.; Ball, K. Three-dimensional trunk kinematics in golf: Between-club differences and relationships to clubhead speed. Sports Biomech. 2013, 12, 108–120. [Google Scholar] [CrossRef] [PubMed]
  3. Miura, K.; Sato, F. The Initial Trajectory Plane After Ball Impact. In Science and Golf III: Poceedings of the 3rd World Scientific Congress of Golf; Farrally, M.R., Cochran, A.J., Eds.; Leeds: Human Kinetics Publishers: Champaign, IL, USA, 1998; pp. 535–542. [Google Scholar]
  4. Neal, R.J.; Lumsden, R.G.; Holland, M.; Mason, B. Body Segment Sequencing and Timing in Golf. Int. J. Sports Sci. Coach. 2007, 2, 25–36. [Google Scholar] [CrossRef]
  5. Smith, A.; Roberts, J.; Wallace, E.; Kong, P.W.; Forrester, S.; Mackenzie, S.; Robertson, S. Golf Coaches’ Perceptions of Key Technical Swing Parameters Compared to Biomechanical Literature. Int. J. Sports Sci. Coach. 2015, 10, 739–756. [Google Scholar] [CrossRef] [Green Version]
  6. McHardy, A.J.; Pollard, H.P. Golf and upper limb injuries: A summary and review of the literature. Chiropr. Osteopat. 2005, 13, 7–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Chan, F.H.Y.; Yang, Y.S.; Lam, F.K.; Zhang, Y.T.; Parker, P.A. Fuzzy EMG classification for prosthesis control. IEEE Trans. Rehabil. Eng. 2000, 8, 305–311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Sebelius, F.C.; Rosén, B.N.; Lundborg, G.N. Refined Myoelectric Control in Below-Elbow Amputees Using Artificial Neural Networks and a Data Glove. J. Hand Surg. 2005, 30, 780–789. [Google Scholar] [CrossRef] [PubMed]
  9. Al-Assaf, Y. Surface Myoelectric Signal Analysis: Dynamic Approaches for Change Detection and Classification. IEEE Trans. Biomed. Eng. 2006, 53, 2248–2256. [Google Scholar] [CrossRef] [PubMed]
  10. Liu, Y.H.; Huan, H.P.; Weng, C.H. Recognition of Electromyographic Signals Using Cascaded Kernel Learning Machine. IEEE/ASME Trans. Mechatron. 2007, 12, 253–264. [Google Scholar] [CrossRef]
  11. Khezri, M.; Jahed, M. An exploratory study to design a novel hand movement identification system. Comput. Biol. Med. 2009, 39, 433–442. [Google Scholar] [CrossRef] [PubMed]
  12. Castellini, C.; Gruppioni, E.; Davalli, A.; Sandini, G. Fine detection of grasp force and posture by amputees via surface electromyography. J. Physiol. 2009, 103, 255–262. [Google Scholar] [CrossRef] [PubMed]
  13. Bu, N.; Okamoto, M.; Tsuji, T. A Hybrid Motion Classification Approach for EMG-Based Human-Robot Interfaces Using Bayesian and Neural Networks. IEEE Trans. Robot. 2005, 25, 502–511. [Google Scholar] [CrossRef]
  14. Khushaba, R.N.; Al-Ani, A. Orthogonal Fuzzy Neighborhood Discriminant Analysis for Multifunction Myoelectric Hand Control. IEEE Trans. Biomed. Eng. 2010, 57, 1410–1419. [Google Scholar] [CrossRef] [PubMed]
  15. Kim, K.S.; Choi, H.H.; Moon, C.S.; Mun, C.W. Comparison of k-nearest neighbor, quadratic discriminant and linear discriminant analysis in classification of electromyogram signals based on the wrist-motion directions. Curr. Appl. Phys. 2011, 11, 740–745. [Google Scholar] [CrossRef]
  16. Khezri, M.; Jahed, M. A Neuro-Fuzzy Inference System for sEMG-Based Identification of Hand Motion Commands. IEEE Trans. Ind. Electron. 2011, 58, 1952–1960. [Google Scholar] [CrossRef]
  17. Naik, G.R.; Kumar, D.K. Identification of Hand and Finger Movements Using Multi Run ICA of Surface Electromyogram. J. Med. Syst. 2012, 36, 841–851. [Google Scholar] [CrossRef] [PubMed]
  18. Sueaseenak, D.; Chanwimalueang, T.; Pintavirooj, C.; Sangworasil, M. An accurate forearm EMG signal classification method using two-channel electrode. IEEJ Trans. Electr. Electron. Eng. 2013, 8, 328–338. [Google Scholar] [CrossRef]
  19. Jiang, N.; Vujaklija, I.; Rehbaum, H.; Graimann, B.; Farina, D. Is Accurate Mapping of EMG Signals on Kinematics Needed for Precise Online Myoelectric Control? IEEE Trans. Neural Syst. Rehabil. Eng. 2014, 22, 549–558. [Google Scholar] [CrossRef] [PubMed]
  20. MacIsaac, D.T.; Parker, P.A.; Englehart, K.B.; Rogers, D.R. Fatigue Estimation with a Multivariable Myoelectric Mapping Function. IEEE Trans. Biomed. Eng. 2006, 53, 694–700. [Google Scholar] [CrossRef] [PubMed]
  21. Cifrek, M.; Medved, V.; Tonkovi´, S.; Ostoji´, S. Surface EMG based muscle fatigue evaluation in biomechanics. Clin. Biomechan. 2009, 24, 327–340. [Google Scholar] [CrossRef] [PubMed]
  22. Xie, H.B.; Guo, J.Y.; Zheng, Y.P. Fuzzy Approximate Entropy Analysis of Chaotic and Natural Complex Systems: Detecting Muscle Fatigue Using Electromyography Signals. Ann. Biomed. Eng. 2010, 38, 1483–1496. [Google Scholar] [CrossRef] [PubMed]
  23. Al-Mulla, M.R.; Sepulveda, F.; Colley, M. Evolved pseudo-wavelet function to optimally decompose sEMG for automated classification of localized muscle fatigue. Med. Eng. Phys. 2011, 33, 411–417. [Google Scholar] [CrossRef] [PubMed]
  24. Al-Mulla, M.R.; Sepulveda, F.; Colley, M. A Review of Non-Invasive Techniques to Detect and Predict Localised Muscle Fatigue. Sensors 2011, 11, 3545–3594. [Google Scholar] [CrossRef] [PubMed]
  25. Moshou, D.; Hostens, I.; Papaioannou, G.; Ramon, H. Dynamic muscle fatigue detection using self-organizing maps. Appl. Soft Comput. 2005, 5, 391–398. [Google Scholar] [CrossRef]
  26. Knaflitz, M.; Molinari, F. Assessment of muscle fatigue during biking. IEEE Trans. Neural Syst. Rehabil. Eng. 2003, 11, 17–23. [Google Scholar] [CrossRef] [PubMed]
  27. Dingwell, J.B.; Joubert, J.E.; Diefenthaeler, F.; Trinity, J.D. Changes in muscle activity and kinematics of highly trained cyclists during fatigue. IEEE Trans. Biomed. Eng. 2008, 55, 2666–2674. [Google Scholar] [CrossRef] [PubMed]
  28. Tenan, M.S.; McMurray, R.G.; Blackburn, T.B.; McGrath, M.; Leppert, K. The relationship between blood potassium, blood lactate, and electromyography signals related to fatigue in a progressive cycling exercise test. J. Electromyogr. Kinesiol. 2011, 21, 25–32. [Google Scholar] [CrossRef] [PubMed]
  29. Kumar, D.K.; Arjunan, S.P.; Naik, G.R. Measuring Increase in Synchronization to Identify Muscle Endurance Limit. IEEE Trans. Neural Syst. Rehabil. Eng. 2011, 19, 578–587. [Google Scholar] [CrossRef] [PubMed]
  30. Katsis, C.D.; Exarchos, T.P.; Papaloukas, C.; Goletsis, Y.; Fotiadis, D.I.; Sarmas, I. A two-stage method for MUAP classification based on EMG decomposition. Comput. Biol. Med. 2007, 37, 1232–1240. [Google Scholar] [CrossRef] [PubMed]
  31. Rasheed, S.; Stashuk, D.; Kamel, M. Fusion of multiple classifiers for motor unit potential sorting. Biomed. Signal Process. Control 2008, 3, 229–243. [Google Scholar] [CrossRef]
  32. Koçer, S. Classification of EMG Signals Using Neuro-Fuzzy System and Diagnosis of Neuromuscular Diseases. J. Med. Syst. 2010, 34, 321–329. [Google Scholar] [CrossRef]
  33. Naik, G.R.; Selvan, S.E.; Nguyen, H.T. Single-Channel EMG Classification with Ensemble-Empirical-Mode-Decomposition-Based ICA for Diagnosing Neuromuscular Disorders. IEEE Trans. Neural Syst. Rehabil. Eng. 2015, 99. [Google Scholar] [CrossRef] [PubMed]
  34. Nair, S.S.; French, R.M.; Laroche, D.; Thomas, E. The Application of Machine Learning Algorithms to the Analysis of Electromyographic Patterns from Arthritic Patients. IEEE Trans. Neural Syst. Rehabil. Eng. 2010, 18, 174–184. [Google Scholar] [CrossRef] [PubMed]
  35. Mitchell, P.; Krotish, D.; Shin, Y.J.; Hirth, V. Cross Time-Frequency Analysis of Gastrocnemius Electromyographic Signals in Hypertensive and Nonhypertensive Subjects. EURASIP J. Adv. Signal Process. 2010, 15. [Google Scholar] [CrossRef]
  36. Arslan, Y.Z.; Adli, M.A.; Akan, A.; Baslo, M.B. Prediction of externally applied forces to human hands using frequency content of surface EMG signals. Comput. Methods Programs Biomed. 2010, 98, 36–44. [Google Scholar] [CrossRef] [PubMed]
  37. Hou, Y.; Zurada, J.M.; Karwowski, W.; Marras, W.S.; Davis, K. Estimation of the Dynamic Spinal Forces Using a Recurrent Fuzzy Neural Network. IEEE Trans. Syst. Man Cybern. Part B Cybern. 2007, 37, 100–109. [Google Scholar] [CrossRef]
  38. Hahn, M.E. Feasibility of estimating isokinetic knee torque using a neural network model. J. Biomech. 2007, 40, 1107–1114. [Google Scholar] [CrossRef] [PubMed]
  39. Naik, G.R.; Nguyen, H.T. Nonnegative Matrix Factorization for the Identification of EMG Finger Movements: Evaluation Using Matrix Analysis. IEEE J. Biomed. Health Inform. 2015, 19, 478–485. [Google Scholar] [CrossRef] [PubMed]
  40. Barniv, Y.; Aguilar, M.; Hasanbelliu, E. Using EMG to anticipate head motion for virtual-environment applications. IEEE Trans. Biomed. Eng. 2005, 52, 1078–1093. [Google Scholar] [CrossRef] [PubMed]
  41. Shrirao, N.A.; Reddy, N.P.; Kosuri, D.R. Neural network committees for finger joint angle estimation from surface EMG signals. BioMed. Eng. OnLine 2009, 8, 1–11. [Google Scholar] [CrossRef] [PubMed]
  42. Chiang, J.; Wang, Z.J.; McKeown, M.J. A Windowed Eigenspectrum Method for Multivariate sEMG Classification During Reaching Movements. IEEE Signal Process. Lett. 2008, 15, 293–296. [Google Scholar] [CrossRef]
  43. Arjunan, S.P.; Kumar, D.K.; Naik, G.R. Independence Between Two Channels of Surface Electromyogram Signal to Measure the Loss of Motor Units. Meas. Sci. Rev. 2015, 15, 152–155. [Google Scholar] [CrossRef]
  44. Choi, C.; Micera, S.; Carpaneto, J.; Kim, J. Development and Quantitative Performance Evaluation of a Noninvasive EMG Computer Interface. IEEE Trans. Biomed. Eng. 2009, 56, 188–191. [Google Scholar] [CrossRef] [PubMed]
  45. Naik, G.R.; Kumar, D.K.; Palaniswami, M. Signal processing evaluation of myoelectric sensor placement in low-level gestures: Sensitivity analysis using independent component analysis. Expert Syst. 2014, 31, 91–99. [Google Scholar] [CrossRef]
  46. Kurta, M.B.; Sezgin, N.; Akin, M.; Kirbas, G.; Bayram, M. The ANN-based computing of drowsy level. Expert Syst. Appl. 2009, 36, 2534–2542. [Google Scholar] [CrossRef]
  47. Palmes, P.; Ang, W.T.; Widjaja, F.; Tan, L.C.S.; Au, W.L. Pattern mining of multichannel sEMG for tremor classification. IEEE Trans. Biomed. Eng. 2010, 57, 2795–2805. [Google Scholar] [CrossRef] [PubMed]
  48. McHardy, A.J.; Pollard, H.P. Muscle activity during the golf swing. Br. J. Sports Med. 2005, 39, 799–804. [Google Scholar] [CrossRef] [PubMed]
  49. Marta, S.; Silva, L.; Castro, M.A.; Pezarat-Correia, P.; Cabri, J. Electromyography variables during the golf swing: A literature review. J. Electromyogr. Kinesiol. 2012, 22, 803–813. [Google Scholar] [CrossRef] [PubMed]
  50. Silva, L.; Marta, S.; Vaz, J.; Fernandes, O.; Castro, M.A.; Pezarat-Correia, P. Trunk muscle activation during golf swing: Baseline and threshold. J. Electromyogr. Kinesiol. 2013, 23, 1174–1182. [Google Scholar] [CrossRef] [PubMed]
  51. Myers, J.B.; Lephart, S.M.; Tsai, Y.S.; Sell, T.C.; Smoliga, J.M.; Jolly, J.T. The role of upper torso and pelvis rotation in driving performance during the golf swing. J. Sports Sci. 2008, 26, 181–188. [Google Scholar] [CrossRef] [PubMed]
  52. Chu, Y.; Sell, T.C.; Lephart, S.M. The relationship between biomechanical variables and driving performance during the golf swing. J. Sports Sci. 2010, 28, 1251–1259. [Google Scholar] [CrossRef] [PubMed]
  53. Okuda, I.; Armstrong, C.W.; Tsunezumi, H.; Yoshiike, H. Biomechanical analysis of a professional golfer’s swing: Hidemichi Tanaka. In Science and Golf IV: Proceedings of the World Scientific Congress of Golf; Thain, E., Ed.; Routledge: London, UK, 2002; pp. 8–16. [Google Scholar]
  54. Healy, A. Identification of the Biomechanical Performance Determining Factors of the 5 Iron Golf Swing When Hitting for Maximum Distance. Master’s Thesis, Dublin City University, Dublin, Ireland, 2009. [Google Scholar]
  55. Bartlett, R. Introduction to Sports Biomechanics: Analysing Human Movement Patterns, 2nd ed.; Taylor & Francis: London, UK, 2007. [Google Scholar]
  56. Luca, C.J.D.; Gilmore, L.D.; Kuznetsov, M.; Roy, S.H. Filtering the surface EMG signal: Movement artifact and baseline noise contamination. J. Biomech. 2010, 43, 1573–1579. [Google Scholar] [CrossRef] [PubMed]
  57. O’Haver, T.C. A Pragmatic Introduction to Signal Processing with Applications in Scientific Measurement, 2015. Software for peak finding and measurement: iPeak. Available online: http://terpconnect.umd.edu/%7Etoh/spectrum/PeakFindingandMeasurement.htm#ipeak (accessed on 19 April 2016).
  58. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  59. Jaiantilal, A. Random Forest (Regression, Classification and Clustering) implementation for Matlab (and Standalone), 2012. Available online: http://code.google.com/p/randomforest-matlab/ (accessed on 19 April 2016).
  60. Verikas, A.; Gelzinis, A.; Bacauskiene, M. Mining data with random forests: A survey and results of new tests. Pattern Recognit. 2011, 44, 330–349. [Google Scholar] [CrossRef]
  61. Ho, T.K. Data Complexity Analysis for Classifier Combination. In Multiple Classifier Systems; Springer Berlin Heidelberg: Berlin/Heidelberg, Germany; pp. 53–67.
  62. Kuncheva, L.I. Combining Pattern Classifiers; Wiley-Blackwell: Oxford, UK, 2004. [Google Scholar]
  63. Pampara, G.; Engelbrecht, A.P.; Franken, N. Binary Differential Evolution. In Proceedings of the IEEE Congress on Evolutionary Computati on (CEC), Vancouver, BC, Canada, 16–21 July 2006; pp. 1873–1879.
  64. Mukhopadhyay, A.; De, M.; Maulik, U. Selection of GO-Based Semantic Similarity Measures through AMDE for Predicting Protein-Protein Interactions. In Proceedings of the 2nd International Conference on Swarm, Evolutionary, and Memetic Computing (SEMCCO) – Part II, Visakhapatnam, India, 19–21 December 2011; Lecture Notes in Computer Science. Panigrahi, B., Suganthan, P., Das, S., Satapathy, S., Eds.; Springer-Verlag: Berlin/Heidelberg, Germany, 2011; Volume 7077, pp. 55–62. [Google Scholar]
  65. Vaiciukynas, E.; Verikas, A.; Gelzinis, A.; Bacauskiene, M.; Kons, Z.; Satt, A.; Hoory, R. Fusion of voice signal information for detection of mild laryngeal pathology. Appl. Soft Comput. 2014, 18, 91–103. [Google Scholar] [CrossRef]
  66. Sáenz-Lechón, N.; Godino-Llorente, J.I.; Osma-Ruiz, V.; Gómez-Vilda, P. Methodological issues in the development of automatic systems for voice pathology detection. Biomed. Signal Process. Control 2006, 1, 120–128. [Google Scholar] [CrossRef]
  67. Brümmer, N.; de Villiers, E. The BOSARIS Toolkit: Theory, Algorithms and Code for Surviving the New DCF. arXiv 2013, 1304, 1–23. [Google Scholar]
Figure 1. Participant body with eight-channel EMG recorder electrodes attached. Electrodes on extensor digitorum communis (EDC) arm muscles can be seen in the front side view (Left); electrodes on flexor carpi radialis (FCR) arm muscles, as well as electrodes on rhomboideus and trapezius muscles can be seen in the back side view (Right). Mapping between a channel number and a muscle: (1) right FCR; (2) right EDC; (3) left FCR; (4) left EDC; (5) right rhomboideus; (6) right trapezius; (7) left rhomboideus; and (8) left trapezius.
Figure 1. Participant body with eight-channel EMG recorder electrodes attached. Electrodes on extensor digitorum communis (EDC) arm muscles can be seen in the front side view (Left); electrodes on flexor carpi radialis (FCR) arm muscles, as well as electrodes on rhomboideus and trapezius muscles can be seen in the back side view (Right). Mapping between a channel number and a muscle: (1) right FCR; (2) right EDC; (3) left FCR; (4) left EDC; (5) right rhomboideus; (6) right trapezius; (7) left rhomboideus; and (8) left trapezius.
Sensors 16 00592 g001
Figure 2. Scatterplot, showing the effectiveness for men (Left) and women (Right) players. Notes: grid lines correspond to sex-specific averages; shots are colored according to the player’s overall effectiveness, where a lighter color means a more effective player, based on average club head speed and average ball carry distance through all of his/her shots.
Figure 2. Scatterplot, showing the effectiveness for men (Left) and women (Right) players. Notes: grid lines correspond to sex-specific averages; shots are colored according to the player’s overall effectiveness, where a lighter color means a more effective player, based on average club head speed and average ball carry distance through all of his/her shots.
Sensors 16 00592 g002
Figure 3. EMG signals recorded from eight muscles of golf players G2 (Top) and G6 (Bottom). Note: position in samples on the X (horizontal) axis. Channel numbers correspond to the following muscles: (1) right FCR; (2) right EDC; (3) left FCR; (4) left EDC; (5) right rhomboideus; (6) right trapezius; (7) left rhomboideus; and (8) left trapezius.
Figure 3. EMG signals recorded from eight muscles of golf players G2 (Top) and G6 (Bottom). Note: position in samples on the X (horizontal) axis. Channel numbers correspond to the following muscles: (1) right FCR; (2) right EDC; (3) left FCR; (4) left EDC; (5) right rhomboideus; (6) right trapezius; (7) left rhomboideus; and (8) left trapezius.
Sensors 16 00592 g003
Figure 4. EMG signal pre-processing: raw signal (Top) and normalized rectified signal (Bottom) with the result of Butterworth filtering (black curve). Note: the position is in seconds on the X (horizontal) axis.
Figure 4. EMG signal pre-processing: raw signal (Top) and normalized rectified signal (Bottom) with the result of Butterworth filtering (black curve). Note: the position is in seconds on the X (horizontal) axis.
Sensors 16 00592 g004
Figure 5. Finding peaks in the filtered EMG signal. Interface window of the ipeak function: entire pre-processed signal (Bottom) with 17 peaks detected and the zoomed-in portion of the signal (Top) containing the two highest peaks (# 10 and # 11). Note: the position is in seconds on the X (horizontal) axis.
Figure 5. Finding peaks in the filtered EMG signal. Interface window of the ipeak function: entire pre-processed signal (Bottom) with 17 peaks detected and the zoomed-in portion of the signal (Top) containing the two highest peaks (# 10 and # 11). Note: the position is in seconds on the X (horizontal) axis.
Sensors 16 00592 g005
Figure 6. Architecture of the random forest model.
Figure 6. Architecture of the random forest model.
Sensors 16 00592 g006
Figure 7. Activation profiles of eight muscles during golf swing for 16 different golf players. Note: channel 8 (left trapezius) was used as the reference for measuring time differences for remaining channels.
Figure 7. Activation profiles of eight muscles during golf swing for 16 different golf players. Note: channel 8 (left trapezius) was used as the reference for measuring time differences for remaining channels.
Sensors 16 00592 g007
Figure 8. Decision-level fusion for detection task: evaluation by the detection error trade-off (DET) (Left) and ROC (Right) curves.
Figure 8. Decision-level fusion for detection task: evaluation by the detection error trade-off (DET) (Left) and ROC (Right) curves.
Sensors 16 00592 g008
Figure 9. RF permutation-based variable importance in decision-level fusion for detection (Left) and regression (Right) tasks.
Figure 9. RF permutation-based variable importance in decision-level fusion for detection (Left) and regression (Right) tasks.
Sensors 16 00592 g009
Figure 10. Decision-level fusion for regression task: results for club head speed (Left) and ball carry distance (Right) prediction. Note: the diagonal dotted line corresponds to an ideal prediction.
Figure 10. Decision-level fusion for regression task: results for club head speed (Left) and ball carry distance (Right) prediction. Note: the diagonal dotted line corresponds to an ideal prediction.
Sensors 16 00592 g010
Table 1. Descriptions of feature sets. Note: the size count does not include sex.
Table 1. Descriptions of feature sets. Note: the size count does not include sex.
#Feature SetSizeShort Description
Channel statistics & correlations
1RightChanStats49Duration of swing and statistics for the right-side channels
2LeftChanStats49Duration of swing and statistics for the left-side channels
3CorrDelta28Difference between Pearson and Spearman correlations
4CorrRatio28Ratio of Pearson correlation to Spearman correlation
5StatsCorrDelta40CorrDelta features and statistics of differences
6StatsCorrRatio40CorrRatio features and statistics of ratios
Properties of the 1st peak ( P 1 )
7PeakP140The main properties of the 1st peak
8PeakP1Sync40Difference between positions and statistics
9PeakP1AbsSync40Absolute difference between positions and statistics
Properties of the 2nd peak ( P 2 )
10PeakP240The main properties of the 2nd peak
11PeakP2Sync40Difference between positions and statistics
12PeakP2AbsSync40Absolute difference between positions and statistics
Properties of both peaks ( P 1 & P 2 )
13PeakDelta40Difference between the main properties of peaks
14PeakRatio40Ratio between the main properties of peaks
15PeakSyncDelta40Differences between positions subtracted and statistics
16PeakSyncSum40Differences between positions summed and statistics
17PeakSyncProd40Differences between positions multiplied and statistics
18PeakSyncDiv40Differences between positions divided and statistics
19PeakAbsSyncDelta40Absolute differences between positions subtracted and statistics
20PeakAbsSyncSum40Absolute differences between positions summed and statistics
21PeakAbsSyncProd40Absolute differences between positions multiplied and statistics
22PeakAbsSyncDiv40Absolute differences between positions divided and statistics
Table 2. RF OOB performance for the detection task.
Table 2. RF OOB performance for the detection task.
#Feature SetClub Head SpeedBall Carry Distance
CllrEER, %AUCRankCllrEER, %AUCRank
1RightChanStats0.95435.710.70241.08847.580.54619
2LeftChanStats1.03540.800.626151.14450.000.50022
3CorrDelta0.91432.040.74211.09349.000.51520
4CorrRatio0.94935.180.71631.04345.590.55612
5StatsCorrDelta0.91831.070.75621.10449.640.50721
6StatsCorrRatio0.97637.060.68261.07148.230.53417
7PeakP11.03646.060.557161.07049.710.50516
8PeakP1Sync1.08648.190.535221.00440.640.6205
9PeakP1AbsSync1.01841.100.611110.98141.010.6382
10PeakP20.96133.310.70851.02144.640.5916
11PeakP2Sync0.99141.170.63771.06343.430.58915
12PeakP2AbsSync1.00441.230.62081.03541.100.64111
13PeakDelta1.04045.070.559170.99736.990.6503
14PeakRatio1.06950.000.500210.99938.200.6474
15PeakSyncDelta1.01642.530.61691.06246.430.56114
16PeakSyncSum1.05344.500.577191.03342.590.6099
17PeakSyncProd1.03443.680.590141.05145.620.57413
18PeakSyncDiv1.01743.800.599101.02242.830.6017
19PeakAbsSyncDelta1.05047.820.535180.97539.050.6611
20PeakAbsSyncSum1.02941.950.617131.03043.210.5918
21PeakAbsSyncProd1.01937.800.648121.03345.190.58710
22PeakAbsSyncDiv1.05749.030.517201.07947.540.54718
Table 3. Decision-level fusion for detection task: summary of results. Note: 95CI stands for 95% confidence interval (mean and corresponding 95% CI are obtained after repeating the task 100 times).
Table 3. Decision-level fusion for detection task: summary of results. Note: 95CI stands for 95% confidence interval (mean and corresponding 95% CI are obtained after repeating the task 100 times).
Target AttributeCllr ± 95CIEER ± 95CI, %AUC ± 95CI
Club head speed0.832 ± 0.00124.32 ± 0.180.825 ± 0.001
Ball carry distance0.895 ± 0.00031.71 ± 0.140.779 ± 0.002
Table 4. RF OOB performance for the regression task.
Table 4. RF OOB performance for the regression task.
#Feature SetClub Head SpeedBall Carry Distance
RMSEMAPE, %Correl.RankRMSEMAPE, %Correl.Rank
1RightChanStats8.395.070.580716.9510.410.51915
2LeftChanStats8.435.170.5751014.728.990.6702
3CorrDelta9.475.780.3942119.6912.060.11722
4CorrRatio9.785.880.3152219.5211.990.17221
5StatsCorrDelta7.854.820.648216.6910.080.53911
6StatsCorrRatio8.114.960.616417.0110.180.51316
7PeakP17.734.700.661116.4310.380.5599
8PeakP1Sync8.805.290.5191617.3210.900.48618
9PeakP1AbsSync8.685.470.5391416.069.750.5866
10PeakP28.675.280.5401315.569.300.6194
11PeakP2Sync9.095.540.4701916.7110.200.53812
12PeakP2AbsSync9.125.410.4662017.9110.940.42820
13PeakDelta8.144.920.612516.1710.050.5787
14PeakRatio8.084.980.621316.199.950.5778
15PeakSyncDelta8.254.850.599615.9110.090.5965
16PeakSyncSum8.805.240.5201517.1310.550.50317
17PeakSyncProd8.935.330.4991716.9110.410.52214
18PeakSyncDiv8.425.020.576914.298.720.6931
19PeakAbsSyncDelta8.515.140.5631116.5010.180.55410
20PeakAbsSyncSum8.655.200.5421216.8110.450.53013
21PeakAbsSyncProd8.985.410.4911817.7210.890.44819
22PeakAbsSyncDiv8.415.050.578815.479.600.6253
Table 5. Decision-level fusion for regression task: summary of results. Note: 95CI stands for 95% confidence interval (mean and corresponding 95% CI are obtained after repeating the task 100 times).
Table 5. Decision-level fusion for regression task: summary of results. Note: 95CI stands for 95% confidence interval (mean and corresponding 95% CI are obtained after repeating the task 100 times).
Target AttributeRMSE ± 95CIMAPE ± 95CI, %Correl. ± 95CI
Club head speed5.70 ± 0.013.22 ± 0.010.813 ± 0.001
Ball carry distance10.70 ± 0.046.38 ± 0.040.842 ± 0.001

Share and Cite

MDPI and ACS Style

Verikas, A.; Vaiciukynas, E.; Gelzinis, A.; Parker, J.; Olsson, M.C. Electromyographic Patterns during Golf Swing: Activation Sequence Profiling and Prediction of Shot Effectiveness. Sensors 2016, 16, 592. https://doi.org/10.3390/s16040592

AMA Style

Verikas A, Vaiciukynas E, Gelzinis A, Parker J, Olsson MC. Electromyographic Patterns during Golf Swing: Activation Sequence Profiling and Prediction of Shot Effectiveness. Sensors. 2016; 16(4):592. https://doi.org/10.3390/s16040592

Chicago/Turabian Style

Verikas, Antanas, Evaldas Vaiciukynas, Adas Gelzinis, James Parker, and M. Charlotte Olsson. 2016. "Electromyographic Patterns during Golf Swing: Activation Sequence Profiling and Prediction of Shot Effectiveness" Sensors 16, no. 4: 592. https://doi.org/10.3390/s16040592

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop