Fly ash-based geopolymer binders have been identified as one of the alternatives to the Ordinary Portland Cement (OPC), which qualify the criteria of green construction material. In the process of enhancement of the properties of geopolymer; engineered geopolymer composite (EGC) is a recent development that is classified as the high-performance fiber reinforced geopolymer matrix. The philosophy of the development of EGC is to achieve high compressive strength and ductility. In this paper, statistical models are developed to predict the mechanical and post-cracking properties of EGC using Response Surface Methodology (RSM). In this regard, effects of three principal variables; molarity of sodium hydroxide, sodium silicate to sodium hydroxide ratio, and curing temperature on the properties of fresh and hardened EGC (setting time, compressive strength, elastic modulus, flexural strength, flexural toughness, ductility index, tensile strength, tensile first crack strength and tensile strain capacity) were investigated. All models are found significant because, and the difference between predicted R2 and adjustable R2 was found less than 0.2. The optimized mixture proportion was assessed using multi-objective optimization technique as discussed in the RSM literature. Finally, an experimental program was performed to validate the optimized mixture proportion. The predicted and experimental results were found in a close agreement because the variation between the predicted and the experimental results was obtained less than 5. The proposed method can be performed for any objective value of the EGC property with desirability almost equal to one, improving the yield, the reliability of the product and the processes. © 2018 Elsevier Ltd