Minimally Invasive Surgery
Volume 2016 (2016), Article ID 5624630, 8 pages
http://dx.doi.org/10.1155/2016/5624630

Research Article

Prediction of Muscle Fatigue during Minimally Invasive Surgery Using Recurrence Quantification Analysis

Ali Keshavarz Panahi and Sohyung Cho

Industrial Engineering, Southern Illinois University Edwardsville, Edwardsville, IL 62026, USA

Received 29 October 2015; Revised 29 February 2016; Accepted 12 April 2016

Academic Editor: Othmar Schöb

Copyright © 2016 Ali Keshavarz Panahi and Sohyung Cho. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Due to its inherent complexity such as limited work volume and degree of freedom, minimally invasive surgery (MIS) is ergonomically challenging to surgeons compared to traditional open surgery. Specifically, MIS can expose performing surgeons to excessive ergonomic risks including muscle fatigue that may lead to critical errors in surgical procedures. Therefore, detecting the vulnerable muscles and time-to-fatigue during MIS is of great importance in order to prevent these errors. The main goal of this study is to propose and test a novel measure that can be efficiently used to detect muscle fatigue. In this study, surface electromyography was used to record muscle activations of five subjects while they performed fifteen various laparoscopic operations. The muscle activation data was then reconstructed using recurrence quantification analysis (RQA) to detect possible signs of muscle fatigue on eight muscle groups (bicep, triceps, deltoid, and trapezius). The results showed that RQA detects the fatigue sign on bilateral trapezius at 47.5 minutes (average) and bilateral deltoid at 57.5 minutes after the start of operations. No sign of fatigue was detected for bicep and triceps muscles of any subject. According to the results, the proposed novel measure can be efficiently used to detect muscle fatigue and eventually improve the quality of MIS procedures with reducing errors that may result from overlooked muscle fatigue.

1. Introduction

Since minimally invasive surgery (MIS) is performed using an endoscope and several thin instruments through small incisions made in patients, it provides patients with some advantages like shorter recovery time, less postoperative pain, earlier resumption of normal activity, and cost savings. However, it is more ergonomically challenging to performing surgeons due to its inherent complexity such as limited work volume and degree of freedom. Specifically, such complexity induces abnormal movements of arm and shoulder while holding certain body postures, for example, head and torso, for prolonged time [1–3]. Importantly, this complexity often exposes performing surgeons to increased ergonomic risks including muscle fatigue that can result in critical errors during surgical procedures [4–6]. The vulnerability to ergonomic risk is well confirmed in a literature stating that “performing laparoscopic surgery is significantly more stressful for the surgeon than open surgery” [7]. In addition, it has been reported that MIS carries more complications than open surgery [8, 9]. Therefore, it is important to quantitatively measure ergonomic risks on performing surgeons, particularly muscle fatigue during MIS procedures. Note that timely intervention with information about muscle fatigue can ensure improved quality of MIS procedures.

Some existing study for fatigue analysis has focused on quantitative measures and their applications for isometric or isotonic contractions during relatively short bouts of high-force activities. Among these studies, frequency banding analysis predicted muscle fatigue when lower frequency bands increase and higher frequency bands decrease [10]. Standard discomfort analysis was also used to detect the level of muscle fatigue on subjects, which may lack objective analysis [11]. In addition, heart rate [12] and tissue oxygen saturation [13] were used as indicators of fatigue development. Importantly, power spectral density (PSD) of transformed EMG data through Fast Fourier Transform approach was widely used to detect muscle fatigue [14–16]. In this approach, frequency and amplitude changes are considered as a result of a reduction in conduction velocity in the muscle fibers and larger motor unit synchronization. Figure 1 schematically depicts these changes in various stages of the muscle fatigue.

Figure 1: Schematic of changes in EMG data during muscle fatigue.

