Development of a Method toMeasure the Quality ofWorking Life Using the Improved Metaheuristic Grasshopper Optimization Algorithm

'is paper provides a method to numerically measure the quality of working life based on the reduction of human resource risks. It is conducted through the improved metaheuristic grasshopper optimization algorithm in two phases. First, a go-to study is carried out to identify the relationship between quality of working life and human resource risks in the capital market and to obtain the factors from quality of working life which reduce the risks.'en, a method is presented for the numerical measurement of these factors using a fuzzy inference system based on an adaptive neural network and a new hybrid method called the improved grasshopper optimization algorithm.'is algorithm consists of the grasshopper optimization algorithm and the bees algorithm. It is found that the newly proposed method performs better and provides more accurate results than the conventional one.


Introduction
Working life keeps changing like any other environmental condition of human life in today's world. To ensure survival in this turbulent environment, organizations and employees have no chance but to adapt to the working environment changes, which will not happen unless the quality of working life (QWL) is improved. Focusing on the improvement of the working-life quality leads to increased employee vitality and satisfaction and, in turn, the success of employees, organizations, and stakeholders [1]. Since the satisfaction of employees affects their physical conditions and mental health, organizations need to keep their employees satisfied in terms of both the work environment and the work itself, thus accomplishing the intended goals [2]. Previous studies have shown that employees with good mental health are more committed and efficient than those in low spirits [3].
erefore, the quality of working life reflects organizational or management culture, based on which employees can experience being self-directed, responsible, and selfesteemed. In this regard, an organization's effectiveness and efficiency can be enhanced if it provides opportunities for the employees to progress [4].
Furthermore, human resources are increasingly viewed as capital [5]. Accordingly, the viability and profitability of organizations in today's turbulent situation largely depend on their ability to cope with the rapid changes in the environment and utilize the created opportunities [6]. is is possible only by the use of enterprise risk management techniques, such as holistic threats and opportunities management. To increase the value for stakeholders in the short or long term, the organization provides funds and supervises its risks by managing the market value, controlling and utilizing any possible source [7]. From an organizational perspective, the accepted definitions of risk often refer to adverse outcomes caused by an unpredictable event in the business environment [8].
ere are other opinions such that risk has variable consequences, i.e., sometimes positive and sometimes negative [9]. In a general definition, human resource risks refer to the set of risks associated with the programs and processes in the human resource of an organization [10]. Ernest & Young Institute considers the risks in employee-related plans and processes as the risks in human resources that, if properly managed, will make the organization one of the market leaders [11].

Literature Review
Previous studies on the quality of working life mainly measured QWL based on the level of need satisfaction [12]. For example, investigations were conducted on the relationship between QWL and job satisfaction [13] and the relationship between QWL and the organizational commitment among the lecturers in a Malaysian public research university [14]. One can also refer to the reports on QWL in child protection in two work environments [15], QWL amongst the workers of state-owned commercial banks in a Bangladesh northern region [16], QWL and the mental health of primary care physicians [17], QWL with an emphasis on human resource interventions in health units [18], QWL and the related factors for migrating workers in Guangdong Province, China [19], the effect of income and the relationship between leadership style and QWL [20] as well as QWL and burnout among nurses in critical care sections [21], quality of work life and contribution to productivity: assessing the moderator effects of burnout syndrome [22], quality of work life among young medical doctors [23], the effect of quality of work life (QWL) on job satisfaction and organization citizenship behavior [24], and determinants of quality of work life among nurses working [25].
ere are also numerous studies on human resource risks such as investigations on substitution plans in New York State [10], risk perspective in human resources management [9], human resource risk management using the normal accident theory [26], and the impact of perceived human resource risk factors on financial problems [27]. According to the literature, the studies ever conducted in the field are either on quality of working life or on human resource risks. ere is little research to examine these two issues together and their mutual effects. Accordingly, the present research seeks to propose a method of evaluating quality of working life on the basis of human resource risk reduction so as to enhance the quality of working life. It is to be noted that this type of study has never been carried out in the capital market. erefore, in this study, a hybrid system is modelled to measure and improve working life quality based on the reduction of human resource risks.

