Optimizing the Processing Factor and Formulation of Oat-Based Cookie Dough for Enhancement in Stickiness and Moisture Content Using Response Surface Methodology and Superimposition

: Despite the utilization of dusting ﬂour and oil to reduce dough stickiness during the production process in food industry, they do not e ﬀ ectively help in eliminating the problem. Stickiness remains the bane of the production of bakery and confectionery products, including cookies. In addition, the high moisture content of cookie dough is unduly important to obtain a high breaking and compression strengths (cookies with high breaking tolerance). This study was conducted in light of insu ﬃ cient research hitherto undertaken on the utilization of response surface methodology and superimposition to enhance the stickiness and moisture content of quick oat-based cookie dough. The study aims at optimizating, validating and superimposing the best combination of factors, to produce the lowest stickiness and highest moisture content in cookie dough. In addition, the e ﬀ ect of ﬂour content and resting time on the stickiness and moisture content of cookie dough was also investigated, and microstructure analysis conducted. The central composite design (CCD) technique was employed and 39 runs were generated by CCD based on two factors with ﬁve levels, which comprised ﬂour content (50, 55, 60, 65, and 70%), resting time (10, 20, 30, 40, and 50 min) and three replications. Results from ANOVA showed that all factors were statistically signiﬁcant at p < 0.05. Flour content between 56% and 62%, and resting time within 27 and 50 min, resulted in dough with high stickiness. High-region moisture content was observed for ﬂour content between 60% and 70%, and within 10 and 15 min of resting time. The optimized values for ﬂour content ( V 1 ) = 67% and resting time ( V 2 ) = 10 min. The predicted model (regression coe ﬃ cient model) was found to be accurate in predicting the optimum value of factors. The experimental validation showed the average relative deviation for stickiness and moisture content was 8.54% and 1.44%, respectively. The superimposition of the contour plots was successfully developed to identify the optimum region for the lowest stickiness and highest moisture content which were at 67–70% ﬂour content and 10–15 min resting time.

. Difference in number of experimental runs between full factorial design and response surface methodology (RSM) [39].

Factors Levels Total Number of Experimental Runs
Full Factorial Design RSM 2 5 25 (5 2 ) 13 3 5 125 (5 3 ) 20 4 5 625 (5 4 ) 31 5 5 3125 (5 6 ) 52 Since response surface methodology (RSM), specifically central composite design (CCD), has been widely used in process optimization [40][41][42][43][44], it will, therefore, be used in this study. The purpose of the present study was: (i) to investigate the effect of factors including flour content and dough resting time on the stickiness and moisture content of oat-based cookie dough with large inclusion; (ii) to perform optimization, validation and superimposition processes to obtain the best combination of factors that produce the lowest stickiness and highest moisture content in cookie dough; and (iii) to conduct a microstructure analysis of cookie dough and evaluate the result.

Design of Experiment Based on RSM
The experiment using RSM was designed based on two factors, namely flour content (V 1 ) and resting time (V 2 ). The responses, i.e., stickiness and moisture content, were assumed to be influenced by these factors. The value for other factors (ingredients) was fixed. Samples with 100% wheat flour acted as control samples, and were compared with the optimized sample after the optimization process. For central composite design (CCD), five levels, and their working ranges were chosen, as shown in Table 2.

Experiment Design
A central composite design (CCD) involving five levels and two factors, with three replications, were applied for all design points, comprising a total of 39 experimental runs. The factors were chosen based on the literature, their significant effect on the responses, and their workability (within the allowed working range). The complete CCD with coded and uncoded levels of the factors are shown in Table 3. The value for the total block is 1, with the experiments carried out in a randomized order. An analysis of variance (ANOVA) was carried out to determine the significance of the main factors and their interactions. The significance level was set at 95% with a p-value of 0.05. The regression coefficient model (mathematical model) generated from the ANOVA table will be used for optimization purpose, depending on the value of the correlation coefficient, R 2 . Experimental data were fitted to the second-order polynomial model to derive the regression coefficient model. The model for the response surface analysis is shown in Equation (1), where Y is the response, β 0 , β i , β ii , and β ij are regression coefficients for the intercept, linear, quadratic, and interaction terms, respectively. X i and X j are coded values for the independent variables [45].

