To establish the near infrared measurement model of nicotine content in e-liquid, 63 samples with different nicotine contents and tastes were collected. After eliminating the abnormal samples via Monte Carlo cross validation, the interval partial least squares (iPLS), and synergy interval partial least squares (SiPLS) methods were used to select the optimal wavelength, and subsequently, the correction model was established. The results show that iPLS identified the characteristic wavelength bands of 1090?1228 nm and 1370?1508 nm. By combining the different bands, SiPLS determined the optimal bands as 1126?1240 nm, 1358?1414 nm, and 1474?1530 nm. When compared with the full-spectrum PLS model, the variables used in the SiPLS-PLS model reduced by two-thirds and the mean square error of prediction value reduced from 1.188 to 0.963; thus, improving the accuracy of the model.