Statistical methodology for the analysis of proteomic mass spectrometry data is proposed using mixed effects models. Each high dimensional spectrum is represented by using a near orthogonal low dimensional representation with a basis of Gaussian mixture functions. Linear mixed effect models are proposed in the lower dimensional space. In particular, differences between groups are investigated by using fixed effect parameters, and individual variability of spectra is modelled by using random effects. A deterministic peak fitting algorithm provides estimates of the near orthogonal Gaussian basis. The mixed effects model is fitted by using restricted maximum likelihood, and a parallel fitting procedure is used for computational convenience. The methodology is applied to proteomic mass spectrometry data from serum samples from melanoma patients who were categorized as stage I or stage IV, and significant locations of peaks are identified. Copyright (c) 2010 Royal Statistical Society.