In this study, the fit to the experimental data of two second order models and two phase model were compared with first order models. The best fit was observed with the two parameter model with rectangular distribution of floatabilities (model II) and the three parameter fast/slow floating particles model (model I). Also, the three parameter gamma distribution model (model VI) showed good fit to the experimental data. The worst fit was observed with first order two stage kinetic model (model V) which indicates that there is no need to divide the flotation rate into two flotation rates. Model VII which includes six parameters gave acceptable fit to experimental data. This result is in contrast with the results obtained by other investigators which states that increase of parameters in model leads to parameter dilution and increase on model error. Hence, it may be concluded that the number of parameters to obtain an adequate model needs more investigation. The confidence intervals for the three best models (model I, model II and model VI) were estimated at 95% confidence level. Also, model II shows more discrete regions with regard to flotation rate and recovery than model VI. This result indicates the ability of such model to express the changes in flotation reagents and/or flotation conditions over other tested models.