As MIS operations can be considered as prolonged light muscle activation, PSD based measure, which is mainly for static activities, may not be efficiently used to predict muscle fatigue during MIS operations. Specifically, when PSD based measures are used for light muscle activations, conflicting results and complicated relationship between subjective and objective fatigue have been reported [17, 18]. Such conflicts result from stationary data requirement for PSD based measures to be successfully used. In fact, EMG data collected from light muscle activations are nonstationary due to changing distance between muscle and EMG sensor, as well as changing muscle length. From the literature, it is noted that fatigue mechanism during these activities is not well understood in comparison to high-force (or intense) muscle activities [19]. Although there are no significant reductions in muscle force level during MIS operations, there are some physiological changes in muscles that can lead to muscle fatigue affecting the capacity of muscles and the performance of subjects [20, 21]. Uhrich et al. [22] assessed the muscle fatigue during a relatively short time during simulated laparoscopic surgery. In their study, the effects of fatigue, monitor placement, and surgical experience have been compared. After obtaining the EMG activity and muscular discomfort scores before and after a fatigue session, it was found that the EMG data and discomfort scores demonstrated a fatigue response in several muscle groups. They found minimal differences between the two monitor positions and less muscle activity and discomfort in the attending surgeons. In another study, Slack et al. [23] studied the effect of operating time on surgeon’s muscular fatigue. The length of the operations varied within 1–10 hr, but only one minute period before and after operations were analyzed using frequency analysis. They found that brachioradialis was being used more and fatigued faster. In addition, the muscle fatigue increases proportional to time. On the other hand, recurrence quantification analysis (RQA) has been tested on biceps brachii and shown to be more sensitive to muscle fatigue than FFT variable spectral center frequency (fC) [24]. In another study, Filligoi and Felici [25] showed that Determinism%, which is one of the RQA variables, is more effective than median frequency to detect EMG signal changes in biceps brachii muscle. For the first time, this novel data analysis method is used in our study to quantify any possible muscle fatigue in real laparoscopic surgery operations which is known as prolonged and low-force muscle activity compared to some intense tasks which involve isometric muscle contraction throughout the whole task.

The main goal of this study is to answer the following questions:

(i)Can objective manifestations of muscle fatigue be detected from EMG data during a laparoscopic surgery as prolonged light muscle activation?(ii)What is time-to-fatigue for the muscles experiencing fatigue during laparoscopic surgery?(iii)Which muscle has the highest possibility level of fatigue in laparoscopic surgery among the tested muscle groups?

2. Methods

As the first step, EMG data for fatigue analysis was collected. Specifically, surface EMG electrodes were attached to upper arm muscles of participants to collect muscle activations while performing various MIS procedures. Next, EMG data was converted into higher dimensional data using time shift and represented by recurrence plots that facilitate recurrence quantification analysis (RQA), particularly computation of determinism values. Then, moving average technique was applied for trajectories of determinism values to detect any changes as the indicator of muscle fatigue.

2.1. EMG Data Collection

Under an IRB-approved protocol, five right-hand-dominant expert laparoscopic surgeons, who have performed more than 100 laparoscopic surgeries, performed fifteen MIS procedures in a local hospital. Table 1 shows more details about those surgical operations.

Table 1: MIS operations from which EMG data was collected.

Subject
Lap #
Completion time (min)
Description

Surgeon 1
Lap 1
96
Ventral hernia repair

Lap 2
113
Cholecystectomy

Lap 3
132
Inguinal hernia repair

Surgeon 2
Lap 1
55
Sleeve gastrectomy

Lap 2
66
Sleeve gastrectomy

Lap 3
128
Sleeve gastrectomy

Surgeon 3
Lap 1
131
Heller myotomy, Dor fundoplication, and liver biopsy

Lap 2
152
Ventral hernia repair

Lap 3
171
Paraesophageal hernia repair

Surgeon 4
Lap 1
122
Sleeve gastrectomy

Lap 2
101
Inguinal hernia repair

Lap 3
88
Sleeve gastrectomy

Surgeon 5
Lap 1
71
Heller myotomy

Lap 2
116
Ventral hernia repair

Lap 3
144
Cholecystectomy

