Volatile compounds from 90 Italian citrus honeys of different varieties (lemon, Citrus limon (L.) Burm., orange, Citrus sinensis (L.) Osbeck, and citrus, Citrus spp.), 30 provided by CRA-API (Consiglio per la Ricerca e la Sperimentazione in Agricoltura-Unita di ricerca di apicoltura e bachicoltura) (Bologna, Italy) and 60 purchased on the markets, were analysed to classify them according to the characteristic aroma profiles. A HS-SPME-GC method previously developed and validated was applied. Factor analysis and general discriminant analysis (GDA) were carried out on peaks and intensity data. Statistical analyses revealed a good separation between samples; the best model based on peaks data shows 66.5% and 33.5% of total variance explained by the first and the second canonical functions and a predictive capacity of 78.9%. Regarding the intensity data, the CFs of the GDA model explain 53.8% and 46.2% of the total variance and the crossvalidation shows a predictive capacity of 96.7%.