It should be mentioned that the data obtained for COD reduction could not be fitted well to the models including linear, two factorial,quadratic and cubic. It seems that using a suitable experimentaldesign or employing other nonlinear modeling techniques such asneural networks, fuzzy logic or a combination of these approaches could be useful for this aim