The extraction of Rhodiola rosea rhizomes using natural deep eutectic solvent (NADES) consisting of lactic acid, glucose, fructose, and water was investigated. A two-level Plackett-Burman design with five variables, followed by the steepest ascent method, was undertaken to determine the optimal extraction conditions. Among the five parameters tested, particle size, extraction modulus, and water content were found to have the highest impact on the extrability of phenyletanes and phenylpropanoids. The concentration of active compounds was analyzed by HPLC. The predicted results showed that the extraction yield of the total phenyletanes and phenylpropanoids (25.62 mg/g) could be obtained under the following conditions: extraction time of 154 min, extraction temperature of 22 • C, extraction modulus of 40, molar water content of 5:1:11 (L-lactic acid:fructose:water, mol/mol), and a particle size of rhizomes of 0.5-1 mm. These predicted values were further verified by validation experiments in predicted conditions. The experimental yields of salidroside, tyrosol, rosavin, rosin, cinnamyl alcohol and total markers (sum of phenyletanes and phenylpropanoids in mg/g) were 11.90 ± 0.02, 0.36 ± 0.02, 12.23 ± 0.21, 1.41 ± 0.01, 0.20 ± 0.01, and 26.10 ± 0.27 mg/g, respectively, which corresponded well with the predicted values from the models.