Karyotypic variations in five lily taxa (Lilium candidum, L. martagon, L. ciliatum, L. carniolicum ssp. ponticum and L. carniolicum var. ponticum) have been investigated by means of numerical methods (Cluster and Principal Component Analyses). By the analysis of this species from six chromosome related characters, it has been determined that relative length is the best character which represents the variations in them. It was also found that the results from numerical analysis of the karyologic characters provide distinct evidences, which are corresponding to the morphology for recognition of the taxa.