Near-infrared (NIR) diffuse reflectance spectra have been measured by use of a rotating drawer for pellets of 12 kinds of ethylene/vinyl acetate (EVA) copolymers with vinyl acetate (VA, the comonomer) varying in the 7-44 wt % range. They are unambiguously discriminated from one another by a score plot of the principal component analysis (PCA) Factor 1 and 2, based upon the NIR spectra pretreated by multiplicative scatter correction (MSC). Principal component (PC) weight loadings for Factor 1 show that the discrimination relies largely upon bands due to the overtone and combination modes arising from the VA unit. We have found one ''outlier'' in the score plot and elucidated its spectral characteristics based upon PC weight loadings for Factor 2. Partial least-squares (PLS) regression has been applied to propose calibration models which predict the VA content in EVA. The models have been prepared for three kinds of pretreatment, the first derivative, the second derivative, and MSC; and four kinds of wavelength regions. The NIR spectra in the 1100-2200 nm region after the MSC treatment has given the best correlation coefficient and standard error of prediction (SEP) of 0.998 and 0.70%, respectively. The calibration models, prepared by NIR diffuse reflectance spectroscopy for the pellet samples, are compared with previously reported models by NIR transmission spectroscopy for the flowing molten samples, and with those by Raman spectroscopy for the pellet samples. PLS regression has also allowed us to predict melting points of the copolymers with the correlation coefficient and SEP of 0.997 and 0.78ะC, respectively.