The stochastic theory of chromatography and an equilibrium based approach were used for the prediction of peak shape and retention data of anions. This attempt incorporating the potential advantages of two different chromatographic phenomena for analytical purposes. It is an integrated method to estimate kinetic and thermodynamic properties for the same chromatographic run of ions. The stochastic parameters of eluted anions, such as the residence time of the molecule on the surface of the stationary phase, and the average number of adsorption steps were determined on the basis of a retention database of organic and inorganic anions (formate, chloride, bromide, nitrate, sulphate, oxalate, phosphate) obtained by using carbonate/bicarbonate eluent system at different pHs (9-11) and concentrations (7-13 mM). In the investigated IC system the residence times are much higher and the average number of sorption steps is somewhat smaller than in RP-HPLC. The simultaneous application of the stochastic and the multispecies eluent/analyte model was utilized to peak shape simulation and the retention controlling of various anions under elution conditions of practical importance. The similarities between the measured and the calculated chromatograms indicates the predictive and simulation power of the combined application of the stochastic theory and the multiple species eluent/analyte retention model.