Ventral hernia repair: the hernia is repaired by mesh or sutures entered through instruments placed into small incisions in the abdomen. Cholecystectomy: to remove the gallbladder using several small incisions. Inguinal hernia repair: to repair a hernia in the abdominal wall of the groin. Sleeve gastrectomy: to remove a large portion of the stomach to help with weight loss. Heller myotomy: to treat the achalasia by cutting the muscles of the cardia to allow food and liquid to pass to the stomach. Dor fundoplication: to prevent reflux from the stomach into the esophagus by partially wrapping the stomach around the esophagus. Paraesophageal hernia repair: the diaphragm at the esophageal hiatus is closed to prevent the stomach from reherniating, and then the fundoplication is performed to keep the stomach from herniating back into the chest cavity.

On each participant, a total of eight surface EMG (sEMG) electrodes were attached to the following four bilateral muscles including bicep, triceps, deltoid, and trapezius. Note that since lower arms are scrubbed from the fingertips to the elbow, electrodes were not placed on these sterilized muscle compartments. Tricep and bicep were selected due to their high activity level during arm flexion and extension. Since the shoulder and neck are the common area of muscle fatigue, trapezius and deltoid were also included in our study. All sEMG data were collected using an 8-channel Bioradio® 150 physiologic data acquisition system (Great Lakes Neurotechnologies, Incorporated®, Cleveland, OH). After wiping skin surfaces overlying target muscle groups with rubbing alcohol and allowing the alcohol to dry, two MVAP-II® electrodes (MVAP Medical Supplies, Incorporated, Newbury Park, CA) were placed over the muscle bellies of each muscle group and connected to the positive and negative input poles for each channel. An electrode was also attached to the right elbow and connected to the ground input on the Bioradio 150 to complete the input circuit. The Biocapture® data acquisition software package (Great Lakes Neurotechnologies, Incorporated, Cleveland, OH) was used to capture and filter sEMG data. More specifically, sEMG data was sampled at frequency of 256 Hz from each channel. Digital signal processing filters were then applied to exclude the noise that are at low (<10 Hz) and high frequency (>127 Hz) signals.

2.2. Recurrence Quantification Analysis (RQA)

RQA is a time series analysis method and detects the deterministic structure of the dynamical systems. It also provides some valuable information from recurrence plots [26]. The principle of recurrence plots analysis can be summarized as follows. Let be the th point on the orbit describing a dynamical system in -dimensional space, for . The recurrence plot is an array of dots in a square, where a dot is placed at whenever is sufficiently close to . In order to obtain a recurrence plot from a time series , the following procedures are required. First, we choose an embedding dimension and construct the -dimensional orbit of by the method of time delays: if and are scalar, . Next, we choose such that the ball of radius centered at in contains a reasonable number of other points of the orbit. Finally, we plot a dot at each point for which is in the ball of radius centered at . We call this picture a recurrence plot (RP) and an example of RPs is shown in Figure 2. Note that and are, in fact, times; therefore, a recurrence plot describes natural time correlation information.

Figure 2: Example of a recurrence plot.

Recurrence plots tend to be fairly symmetric with respect to the diagonal because if is close to , then is close to . There is, however, no complete symmetry because we do not require [27]. When is more than 2 dimensions in a phase space, projection is performed. Starting from the time series , the attractor of the underlying dynamics is reconstructed in a phase space by applying the time-delay vector method by Takens [28]. The reconstructed trajectory can be expressed as a matrix where each row is a phase space vector:

where , , is the embedding dimension, and is the delay time. The recurrence plot is a tool that can be used to investigate higher dimensional dynamics through a two-dimensional binary plot of its recurrences. Any recurrence of state with state is pictured on a Boolean matrix expressed by

where are the embedded vectors, is the Heaviside step function, and is an arbitrary threshold. In the graphical representation, each nonzero entry of is marked by a black dot in the position . Since any state is recurrent with itself, the recurrent plot (RP) matrix fulfills which hence contains the diagonal line of identity (LOI).

