The use of experimental design and in particular the response surface methodology (RSM) allowed the determination of the influence of the simultaneous effects and the interaction of the operating parameters on the methylene blue removal efficiency. The parameters studied were the initial concentration of the adsorbate, the stirring speed, mass and particle size of the adsorbent. The results show that the application of RSM allows describing the influence of these four experimental parameters on the treatment effectiveness. The second-order model obtained, for the Methylene Blue (MB) removal efficiency was validated by using different statistical approaches. The use of the ANOVA showed that the model is significant and in functional adequacy with the experimental results.