Development of modified rapid entire body assessment (MOREBA) method for predicting the risk of musculoskeletal disorders in the workplaces

Background Rapid Entire Body Assessment (REBA) technique is one of the tools developed for predicting the risk of musculoskeletal disorders based on the effective risk factors. This method has several limitations. The present study was aimed to develop the Modified Rapid Entire Body Assessment (MOREBA) method to more accurately predict the risk of musculoskeletal disorders. Materials This cross-sectional study was performed on 300 male workers of a steel factory with a variety of job tasks in Iran. Then, the information related to the various physical risk factors was extracted through observation of their duties and conversation with them. Also, the subjects were asked to complete the Persian version of Cornell musculoskeletal discomfort questionnaires (CMDQ). Then, a theoretical model was drawn in AMOS software. Computed coefficients were used to develop the MOREBA equation. In the end, the final scores were categorized by ROC curves, and the validation of the novel method was investigated using linear regression analysis. Results The parameters evaluated in the MOREBA method included contact stress, rapid and sudden movement, throwing motion, hand-arm vibration, whole-body vibration, temperature, and work-rest cycle in addition to the parameters of the REBA method, including awkward posture, coupling, force, load, static activity, and repetitive activity. The results showed that the strain produced by the physical risk factors with the total effect coefficient of 0.783 could significantly affect the musculoskeletal symptoms. The computed coefficients of the risk factors were applied to develop a novel index. The final score of the MOREBA method was categorized into four levels by optimal cut-off points of 12.37, 16.51, and 24.35. Based on the results, the MOREBA and REBA methods could justify 67 and 55% of the variations of musculoskeletal symptoms, respectively. Conclusions The results revealed that modifications conducted in the REBA method were effective, and the MOREBA method can provide a more accurate prediction of the risk of musculoskeletal disorders.


Introduction
Musculoskeletal disorders (MSDs) are one of the major problems in various workplaces, such as office and operational environments. Musculoskeletal disorders are defined as a group of preventable disorders affecting nerves, tendons, muscles, and supporting structures posture, and contact stress are the main ergonomic factors associated with upper extremities musculoskeletal disorders [19]. Various tools have been developed for predicting the risk of musculoskeletal disorders based on the effective risk factors. Some of the observational methods include Rapid Upper Limb Assessment (RULA), Rapid Entire Body Assessment (REBA), Quick Exposure Check (QEC), Ovako Working Posture Analysis System (OWAS), and Novel Ergonomic Postural Assessment Method (NERPA). Among them, REBA is one of the known methods to evaluate different body parts, including upper limbs (arm, forearm, and wrist), lower extremities, trunk, and neck [20]. However, this tool has several limitations. REBA assesses only some of the mentioned risk factors, consisting of awkward posture, load/force, coupling, and repetitive and static activities [20]. Also, several assessments are required for different tasks of one person [21]. Moreover, REBA can overestimate the risk of MSDs in the employees. Sabini et al. observed that there is an overestimation of 45% in the cases assessed with REBA [22]. Therefore, this method can more accurately predict the risk of musculoskeletal disorders if these limitations are corrected. For this reason, the present study was aimed to develop the Modified Rapid Entire Body Assessment (MOREBA) method and evaluate its validation.