The percentage of determinism (%DET) derived from diagonal lines and related to predictability of the system is an important parameter in RQA which quantifies the ratio of the recurrence points. The following formula is typically used to calculate %DET from a recurrence plot [26]:

where is histogram or frequency distribution of diagonal line lengths, is the length of a data series, and is predefined minimal length of a diagonal line. In this study, the above %DET for EMG data was used as a muscle fatigue indicator. Cross recurrence plot (CRP) toolbox available in MATLAB was used for RQA analysis [26]. In order to ensure optimized results, RQA parameters including embedding dimension (), time delay (), and threshold () should be selected carefully. In the CRP toolbox, one can set the parameters of optimal embedding dimension and time delay to the obtained values from false nearest neighbors (FNN) and average mutual information (AMI) function which may lead to the optimal recurrence plots as shown in Figure 3.

(a)

(b)

Figure 3: Optimal embedding dimension using FNN (a) and time delay using AMI (b).

Using these functions, embedding dimension and time-delay parameters were set to 10 and 4, respectively. To facilitate optimal threshold value, it has been suggested to consider the threshold value only a few percent of the maximum phase space diameter [29]. Therefore, the threshold value was chosen according to the 10% of the minimum value of maximum phase space diameter of the data. In this study, the minimum number of points constructing a deterministic diagonal line was set to two.

The whole duration of MIS operations summarized in Table 1 was analyzed by applying RQA for each minute of the operation with a window length of 15,360 data points. Then, %DET value was derived from RQA analysis and the results were plotted against operation time. Finally, a moving average analysis with the interval (or window size) of ten was applied to the %DET values to determine time-to-fatigue.

3. Results

Figure 4 shows %DET values of all eight muscles for surgeon 1 performing case  3 for 132 minutes as an example. It was observed from Figures 4(g) and 4(h) that moving average of %DET of bilateral trapezius muscles increased after 45–50 minutes of operation. In addition, moving average of %DET of bilateral deltoid increased after 55–60 minutes of operation as shown in Figures 4(e) and 4(f). Between deltoid and trapezius, bilateral trapezius became more deterministic with %DET close to 75% at the end of operation while bilateral deltoid was less deterministic with %DET value close to 50%. In this study, the increase in the moving average of %DET values is considered as the development of muscle fatigue. Thus, it is conjectured from these observations that trapezius became more fatigued than deltoid for surgeon 1 performing case  3. However, there were no changes in moving averages of %DET of bilateral bicep and triceps throughout the entire operation as shown in Figure 4.

(a) Left bicep

(b) Right bicep

(c) Left triceps

(d) Right triceps

(e) Left deltoid

(f) Right deltoid

(g) Left trapezius

(h) Right trapezius

Figure 4: %DET value for all the muscle groups of surgeon 1, case 3. The vertical axis is Determinism % and the horizontal axis is time (min). (a) left bicep, (b) right bicep, (c) left triceps, (d) right triceps, (e) left deltoid, (f) right deltoid, (g) left trapezius, and (h) right trapezius.

As the next step, %DET was tested for all the MIS operations and results are summarized in Table 2. DET results did not detect any changes in any of the muscles for surgeon 2 in case number 1 which had the least completion time of 55 minutes. Also, surgeon 4 in case number 2 as well as surgeon 5 in case number 3 did not experience any fatigue in their deltoid muscle. Specifically, no fatigue signs were detected on bilateral bicep and triceps in any of the subjects. Figure 5 shows the mean of moving average values for DET% of all the operations and subjects. Table 2 also shows that fatigue signs are developed at least 45 minutes after operations begin and trapezius gets fatigue sign earlier than deltoid.

Table 2: Time-to-fatigue for all the muscle groups (minutes).

Surgeon
Case
L-bicep
R-bicep
L-tricep
R-tricep
L-deltoid
R-deltoid
L-trap
R-trap

1
1
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(45–50)
(45–50)

2
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(45–50)
(45–50)

3
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(45–50)
(45–50)