Raw Materials
The raw materials were bought from local manufacturers and certified halal, according to the Malaysian Standard on Halal Food (MS 1500:2009). The halal certificate, issued by Malaysia's Department Islamic Development (JAKIM), is a document that certifies the products that have followed and cover the Islamic dietary guidelines [46]. The superfine wheat flour (9 g/100 g protein) used in this study was produced by Prestasi Flour Mill (M) Sdn. Bhd, Malaysia. The quick oats were obtained from Cimei Food Ingredients Sdn. Bhd, Malaysia. The characteristics of wheat flour and quick oats are shown in Table 4. In this study, only wheat flour was used. Instead of oat flour, quick oat was incorporated in the dough at the designated ratio. Quick oats are flakes with thickness between 0.356 and 0.457 mm. It is produced from steel-cut groats, which are usually produced by sectioning groats into three or four pieces, before the final steaming and flaking process [47]. Quick oats were used in this study, due to the health benefit of oats and increasing interest among researchers and manufacturer.
Other ingredients used in the cookie recipe were coarse grain sugar and soft brown sugar (

Cookie Dough Preparation
Except for flour and oats, the formulations for all samples were prepared using the same amount of ingredients, as shown in Table 5. The oat to flour composition was measured based on cup percentage. The weight of these ingredients was based on the American Standard Matrix Measurement and Cups Conversion [48]. One cup of oat and flour was 90 g and 120 g, respectively. The following ingredients (as g/345 g on dough basis) were used: oat (27-54 g), flour (60-84 g), brown sugar (25 g), coarse grain sugar (25 g), chocolate chips (30 g), cashew nuts (7 g), almond slices (10 g), egg (15 g), butter (125 g), and baking powder (3 g). Butter was sliced into small pieces and creamed with brown and coarse grain sugar at speed of 3 for 240 s in a heavy-duty mixer (KitchenAid, Model 5K5SS, MI, USA), using a flat beater and with scraping down every 60 s. Egg was then added, and the mixing continued at speed 4 for 30 s, and followed at speed 5 for 30 s, to form a soft and thick batter. The water contained in egg and butter is important in hydrating the starch particles, and can be absorbed by the protein to form the gluten of the dough. Egg is used as a source of liquid in the dough preparation, since water constitutes the major part of an egg [49].
In another mixing bowl, flour, oat, and baking powder were manually mixed. The amount of flour and oat required was based on the experimental design (Table 3) and mixture composition ( Table 6). Large inclusions were then added, comprising chocolate chips, cashew nuts, and almond slices, and homogeneously mixed with the flour-oat-baking powder mixture. The batter was scraped down and poured into the mixing bowl, followed by manual mixing for 120 s. The dough was then stored in a closed container prior to testing.

Cookie Dough Formulation
The amount of oat and flour required in preparing the cookie dough was calculated based on the cup percentage, as shown in Table 6. The oat to flour composition of 50O/50F indicated that the dough comprised 50% of a cup of oat and similarly 50% of flour. This indicated that half cup of oat (45 g) and half cup of flour (60 g) were required for the 50O/50F composition. Likewise, 30% cup of oat and 70% cup of flour was designed for the 30O/70F composition. The total number of cups for oat and flour were 0.3 cup (27 g) and 0.7 cups (84 g), respectively. Samples were prepared according to the experimental design in Table 3.

Stickiness Analysis
The value of stickiness was determined using a texture analyzer (TA.XT PLUS, Stable Micro Systems, Surrey, U.K.), with a 75 mm diameter cylinder probe (P/75P), and at pre-test speed of 0.5 mm/s, test speed of 0.5 mm/s, post-test speed of 10.0 mm/s, return distance of 5 mm, applied force of 5 g, contact time of 0.1 s, and trigger type of Auto-5 g [50]. Measurements were made on triplicate from each sample composition, with three measurements per replicate. The schematic diagram of the apparatus used to conduct the stickiness test is shown in Figure 1. from each sample composition, with three measurements per replicate. The schematic diagram of the apparatus used to conduct the stickiness test is shown in Figure 1. The TestXpert ® 10.11 software (Zwick GmbH & Co. KG, Ulm, Germany) was employed to calculate stickiness, as the energy needed to remove the working plate from the dough surface after compression. The dough stickiness value was obtained based on the maximum positive force on the force versus time curve. According to Agrahar-Murugkar and Dixit-Bajpai [51], the maximum force reading, which was the highest positive peak from the positive region, indicates the stickiness value.

Moisture Content Analysis
The moisture content of the cookie dough was examined using Moisture Analyzer (Model MB45, Othaus, US). A total of 2 g of sample was placed in a tray. The starting and end temperature was set at 105 ℃ and 120 ℃, respectively. Measurement was made on triplicate from each sample composition with three measurements per replicate.

Microstructure Analysis
The microstructure of sample which exhibited the highest and lowest value of stickiness and moisture content was studied, using the Hitachi S-3400N variable scanning electron microscope (SEM). The samples were mounted on an 8mm diameter stub, with a carbon conductive tape and sputter-coated with gold in a vacuum JSM-IT100 (InTouchScope, Hitachi, Japan). All samples were observed at 1.0 kV for 50x and 100x magnifications.

Results and Discussion
Data obtained from the experiments were analyzed using MINITAB software and results on the responses (stickiness and moisture content) are given in Table 7. Flour content and resting time are represented by V1 and V2, respectively. The TestXpert ® 10.11 software (Zwick GmbH & Co. KG, Ulm, Germany) was employed to calculate stickiness, as the energy needed to remove the working plate from the dough surface after compression. The dough stickiness value was obtained based on the maximum positive force on the force versus time curve. According to Agrahar-Murugkar and Dixit-Bajpai [51], the maximum force reading, which was the highest positive peak from the positive region, indicates the stickiness value.

Moisture Content Analysis
The moisture content of the cookie dough was examined using Moisture Analyzer (Model MB45, Othaus, US). A total of 2 g of sample was placed in a tray. The starting and end temperature was set at 105 • C and 120 • C, respectively. Measurement was made on triplicate from each sample composition with three measurements per replicate.

Microstructure Analysis
The microstructure of sample which exhibited the highest and lowest value of stickiness and moisture content was studied, using the Hitachi S-3400N variable scanning electron microscope (SEM). The samples were mounted on an 8 mm diameter stub, with a carbon conductive tape and sputter-coated with gold in a vacuum JSM-IT100 (InTouchScope, Hitachi, Japan). All samples were observed at 1.0 kV for 50x and 100x magnifications.

Results and Discussion
Data obtained from the experiments were analyzed using MINITAB software and results on the responses (stickiness and moisture content) are given in Table 7. Flour content and resting time are represented by V 1 and V 2 , respectively. Table 7. Design matrix and response value for the stickiness and moisture content test.

Sample
Coded Factor Uncoded Factor Response

Statistical Analysis of Stickiness
The ANOVA table is shown in Table 8. The results showed that the squared terms for both factors are highly significant at p-value of 0.000. However, the interaction term (interaction between different factors, V 1 *V 2 ) showed a p-value of 0.667 which indicates that there is no significant relationship between the factors.
The coefficient of determination (R 2 ) is 0.9644, which indicates that 96.44 percent of the sample variation in the stickiness was attributed to the factors V 1 and V 2 . The R 2 of 1.0000 means that the regression coefficient model is capable in predicting the optimum value with high accuracy. All main factors (V 1 and V 2 ) and their squared terms (V 1 *V 1 and V 2 *V 2 ) are highly significant at 0.000. The squared term of V 1 *V 1 was found to has the strongest effect on the response, due to the highest absolute coefficient value of 1073.78. The squared term, which is significant (p-value < 0.005), indicates that the relationship between the factor and the response follows a curved line. The regression coefficient model for the factors on the stickiness of cookie dough is given in Equation (2).
where Y s represents the stickiness (response) while V 1 and V 2 are the flour content and resting time, respectively. This regression coefficient model (mathematical model) can be used to calculate and analyze the effect of factors on the stickiness of the cookie dough.

Effect of Factors on Stickiness
From the analysis, ANOVA and regression coefficient model were used in analyzing the effect of factors on stickiness. Contour and surface plots, which showed the effect of flour content and resting time on the stickiness of cookie dough, were used for better illustration, and are shown in (b).
It was observed in Figure 2b that at all resting times, the stickiness value for the sample containing 70% of flour was lower than that of other samples. High stickiness was observed for flour content between 56% and 62%, and within 27 and 50 min for resting time. The results agreed with those reported by Yildiz et al. [52], who showed that low stickiness values were recorded under high flour content. As shown in Figure 2a, at the resting time of 40 min, the stickiness value increases when the flour content increases from 50% to 60%. Dough stickiness is generally associated with the development of gluten, and the interaction with other ingredients in the formula including sugar, fat, and water [53]. Thus, low flour content (conversely high oat content) indicates a low amount of gluten in the dough matrix, which resulted in a low stickiness value.
The further addition of flour from 60% to 70% resulted in a decrease in stickiness value from approximately 6814 N to 1524 N, respectively. The trend is consistent with those for other resting times. The results agreed with those of Yildiz et al. [52], where the dough was stickier at a longer resting time. Other studies likewise reported that the stickiness value increases with resting time. However, no systematic trend of stickiness values on the resting time was recorded in these studies [26]. During the production process, a longer contact time between dough and its contact surface may result in an increase in the stickiness value, due to an increase in the contact area [54]. The dough will maintain contact with the contact surfaces, create a bigger contact area, and subsequently increase surface wettability [55]. The further addition of flour from 60% to 70% resulted in a decrease in stickiness value from approximately 6814 N to 1524 N, respectively. The trend is consistent with those for other resting times. The results agreed with those of Yildiz et al. [52], where the dough was stickier at a longer resting time. Other studies likewise reported that the stickiness value increases with resting time. However, no systematic trend of stickiness values on the resting time was recorded in these studies [26]. During the production process, a longer contact time between dough and its contact surface may result in an increase in the stickiness value, due to an increase in the contact area [54]. The dough will maintain contact with the contact surfaces, create a bigger contact area, and subsequently increase surface wettability [55].
In Figure 2a, the highest stickiness value of flour content with 58% at all resting times. This could be associated with the presence of excess water that is not bound by proteins [26]. The gluten present in the flour absorbed water to induce protein interactions that play a critical role in dough stickiness [56]. For samples with flour content of 50%, 55%, and 60%, similar results were obtained by Tang and Liu [57]. Thus, they concluded that as the protein content increased, its hydration increased surface adhesion, and consequently effected higher dough stickiness. The surface tension between the dough and contact surface is developed by the mixture of water and water-soluble materials that dissolved [58]. It was also suggested that doughs made through mixing non-sticky dough flour with a sufficient amount of water-soluble fiber would increase their stickiness [50].
The decrease in the stickiness value of the sample with 65% to 70% of flour content, at all resting times, could be associated with the limited water availability in the protein complex. The hydration In Figure 2a, the highest stickiness value of flour content with 58% at all resting times. This could be associated with the presence of excess water that is not bound by proteins [26]. The gluten present in the flour absorbed water to induce protein interactions that play a critical role in dough stickiness [56]. For samples with flour content of 50%, 55%, and 60%, similar results were obtained by Tang and Liu [57]. Thus, they concluded that as the protein content increased, its hydration increased surface adhesion, and consequently effected higher dough stickiness. The surface tension between the dough and contact surface is developed by the mixture of water and water-soluble materials that dissolved [58]. It was also suggested that doughs made through mixing non-sticky dough flour with a sufficient amount of water-soluble fiber would increase their stickiness [50].
The decrease in the stickiness value of the sample with 65% to 70% of flour content, at all resting times, could be associated with the limited water availability in the protein complex. The hydration capability and complexing with ingredients play an important role in imparting stickiness to the dough [26]. The changes in dough stickiness over resting time and flour content may probably be due to the structural relaxation [59] and network degradation caused by enzyme activity after a longer resting time [60]. Figure 3 shows the micrograph images at 100x magnification after 30 min of resting time for sample S4 and S8, which produced the highest (6918 N) and lowest (1528 N) stickiness values, respectively.
Flour content for sample S4 and S8 were 60% and 70%, respectively. Oat starch granules can be determined by a much smaller size than wheat granules, and are irregular and polyhedral [61]. It can be observed that wheat and oat starch granules were distributed in the protein matrix. These results were consistent with those observed by Dachana et al. [62], which stated that a thin sheet representing the protein matrix, along with small and large embedded starch granules, can be observed in cookie dough. The starch granules, which were almost entirely covered by thin layers of protein film, were melted and re-solidified fat [63]. It can be seen in Figure 3b that the area covered with the protein matrix is larger than that shown in Figure 3a. In Figure 3a, the starch granules are partially embedded in the protein matrix, due to the surface adhesion induced by the protein. Hollows and ditches were observed on the dough surface, which indicates that the continuity of the gluten matrix has been disrupted by the protein [57]. Figure 3a shows sample with higher number of holes, and with a more uneven structure, as compared to that shown in Figure 3b. This observation is in agreement with results by Majzoobi et al. [64], which concluded that a higher percentage of oat content (lower content of flour) resulted in a more uneven surface structure, and with an increased number of holes. Figure 3b shows the micrograph image of the sample with the lowest stickiness value. The starch granules were mostly covered with a thin protein film, which was responsible for network formation, resulting in low stickiness properties and a good handling of dough [62,65].

Statistical Analysis of Moisture Content
The ANOVA table is shown in Table 9. The results showed that the main factors, which are flour content and resting time, are significant at p-value of 0.005 and 0.000, respectively. The interaction term (interaction between different factors, V 1 *V 2 ) showed a p-value of 0.000, which indicates that there was a very strong relationship between the factors.
The coefficient of determination (R 2 ) was 0.9891, which indicates that 98.91 percent of the sample variation in the stickiness, was attributed to the factors V 1 and V 2 . All terms were significant at p-value < 0.050, with the squared term of V 2 *V 2 being the least significant term. Resting time (V 2 ) was found to have the strongest effect on the response, due to the highest absolute coefficient value of 0.3081. All terms have negative effect on the moisture content, as indicated by the negative sign of the coefficient value. The squared term which is significant (p-value < 0.005) indicates that the relationship between the factor and the response follows a curved line. The regression coefficient model for the factors on the stickiness of cookie dough is given in Equation (3).
where Y m represents the moisture content (response), while V 1 and V 2 are the flour content and resting time, respectively. This regression coefficient model (mathematical model) can be used to calculate and analyze the effect of factors on the moisture content of the cookie dough.

Effect of Factors on Moisture Content
From the analysis, ANOVA and a regression coefficient model were used in analyzing the effect of factors on moisture content. Contour and surface plots, which showed the effect of flour content and resting time on the moisture content of cookie dough, were used for better illustration, and are shown in Figure 4.
Except for the sample with 50% flour content, the moisture content decreases when the resting time increases, as shown in Figure 4b. Similarly, a high moisture content was recorded at 10 min resting time, and flour content between 55% and 70%. The high region moisture content was observed for flour content between 60% and 70%, and within 10 and 15 min for resting time. Moisture reduction is probably due to evaporation from the cookie dough, which is affected by the environment. Although water is a minor component in cookie batter formula, its rheological behavior and machinability are largely influenced by the water content, and its distribution within the batter [66,67].
The effect of resting time and dough composition on the moisture content can be analyzed from the surface plot, as shown in Figure 4a. At 10 min resting time, the sample with 60% flour content (sample S3), which had the lowest stickiness value (Figure 2b), recorded the highest moisture content. At flour content of 70%, the moisture content decreased with an increase in resting time. Sample S13 recorded the lowest moisture content (10.69%) at 50 min of resting time. The high content of flour resulted in a dough with large surface area, which absorbed more moisture within the first 20 min of resting time. An increase in the surface area of starch led to higher water absorption [68]. Bushuk [69] showed that moisture was bounded variously to the flour constituents; namely 46% to starch, 31% to proteins, and 23% to pentosans.
At 30, 40, and 50 min of resting time, the samples with a high stickiness value (those with 60% flour content) exhibited comparably high moisture content. Based on this trend, it can be said that the lower the moisture content in the cookie dough, the lower the stickiness value. Similar results were also obtained by Milašinović Šeremešić et al. [70]. Higher water absorption provided better wetting properties, and the dough is in better contact with the solid contact surface. Thus, it gave higher surface adhesion, and subsequently increased dough stickiness [71,72]. The moisture absorption ability of the flour was found to be directly proportional to the protein level. An increase in protein level increased water absorption and stickiness properties [56].

Effect of Factors on Moisture Content
From the analysis, ANOVA and a regression coefficient model were used in analyzing the effect of factors on moisture content. Contour and surface plots, which showed the effect of flour content and resting time on the moisture content of cookie dough, were used for better illustration, and are shown in Figure 4.  Except for the sample with 50% flour content, the moisture content decreases when the resting time increases, as shown in Figure 4b. Similarly, a high moisture content was recorded at 10 min resting time, and flour content between 55% and 70%. The high region moisture content was observed for flour content between 60% and 70%, and within 10 and 15 min for resting time. Moisture reduction is probably due to evaporation from the cookie dough, which is affected by the environment. Although water is a minor component in cookie batter formula, its rheological behavior and machinability are largely influenced by the water content, and its distribution within the batter [66,67].
The effect of resting time and dough composition on the moisture content can be analyzed from the surface plot, as shown in Figure 4a. At 10 min resting time, the sample with 60% flour content (sample S3), which had the lowest stickiness value (Figure 2b), recorded the highest moisture content. At flour content of 70%, the moisture content decreased with an increase in resting time. Sample S13 recorded the lowest moisture content (10.69%) at 50 min of resting time. The high content of flour resulted in a dough with large surface area, which absorbed more moisture within the first 20 min of resting time. An increase in the surface area of starch led to higher water absorption [68]. Bushuk [69] showed that moisture was bounded variously to the flour constituents; namely 46% to starch, 31% to proteins, and 23% to pentosans.
At 30, 40, and 50 min of resting time, the samples with a high stickiness value (those with 60% flour content) exhibited comparably high moisture content. Based on this trend, it can be said that the lower the moisture content in the cookie dough, the lower the stickiness value. Similar results were also obtained by Milašinović Šeremešić et al. [70]. Higher water absorption provided better wetting properties, and the dough is in better contact with the solid contact surface. Thus, it gave  Figure 5 shows the microscopic images of sample S3 and S13 at 10 and 50 min of resting time, respectively. Both samples contained 60% of flour. The starch granules in Figure 5a are completely covered with the protein matrix, due to high moisture content. These results concurred with observations made by Létang et al. [73], which showed that starch granules were less visible in highly hydrated dough, due to a covering of continuous film. The dough with high moisture content, as shown in Figure 5a, entailed a highly aggregated protein phase, which is scattered, clustered, and less interconnected. In Figure 5b, the sample with the lowest moisture content possesses a slightly clustered, yet interconnected, network of the protein matrix. Fewer empty areas can be seen between the protein matrix, as compared with the more compact structure, as shown in Figure 5a.

Optimization of the Responses
The optimization plot is shown in Figure 6. The purpose of conducting optimization is to obtain the best combination of flour content and resting time in producing the cookie dough with the minimum value of stickiness and the maximum value moisture content. The composite desirability value, D, was calculated to be close to 1, therefore the factors were within the working range. The optimized values of factors were flour content (V 1 ) = 67% and resting time (V 2 ) = 10 min. As compared to the control sample with 100% wheat flour, the control sample recorded low stickiness and moisture content at 668N and 11.51%, respectively. These results are in agreement with other studies, which obtained low stickiness [34] and low moisture content [74] for a sample with 100% wheat flour.
clustered, yet interconnected, network of the protein matrix. Fewer empty areas can be seen between the protein matrix, as compared with the more compact structure, as shown in Figure 5a.

Optimization of the Responses
The optimization plot is shown in Figure 6. The purpose of conducting optimization is to obtain the best combination of flour content and resting time in producing the cookie dough with the minimum value of stickiness and the maximum value moisture content. The composite desirability value, D, was calculated to be close to 1, therefore the factors were within the working range. The optimized values of factors were flour content (V1) = 67% and resting time (V2) = 10 min. As compared to the control sample with 100% wheat flour, the control sample recorded low stickiness and moisture content at 668N and 11.51%, respectively. These results are in agreement with other studies, which obtained low stickiness [34] and low moisture content [74] for a sample with 100% wheat flour.

Experimental Validation
Experimental validation is the final step in the modelling process, and was used to verify the accuracy of the predicted model (regression coefficient model) [75]. A validation experiment was

Experimental Validation
Experimental validation is the final step in the modelling process, and was used to verify the accuracy of the predicted model (regression coefficient model) [75]. A validation experiment was carried out for three validation samples (SV1, SV2 and SV3) under the optimal conditions obtained from the optimization plot (see Figure 6). As shown in Table 10, the average relative deviation for stickiness and moisture content was 8.54% and 1.44%, respectively. This verified the predictability of the model with a comparison of the experimental (actual) values against the predicted figures (1720 N and 12.25%), implying that the RSM-based empirical model can adequately describe the relationship between the independent variables and the target response and, therefore, successfully reveal the optimum process conditions.

Contour Plots Superimposition
A method used to plot the overlaid graphs for various response surfaces is through the superimposition of the contour plots. This method is superior to the classical one-factor-at-a-time (OFAT) approach, which lacks the interaction of selected variables, and with experimental runs that are cumbersome [76]. The superimposed contour plots further serve as a reference template, which is easy to use in estimating the response for any given factor value within the range available. Figure 7 shows a feasible representation of the optimum range for the processing condition and formulation of cookie dough. The solid and dotted lines represent the contour line for stickiness and moisture content, respectively. Based on the overlaid contour plots, the optimum range for the lowest stickiness and highest moisture content (represented by the grey region) were found to be 67-70% flour content and 10-15 min resting time. carried out for three validation samples (SV1, SV2 and SV3) under the optimal conditions obtained from the optimization plot (see Figure 6). As shown in Table 10, the average relative deviation for stickiness and moisture content was 8.54% and 1.44%, respectively. This verified the predictability of the model with a comparison of the experimental (actual) values against the predicted figures (1720 N and 12.25%), implying that the RSM-based empirical model can adequately describe the relationship between the independent variables and the target response and, therefore, successfully reveal the optimum process conditions.

Contour Plots Superimposition
A method used to plot the overlaid graphs for various response surfaces is through the superimposition of the contour plots. This method is superior to the classical one-factor-at-a-time (OFAT) approach, which lacks the interaction of selected variables, and with experimental runs that are cumbersome [76]. The superimposed contour plots further serve as a reference template, which is easy to use in estimating the response for any given factor value within the range available. Figure  7 shows a feasible representation of the optimum range for the processing condition and formulation of cookie dough. The solid and dotted lines represent the contour line for stickiness and moisture content, respectively. Based on the overlaid contour plots, the optimum range for the lowest stickiness and highest moisture content (represented by the grey region) were found to be 67-70% flour content and 10-15 min resting time.

Conclusions
The effects of flour content and dough resting time on stickiness and moisture content were studied, and the statistical analysis using RSM and CCD showed that flour content and resting time

Conclusions
The effects of flour content and dough resting time on stickiness and moisture content were studied, and the statistical analysis using RSM and CCD showed that flour content and resting time have a significant effect on the stickiness and moisture content of the cookie dough. High stickiness was observed for flour content between 56% and 62%, and within 27 and 50 min for resting time. The lowest stickiness is associated with limited water available in the enlarged area of the protein matrix. Except for the sample with 50% flour content, the moisture content decreases when the resting time increases. Similarly, high moisture content was recorded at 10 min resting time and flour content between 55% and 70%. The high region moisture content was observed for flour content between 60% and 70%, and within 10 and 15 min of resting time. The sample with the lowest moisture content possessed a slightly clustered yet interconnected network of the protein matrix. The optimized values for factors were flour content (V 1 ) = 67% and resting time (V 2 ) = 10 min. The validation results showed the average relative deviation for stickiness and moisture content were 8.54% and 1.44%, respectively. The predicted model (regression coefficient model) thus developed was verified as highly accurate. The superimposition of the contour plots was found to successfully identify the optimum range for the lowest stickiness and highest moisture content, which were identified at 67-70% flour content and 10-15 min resting time. Several potential studies have been identified, including sensory evaluation and customer satisfaction. In addition, a study on different types of wheat flours may potentially improve the properties of cookie dough.