Evaluation of model adequacy There are many statistical techniques for the evaluation of model adequacy, but graphical residual analysis is the primary statistical method for assessment of model adequacy [59]. The normal probability plot indicates that the points on this plot are formed a nearly linear pattern (Fig. 2 (a)). Therefore, the normal distribution is a good model for this data set. Random scattering of the points of internally studentized residual (the residual divided by the estimated standard deviation of that residual) versus predicted values between −3 and +3 emphasizes highly accurate prediction of the experimental data through the derived quadratic model (Fig. 2 (b)). Fig 2 The plot of: (a) Normal Plot of Residuals; (b) residuals vs Predicted Response; (c) predicted vs Actual values; (d) Residuals vs Run Order The plot of predicted vs Actual values (Fig. 2 (c)) indicate a higher correlation and low differences between actual and predicted values. Hence, the predictions of the experimental data by developed quadratic models for the TC degradation is perfectly acceptable and this model fits the data better. Also, the random spread of the residuals across the range of the data between −3 and +3 implies that there are no evident drift in this process and the model was a goodness fit (Fig. 2 (d)). The Box-cox plot is used for determine the suitability of a power low transformation for the selected data (Fig. 3 (a)). In this study, the best lambda values of 0.92 was obtained with low and high confidence interval 0.73 and 1.11, respectively. Therefore, recommend the standard transformation by the software is 'None'. The plot of points Leverage vs Run order is shown in Fig. 3 (b). The factorial and axial points have the most influence with a leverage of approximately 0.59, while the center points have the least effect with a leverage of 0.16. Fig 3 The plot of: (a) Box-Cox Plot for Power Transforms; (b) The points Leverage vs Run order for the CCD Design