2
1
N/F
N/F
N/F
N/F
N/F
N/F
N/F
N/F

2
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(45–50)
(45–50)

3
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(45–50)
(45–50)

3
1
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(45–50)
(45–50)

2
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(45–50)
(45–50)

3
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(45–50)
(45–50)

4
1
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(45–50)
(45–50)

2
N/F
N/F
N/F
N/F
N/F
N/F
(45–50)
(45–50)

3
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(55–60)
(55–60)

5
1
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(55–60)
(55–60)

2
N/F
N/F
N/F
N/F
(55–60)
(55–60)
(55–60)
(55–60)

3
N/F
N/F
N/F
N/F
N/F
N/F
(55–60)
(55–60)

Note: N/F: no fatigue.

Figure 5: Average of moving average values for DET% of all the operations and subjects for the beginning and end of operations.

4. Discussion

The moving average of %DET for bilateral trapezius deltoid increased after 45–50 and 55–60 minutes of operation, respectively. This might be caused due to the following nature of MIS procedures: (i) during MIS, surgeons should maintain their head positions at certain orientation to keep watching the monitor and this may impose significant stress on trapezius muscle, and (ii) as the surgeon operates thin and long instrument through a small incision, they often use excessive shoulder movements to overcome the limited degree of freedom and this may impose fatigue on deltoid muscle. In this study, higher %DET value implies that the EMG data are getting more deterministic and periodic. This deterministic and periodic pattern is the result of synchronization of motor units in fatiguing stage [30]. In order to supply necessary force to continue certain tasks, more motor units are involved during muscle fatigue. If almost all motor units have already been involved, motor unit synchronization takes place to continue task operation and prevent the failure. The fact that more motor units are recruited and synchronized during the course of muscle fatigue has been shown in literature [31]. For conventional occupational tasks, increase in synchronization of motor units has been successfully detected using %DET values through a computer simulation model [32].

On the other hand, shoulder disorder is considered as one of most important musculoskeletal disorders, especially in prolonged repetitive activations with high precision [33–35]. An increase in interstitial potassium concentrations in trapezius muscle can be the cause of decreased conduction velocity [36]. Also, the intra- and extracellular sodium and potassium concentration changes can be the main reason for changes in EMG data which are related to muscle fatigue [37]. Although the recovery of these metabolic changes is quick, in prolonged light muscle activation, muscles need longer recovery times in order to regain the primary force capacity [37, 38]. Because of orderly recruitment of motor units, low threshold fibers are vulnerable in prolonged muscle activation, which necessitates the importance of recovery time to avoid myalgic disorders in trapezius [39, 40]. It can be summarized from the literature review that fatigue analysis is highly difficult for prolonged light muscle activation tasks such as MIS procedures.

This study shows that trapezius was the first muscle to show fatigue sign in prolonged light muscle activation during MIS procedures. It is interesting to note that similar results were reported in a light assembly task in industrial operations [41]. While the existing study also suggested recovery/break time for workers at 90 minutes after the task begins, the results in this study suggest recovery/break time after 45–50 minutes after MIS operations begin. More likely, this difference results from intensity and precision level of two tasks: industrial and surgical. Bicep and triceps did not experience any muscle fatigue. This might be because of less muscle activation and less muscle contraction in laparoscopic surgery for these two muscles.

As muscles are vulnerable to fatigue during MIS operations which have prolonged and intense nature, detecting the muscle fatigue and estimation of time-to-fatigue are of great importance. To meet this need, this study proposed and tested a novel measure that can be efficiently used for detection of muscle fatigue and time-to-fatigue from EMG data. In the future study, correlation between objective and subjective results would be important. Here, objective results mean fatigue analysis using the proposed measure of this study and subjective results correspond to survey analysis to subjects asking whether they feel fatigue symptoms during MIS procedures. Also, it would be interesting to compare performance analysis with fatigue analysis in the future, since the final goal for muscle fatigue analysis is to detect any possible effect on surgeon’s performance. This will emphasize the importance of quantitative muscle fatigue analysis, although it might be challenging to apply performance analysis to real surgical operations. However, performance analysis can be easily applied to dry lab experiments such as Fundamentals of Laparoscopic Surgery (FLS) tasks, since there are validated methods to analyze the performance of subjects. Therefore, the combination of RQA, as a quantitative muscle fatigue analysis method, and standard performance analysis of FLS tasks will open up new doors for future studies related to muscle fatigue and performance analysis of real MIS operations which would be very beneficial for surgical community. Also, future studies should include larger number of subjects and more surgical operations in order to confirm findings of this study.

