This article is part of the series Proceedings of the International Congress in Honour of Professor Hari M. Srivastava.

Open Access Research

Comparison of hierarchical cluster analysis methods by cophenetic correlation

Sinan Saraçli1*, Nurhan Doğan2 and İsmet Doğan2

Author Affiliations

1 Department of Statistics, Faculty of Arts and Sciences, Afyon Kocatepe University, Afyonkarahisar, 03200, Turkey

2 Department of Biostatistics, Faculty of Medicine, Afyon Kocatepe University, Afyonkarahisar, 03200, Turkey

For all author emails, please log on.

Journal of Inequalities and Applications 2013, 2013:203  doi:10.1186/1029-242X-2013-203


The electronic version of this article is the complete one and can be found online at: http://www.journalofinequalitiesandapplications.com/content/2013/1/203


Received:31 December 2012
Accepted:10 April 2013
Published:23 April 2013

© 2013 Saraçli et al.; licensee Springer

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M1">View MathML</a> and the second one is with outliers (5%) for <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M2">View MathML</a>. The proposed method is applied to simulated multivariate normal data via MATLAB software.

Results

According the results of simulation the Average (especially for <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M3">View MathML</a>) and Centroid (especially for <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M4">View MathML</a> and <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M5">View MathML</a>) methods are recommended at both conditions.

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 measures

1 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, city-block (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 pair-group average (UPGMA), Weighted pair-group average (WPGMA), Unweighted pair-group centroid (UPGMC), Weighted pair-group centroid (median), Ward’s method, etc.) are available; many of these begin with an n-dimensional 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 <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M6">View MathML</a> have been modeled using a cluster method to produce a dendrogram <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M7">View MathML</a>; that is, a simplified model in which data that are ‘close’ have been grouped into a hierarchical tree. Define the following distance measures. <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M8">View MathML</a>, the ordinary Euclidean distance between the ith and jth observations. <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M9">View MathML</a> the dendrogrammatic distance between the model points <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M10">View MathML</a> and <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M11">View MathML</a>. This distance is the height of the node at which these two points are first joined together. Then, letting x be the average of the <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M12">View MathML</a>, and letting t be the average of the <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M13">View MathML</a>, the cophenetic correlation coefficient c is defined as in (1) [9].

<a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M14">View MathML</a>

(1)

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 cluster-based 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 Goodman-Kruskal 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 Goodman-Kruskal 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: k-means, SOM, EM-based clustering, classification EM clustering, fuzzy k-means, 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 (<a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M3">View MathML</a>, <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M4">View MathML</a> and <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M5">View MathML</a>), variables number (<a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M18">View MathML</a>, <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M19">View MathML</a> and <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M20">View MathML</a>) and distance measures via a simulation study. The simulation program is developed in a MATLAB software development environment by the authors. We have 567 different simulation scenarios and 100,000/n replications for each scenario. The performance is monitored by two different conditions that are mentioned in Table 1 and Table 2 with 7 cluster methods, 9 distance measures by cophenetic correlation coefficient in various settings of subgroup means, variances, sample size and variable numbers simultaneously.

Table 1. The cophenetic correlation coefficient values for<a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M21">View MathML</a>,<a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M22">View MathML</a>(without outliers)

Table 2. The cophenetic correlation coefficient values for<a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M21">View MathML</a>,<a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M36">View MathML</a>(with outliers)

For 567 different simulation scenarios, the data was derived from multivariate normal distribution for <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M49">View MathML</a>, <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M50">View MathML</a> with and without outliers, respectively. The data set for outliers is obtained according to Dixon’s [22] ‘Outlier Model’ like <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M51">View MathML</a>. In this study, <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M52">View MathML</a> means that while 95% of the data set does not include any outliers, 5% of the data set includes outliers.

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 <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M3">View MathML</a>. Moreover, increasing the sample size to <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M4">View MathML</a> and <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M5">View MathML</a> favors the complete, weighted, and centroid methods for all measures. However, the cophenetic correlation coefficient for the Mahalanobis measure cannot be calculated in both distributions when there are 10 variables with sample size <a onClick="popup('http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.journalofinequalitiesandapplications.com/content/2013/1/203/mathml/M3">View MathML</a>, whereas there is not any meaningful explanation for this unexpected result, we still could not find the main reason for this situation, but the same result is obtained for more than three times run of the simulation program.

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 decision-makers 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

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

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

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

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

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

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

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

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

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

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

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

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

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

  14. Hubert, L: Approximate evaluation techniques for the single-link and complete-link hierarchical clustering procedures. J. Am. Stat. Assoc.. 69, 698–704 (1974). Publisher Full Text OpenURL

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

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

  17. 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 OpenURL

  18. 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 OpenURL

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

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

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

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

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