The chemometrics approach was applied for the separation optimization of flavonoid markers (quercetin, hesperetin and chrysin) in honey using micellar liquid chromatography (MLC). The investigated method combines SPE of flavonoids from honey using C(18) cartridge and their separation and quantification by micellar liquid chromatography. A two level full factorial design was carried out to evaluate the effect of four experimental factors including concentration of SDS, alkyl chain length of the alcohol used as the organic modifier (N), volume percentage of the organic modifier (V(m)) and volume percentage of acetic acid (AcOH) in mobile phase on analytes retention times. Experiments for analytes retention times modeling and optimization of separation were performed according to central composite design. Multiple linear regression method was used for the construction of the best model based on experimental retention times. Pareto optimal method was used to find suitable compatibility between resolution and analysis time of analytes in honey. The optimum mobile phase composition for separation and determination of analytes in honey were [SDS]=0.124 mol/L; 7.8% v/v ethanol and 5.0% v/v AcOH. Limits of detection and linear range of flavonoid markers were 0.0079-0.0126, 0.05-50.0 mg/L, respectively.