Abstract
Affective responses concern customers’ affective needs and have received increasing attention in consumer-focused research. To design a product that appeals to consumers, designers should consider multiple affective responses (MARs). Designing products capable of satisfying MARs falls into the category of multiobjective optimization (MOO). However, when exploring optimal product form design, most relevant studies have transformed multiple objectives into a single objective, which limits their usefulness to designers and consumers. To optimize product form design for MARs, this paper proposes an integrated model based on MOO and multicriteria decision-making (MCDM). First, design analysis is applied to identify design variables and MARs; quantification theory type I is then employed to build the relationship models between them; on the basis of these models, an MOO model for optimization of product form design is constructed. Next, we use nondominated sorting genetic algorithm-II (NSGA-II) as a multiobjective evolutionary algorithm (MOEA) to solve the MOO model and thereby derive Pareto optimal solutions. Finally, we adopt the fuzzy analytic hierarchy process (FAHP) to obtain the optimal design from the Pareto solutions. A case study of car form design is conducted to demonstrate the proposed approach. The results suggest that this approach is feasible and effective in obtaining optimal designs and can provide great insight for product form design.
1. Introduction
Because of the intense competition in contemporary markets, it is essential that companies produce products that meet the needs of consumers. Affective responses have attracted increasing attention as indicators of consumers’ affective needs. Kansei engineering (also known as affective engineering) is a method for translating affective responses into design variables [1]. It has achieved considerable success in the field of product design [2, 3].
Affective responses to products have various aspects. Thus, the process of optimizing the form of a product capable of satisfying multiple affective responses (MARs) can be considered as a multiobjective optimization (MOO) problem, which is an important issue of Kansei engineering. Deriving an optimal product design is a complicated process in which several methods need to be combined. For example, Hsiao et al. [4] used quantification theory type I (QT1) to build the predictive models for MARs and then adopted analytical hierarchy process (AHP) to obtain the weights for MARs; on the basis of these weights, the values of the affective responses were integrated into a single objective value for product form optimization. Guo et al. [5] integrated the values of MARs, which were predicted using a back propagation neural network (BPNN) into a single objective value for design optimization by using a genetic algorithm (GA). These studies have addressed MARs; however, they transformed the MOO problem into a single objective optimization (SOO) problem in which only one solution is provided in each simulation run. This approach, which has the advantage of simplicity, can be considered a classical MOO [6]. However, it always requires advance preference information prior to optimization in circumstances with no knowledge of the possible consequences; furthermore, it precludes the consideration of all possible solutions, which provides the decision-makers with no alternative options because only one optimal solution is obtained after optimization. Thus, the classical MOO approach is of limited usefulness to designers and consumers.
In recent years, it has been recognized that MOO and multicriteria decision-making (MCDM) are mutually enhancing; considerable attention has been focused on combining MOO with MCDM to assist decision-makers in solving problems in real-world applications [7–9]. One essential approach for combining MOO with MCDM is the use of MOO to derive numerous Pareto solutions and then apply MCDM to select the best solution from the Pareto set [10]. This approach does not require advance information related to the preferences of the decision-maker. Moreover, it can provide a number of Pareto solutions in a single run in the form of a set of nondominated solutions rather than a single best solution. On the basis of these Pareto solutions, the decision-maker can make a trade-off decision. Therefore, this approach is more appropriate for addressing multiple objectives simultaneously [11] and tends to be more practical and less subjective [6].
Recently, several researchers in the product design field have investigated the application of MOO for obtaining Pareto optimal solutions. For example, Yang [12] used support vector regression (SVR) and nondominated sorting genetic algorithm-II (NSGA-II) in the construction of a hybrid Kansei engineering system (KES) for product form design. Jiang et al. [13] proposed a chaos-based NSGA-II for solving an MOO model in product design. However, these studies have focused mainly on obtaining Pareto solutions, not on further combining MCDM methods to find optimal design solution according to the consumers’ preference. Few studies have investigated using a combined MOO-MCDM approach to realize optimal product form design.
This paper proposes a new design model which combines MOO with MCDM for addressing the problem of optimizing product form design. We employ NSGA-II as an MOO algorithm to generate Pareto optimal solutions. The fuzzy analytic hierarchy process (FAHP) is then adopted as the MCDM method to derive the optimal product form design. The rest of this paper is organized as follows: Section 2 briefly reviews the research context. Section 3 presents our proposed methodology. Section 4 demonstrates a case study with the proposed method, the results of which are discussed in Section 5. Finally, the conclusion is presented in Section 6.
2. Literature Review
2.1. Kansei Engineering System
Kansei engineering was founded at Hiroshima University about 30 years ago [14]. It focuses on the implementation of the consumer’s affective response into product design. The Kansei implies the consumer’s psychological feeling as well as embracing physiological issues. Kansei engineering is defined as translating the consumer’s Kansei into design. The KES means a computer assisted design system of Kansei engineering which is usually based on artificial intelligence technology [15]. It can be used to meet the needs of consumers and facilitate more efficient product development. When implementing KES studies, the products are usually decomposed into design variables [16, 17]. The design variables can be discrete variables such as the type of design components [18, 19] and can also be continuous variables such as the coordinates of the feature points on the contour [20, 21]. The result of KES study can provide deep insight for design.
Obtaining optimal designs that satisfy various consumers’ affective responses, which falls into the category of MOO, has been a crucial subject in KES studies. Essentially, two approaches to solving MOO problems can be applied. One is the classical MOO method that converts multiple objectives into a single objective. This approach has been extensively used in KES [4, 5, 22] and begins with the determination of a preference weight vector for use in the construction of a composite fitness function. Using this approach to obtain the most favourable product design is simple but somewhat subjective [6]. The other approach that involves combining MOO with MCDM has been used in engineering in recent years [8, 9, 23]. However, few studies have investigated using KES for product form design. The KES approach integrating MOO with MCDM does not require advance information on preferences and provides numerous Pareto solutions in a single run. Designers and consumers can evaluate the advantages and disadvantages of each of these Pareto solutions by using the MCDM method based on preference information and comparing them to derive the optimal product design.
These two approaches differ fundamentally with regard to the use of preference information. The classical MOO approach requires that the preference vector be supplied without any knowledge of the possible consequences. By contrast, combining MOO with MCDM uses preference information for selection of a single solution from the obtained Pareto set. Hence, this approach is more practical and less subjective. Therefore, the combined MOO-MCDM approach more effectively facilitates the product design process and is more appropriate to apply in KES. In this study, we apply this approach to KES for the purpose of obtaining optimal product form design.
2.2. Multiobjective Optimization
MOO addresses the multiple, often conflicting objectives that occur in real-world scenarios. An MOO problem can be formally defined as follows. Find the decision vector , which satisfies inequality constraints and equality constraints and optimizes the vector of the objective function [24]. It is not possible to optimize all objective functions that are in conflict with each other simultaneously by using one solution vector . Therefore, a logical approach to MOO is to obtain a set of nondominated solutions which are known as Pareto optimal solutions. A solution is called “Pareto optimal” if and only if there is no for which dominates . In other words, none of the solutions can be deemed better than the others in relation to all the objectives, and each ideally represents an optimal compromise between the objectives [25].
Evolutionary algorithms are generic population-based metaheuristic optimization algorithms that mimic evolutionary principles. They are often used to solve optimization problems with multiple objectives, and are then called multiobjective evolutionary algorithms (MOEAs). MOEAs have the advantage of entailing no assumptions about the shape of Pareto fronts and being able to simultaneously search for a set of approximating Pareto fronts. Hence, they are recognized as powerful and versatile tools and have been used for solving MOO for over two decades. The nondominated sorting genetic algorithm (NSGA) proposed by Srinivas and Deb [26] is an essential type of MOEA. It has been relatively successful for several years. However, it is an inefficient algorithm mainly because of its method for classifying individuals. Deb et al. [27] proposed an improved version of NSGA, named NSGA-II. NSGA-II is currently the most commonly used MOO algorithm and has been used as a foundation for designing other algorithms. It has less time complexity and uses elitism and crowded comparison operator to ensure diversity that is very important to derive representative design solutions. Consequently, we adopt NSGA-II for Pareto set generation.
2.3. Multicriteria Decision-Making
MCDM is a discipline aimed at supporting decision-makers who are faced with numerous and conflicting alternatives in reaching an optimal decision [28]. MCDM has been successfully applied in practice and investigated by many researchers and scientists. The analytic hierarchy process (AHP) is an MCDM method that was originally developed by Saaty [29]. It is a structured qualitative and quantitative technique for addressing multiple criteria decision problems and has been successfully applied in many fields. However, when using AHP, the result of pairwise comparison must be expressed with crisp values, which does not conform to the fact that real-world judgement is subjective, uncertain, and fuzzy. van Laarhoven and Pedrycz [30] first added the fuzzy concept to AHP. Triangular fuzzy numbers were used to express judgements of the relative importance of two criteria. Because of insufficient use of the AHP method, Buckley [31] improved Saaty’s traditional AHP and developed FAHP to incorporate fuzzy numbers into pairwise comparison. FAHP is very suitable for the user to express subjective preference related to affective responses. Therefore, we adopt FAHP as the MCDM method to derive the optimal design solution from the Pareto set.
3. Proposed Methodology
This paper proposes an integrated product form design model based on MOO and MCDM. To facilitate comprehension, assume that a consumer requires a product, such as a car, and has preferences that include “Modern,” “Comfortable,” “Rounded,” and “Simple.” Thus, the research question concerns how to derive an optimal car form design which corresponds to a specific consumer’s preference weights for “Modern,” “Comfortable,” “Rounded,” and “Simple.” The proposed approach entails first obtaining many cars which satisfy Pareto optimality and then selecting the optimal car corresponding to the preference weights from these obtained cars.
Figure 1 shows the framework of the proposed methodology. QT1 is used to build the predictive models for affective responses. Based on these predictive models, an MOO model is constructed. NSGA-II is then used to search for the Pareto optimal design solutions. Subsequently, FAHP is employed to determine the optimal product form design solution.
3.1. Construction of Multiobjective Optimization Model
In the proposed methodology, design analysis that includes product form analysis and affective response analysis should first be conducted. Product form analysis can facilitate clarifying the character of the form for ascertaining the design variables. Our method uses morphological analysis, which can decompose a product into a set of design variables; optimal product form can then be identified by combining these design variables. The number of possible combinations is typically extremely high, and they entail not only existing and conventional solutions but also a wide range of variations and completely novel solutions [32]. To derive the experimental samples, we adopt orthogonal arrays that are a set of tables devised by Taguchi and are used to specify the minimum number of experiments needed. With regard to analysis of affective responses, Kansei words are adopted to describe people’s responses to a product and factor analysis is then applied to identify the affective factors underlying these words. According to the affective response factors, the objectives for MOO and the criteria for making trade-off decision can be determined.
To construct the MOO model, it is necessary to first build models that use design variables to predict affective responses. In this study, QT1 is adopted, which has most often been used to build the predictive models for affective responses [1, 4]. QT1 is an extension of multiple linear regression analysis by using dummy coding. It can calculate nominal scaled explanatory variables and interval or proportionally scaled objective variables. Furthermore, it has the advantage of showing the contribution of each design variable for affective response quantitatively [33]. The degree of influence of each level of design variables on the affective responses is calculated as follows:where is the affective response, is the level of design variable, is the weight for the level, is the number of design variables, is the number of levels, and is the random variable.
The MOO problem in this study involves maximizing or minimizing all of the MRAs. For the maximization problem, the MOO model can be constructed as follows:where is the affective response predicted using QT1.
3.2. Generation of Pareto Optimal Solutions by Using NSGA-II
The nondominated solutions can be obtained by solving the MOO model. Figure 2 illustrates the nondominated solutions for an optimization problem with two affective responses. The feasible area is inside (and the infeasible area is outside) the region surrounded by the curve. Each point in the feasible area represents a solution. The points located on the bold curve, which denotes the Pareto optimal frontier, are Pareto optimal solutions. Thus, the points , , , , , and are the Pareto optimal solutions. Of any two such solutions, one is superior with respect to a specific affective response, but this superiority occurs only from sacrificing one of the other affective responses. Consequently, a decision-making process is necessary to determine which of them is the optimal product form design that most satisfies the affective needs of the consumer.
To generate the nondominated solutions, NSGA-II is adopted. Figure 3 shows the schematic of the NSGA-II algorithm [6]. This algorithm combines a parent population of size and an offspring population of size to generate population with a size . The ranks of individuals in a population are then calculated through nondominated sorting. Next, classify all individuals according to their ranks, and let be a subpopulation of individuals with the same rank. The next generation is filled according to the order of the priority of these ranks, starting with the best nondomination ranks, followed by the second nondomination ranks and so on, until is filled. Regarding the final rank, the solutions which have a larger crowding distance are prioritized for filling . Finally, the population of offspring population is generated from the parent population by using tournament selection, crossover, and mutation. NSGA-II uses crowding distance to maintain diversity among nondominated individuals. Thus, it provides well-distributed individuals in the objective space, which is crucial for generating representative product design solutions.
3.3. Selection of Best Solution by Using FAHP
After obtaining Pareto optimal solutions, a decision-maker can use the MCDM method to determine which of the solutions is most suitable for the specific affective responses preference. In this research, FAHP is employed to determine the best design solution.
To represent subjective pairwise comparisons of customer requirements, triangular fuzzy numbers are used, which have the advantages of computational efficiency and ease of data acquisition. A triangular fuzzy number has the following membership function:where and denote the lower and upper bounds of the fuzzy number , respectively, and denotes the modal value (Figure 4). The triangular fuzzy number can be denoted by .
-cut is a means for converting a fuzzy set into a crisp set, which is very useful for defining the arithmetic operations on fuzzy numbers [34]. The -cut of a fuzzy set is a crisp set that contains all the elements that have a membership grade in fuzzy set greater than or equal to the specified value of (Figure 4). By defining the level of , the triangular fuzzy number can be characterised as
To obtain the importance weights for multiple affective responses using FAHP has the following steps.
Step 1 (construct the fuzzy comparison matrix for MARs). Triangular fuzzy numbers , , , , and indicate that criterion is of relative importance to criterion and then , , , , and indicate that criterion is of relative importance to . The linguistic terms, corresponding fuzzy numbers, and membership functions are listed in Table 1, which were defined by Kwong and Bai [35].
The fuzzy pairwise comparison matrix for all the criteria is constructed as follows:where
Step 2 (determine the importance weight of MARs). The lower limit and upper limit of the fuzzy numbers with respect to level can be defined as follows:Using the interval arithmetic and cuts, can be expressed as follows: cuts of fuzzy matrix can then be denoted as follows:The degree of satisfaction for the judgement matrix is estimated by the index of optimism . A higher value of indicates a higher degree of optimism. In practical applications, , , and can be used to indicate that the decision-maker has an optimistic, moderate, or pessimistic view, respectively [36]. The index of optimism is a linear convex combination defined as [37]The following matrix can be obtained after setting the index of optimism :The eigenvalues of matrix can be calculated by solvingThe eigenvector with respect to the highest eigenvalue can be calculated usingThe weighting vector for the criteria can be determined after normalizing .
Step 3 (verify the consistency of the evaluation). To ensure that the evaluations are consistent when the user conducts pairwise comparisons, a consistency index (C.I.) and consistency ratio (C.R.) are adopted. If C.I. and C.R. are both less than 0.1, this indicates that the pairwise comparison is consistent [38]. The equations of C.I. and C.R. are as follows:where is the highest eigenvalue of , is the number of evaluation criteria, and R.I. refers to a random consistency index whose value increases if the number of criteria increases, as shown in Table 2 [29].
Step 4 (select the best design solution). The decision matrix for criteria and alternatives can be expressed asThus, the performance value of the th design solution can be calculated aswhere denotes the normalized weight of the th criterion, which is calculated using FAHP. Finally, rank the design solutions according to the performance values. The solution ranking 1st is the optimal design.
4. Implementation Procedures
4.1. MOO Model Construction
4.1.1. Design Analysis
A car is a representative example of a mature product. Car profile, which is a concentrated expression of the car form and provides a rich form language of car design, is receiving increasing emphasis in consumer car choice; considerable research has been conducted on this topic [18–21]. Therefore, this study adopts a car profile case to demonstrate the proposed method. A total of 120 profile images of cars currently on the market were collected. These images were reviewed by seven experts in the field of car design who used morphological analysis to construct the design variables that influence consumer feeling. These design variables can be divided into two categories: one is relevant to the form of different profile parts, as illustrated in Table 3; the other is relevant to the ratio of different profile parts, which also contributes greatly to the affective responses of product form [16], as shown in Table 4. The design variables of 6, 8, and 12 and the corresponding levels were adopted from Lai et al. [19].
The possible combinations of these design variables can generate 1,594,323 ( = 1,594,323) design solutions. To derive the experimental samples, Taguchi’s orthogonal array was adopted, and 27 combinative samples were obtained, as shown in Figure 5.
In Kansei engineering studies, Kansei words have been used to characterise the affective responses of consumers to products. A total of 258 affective descriptors were collected from car magazines and car company websites. From them, 20 pairwise Kansei words were selected as most relevant to car design. We employed factor analysis to investigate the factor structure of these words, wherein a minimum eigenvalue of 1 was adopted as the cutoff value. Table 5 lists the results of the factor analysis. The 20 pairwise Kansei words were grouped into four main factors, which account for 92.752% of the cumulative variance. When naming the factors, we employed pairwise Kansei words with maximum loading coefficients. Therefore, the four factors were named as “Traditional–Modern,” “Uncomfortable–Comfortable,” “Orthogonal–Rounded,” and “Complex–Simple.” On the basis of these factors, we can establish the objectives for MOO and also the criteria for MCDM.
4.1.2. Model Construction
To construct an MOO model, experimental data must first be obtained. A questionnaire was designed by combining the 27 representative samples with the four affective factors by using a 7-point semantic differential method (1–7), in which, for example, 1 is the most traditional and 7 is the most modern for the “Traditional–Modern” Kansei word pair. We recruited 60 participants (30 females and 30 males, aged 22–30 years) to fill out the questionnaire. The results are shown in Table 6.
We applied the obtained experimental data to construct predictive models for affective responses using QT1; the affective response values for the four factors could thereby be predicted using the design variables, and the contribution of every design variable to affective response can also be obtained from the partial correlation coefficient. For example, the predictive model for “Traditional–Modern” is presented in Table 7. From this model, the affective response can be predicted, and it is clear that the car’s height to length ratio (6) is the most crucial design variable affecting “Traditional–Modern,” which has the highest partial correlation coefficient (0.468), followed by the shape of the top part (3) and the ratio of rear part’s length to the car’s length (11); the ratio of the rear overhang’s length to the car’s length (13) is the least crucial design variable, which has the lowest partial correlation coefficient (0.035).
For the semantic differential scale, to derive a design having a right-side affective response, it is necessary to maximize the scale; to derive a design having a left-side affective response, it is necessary to minimize the scale. In this study, it is assumed that the consumer requires a product which appears “Modern,” “Comfortable,” “Rounded,” and “Simple,” all of which are factors that lie on the right side of the semantic differential scale. Therefore, the MOO problem is the maximization of the four scales and the MOO model can be presented as follows:
4.2. Pareto Optimal Solutions Generation
NSGA-II was employed to solve the MOO model. Simulated binary crossover (SBX) and polynomial mutation were adopted in NSGA-II. The SBX and polynomial mutation are widely used in MOEA research [25]. The tuned parameters for optimization were specified as follows: 25,000 evaluations; a population size of 100; a distribution index of 15 for SBX; a crossover rate of 0.8; a distribution index of 20 for mutation; and a mutation rate of , where indicates the number of variables. A total of 100 Pareto solutions were obtained, as shown in Table 8. The Pareto optimal frontier in objectives’ space is illustrated in Figure 6.
(a)
(b)
4.3. Best Solution Selection
In this study, we used FAHP to determine the best design solution. The computational procedure is as follows.
Step 1. Assume that a consumer requires a car and has affective preferences that are expressed using pairwise comparisons, as shown in Table 9.
On the basis of the pairwise comparison, the fuzzy comparison matrix can be obtained as follows:
Step 2. For a moderate decision-maker, set and [39]. Combining (7), (8), (9), (10), and (11), can be derived as follows:Substituting into (12) yields the highest eigenvalue . Substituting into (13) provides the corresponding eigenvectors. Having been normalized, the importance weights of the affective responses can be determined as follows:
Step 3. In conducting the conformance test, we can derive that the C.I. is 0.044 and the C.R. is 0.049 by using (14). Because the values of C.I. and C.R. are both less than 0.10, we can conclude that the consistency of the judgement matrix is satisfied.
Step 4. On the basis of Table 8, the decision matrix for the four affective responses and the 100 alternatives can be calculated as follows:According to the importance weights of the affective responses, we can compute the performance value matrix by using (16) as follows:Rank the design alternatives by the performance values asThe ranking of the Pareto solutions is depicted in Figure 7. Solution number 37 ranks 1st. Consequently, we can conclude that the optimal product form design solution is solution number 37.
4.4. Verification Result
To verify the result, typical design solutions were selected which corresponded to the indicative points in the ranking shown in Figure 7, such as the starting point, the end point, and the jumping points. The design variables, affective responses, and design sketches of these solutions are shown in Table 10.
We combined the design sketches with the four affective responses to construct the questionnaire, which was filled in by a total of 60 respondents. The statistical scores are shown in Table 11, where the performance values were calculated according to the importance weights. It is evident that solution number 37 solution number 27 solution number 2 solution number 98; this sequence indicates that the results of the evaluation for the Pareto set are the same as the statistics of the affective responses of actual respondents. Therefore, this method can be used by designers to derive various Pareto optimal design solutions corresponding to consumers’ multiple affective responses and can also obtain the optimal design solution on the basis of consumers’ preference weights for these responses.
5. Discussion
This paper proposes an integrated product form design method in which the MOO process is applied to obtain a Pareto set, and then the MCDM process is performed to derive the optimal design according to the preference weights. The MOO is a type of search process and the MCDM is a type of decision-making process. To obtain an optimal design solution, the search process and the decision process need to be integrated. In the proposed product form design model, the decision process is posterior to the search process; hence, it can be categorised as a posterior preference articulation approach [9]. The decision-maker directly applies higher-level information related to affective preference to a set of preestablished Pareto optimal designs and thus enables clearer insight into design decisions. By contrast, for a traditional prior preference articulation approach which transforms MOO into SOO, the decision-maker must quantify opinions before viewing the potential design solutions. Therefore, the final solution with the posterior preference articulation approach can reflect the decision-maker’s preference more accurately [40]. Furthermore, this design method is particularly useful in cases where the user has no prior knowledge through which to ascertain the preference weights for MARs [41].
In this integrated design model, when the user wants to change the importance weights for multiple affective responses, the optimal design can be chosen from the already obtained Pareto set merely by executing the decision-making process; the search process, which generally consumes considerable time by using GA [4, 5], does not require execution. Thus, the proposed model is more convenient and efficient.
We use QT1 to construct the predictive models for affective responses, on the basis of which the MOO model is established and the Pareto solutions are then obtained. QT1 provides an indication of the influence of various design variables on affective responses through the partial correlation coefficients, thereby facilitating the identification of crucial design variables. Nevertheless, the prediction accuracy tends to be low. If a more accurate predictive model is required, BPNN [5] or SVR [42] can be adopted. In addition, there are only 27 representative samples used in the experiment for prediction. The prediction reliability may be low because of the small sample size. Thus, more representative samples could be added to improve the prediction reliability.
The execution of the proposed method is demonstrated using a case study involving the design of a car profile. It should be mentioned that the design variables and corresponding levels can be changed according to the cars available on the market and the research focus. This approach can be applied to other products with various design elements, in which case design analysis should be performed to ascertain both the design variables and affective responses. Subsequently, the MOO and the MCDM can be executed using NSGA-II and FAHP, respectively.
6. Conclusion
Product form design has been recognized to contribute critically to a corporation’s performance in the highly competitive contemporary global market. This paper has presented an integrated product form design model based on MOO and MCDM. The car profile design was used to illustrate the methodology, and the Pareto solution set satisfying MARs was derived and the optimal design solution most corresponding to the preference of consumer was then obtained. The results indicate that the proposed product form design model can address multiple objectives for MARs simultaneously and ensure that the design decision process is less subjective and more practical. Furthermore, this approach can provide great insight for design and is particularly helpful for users without prior knowledge in specifying the relative importance of criteria.
It should be noted that in this study all possible product form designs are generated from the combination of the defined design variables. Therefore, the obtained new forms are always from the known specified patterns. How to develop a new form design different from the collected forms is important for the realistic design. Further research would focus on this topic. In addition, there are many types of MOO methods, for example, NSGA-II, PESA-II, and SPEA2, as well as many types of MCDM methods, for example, AHP, FAHP, and TOPSIS. In the future, comparison studies could be carried out to make sure which type of MOO or MCDM method is most suitable in the context of Kansei engineering.
Competing Interests
The authors declare that there are no competing interests regarding the publication of this paper.
Acknowledgments
This work was supported by the Ministry of Science and Technology, Taiwan (Project no. MOST 102-2221-E-006-213). The authors would like to express their sincere thanks to the participants for their assistance in this study.