5. Conclusion

In this study, recurrence quantification analysis (RQA) was applied to EMG data recorded from eight muscle groups of five surgeons while doing fifteen MIS operations. The results showed that this novel measure could detect the sign of muscle fatigue on bilateral deltoid and trapezius at 45–55 minutes after operations began, and no sign of fatigue was found on other muscles. It was found from this study that trapezius and deltoid were the most vulnerable muscles among all eight muscle groups tested. Here, it is worthwhile to note the nature of MIS procedures such that surgeons manipulate long and thin instruments for the frequent changes of their orientation and position mainly using bicep and triceps and maintain their arm posture using deltoid while gazing through monitor using trapezius muscle for prolonged time. Considering this nature, deltoid and trapezius may show clear sign of fatigue while other muscles do not show any sign of fatigue as tested in this study. Based on the results, it could be suggested that surgeons need to take break time at 45–50 minutes after MIS operations begin in order to minimize muscle fatigue and other possible muscle disorders. To have a better understanding about the effect of recovery time on surgeons doing such a complicated laparoscopic operation, future research might test different recovery times to find an optimum amount of time for muscle relaxation before muscles fatigue. In addition, ergonomic improvement of surgical equipment to reduce the intensity level might be efficient in reducing muscle fatigue or at least increasing time-to-fatigue.

Competing Interests

The authors declare that they have no competing interests.

References

R. Berguer, J. Chen, and W. D. Smith, “A comparison of the physical effort required for laparoscopic and open surgical techniques,” Archives of Surgery, vol. 138, no. 9, pp. 967–970, 2003.

R. Berguer, D. L. Forkey, and W. D. Smith, “Ergonomic problems associated with laparoscopic surgery,” Surgical Endoscopy, vol. 13, no. 5, pp. 466–468, 1999.

U. Matern, P. Waller, C. Giebmeyer, K. D. Rückauer, and E. H. Farthmann, “Ergonomics: requirements for adjusting the height of laparoscopic operating tables,” Journal of the Society of Laparoendoscopic Surgeons, vol. 5, no. 1, pp. 7–12, 2001.

T. M. Kowalewski, Real-time quantitative assessment of surgical skill [Ph.D. dissertation], University of Washington Department of Electrical Engineering, 2012.

M. Riojas, C. Feng, A. Hamilton, and J. Rozenblit, “Knowledge elicitation for performance assessment in a computerized surgical training system,” Applied Soft Computing, vol. 11, no. 4, pp. 3697–3708, 2011.

E. Domany, “Protein folding in contact map space,” Physica A: Statistical Mechanics and Its Applications, vol. 288, no. 1–4, pp. 1–9, 2000.

R. Berguer, W. D. Smith, and Y. H. Chung, “Performing laparoscopic surgery is significantly more stressful for the surgeon than open surgery,” Surgical Endoscopy, vol. 15, no. 10, pp. 1204–1207, 2001.

A. G. Gallagher and R. M. Satava, “Virtual reality as a metric for the assessment of laparoscopic psychomotor skills,” Surgical Endoscopy and Other Interventional Techniques, vol. 16, no. 12, pp. 1746–1752, 2002.

J. H. Peters, G. M. Fried, L. L. Swanstrom et al., “Development and validation of a comprehensive program of education and assessment of the basic fundamentals of laparoscopic surgery,” Surgery, vol. 135, no. 1, pp. 21–27, 2004.

