The calyxes of Hibiscus sabdariffa are used in traditional medicine around the world. However, quality assurance protocols and chemical variability have not been previously analyzed. In the present study, chemical characterization of a set of samples of H. sabdariffa calyxes commercialized in Colombia was accomplished with the aim to explore the chemical variability among them. Chemometrics-based analyses on the data obtained from the HPLC-UV-DAD-derived profiles were then performed. Thus, the pre-processed single-wavelength data were subjected to principal component analysis (PCA). The PCA-derived results evidenced different groups which were well-correlated to the corresponding total phenolic and total anthocyanin contents. Multi-wavelength chromatographic (HPLC-UV-DAD surfaces) data were additionally examined via parallel factor analysis (PARAFAC) as data reduction method and the obtained loadings were subsequently submitted to PCA and orthogonal partial least squares discriminant analysis (OPLS-DA). Results were thus consistent with those from single-wavelength data. PCA loadings were employed to determine those chemical components responsible for the data variance and OPLS-DA model, constructed from PARAFAC loadings, and indicated differentiation according total anthocyanin contents among samples. The present chemometric analysis therefore demonstrated to be an excellent tool for differentiation of H. sabdariffa calyxes according to their chemical composition.