Abstract
Purpose
This study proposes the best clustering method(s) for different distance measures under two different conditions using the cophenetic correlation coefficient.
Methods
In the first one, the data has multivariate standard normal distribution without outliers
for
Results
According the results of simulation the Average (especially for
Conclusions
This study hopes to contribute to literature for making better decisions on selection of appropriate cluster methods by using subgroup sizes, variable numbers, subgroup means and variances.
Keywords:
cophenetic correlation; hierarchical clustering methods; distance measures1 Introduction
Classification, in its widest sense, has to do with forms of the relatedness and with the organization and display of the relations in a useful manner. The items to be studied could be anything: people, bacteria, religions, books, etc. The attributes in each case would be those features of the items that are of interest for the purpose of the study [1]. Classifications are generally pictured in the form of hierarchical trees, also called a dendrogram. A dendrogram is the graphical representation of an ultrametric (= cophenetic) matrix; so dendrograms can be compared to one another by comparing their cophenetic matrices [2].
Cluster Analysis (CA), Principal Components Analysis (PCA) and Discriminant Analysis (DA) are three of the primary methods of modern multivariate analysis. Because of its utility, clustering has emerged as one of the leading methods of multivariate analysis [3].
Cluster analysis is a multivariate statistical technique which was originally developed for biological classification. Biologists Robert Soka1 and Peter Sneath published their seminal text ‘Principles of Numerical Taxonomy’ in 1963. Sokal and Sneath demonstrated that cluster analysis could be utilized to efficiently classification a data set which contained all relevant characteristics of an organism. When the organisms had been classified based on these characteristics, it could be determined in which way they differed, and if they belonged to different species. In this way, Sokal and Sneath asserted, researchers could trace the path of evolution from one species to another [4].
In this study for clustering, two measures of cluster ‘goodness’ or quality are used. One type of measure allows us to compare different sets of clusters without reference to external knowledge and is called an internal quality which is used as a measure of ‘overall similarity’ based on the pairwise similarity of documents in a cluster. The other type of measures allows evaluating how well the clustering is working by comparing the groups produced by the clustering techniques to known classes. This type of measure is called an external quality measure, which is not scope of this study [5].
The joining or tree clustering method uses the dissimilarities (similarities) or distances (Euclidean distance, squared Euclidean distance, cityblock (Manhattan) distance, Chebychev distance, power distance, Mahalanobis distance, etc.) between objects when forming the clusters. Similarities are a set of rules that serve as criteria for grouping or separating items. These distances (similarities) can be based on a single dimension or multiple dimensions, with each dimension representing a rule or condition for grouping objects. The joining algorithm does not ‘care’ whether the distances that are ‘fed’ to it are actual real distances, or some other derived measure of distance that is more meaningful to the researcher; and it is up to the researcher to select the right method for his/her specific application [6].
The next step is to identify how one can find the natural clusters among items characterized by many attributes. A number of cluster analysis procedures (single linkage (nearest neighbor), Complete linkage (furthest neighbor), Unweighted pairgroup average (UPGMA), Weighted pairgroup average (WPGMA), Unweighted pairgroup centroid (UPGMC), Weighted pairgroup centroid (median), Ward’s method, etc.) are available; many of these begin with an ndimensional space in which each entity is represented by a single point. The dimensions in the space represent the characteristics upon which the entities are to be compared. Similarity between entities can be measured by: (1) the correlation of entities’ scores on the dimensions (cophenetic correlation) or (2) the distance between points in the space (points closest to each other are most similar) [7,8].
Suppose that the original data
Since its introduction by Sokal and Rohlf [10], the cophenetic correlation coefficient has been widely used in numerical phenetic studies, both as a measure of degree of fit of a classification to a set of data and as a criterion for evaluating the efficiency of various clustering techniques [11]. In statistics, and especially in biostatistics, cophenetic correlation (more precisely, the cophenetic correlation coefficient) is a measure of how faithfully a dendrogram preserves the pairwise distances between the original unmodeled data points. Although it has been most widely applied in the field of biostatistics (typically to assess clusterbased models of DNA sequences, or other taxonomic models), it can also be used in other fields of inquiry where raw data tend to occur in clumps, or clusters. This coefficient has also been proposed for use as a test for nested clusters [12].
The problem of comparing classifications with numerical methods is not new; the first effective numerical method known to us is the ‘cophenetic correlation’ technique of Sokal and Rohlf [10]. Beginning with the development of cophenetic correlations methods for comparison of dendrograms have recently been the object of strong interest. Baker [13] investigated the impact of observational errors on the dendrograms produced by the complete linkage and single linkage hierarchical grouping techniques. The goodness of fit of the dendrograms was measured by means of the GoodmanKruskal gamma coefficient. The gamma coefficients indicated that the single linkage grouping technique was more sensitive to the type of data errors employed than the complete linkage technique. Hubert [14] compared two rank orderings of the object pairs. He tested hypothesis that the given set of proximity values have been assigned randomly by referring the GoodmanKruskal rank correlation γ statistic to an approximate permutation distribution. Kuiper and Fisher [15] compared six hierarchical clustering procedures (single linkage, complete linkage, median, average linkage, centroid and Ward’s method) for multivariate normal data, assuming that the true number of clusters was known. The authors used the Rand index, which gives a proportion of correct groupings, to compare the clustering methods. In their study for clusters of equal sizes, Ward’s method and complete linkage method, with very unequal cluster sizes centroid and average linkage method found best, respectively. Blashfield [16] compared four types of hierarchical clustering methods (single linkage, complete linkage, average linkage and Ward’s method) for accuracy in recovery of original population clusters. He used Cohen’s statistic to measure the accuracy of the clustering methods. According to his results, Ward’s method performed significantly better than the other clustering procedures and average linkage gave relatively poor results. According to Milligan [17], complete linkage and Ward’s method reacted badly when outliers were introduced into the simulated data.
Hands and Everitt [18] compared five hierarchical clustering techniques (single linkage, complete linkage, average, centroid, and Ward’s method) on multivariate binary data. They found that Ward’s method was the best overall than other hierarchical methods. Yao [19] discussed six classical clustering algorithms: kmeans, SOM, EMbased clustering, classification EM clustering, fuzzy kmeans, leader clustering and different combination scenarios of these algorithms. He used a count of cluster categories, classification accuracy and cluster entropy. Ferreira and Hitchcock [20] compared the performance of four major hierarchical methods (single linkage, complete linkage, average linkage and Ward’s method) for clustering functional data. They used the Rand index to compare the performance of each clustering method. According to their study, Ward’s method was usually the best, while average linkage performed best in some special situations, in particular, when the number of clusters is over specified. Milligan and Cooper [21] used four agglomerative hierarchical clustering methods to generate partition solutions and formed one factor in the overall design. These were the single link, complete link, group average (UPGMA) and Ward’s minimum variance methods. As a result, they found that the single link technique was least effective while the group average and Ward’s methods gave the best overall recovery.
Consider the studies in the literature and the importance of using the most convenient cluster method under different conditions (sample size, variables number and distance measures), a detailed simulation study is undertaken. This study gives more insight into the functioning of the cluster method under different conditions. The purpose of this research is to investigate the best clustering method under different conditions.
2 Method
In this study, seven cluster analysis methods are compared by the cophenetic correlation
coefficient computed according to different clustering methods with a sample size
(
Table 1. The cophenetic correlation coefficient values for
Table 2. The cophenetic correlation coefficient values for
For 567 different simulation scenarios, the data was derived from multivariate normal
distribution for
3 Results and discussion
All numerical results, obtained by running the simulation program, are given in Table
1 and Table 2. According to Table 1 and Table 2, the average method gives the best results at all measures and at all variable numbers
for both distributions with sample size
4 Conclusion
In general, researchers especially nonstatisticians use cluster analysis methods and distance measures in different conditions. In addition, they choose to use the most famous cluster analysis methods and distance measures, which are available in statistical packages, without evaluating the validity of different conditions. When the different conditions are considered, drawn inferences are dubious, and may lead the decisionmakers to incorrect decisions. It is noted that, with respect to the selection of a distance measures, the researcher must be aware that their choice can often significantly affect the results of the clustering. For example, some distance measures are inappropriate when different conditions of the variables are not met. On this point, the determination of the correct distance measures to use under various cases is the main motivation of researchers working on this subject to determine which distance measures should be used in case of different conditions.
One may conclude that the results of this study, which is similar to findings of Johnson and Wichern [23], indicate the data set with outliers have higher cophenetic correlation values than the data set without outliers.
This study hopes to contribute to literature for making better decisions on selection of appropriate cluster methods by using subgroup sizes, variable numbers, subgroup means and variances.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
SS has made intellectual contributions in order to carry out this study and also has carried out the simulation study. ND has determined the research design as well as has coordinated the whole process. İD has made theoretical contributions and has performed statistical analysis of the study. All authors read and approved the final manuscript.
Acknowledgements
Dedicated to Professor Hari M Srivastava.
The authors would like to thank Rıdvan ÜNAL for support of technical help. He is a lecturer at the Afyon Kocatepe University, Faculty of Science, Department of Physics, Afyonkarahisar/Turkey.
References

Carmichael, JW, George, JA, Julius, RS: Finding natural clusters. Syst. Zool.. 17(2), 144–150 (1968). Publisher Full Text

Lapointe, FJ, Legendre, P: Comparison tests for dendrograms: a comparative evaluation. J. Classif.. 12, 265–282 (1995). Publisher Full Text

Kettenring, JR: The practice of cluster analysis. J. Classif.. 23, 3–30 (2006). Publisher Full Text

Gunnarsson, J: PortfolioBased Segmentation and Consumer Behaviour: Empirical Evidence and Methodological Issues. Ph.D. Dissertation, Stockholm School of Economics, The Economic Research Institute, p. 274 (1999)

Steinbach, M, Karypis, G, Kumar, V: A comparison of document clustering techniques. Text mining workshop. Boston, MA. (2000)

Hill, T, Lewicki, P: STATISTICS: Methods and Applications, StatSoft, Tulsa (2007)

Lessig, VP: Comparing cluster analyses with cophenetic correlation. J. Mark. Res.. 9(1), 82–84 (1972). Publisher Full Text

Sneath, HA, Sokal, RR: Numerical Taxonomy: The Principles and Practice of Numerical Classification,p. 573. Freeman, San Francisco (1973)

Mathworks statistics toolbox: http://www.mathworks.com/help/stats/cophenet.html (2012)

Sokal, RR, Rohlf, FJ: The comparison of dendrograms by objective methods. Taxon. 11, 33–40 (1962). Publisher Full Text

Farris, JS: On the cophenetic correlation coefficient. Syst. Zool.. 18(3), 279–285 (1969). Publisher Full Text

Rohlf, FJ, David, LF: Test for hierarchical structure in random data sets. Syst. Zool.. 17, 407–412 (1968). Publisher Full Text

Baker, FB: Stability of two hierarchical grouping techniques  case I: sensitivity to data errors. J. Am. Stat. Assoc.. 69, 440–445 (1974)

Hubert, L: Approximate evaluation techniques for the singlelink and completelink hierarchical clustering procedures. J. Am. Stat. Assoc.. 69, 698–704 (1974). Publisher Full Text

Kuiper, FK, Fisher, LA: A Monte Carlo comparison of six clustering procedures. Biometrics. 31, 777–783 (1975). Publisher Full Text

Blashfield, RK: Mixture model tests of cluster analysis: accuracy of four agglomerative hierarchical methods. Psychol. Bull.. 83, 377–388 (1976)

Milligan, GW: An examination of the effect of six types of error perturbation on fifteen clustering algorithms. Psychometrika. 45, 325–342 (1980). Publisher Full Text

Hands, S, Everitt, B: A Monte Carlo study of the recovery of cluster structure in binary data by hierarchical clustering techniques. Multivar. Behav. Res.. 22, 235–243 (1987). Publisher Full Text

Yao, KB: A comparison of clustering methods for unsupervised anomaly detection in network traffic. Ph.D. Thesis, University of Copenhagen (2006)

Ferreira, L, Hitchcock, DB: A comparison of hierarchical methods for clustering functional data. Commun. Stat., Simul. Comput.. 38, 1925–1949 (2009). Publisher Full Text

Milligan, GW, Cooper, MC: A study of standardization of variables in cluster analysis. J. Classif.. 5, 181–204 (1988). Publisher Full Text

Dixon, WJ: Analysis of extreme value. Ann. Math. Stat.. 21, 488–506 (1950). Publisher Full Text

Johnson, RA, Wichern, DW: Applied Multivariate Statistical Analysis, Prentice Hall, New York (2002)