Response surface methodology for split-plot designs with categorical factors