An off-line system is proposed consisting of ion-pair reversed-phase liquid chromatography, collections of fractions at the outflow of the column and furnace atomic absorption spectrometry. The so-called system allowed determination of both arsenic and selenium species mainly found in the environment and in mammals (arsenite, arsenate, monomethylarsonate, dimethylarsinate, selenite, selenate, selenocystamine, selenocystine, selenomethionine and selenoethionine). In order to study the retention behaviour of these compounds and to estimate the optimal conditions for the chromatographic separation, central composite designs were used to evaluate the influence of the eluent parameters such as pH, tetrabutylammonium phosphate (TBA) concentration and sodium hydrogenphosphate amounts. The retention factors of each species and the selectivity were established as response criteria. Response surfaces and isoresponse curves were drawn from the mathematical models and enabled one to determine the optimal conditions and to visualise the method robustness. The predicted optimal zone was situated at pH 5.5-6.5, 4.0 mM Na2HPO4 and 3.0-4.0 mM TBA. Regression models suggested linearity for the studied compounds in the range 25-200 microg selenium and arsenic per litre investigated.