A response surface optimized USAEME followed by GC–MS for analysis of PAHs.
A quadratic model established to predict efficiency as a function of 3 variables.
GC–MS data analysis by MCR-ALS allowed comprehensive determination of PAHs.
Sensitivity and selectivity of analysis were improved using second-order calibration.
AFOM obtained by this method was superior over that of univariate calibration.