Owing to the increasing number of erythropoietin biosimilars being approved, the comparison of different erythropoietin preparations in the pharmaceutical area is gaining in importance. Erythropoietin has a distinct natural heterogeneity arising from its glycosylation, which shows strong composition variations. This heterogeneity increases the complexity of the analysis of erythropoietin considerably, but may also be used to distinguish different preparations. Here, a method is presented for the differentiation of various erythropoietin preparations by capillary electrophoresis-mass spectrometry and the subsequent application of multivariate statistics. Relative peak areas of selected intact erythropoietin isoforms were used as variables in principal component analysis and hierarchical agglomerative clustering. Both of these strategies were suited for the clear differentiation of all erythropoietin preparations, including marketed products and preproduction preparations, which differ in the manufacturer, the production cell line, and the batch number. By this means, even closely related preparations were distinguished on the basis of the combined information on the antennarity, the sialoform, and the acetylation of the observed isoforms.