Participants
This cross-sectional study was performed on 300 male workers of a steel factory with a variety of job tasks in Iran. For selecting the desired subjects, various parts of the stated industry including steelmaking, forging, machining, supporting, and administrative units carefully visited, and several occupations were chosen so that the extensive ranges of values related to the risk factors are collected. Then, 410 subjects employed in these job tasks were invited to the study and evaluated in terms of inclusion and exclusion criteria, of which 300 persons participated in the study. Data collected of these 300 subjects were used to develop and validate the novel method. Inclusion criteria included age range from 19 to 55 years, having work experience more than one year, not having a second job, and not having heavy physical activity in leisure time. Exclusion criteria also consisted of having a history of major trauma )such as driving, sports, and occupational accidents(, rheumatic diseases, spine surgeries, and large joint surgeries, and having musculoskeletal structural deformities )such as spinal abnormalities and genu varum(. Other exclusion criteria were longterm intake of corticosteroids and immunosuppressants and non-cooperation during the study.

Sample size calculation
Given the aim of the present study, the minimum probable correlation between the score of the MOREBA and the score of musculoskeletal symptoms was considered by 0.2 and the sample size with a confidence level of 95% and a test power of 90% was computed as follow: Where Z 1− α 2 is equal to 1.96 based on a confidence level of 95%, Z 1−β is equal to 1.29 based on a test power of 90%, and W is equal to 0.203 based on a minimum correlation coefficient of 0.2. Therefore, the lowest sample size in this study was 259 persons.

Data collection
The protocol of this study was reviewed and approved by the Medical Ethics Committee of Baqiyatallah University of Medical Sciences (IR.BMSU.REC.1399.366). The demographic characteristics of the participants, including age, work experience, height, weight, and physical activity, were collected. Then, the information related to the risk factors of awkward posture, coupling, contact stress, static activity, repetitive activity, rapid and sudden movement, and throwing motion was recorded through observing the duties of the subjects and the data related to the risk factors of load, force, hand-arm vibration, whole-body vibration, temperature, and the work-rest cycle was gathered by an interview with them. Two persons of the researchers, who were well versed with the method, performed the assessments. Table 1 describes  the guidance of scoring the factors. Based on this table, for evaluating the risk factors, the structured answers were used in steps of observation and interview. These risk factors were determined via a literature review in the valid databases and consultants with the specialists of occupational health and occupational medicine. The most worst and frequent positions related to each of body organs during work time were evaluated, and final scores related to each group (A and B) were estimated by tables of the REBA method [23]. The scores related to other factors were obtained based on the five-point Likert scales from 0 to 4 during work time. The final score related to force was calculated by scores of maximum load weight and load carrying time and the final score related to force was computed by scores of maximum force value and work time. Moreover, the most worst and frequent posture of each person was evaluated by the REBA method, and the final score was calculated using the assessment of factors of awkward posture, coupling, force/load, and muscle use [24]. Also, the subjects were asked to complete the Persian version of the Cornell musculoskeletal discomfort questionnaire (CMDQ) developed by Hedge et al. during rest periods. This questionnaire with a body map evaluates the frequency (never, 1-2 times last week, 3-4 times last week, once every day, and several times every day), intensity (slightly uncomfortable, moderately uncomfortable, and very uncomfortable), and interference (not at all, slightly interfered, and Substantially interfered) related to musculoskeletal symptoms in 12 regions of the human body at last week [25]. To calculate the total score of discomfort caused by musculoskeletal disorders, scores of frequency (0, 1.5, 3.5, 5, and 10), severity (1, 2, and 3), and interference (1, 2, and 3) at each of body regions was multiplied by each other and then obtained values were summed together [26]. The validity and reliability of the Persian version of this questionnaire were investigated by Afifzadeh et al. and confirmed by them in a steel factory. In this study, Cronbach's alpha coefficient of frequency, intensity, and interference dimensions were computed by 0.955, 0.961, and 0.969 [27]. Finally, The accuracy of the answers in the completed questionnaires was checked by reviewing medical records and performing medical examinations.

Predictive model
Structural equation modeling (SEM) was applied to predict the occurrence of musculoskeletal symptoms in the present study. For this purpose, a theoretical model was drawn, in which the total score of evaluated musculoskeletal discomforts was considered as a criterion variable, and the effect of strain resulting from the defined risk factors on it was investigated. Minimum acceptable factor loading was equal to 0.3 for retaining the item in the model [28]. The values of the Cronbach's alpha coefficient, average variance extracted (AVE), and composite Reliability were also estimated. The minimum acceptable limits of these values were 0.7, 0.5, and 0.7, respectively [29,30]. Finally, the indirect effect coefficients of the risk factors were extracted from the model and applied to develop the MOREBA equation. For calculating the indirect effect coefficient of each item, its direct effect coefficient was multiplied by the effect coefficient of the strain on musculoskeletal symptoms.

Statistical analyses
Data were recorded into the statistical package for the social sciences (SPSS) version 18. For handling the missing data, the method of the mean of nearby points was applied to impute the missing data. The values out of the logical range related to the variables (outliers) were identified and corrected. The normality of variables was examined using Skew and kurtosis curves. Given that the statistical distribution of all variables was normal, the Pearson test was used to calculate the correlation coefficients between them. Then, a theoretical model was drawn in AMOS software. The fitness of this model was evaluated using absolute, comparative, and normed fit indices. After that, the MOREBA equation was written by indirect effect coefficients of the risk factors on the occurrence of musculoskeletal symptoms. Finally, the MOREBA score was categorized into four levels using receiver operator curves (ROC) analysis. The total scores of 450, 900, and 1350 related to the musculoskeletal discomforts were considered as boundaries of risk levels [31]. Nearest points to the ideal state in ROC curves were adopted as optimal cut-off points in the MOREBA index. The validity of the developed method was also examined by linear regression analysis. Moreover, the frequency distributions of the risk levels estimated by CMDQ, REBA, and MOREBA were computed. Table 2 represents the statistical distribution of demographic characteristics and studied variables in the participants. The results of examining the Skew and kurtosis curves showed that the statistical distribution of all variables was normal. The CMDQ score related to musculoskeletal symptoms was variable from 0 to 1674. Table 3 reports the correlation matrix of the studied variables. The results indicated that there were meaningful correlations between all studied risk factors and CMDQ score (P<0.01). The highest correlation coefficients were related  (4) Posture group B Assessing the most worst and frequent positions related to upper arm, lower arm, and wrist during work time and calculating score B using the table of REBA method.

Results
Repetitive activity -Never (0) -Little (1) -Sometimes (2) -Much (3) -Very much (4) Coupling status -Very good (0) -Good (1) -Acceptable (2) -Poor (3) -Very poor (4) Rapid and sudden movement - (4) Throwing motion (such as hitting with a hammer or hand) (4) Load Maximum load weight -Less than 5 kg (0) -5 to 10 kg (1) Figure 1 depicts the theoretical model for predicting the occurrence of musculoskeletal symptoms due to occupational conditions. The results showed that the strain produced by the risk factors with the total effect coefficient of 0.783 could significantly affect the musculoskeletal symptoms. Table 4 describes the effect coefficients of the risk factors on the variations of musculoskeletal symptoms. The factor loading (direct effect coefficient) of any item was not less than 0.3. Therefore, all items retain in the model. Posture group A (0.734) and B (0.714) had the greatest indirect effect coefficients on musculoskeletal symptoms, respectively. The lowest indirect coefficients also were related to variables of whole-body vibration (0.257) and contact stress (0.272), respectively. Table 5 reports the goodness-of-fit indices of the drawn model. Based on the results, the goodnessof-fit of this model was confirmed. The values of the average variance extracted (AVE), Cronbach's alpha coefficient, and composite Reliability were computed by 0.50, 0.93, and 0.91, respectively. The results showed that these values were higher than the minimum acceptable limits.
The indirect effect coefficients of the risk factors were used to develop the equation of Modified Rapid Entire Body Assessment (MOREBA) as follow: Where P A is the score of posture group A, P B is the score of posture group B, C is the score of coupling, CS is the score of contact stress, L is the score of load, F is score of force, SA is the score of static activity, RA is the score of repetitive activity, RM is the score of rapid movements, TM is the score of throwing motions, HAV is the score of hand-arm vibration, WBV is the score of whole-body vibration, T is the score of air temperature, and WRC is the score of work-rest cycle.
The scores of load and force parameters also are calculated as follow: Where t L is the score of load-carrying time and W L is the score of maximum load weight.    Where t w is the score of work time and V F is the score of maximum force value. Figure 2 shows the receiver operating characteristic (ROC) curves. The results revealed that the optimal cut-off points between low and moderate risk, moderate and high risk, and high and very high risk  Table 6 represents the risk levels and equivalent MOREBA scores. Figures 3 and 4 show the linear regression curve between the CMDQ score and MOREBA score and REBA score, respectively. Based on the results, the MOREBA and REBA methods could justify 67 and 55% of the musculoskeletal symptoms, respectively. Moreover, Fig. 5 indicates the frequency  distribution of the risk levels estimated by CMDQ, REBA, and MOREBA. As a result, the REBA method overestimates the risk level of musculoskeletal disorders while the relative frequency of different risk levels in MOREBA method was nearly close to those in CMDQ evaluation.

Discussion
The results of the present study showed that the extensive range of values with normal distributions was collected. It indicates that a variety of situations have been studied, and the novel method can be used in workplaces with different physical risk factors. Based on the results, these various situations have been led to wide variations in the occurrence of musculoskeletal disorders. Moreover, it was observed that strain produced by the physical risk factors with the effect coefficient of 0.783 could significantly affect musculoskeletal disorders in people. Therefore, the MOREBA method can significantly predict the occurrence of musculoskeletal disorders. It may be due to the fact that almost all major physical risk factors were exploited in developing the novel method. The results of a study performed by Maakip et al. also showed that physical demands with 61% had the highest relative importance among various groups of effective factors for predicting musculoskeletal discomforts in Malaysian office workers [32]. However, the results of previous studies show that other groups of the risk factors such as individual and psychosocial items can impress on the occurrence of musculoskeletal disorders [32,33]. So, if higher numbers of the risk factors are entered into the model, a more accurate prediction may be obtained. Based on the results, the highest indirect effect coefficients on the occurrence of musculoskeletal symptoms were related to the factors of posture group A, posture group B, load, force, throwing motion, repetitive activity, static activity, and coupling, respectively. Most of these factors are essential items in other famous ergonomics risk assessment methods. In the methods of Rapid Upper Limb Assessment (RULA) and Novel Ergonomic Postural Assessment (NERPA), the factors of awkward posture, force/load, static activity, and repetitive activity are evaluated. The REBA method assesses the coupling in addition to factors of the RULA method [34]. In the occupational repetitive actions (OCRA) method, the factors of posture, force, frequency, recovery, duration, and some other items are examined [35]. In the Loading on the Upper Body Assessment (LUBA) technique, the body postures are appraised and other important factors are not considered [36]. In Job Strain Index (JSI), the factors of hand and wrist posture, force, duration, frequency, and work speed are evaluated [37]. Chander et al. also developed an observational method for Postural Ergonomic Risk Assessment (PERA) using three factors of posture, force, and work duration [38]. In the mentioned methods, posture and load/force variables have the most effect in calculating the final score of the risk assessment, respectively. Furthermore, Cheong et al. concluded that physical factors of posture and lifting/ pulling had the greatest impact coefficients on workrelated upper extremity musculoskeletal disorders in male workers [39]. the results of a study conducted by Mehralizadeh et al. indicated that the highest total effect coefficient on musculoskeletal complaints in the hospital nurses was related to the factor of body posture [40]. Therefore, the results of the present study are consistent with the results of previous studies. Moreover, the results showed that the factors of whole-body vibration, contact stress, temperature, and hand-arm vibration had the lowest indirect effect coefficients on the occurrence of musculoskeletal symptoms, respectively. It may be argued that the workers had less exposure to those compared to other risk factors in the workplaces and each of them affects fewer body regions. However, these factors also have significant effects on the occurrence of musculoskeletal disorders. Charles et al. concluded that occupational exposures to the whole body or hand-arm vibration resulted in musculoskeletal disorders in the shoulder and neck regions [41]. Based on the results of a study conducted by Pullopdissakul et al., factor of contact stress is associated with upper extremities musculoskeletal disorders [19]. Magnavita et al. found a significant interaction between temperature complaints and strain for upper limb disorders [42]. Therefore, the use of those in assessing the risk of musculoskeletal disorders can result in a more accurate prediction. However, a low number of ergonomic risk assessment methods, such as quick exposure check (QEC), has involved some of these factors.
The area under ROC (AUC) related to different risk levels was greater than 0.80 which indicates good    the majority of workstations (70%) in the high-risk category [44]. These findings are consistent with the results of the present study. As one of the limitations in this study, all participants were male and the female workers were not studied. Moreover, objective examinations were not conducted for the diagnosis of musculoskeletal disorders. However, to ensure the correctness of the musculoskeletal discomforts reported by the participants in the CMDQ questionnaire, their medical records were reviewed and medical examinations were performed. It should be noted that the Iranian workers are annually examined by occupational medicine specialists in terms of various diseases, such as musculoskeletal disorders and their medical records are kept.

Conclusions
The results demonstrated that modifications conducted in the REBA method were effective, and the novel method of MOREBA can provide a more accurate prediction of the risk levels of musculoskeletal disorders. The parameters of awkward posture, load, force, throwing motion, repetitive activity, static activity, and coupling in this method had the greatest effect coefficients, respectively. Moreover, the novel method has nearly rectified the problem of overestimating the risk of musculoskeletal disorders in the REBA method. Therefore, this method can be reliably used to predict the risk of musculoskeletal disorders produced by physical factors in the workplace. However, it is suggested that the MOREBA method be more validated by other assessors in other industries in future studies.