P. Dolan, A. F. Mannion, and M. A. Adams, “Fatigue of the erector spinae muscles-a quantitative assessment using ‘frequency banding’ of the surface electromyography signa,” Spine, vol. 20, no. 2, pp. 149–159, 1995.

M. P. Van der Grinten, “Test-retest reliability of a practical method for measuring body part discomfort,” in Designing for Everyone, Y. Quéinnec and R. Daniellou, Eds., pp. 54–57, Taylor & Francis, 1991.

N. Hubert, M. Gilles, K. Desbrosses, J. P. Meyer, J. Felblinger, and J. Hubert, “Ergonomic assessment of the surgeon's physical workload during standard and robotic assisted laparoscopic procedures,” The International Journal of Medical Robotics and Computer Assisted Surgery, vol. 9, no. 2, pp. 142–147, 2013.

M. Heiden, E. Lyskov, M. Djupsjöbacka, F. Hellström, and A. G. Crenshaw, “Effects of time pressure and precision demands during computer mouse work on muscle oxygenation and position sense,” European Journal of Applied Physiology, vol. 94, no. 1-2, pp. 97–106, 2005.

J. V. Basmajian and C. J. De Luca, Muscles Alive: Their Functions Revealed by Electromyography, Lippincott Williams & Wilkins, Baltimore, Md, USA, 5th edition, 1962.

C. J. De Luca, “Myoelectrical manifestations of localized muscular fatigue in humans,” Critical Reviews in Biomedical Engineering, vol. 11, no. 4, pp. 251–279, 1984.

O. Lippold, “The tremor in fatigue,” Ciba Foundation Symposium, vol. 82, pp. 234–248, 1981.

T. Bosch, M. P. de Looze, and J. H. van Dieën, “Development of fatigue and discomfort in the upper trapezius muscle during light manual work,” Ergonomics, vol. 50, no. 2, pp. 161–177, 2007.

G. P. Y. Szeto, L. M. Straker, and P. B. O'Sullivan, “A comparison of symptomatic and asymptomatic office workers performing monotonous keyboard work–1: neck and shoulder muscle recruitment patterns,” Manual Therapy, vol. 10, no. 4, pp. 270–280, 2005.

J. T. Dennerlein, V. M. Ciriello, K. J. Kerin, and P. W. Johnson, “Fatigue in the forearm resulting from low-level repetitive ulnar deviation,” American Industrial Hygiene Association Journal, vol. 64, no. 6, pp. 799–805, 2003.

K. R. Sahlin and J. M. Ren, “Relationship of contraction capacity to metabolic changes during recovery from a fatiguing contraction,” Journal of Applied Physiology, vol. 67, no. 2, pp. 648–654, 1989.

E.-P. Takala, “Static muscular load, an increasing hazard in modern information technology,” Scandinavian Journal of Work, Environment and Health, vol. 28, no. 4, pp. 211–213, 2002.

M. L. Uhrich, R. A. Underwood, J. W. Standeven, N. J. Soper, and J. R. Engsberg, “Assessment of fatigue, monitor placement, and surgical experience during simulated laparoscopic surgery,” Surgical Endoscopy and Other Interventional Techniques, vol. 16, no. 4, pp. 635–639, 2002.

P. S. Slack, C. J. Coulson, X. Ma, K. Webster, and D. W. Proops, “The effect of operating time on surgeons' muscular fatigue,” Annals of the Royal College of Surgeons of England, vol. 90, no. 8, pp. 651–657, 2008.

C. L. Webber Jr., M. A. Schmidt, and J. M. Walsh, “Influence of isometric loading on biceps EMG dynamics as assessed by linear and nonlinear tools,” Journal of Applied Physiology, vol. 78, no. 3, pp. 814–822, 1995.

G. Filligoi and F. Felici, “Detection of hidden rhythms in surface EMG signals with a non-linear time-series tool,” Medical Engineering and Physics, vol. 21, no. 6-7, pp. 439–448, 1999.