Methodology
3.1. Framework. Using the tenets of previous studies, primary conceptual model has been made.
As shown in Figure 1, this model presents the various indicators widely reported, as in [28][29][30][31], for human resource risks. ere are 11 categories of indicators identified to cover almost all human resource risks. Also to evaluate the quality of working life, Walton model [32] serves as a tool.
First of all, the conceptual model of the research (Figure 1) served as a basis to identify and investigate the factors that affect working life quality and human resource risks. en, the effect of that quality on those risks was evaluated and the factors which could reduce the risks were identified. Finally, a method was proposed to numerically measure these factors. e data used in this study were divided into training and testing categories. For the training part, questionnaires were collected from the experts, and the data were checked by the specialists. e training phase began according to the scores given by the experts and specialists. To improve the QWL model and analyze it, the adaptive neurofuzzy inference system (ANFIS) method was used. e rules of a fuzzy inference system are based on the "if-then" set, which forms the core of the system. Once the rules are put into practice, the relationship between input and output variables can be identified. is system can be used as a forecasting model in conditions where the input or output data are highly uncertain. e basis of this approach is a set of fuzzy inference system (FIS) input and output [33].
e statistical population of this study comprised the employees of the financial institutes in the Iranian capital market. In this regard, references were made to investing companies, investment funds, investment consulting centers, and the share market for securities in Tehran. Some of the experts in these institutes were randomly selected to provide the required data. e steps of this research, as shown in Figure 2, are as follows: (i) Conduction of measurements and the achievement of more precise and reliable results using the improved grasshopper optimization algorithm (ii) Evaluation of the results and the use of the grasshopper optimization algorithm combined with the bees algorithm (iii) Introduction of a measurement system by means of the neurofuzzy inference system (iv) Identification of the factors affecting working life quality and human resource risks through a questionnaire and confirmatory factor analysis (v) Exploration of the relationship between working life quality and human resource risks through a questionnaire and the structural equation modelling technique 3.2. Adaptive Neurofuzzy Inference System. In this paper, the rules for training a triangular membership network are examined. e network has three membership functions, and there are two membership functions for each input. In order to model the neural-fuzzy system, the data were segmented. e segmentation was, indeed, done to train and validate the system. Eighty percent of the data was taken into account as the training data, while the remaining twenty percent was used for testing and validation. e program was then executed based on the code generated in the MATLAB software and through the method already mentioned.
As shown in Figure 3, since the factors in the chain continuity framework were composed of one main layer and several sublayers, the neural-fuzzy inference system was modelled according to the number of the layers. At the end of the modelling phase, there were 8 neural-fuzzy systems.

Grasshopper Optimization Algorithm (GOA).
Grasshopper optimization algorithm (GOA) is a metaheuristic algorithm inspired by nature. It mimics and simulates grasshoppers' behavior and their group movement toward food sources in nature. In order to solve an optimization problem, a GOA mathematical model is of benefit with its ability to mimic grasshoppers' behavior in nature [34]. Simulations have shown that, in comparison with the latest well-known algorithms, GOA can provide better results. Putting simulation results into practice has proven that GOA is capable of solving real problems with unknown space. e mathematical model used to simulate grasshoppers' behavior is as follows [35]: where X i is the position of grasshopper number i, S i is the social interaction, G i is the gravitational force, and A i is the wind flow rate.
To present a random behavior, the equation can be written as follows [36]: where r 1 , r 2 , r 3 are the random numbers in [0, 1].
In the above equation, d ij is the distance between grasshoppers i and j, which can be calculated as d ij � |X i − X j |. Moreover, S function, which stands for social interaction, can be obtained using the following equation: where G i can also be calculated as G i � −ge ∧ g , in which e ∧ g is a unit vector that tends to the center of the earth.
Finally, A i follows the equation below, where u is a fixed drift and e ∧ w is a wind flow unit vector: In general, according to the equations presented, the mathematical model used to simulate the grasshoppers' behavior is as follows: In order to solve optimization problems, a modified version of this equation is presented as follows [37]:

Bees Algorithm.
In the bees algorithm, half of the bees' population is considered as workers and the other half as nonworkers. Since there is only one worker bee for each food source, the number of worker bees is equal to that of food sources around the hive. e working bees whose resources are provided by the supervising bees and the other worker bees transform into scout bees [38]. A supervisor bee selects a food source according to the probability value associated with that food source. e value is obtained from the following equation: where p i is the probability of selecting i bee, fit i is the eligibility rate of i solution, and SN is the number of food sources.
In this algorithm, the worker bees share their information with the supervisor bees, and then the supervisor bees use this information to select a new food source around the previous one. To make this selection, the following equation is used [39]. In this equation, X i is the answer to i optimization problem and v i is the new answer around the answer x i . Figure 4, one of the disadvantages of using the grasshopper optimization algorithm is the convergence of grasshoppers in multiple repetitions [37]. In this study, the GOA and bees algorithm were combined to eliminate the weak point of GOA (i.e., grasshoppers' convergence in large numbers).

Improved Grasshopper Optimization Algorithm (the Recommended Method in is Research). As shown in
e resulting algorithm was used as an improved grasshopper optimization algorithm. In other words, as GOA and bees algorithm were combined, a new grasshopper was added in each repetition, which prevented convergence. Ultimately, this new algorithm led to better, more accurate and more reliable results.

Phase I (Qualitative Phase): Descriptive Results.
To check whether or not the distribution of the research variables was normal, the Kolmogorov-Smirnov test was used. As a statistical hypothesis, if the test significance level was greater than 0.05, the null hypothesis cannot be rejected. It was concluded that the distribution of the variables was not significantly different from a normal distribution. e test results are presented in Table 1.
In the next step, to identify the factors affecting the quality of working life and the human resource risks in the capital market, a confirmatory factor analysis model was employed. is model is made through a structural equation and the Smart-PLS software, which makes use of partial least squares. e data were collected through a questionnaire distributed among 38 capital market specialists selected from a statistical sample of 128 capital market employees. Considering the variables and the predicted relations, the data were fed into the Smart-PLS software for a hypothetical model. After the deletion of the unknown variables, the final model emerged as shown in Figure 5. Figure 5 shows the final factors of QWL and human resource risks. e presented figures are the absolute values of the t scores above 1.96; those at the 95% confidence level were considered significant. To check the reliability of the model (i.e., the internal similarity between the observed variables), two indicators were used, a traditional indicator called Cronbach's alpha and a composite reliability indicator. e acceptable criterion for the internal similarity of the indicators was 0.7, which was confirmed for both models. In order to check the validity, the AVE square root of each structure was compared with the correlation values between the structures. Since the coefficients of the correlation between each structure and the others were less than the AVE root square of that structure, the validity of both models was confirmed. e transverse load test and Fornell-Larcker test were also used to investigate the divergent validity; the output of the matrices proved the existence of divergent validity. To evaluate and prioritize the factors affecting human resource risks, it was necessary to identify them in the capital market. To this end, expert opinions were obtained by a two-section questionnaire. In the first section, the risks were scored, and the second was given to the pairwise comparisons of the RPN criteria. After the risk levels of the RPN factors were determined, the risk priority number (RPN) was used to determine the risk level of each   (RPN � Severity * Occurance * De tection).
As mentioned earlier, the failure modes and effects analysis (FMEA) method has several limitations for risk assessment. In this study, Fuzzy Analytical Hierarchy Process (FAHP) was implemented to address one of those limitations (considering the same importance for the RPN risk factors). To this end, in the second section of the questionnaire, the experts were asked to compare the matrices of the risk factor RPNs in a pairwise manner and state their preferences for the factors according to a fuzzy spectrum table. en, in order to obtain the preferences of the decision-maker, a pairwise comparison of the matrices was conducted using fuzzy numbers (l. m. u). ereafter, the obtained fuzzy values were diffused and normalized by the method of the center of gravity, and the risk factors were obtained for the RPNs.
According to Table 2, the risks of "conflict of interests in the organization, unethical behavior, and existence of lobbying culture in the organization" account for the highest RPNs, which means they are the most critical human resource risks in the capital market. On the other hand, the risks of "incompatibility of tasks with key skills and abilities of the organization's employees, failure to evaluate performance by specific criteria, and lack of proper feedback to employees" have the lowest RPNs, which means these risks are not of much importance for the capital market. To determine the RPNs of human resource risks in the capital market, the critical risks of this market need to be identified. To this end, the Pareto principle or 80-20 rule is used. e rule asserts that 80% of outcomes (or outputs) result from 20% of all the causes (or inputs) for any given event. In other words, 20% of human resource risks cause 80 problems. erefore, it is not necessary to examine all the risks to reduce the adverse risks. According to the Pareto principle, human resource risks are as shown in Figure 6.
As Figure 6 suggests, the risks of r8 and those on its left are the ones that cause 80% of the problems. In other words, the risks that are considered critical in the Pareto principle include lack of cooperation and teamwork culture, lack of law enforcement oversight, incompatibility of training with human resources requirements, lack of up-to-date knowledge and skills in employees, wrong rules and policies, recruitment without rules or without considering the employees' qualifications, lack of human skills, poor leadership, lack of employee commitment and job abandonment, lack of managers' perceptual skills, lack of proper successor training process, low tendency to moral principles at work, lobbying culture in the organization, and unethical behavior and conflict of interests in the organization. In this study, since the Pareto principle is not proved yet, the specialists' opinions in the questionnaire were utilized to validate the results. To this end, in order to determine the risk priority number, the questionnaire also asked whether or not human resource risks were critical. In this case, the specialists' positive answers beyond the third quartile were considered as critical risks. ese critical risks are shown in Table 3. erefore, through a correspondence made between the results obtained from the Pareto principle and those from the specialists, the commonality of the two sets of results was identified as the critical risks of human resources in the capital market. e risks are as follows:       reduction of the critical risks in human resources. Once the normality, validity, and reliability of the model were confirmed, at first, the relationship of QWL with human resource risks was investigated, as shown in Figure 7. After the general model was evaluated and fitted, the quality of working life as related to human resource risks was interpreted. Based on the results gained from the test of the model, it is concluded that the model has a perfect fit. e results of the structural equation modelling, i.e., the calculated characteristics of each model path including path coefficients and t values, were considered at a 5% error level, as presented in Table 4.
In the first step, to qualify the working life model on the basis of reducing the risk of human resources and investigate their relationships, a path model was implemented with the structural equation modelling technique and the Smart-PLS software. e software uses partial least squares (PLS) to provide structural equation models. is approach is less dependent on sample size than the covariance-based approach is; therefore, it can be used for small-size samples. On this basis, a hypothetical research model was introduced to the Smart-PLS software according to the predicted variables and relations. e examined model is displayed in Figure 8. Figure 8 shows the calculated QWL model values based on reducing the risks of human resources significantly. e absolute values of the calculated quantities, if above 1.96, were considered significant at the confidence level of 95%. Based on the above integrated model, the second phase of the study could start. As it can be seen in the figure above, each factor of the quality of working life had at least one significant inverse relationship with one of the critical risk factors of human resources.

Phase II (Quantitative Phase).
In this section, a method was presented to determine the quality of working life based on human resource risks. First, the experts scored the items based on an available scale. e items included providing opportunities for growth and continuous security, fair payment, safe and healthy work environment, integration and cohesion in the organization, legalism in the organization, general living space, and developing human capabilities. en, they scored the quality of working life based on the scores given to those items. After the scoring, to perform the necessary fuzzy inference process, the questionnaire data were converted from qualitative to quantitative and entered into the system. Finally, the fuzzy values had to be converted to quantitative values and go into the defuzzification process. Individual symbols were also assigned to each item for convenience, which are shown in Table 5.
e Gaussian membership functions for the POGS item and the improved grasshopper optimization algorithm (for items 1, 2, and 3 from left to right) are as follows.
As shown in Figure 9, for each of the items, 13 figures were drawn in the same way as the above figure and according to the gradient descent and improved grasshopper algorithms. Figure 9 is just one example from those 13 figures. To create a managerial perspective and quickly evaluate the impacts of the items on the categories of the quality of working life model based on the reduction of human resource risks, the influence of each item is reported in comparison to its subitems; each item has two important subitems. POGS and FP, as two example items, are presented in Figure 10.
In the next step, an adaptive fuzzy neural inference system was used to implement a decision-making assistant system for institutions and organizations operating in the capital market.
e aim was to measure the quality of working life, identify the best option to improve it, and quantify its effect. e inference engine of the adaptive neurofuzzy inference system was TSK level one, the default inference engine. e number of the rules in the adaptive neurofuzzy inference system was different for each item, which can be seen in Table 6.
In addition to the Gaussian membership function, the triangular membership function was also examined. Besides, an adaptive neurofuzzy inference system was performed based once on the gradient descent training algorithm (classic method) and once on the improved grasshopper optimization training algorithm. e results of these conductions were compared later. In order to evaluate the performance error of the employed methods, MAE, RMSE, and R 2 were used as three criteria. ey are formulated as follows: In these three formulations, N represents the total number of the observations done in the predicted period, y t and y ∧ t denote the predicted value and the actual value at time t , respectively, and X K , Y K, and K stand for the observational values, number of data, and estimated values, respectively.
As presented in Table 7, the adaptive neurofuzzy inference system was first implemented based on the Gaussian membership function and the gradient descent training algorithm with two and three membership functions. According to Table 8, the two criteria of MAE and RMSE from the adaptive neurofuzzy inference system with two membership functions at the values of 33.4 and 60.5, respectively, had much poorer performance than that with three membership functions and the values of 71.3 and 63.4. In addition, the R2 index of the adaptive neurofuzzy inference system with two membership functions and the value of 93.0 was smaller than that with three membership functions, which indicates the better performance of the adaptive neurofuzzy inference system with two membership functions. In general, it can be concluded that an adaptive neurofuzzy inference system with three membership 8 Mathematical Problems in Engineering       functions performs better than the one with two membership functions. Accordingly, instead of the conventional method (i.e., gradient descent method), the improved grasshopper optimization algorithm (which is a combination of GOA and bees algorithms) was put into practice. e results gained from this improved algorithm are presented in Table 8. According to Tables 7 and 8, an ANFIS based on the Gaussian membership function and the improved grasshopper optimization training algorithm with three membership functions yields better results than those based on the Gaussian membership function and the gradient descent algorithm with three membership functions. It is to be noted that when the ANFIS was implemented through the improved grasshopper optimization training algorithm with a triangular membership function, as expected, it had a high error. Other membership functions, such as trapezoidal function, were abandoned due to nondifferentiability and lack of flexibility.

Mathematical Problems in Engineering
In general, according to the results of the second phase of the study, the use of the grasshopper optimization algorithm and its combination with the bees algorithm leads to more accurate results. Indeed, since grasshoppers' convergence in numerous repetitions is a disadvantage of using the GOA method, combining it with the bees algorithm provides more accurate and reliable results; it has been known that adding a new grasshopper at every repetition stage helps to avoid convergence.

Conclusion
Considering the gap in the literature on QWL and human resource risks and their evaluation, this study was conducted to devise a method for measuring QWL so as to reduce human resource risks in the capital market. e research was done in two phases. In the first phase, the factors affecting QWL and human resource risks as well as their mutual  impacts were identified and studied. In the second phase, a method was provided to measure the characteristics obtained in the first phase. As the results of the first phase showed, all the identified QWL factors reduced at least one of the human resource risks. In this regard, the quality of working life was found to have the most significant impact on the reduction of the 'unethical behavior' as a risk. e unethical behavior of employees can indeed cause many problems. It can even push a business to failure. One way to cope with this problem is to pay special attention to the quality of working life and the reduction of risks in human resources. e best way to monitor a factor is to measure and display it numerically. is makes it possible to improve the plans for operation and the efficacy of the measures taken. e second phase of the study is dedicated to the method proposed for the measurement of the quality of working life in link with the risks involved in human resources. e measurements were made through the neurofuzzy inference system by means of the improved grasshopper optimization algorithm. is algorithm was a combination of the grasshopper optimization algorithm and the bees algorithm. In this respect, the results of the second phase of the research show that the proposed method to measure the quality of working life based on human resource risks is more accurate and reliable than the conventional method in the adaptive neural network (i.e., reduction gradient method). e proposed method serves to create a decision-assistant system for managers and business owners to measure the quality of working life, select the best option to improve it, and ultimately measure the feedbacks on its effectiveness numerically. Inspired by this method, capital market managers can obtain the QWL of their employees and take appropriate measures to improve it, which, in turn, leads to the reduction of risks in human resources and the growth of their business. According to the account provided in this research for the capital market, it is suggested that the issue is examined in other fields and occupations such as banking and insurance and that the results are compared with those of this research.

Data Availability
Some of the numerical data and software outputs data used to support the findings of this study are included within the article. Also this data in full are currently under embargo. But requests for data, after publication of this article, will be considered by the corresponding author.

Conflicts of Interest
e authors declare that they have no conflicts of interest.