N. Marwan, M. C. Romano, M. Thiel, and J. Kurths, “Recurrence plots for the analysis of complex systems,” Physics Reports, vol. 438, no. 5-6, pp. 237–329, 2007.

J.-P. Eckmann, S. Oliffson Kamphorst, and D. Ruelle, “Recurrence plots of dynamical systems,” EPL, vol. 4, no. 9, pp. 973–977, 1987.

F. Takens, “Detecting strange attractors in turbulence,” in Dynamical Systems and Turbulence, Warwick 1980, D. Rand and L.-S. Young, Eds., vol. 898 of Lecture Notes in Mathematics, pp. 366–381, Springer, New York, NY, USA, 1981.

G. B. Mindlin and R. Gilmore, “Topological analysis and synthesis of chaotic time series,” Physica D: Nonlinear Phenomena, vol. 58, no. 1–4, pp. 229–242, 1992.

Y. Liu, M. Kankaanpää, J. P. Zbilut, and C. L. Webber Jr., “EMG recurrence quantifications in dynamic exercise,” Biological Cybernetics, vol. 90, no. 5, pp. 337–348, 2004.

B. R. Jensen, M. Pilegaard, and G. Sjøgaard, “Motor unit recruitment and rate coding in response to fatiguing shoulder abductions and subsequent recovery,” European Journal of Applied Physiology, vol. 83, no. 2-3, pp. 190–199, 2000.

D. Farina, L. Fattorini, F. Felici, and G. Filligoi, “Nonlinear surface EMG analysis to detect changes of motor unit conduction velocity and synchronization,” Journal of Applied Physiology, vol. 93, no. 5, pp. 1753–1763, 2002.

K. Walker-Bone and C. Cooper, “Hard work never hurt anyone: or did it? A review of occupational associations with soft tissue musculoskeletal disorders of the neck and upper limb,” Annals of the Rheumatic Diseases, vol. 64, no. 10, pp. 1391–1396, 2005.

S. Kvarnström, “Occurrence of musculoskeletal, disorders in a manufacturing industry with special attention to occupational shoulder disorders,” Scandinavian Journal of Rehabilitation Medicine, Supplement, vol. 8, no. 1, pp. 1–114, 1983.

H. Christensen, “Muscle activity and fatigue in the shoulder muscles during repetitive work—an electromyographic study,” European Journal of Applied Physiology and Occupational Physiology, vol. 54, no. 6, pp. 596–601, 1986.

L. Rosendal, A. K. Blangsted, J. Kristiansen et al., “Interstitial muscle lactate, pyruvate and potassium dynamics in the trapezius muscle during repetitive low-force arm movements, measured with microdialysis,” Acta Physiologica Scandinavica, vol. 182, no. 4, pp. 379–388, 2004.

G. Sjøgaard, G. Savard, and C. Juel, “Muscle blood flow during isometric activity and its relation to muscle fatigue,” European Journal of Applied Physiology and Occupational Physiology, vol. 57, no. 3, pp. 327–335, 1988.

G. W. Kroon and M. Naeije, “Recovery of the human biceps electromyogram after heavy eccentric, concentric or isometric exercise,” European Journal of Applied Physiology and Occupational Physiology, vol. 63, no. 6, pp. 444–448, 1991.

E. Henneman, G. Somjen, and D. O. Carpenter, “Excitability and inhibitability of motoneurons of different sizes,” Journal of Neurophysiology, vol. 28, no. 3, pp. 599–620, 1965.

B. Rolander, D. Jonker, A. Karsznia, and T. Öberg, “Evaluation of muscular activity, local muscular fatigue, and muscular rest patterns among dentists,” Acta Odontologica Scandinavica, vol. 63, no. 4, pp. 189–195, 2005.

Santy and S. Z. Dawal, “Investigation on time to fatigue for upper limb muscle during a repetitive light assembly task,” in Proceedings of the International MultiConference of Engineers and Computer Scientists (IMECS '10), vol. 3, pp. 1584–1588, March 2010.