 Research
 Open Access
 Published:
Divisive hierarchical maximum likelihood clustering
BMC Bioinformatics volume 18, Article number: 546 (2017)
Abstract
Background
Biological data comprises various topologies or a mixture of forms, which makes its analysis extremely complicated. With this data increasing in a daily basis, the design and development of efficient and accurate statistical methods has become absolutely necessary. Specific analyses, such as those related to genomewide association studies and multiomics information, are often aimed at clustering subconditions of cancers and other diseases. Hierarchical clustering methods, which can be categorized into agglomerative and divisive, have been widely used in such situations. However, unlike agglomerative methods divisive clustering approaches have consistently proved to be computationally expensive.
Results
The proposed clustering algorithm (DRAGON) was verified on mutation and microarray data, and was gauged against standard clustering methods in the literature. Its validation included synthetic and significant biological data. When validated on mixedlineage leukemia data, DRAGON achieved the highest clustering accuracy with data of four different dimensions. Consequently, DRAGON outperformed previous methods with 3,4 and 5dimensional acute leukemia data. When tested on mutation data, DRAGON achieved the best performance with 2dimensional information.
Conclusions
This work proposes a computationally efficient divisive hierarchical clustering method, which can compete equally with agglomerative approaches. The proposed method turned out to correctly cluster data with distinct topologies. A MATLAB implementation can be extraced from http://www.riken.jp/en/research/labs/ims/med_sci_math/ or http://www.alokailab.com
Background
In unsupervised clustering algorithms, the class label or the state of nature of a sample is unknown. The partitioning of data is then driven by considering similarity or distance measures. In some applications (e.g. genomewide association studies, multiomics data analyses), the number of clusters also remains unknown. Because such biological information usually tends to follow a normal distribution, the distribution of samples of each cluster can be assumed to be Gaussian.
Hierarchical clustering methods, which can be mainly categorized into agglomerative (bottomup) and divisive (topdown) procedures, are well known [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]. In agglomerative procedures, each sample is initially assumed to be a cluster. The two nearest clusters (based on a distance measure or criterion function) are then merged at a time. This merger continues until all the samples are clustered into one group. Consequently, a tree like structure known as dendrogram is yielded. If the number of clusters is provided, the process of amalgamation of clusters can be terminated when the desired number of clusters is obtained. The first step of an agglomerative procedure considers all the possible mergers of two samples, which requires n(n − 1)/2 combinations (where n depicts the number of samples). Divisive procedures, on the other hand, perform clustering in an inverse way as compared to their agglomerative counterparts. They begin by considering a group (having all the samples) and divide it into two groups at each stage until all the groups comprise of only a single sample [21, 22]. In the first step of a divisive procedure all the partitions of a sample set are considered, which amounts to 2^{n − 1} − 1 combinations. This number of combinations grows exponentially and practically makes divisive clustering a difficult procedure to implement. However, there are a few divisive approaches which do not necessarily consider all the divisions [21]. In hierarchical classifications, each subcluster can be formed from one larger cluster split into two, or the union of two smaller clusters. In either case, false decisions made in early stages cannot be corrected later on. For this reason, divisive procedures, which start with the entire dataset, are in general considered safer than agglomerative approaches [21, 23]. Therefore, the accuracy of a divisive procedure is envisaged to be higher than that of an agglomerative procedure [24]. However, the high computational demand (O(2^{n})~O(n ^{5})) of divisive procedures has severely restricted their usage [24, 25] (though for special cases the complexity can be further reduced [26]). Therefore, the divisive procedure has not been generally used for hierarchical clustering, remaining largely ignored in the literature.
Hierarchical approaches do not require initial parameter settings and generally employed either linear or nonlinear regression models [27,28,29]. Over the last few decades, a number of hierarchical approaches have been proposed. Some of these popular schemes are summarized below. The single linkage or link agglomerative hierarchical approach (SLink) [30] merges two adjacent neighbour groups. Euclidean distance for computing the proximity between two clusters. SLink is very sensitive to data location and occasionally generates groups in a long chain (called as chaining effect). This chaining effect can be reduced developing a method based on farthest distance. This was achieved by the complete linkage (CLink) hierarchical approach [2]. Nevertheless, CLink is also sensitive to outliers. Sensitiveness could be further decreased by the average linkage (ALink) hierarchical approach [31, 32]. ALink implements linking by using the average distance between two groups. In a similar way, the median linkage (MLink) hierarchical approach [33] regards median distance for linking. In Ward’s linkage (WaLink), clusters are merged based on the optimal value of an objective function [34]. In weighted average distance linkage (WtLink) hierarchical clustering [35, 36], the group sizes are not considered when computing average distances. Consequently, smaller groups will be assigned larger weights during the clustering process [35]. Similarly, modelbased hierarchical clustering [4, 8] uses an objective function. Whereas the method in [8] follows a Bayesian analysis and uses both Dirichlet priors and multinomial likelihood function, the approach in [4] optimizes the distance between two GMMs. The number of group is previously defined. Most of these approaches are constructed using the agglomerative procedure, though their construction (with higher computational demand) is equally possible using the divisive procedure. Although divisive clustering is generally disregarded, some approaches like DIANA (DIvisive ANAlysis) program has been recently established [21]. In spite of wellestablished methods (i.e. EM algorithm [37, 38]) for estimating the parameters of a Gaussian mixture model, it is worth noting that hierarchical and expectationmaximization (EM) algorithms are very different in nature. The EM algorithm is an iterative optimization method, which requires prior knowledge of the number of clusters. It begins with a random choice of cluster centers and therefore returns different sets of clusters for distinct runs of the algorithm. Hierarchical clustering algorithms, on the other hand, do not require such prior knowledge and return a unique set of clusters. These advantages often make hierarchical clustering methods preferable to the EM algorithm for dealing with biological datasets where unique solutions are of utmost importance [39,40,41].
In this work, we described a new Divisive hieRArchical maximum likelihOod clusteriNg approach, abbreviated as DRAGON hereafter. This is a topdown procedure which does not find pairs. Instead, it takes out one sample at a time, maximally increasing the likelihood function. This process continues until the first cluster is obtained. This cluster is not further subdivided but removed from the sample set. In the remaining sample set, the same procedure is repeated for obtaining all the possible clusters. The removal of one sample out of n samples requires n search. This reduces the total search complexity to O(n ^{2} c) (where c is the number of clusters), which represents a significant reduction of the topdown procedure. The following sections present the mathematical derivation of the proposed model, and the analyses carried over on synthetic as well as on biological data to illustrate its usefulness.
Methods
Maximum likelihood clustering: an overview
This section summarizes an overview of the maximum likelihood method for clustering [22]. Here we are not introducing our method, instead we are proving a brief description of conventional maximum likelihood approach for clustering applications. It is possible to learn from an unlabeled data if some assumptions are taken. We will begin the section with an assumption that probability densities are known and it is required to estimate unknown parametric vector θ. The solution comes out to be similar to supervised learning case of maximum likelihood estimation. However, in the supervised learning case, the topology of groups of data is known. But in an unsupervised learning case one has to assume parametric form of data to reach to the solution. Here we describe how to estimate of maximum likelihood of clusters of a given sample set χ. The label of cluster of the sample sets is defined as ω. Assuming there are c clusters in the sample set (c ≥ 1), we define Ω = {ω _{ j }} (for j = 1, 2, …, c) as the cluster label for jth cluster χ _{ j } (In many clustering problems, the number of c is unknown, this issue we will deal in detail in later section and in Additional file 1). In this paper, we followed the notations from Duda et al. [22] for the convenience of readers. Let a sample set χ = {x _{1}, x _{2}, …, x _{ n }} be defined in a ddimensional space (It is assumed that d < n. For d ≫ n, dimensionality reduction techniques can be first applied for supervised or unsupervised learning tasks [42,43,44,45,46,47]). Let an unknown parameter vector be θ consisting of mean μ as well as covariance Σ. This will specify the mixture density as
where p(x _{ k }ω _{ j }, θ _{ j }) (for j = 1, …, c) is the conditional density, P(ω _{ j }) is the a priori probability and θ = {θ _{ j }}. The joint density is further defined using the log likelihood as
Assuming the joint density p(χθ) is differentiable w.r.t to θ then from Eqs. (1) to (2)
where \( {\nabla}_{{\boldsymbol{\uptheta}}_i}L \) is the gradient of L w.r.t. θ _{ i }. Assuming θ _{ i } and θ _{ j } are independent, and supposing a posteriori probability is
then from Eq. (4) we can observe that \( \frac{1}{p\left(\left.{\mathbf{x}}_k\right\boldsymbol{\uptheta} \right)}=\frac{P\left(\left.{\omega}_i\right{\mathbf{x}}_k,\boldsymbol{\uptheta} \right)}{p\left(\left.{\mathbf{x}}_k\right{\omega}_i,{\boldsymbol{\uptheta}}_i\right)P\left({\omega}_i\right)} \). Substituting this value in Eq. (3), we obtain
Note that in Eq. (5), (1/f(z))∇ _{ z } f(z) is arranged as ∇ _{ z } log f(z). Equation (5) can be equated to 0 (\( {\nabla}_{{\boldsymbol{\uptheta}}_i}L=0 \)) for obtaining maximum likelihood estimate \( {\widehat{\boldsymbol{\uptheta}}}_i \). This will give us the solution as (interested readers may refer to Duda et al. [22] for further details.)
In the case of normal distribution, the unknown mean and covariance {μ, Σ} parameters are replaced in θ in the Eqs. 6, 7 and 8 for yielding maximum likelihood estimates. The parameter θ is usually updated in an iterative fashion to attain \( \widehat{\boldsymbol{\uptheta}} \) by EM algorithms or hill climbing schemes.
DRAGON method: concept
Here we illustrate the clustering method DRAGON. In brief, the proposed procedure is topdown in nature. It initially considers the sample set as one cluster from which one sample is removed at a time. This increases the likelihood function and continues until the maximum likelihood is reached as depicted in Fig. 1 (where L _{1} is the cluster likelihood at the beginning of the process and L _{3} is the maximum likelihood after removing two samples). Once the first cluster is obtained it is removed from the sample set and the procedure is then repeated for attaining the subsequent clusters. Consequently, only one cluster will be retrieved from a sample set. It is assumed that samples are multinomial distributed, however, the number of clusters is not known at the beginning of the process.
To establish the maximum likelihood estimate in the divisive hierarchical context, we investigate the criterion function and the distance measure that satisfy it.
DRAGON method: algorithm
To find the distance measure, we first define the loglikelihood function of a cluster χ _{ s }, where χ _{ s } is a subset of χ. At the beginning, χ _{ s } is the same as χ, however, in every subsequent iteration a sample \( \widehat{\mathbf{x}} \) is removed from χ _{ s } such that the likelihood function
is maximized.
Since we are finding only one cluster in the sample set χ, a priori probability P(ω) can be ignored. We would like to explore how function L changes when a sample \( \widehat{\mathbf{x}} \) is taken out. Let us suppose centroid μ and covariance Σ of χ _{ s } are defined as
where the number of samples in χ _{ s } is depicted as n. Assuming that the component density is normal then Eq. (9) can be simplified as
where trace function is denoted by tr(). Since \( tr\left[{\varSigma}^{1}{\sum}_{\mathbf{x}\in {\chi}_s}\left(\mathbf{x}\boldsymbol{\upmu} \right){\left(\mathbf{x}\boldsymbol{\upmu} \right)}^{\mathrm{T}}\right]= tr\left({nI}_{d\times d}\right)= nd \), we can write L as
If a sample \( \widehat{\mathbf{x}} \) is removed from χ _{ s } then centroid and covariance (Eqs. (10) and (11)) will change as follows
In order to observe the alteration in the likelihood function (of Eq. (12)), we provide the following Lemma.
Lemma 1 Assume point \( \widehat{\boldsymbol{x}} \) is taken out of a set χ _{ s } and this changes the centroid and covariance (as Eqs. ( 13 ) and ( 14 ) described). Thereby the determinant of Σ ^{*} is defined as
Proof From Eq. (14), the determinant of Σ ^{*} will be
For any square matrix of size m × m, we can write AB = AB, and cA = c ^{m}A where c is any scalar, This would enable us to write Eq. (L1) in the following manner
\( \left{I}_{m\times m}+ AB\right \) can be proved to be I _{ n × n } + BA by Sylvester’s determinant theorem (where A ∈ ℝ^{m × n} and B ∈ ℝ^{n × m} are rectangular matrices). This would allow us to write
For any scalar c = c, the Lemma is then proved by substituting this term in Eq. (L1). □
It is now possible to define the change in L as
where ΔL is defined as
and P is expressed as
It can be observed from Eqs. (15) to (17) that when a sample \( \widehat{\mathbf{x}} \) is taken out of cluster χ _{ s }, the change in L mainly depends on the term P as all the other terms are not changing. If we want to select \( \widehat{\mathbf{x}} \) such that L ^{*} > L, this requires to solve the following maximization problem
Therefore, by removing \( \widehat{\mathbf{x}} \) the likelihood should increase until the maximum value is reached.
This procedure can track the location of the cluster having the highest density or likelihood. Because one sample is taken out at a time, it could be sensitive to data positions around the center of the cluster. Thereby, it is possible to locate the center of the cluster whereas its complete topology can be missed. In order to reduce such sensitiveness additional processing for tuning the cluster would be useful.
By taking out one sample at a time, we can obtain a cluster χ _{ s } that provides maximum likelihood. All the samples taken out can be collated in a set defined as \( {\chi}_s^{\ast } \), where \( {\chi}_s\cup {\chi}_s^{\ast }=\chi \). The centroid of cluster χ _{ s } can be obtained by μ _{ s } = E[χ _{ s }]. It can be then employed to compute the distance of all the samples from μ _{ s }; i.e. d _{ x } = δ(x, μ _{ s }) ∀ x ∈ χ, where δ denotes a distance measure (in this case the Euclidean metric) and d _{ x } is a 1dimensional sample or point corresponding to x. Thereafter, a centroidbased clustering scheme can be applied on this distance metric or inner product space (by considering 1dimensional data and partitioning it into 2 groups) to readjust the cluster χ _{ s }. Here clustering is applied on a distance metric, which can be either the Euclidean norm or any form of kernel (as it is derived from dot product). This procedure can be repeated if μ _{ s } is changing dramatically. The overall method is summarized in Table 1 and illustrated in Additional file 1 (slides 1–6).
The next issue is the estimation of the number of clusters (c). If the value of c were given, it is then easier to find the locations. However, in some applications, c is unknown. In such situations, a range of values of c can be inserted to the procedure so that the best value in the range can be estimated. If no clue about c were given, the maximum number of possible clusters C can be investigated and the best among them (c ≤ C) can be chosen. In order to estimate c, we first define the total likelihood function as
where L _{ i } is the likelihood of ith cluster. The total likelihood function L _{ t } can be computed for different values of c. If for a particular number of clusters (k), the variation between two successive total likelihood functions were not significant, then we can estimate c to be k. Let the difference between two successive total likelihoods be δL _{ t }(k) = L _{ t }(k + 1) − L _{ t }(k), this quantity can be normalized as
where δL _{ t } is normalized over all the possible values of k.
Figure S5 in Additional file 1 illustrates the above explanation with a dataset of 4 clusters.
DRAGON method: search complexity
In this section, we briefly discuss the search complexity of the DRAGON method. As explained above the proposed method begins by taking one sample out of the sample set, which increases the likelihood function and requires n search. However, in the second iteration the search reduces to n − 1. Finding a cluster having n _{1} samples requires (1/2)(n − n _{1})(n + n _{1} + 1) total search (see Additional file 2 for details). Therefore, the search for c clusters results O(n ^{2} c). It should be noted here that this search in conventional divisive hierarchical approaches is quite expensive, in the order of O(2^{n}). However, the search of DRAGON is in the order of O(n ^{2} c), which indicates a considerable reduction. Furthermore, DRAGON employs the kmeans clustering algorithm in the dot product space as an intermediate step. The computational complexity of kmeans is considered to be linear (e.g. using Lloyd’s algorithm this is O(2nt) because dimensionality is 1 in an intermediate step, and the number of classes is 2. Here, t represents the number of iterations).
Results and discussion
To validate the DRAGON method, we performed analyses using synthetic and biological data. We further compared its clustering accuracy with that of existing hierarchical methods.
Analysis on synthetic data
For the analysis on synthetic data, we generated Gaussian data of dimensionality d with 4 clusters. This data consisted of 400 samples with similar topology to that described in Figure S1 of Additional file 1). With the help of different random seeds we produced the data 20 times, and for on each occasion we calculated the clustering accuracy. We then computed the average or mean of clustering accuracy over 20 attempts to have a statistically stable value. The dimension of the generated data was increased from 2 to 30. For evaluation purposes, we also used other agglomerative hierarchical methods such as SLink, CLink, MLink, ALink, WaLink, and WtLink. The average clustering accuracies of these approaches over dimensionality d are shown in Fig. 2a. For all the above methods, we provided the number of clusters; i.e. c = 4. For the DRAGON method (of Table 1), we iterated two times (step 5 of Table 1) in all the experiments. Additionally, we assessed seven previously compared divisive hierarchical methods [24]: ALink (average link), SLink (single link), CLink (complete link), DunnsOrg (Dunn’s original), DunnsVar (Dunn’s variant), MacSmith (MacnaughtonSmith) and PDDP (Principal Direction). The average clustering accuracies of these divisive methods are summarized in Fig. 2b. As Fig. 2b shows when the dimensionality of data increases, MacSmith performs poorly whereas ALink, SLink and DunnsOrg slightly improve their performances. The data dimension did not affect the accuracy of DunnsVar whatsoever, which remains around 50%. The highest accuracy is achieved by the PDDP method (roughly 80%), however, this accuracy is still lower than that of DRAGON.
It can be observed from Fig. 2 that on Gaussian data, the DRAGON method provides promising results over other hierarchical methods either agglomerative or divisive. The WaLink hierarchical method also shows good results after the DRAGON method over dimensionality d = 2, …, 30.
To avoid limiting our validation to Gaussian datasets, we carried out additional analyses on synthetic data that included Pathbased [48], Flame [49] and Aggregation [50] datasets. The results are summarized in Tables S1.1, S1.2 and S1.3 of Additional file 1.
Analysis on biological data
We also utilized biological datasets, namely acute leukemia [51], mixedlineage leukemia (MLL) [52] and mutation data from The Cancer Genome Atlas for assessing clustering accuracy of several hierarchical approaches studied in this paper. A summary of datasets can be found below:
Acute leukemia dataset –comprises DNA microarray gene expressions. The samples belong to acute leukemias of humans. Two kinds are available: acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL). The dataset has 25 AML and 47 ALL bone marrow samples with a dimension of 7129.
MLL dataset –has three ALL, MLL and AML classes. It contains 72 leukemia samples where 24 belong to ALL, 20 belong to MLL and 28 belong to AML with dimensionality 12,582.
Mutation dataset – this dataset is derived from The Cancer Genome Atlas project (https://tcgadata.nci.nih.gov/docs/publications/tcga/?). It includes mutation data for breast cancer, glioblastoma, kidney cancer and ovarian cancer. The data is divided into two groups of 416 samples and 90 samples, which contain 1636 genes.
To vary its number of features or dimensions, we employed the Chisquared method for ranking the genes (the InfoGain feature selection method was also employed, see Additional file 3). We then performed clustering on samples and computed the clustering accuracy on dimensionality d = 2, …, 5. The clustering accuracies (was measured with the package AccMeasure 2011: http://www.mathworks.com/matlabcentral/fileexchange/32197clusteringresultsmeasurement/content/AccMeasure.m) on acute leukemia, MLL and mutation datasets are depicted in Tables 2, 3 and 4, respectively. The best outcomes are highlighted in bold faces.
It is noticed from Table 2 that WaLink, WtLink and MLink achieve the highest performance when d = 2. For the other dimensions (d = 3, 4 and 5), DRAGON shows the highest performance. From Table 3, it can be seen that DRAGON achieves reasonable performance for all the dimensions. On mutation data (Table 4), DRAGON is showing the highest performance for d = 2 and slightly low performance (82.2%) for the other dimensions. Dunn’s original and MacnaughtonSmith provide the highest performance when d = 3. Despite the reasonable performance of divisive clustering methods in Table 4, it is worth noting that their running times were extremely slow, specially when the number of samples increases as it is the case with mutation data. In general, it can be summarized that the DRAGON method exhibited promising results in terms of clustering accuracy over other hierarchical methods. Also its search complexity is O(n ^{2} c), which is significantly lower than that of conventional divisive approaches.
Conclusions
In this work, we proposed a divisive hierarchical maximum likelihood clustering method whose search complexity was reduced to O(n ^{2} c). Its overall clustering accuracy showed a significant improvement over that of agglomerative hierarchical clustering methods when compared on both synthetic and biological datasets.
Abbreviations
 ALL:

Acute lymphoblastic leukemia
 AML:

Acute myeloid leukemia
 CLink:

Complete linkage
 DIANA:

Divisive analysis
 DRAGON:

Divisive hierarchical maximum likelihood clustering
 EM:

Expectationmaximization
 GMM:

Gaussian mixture model
 ML:

Maximum likelihood
 MLink:

Median linkage
 MLL:

Mixed lineage leukemia
 SLink:

Single linkage
 WaLink:

Ward’s linkage
 WLink:

Weighted average distance linkage
 WtLink:

Weighted linkage
References
 1.
Castro RM, Coates MJ, Nowak RD. Likelihood based hierarchical clustering. IEEE Trans Signal Process. 2004;52(8):2308–21.
 2.
Defays D. An efficient algorithm for a complete link method. Comput J. 1977;20(4):364–6.
 3.
Fisher D. Iterative optimization and simplification of hierarchical clusterings. J Artif Intell Res. 1996;4(1):147–79.
 4.
Goldberger J, Roweis S. Hierarchical clustering of a mixture model. In: Neural information processing systems. Cambridge: MIT Press; 2005. p. 505–12.
 5.
Heller KA, Ghahramani Z. Bayesian hierarchical clustering. In: 22nd international conference on machine learning. ACM: New York, Bonn; 2005. p. 297–304.
 6.
Horng SJ, MY S, Chen YH, Kao TW, Chen RJ, Lai JL, Perkasa CD. A novel intrusion detection system based on hierarchical clustering and support vector machines. Expert Syst Appl. 2011;38(1):306–13.
 7.
Jain AK, Murty MN, Flynn PJ. Data clustering: a review. ACM Comput Surv. 1999;31(3):264–323.
 8.
Vaithyanathan S, Dom B. Modelbased hierarchical clustering. In: 16th conference in uncertainty in artificial intelligence; 2000. p. 599–608.
 9.
Sharma A, Boroevich KA, Shigemizu D, Kamatani Y, Kubo M, Tsunoda T. Hierarchical maximum likelihood clustering approach. IEEE Trans Biomed Eng. 2017;64(1):112–22.
 10.
Zhang W, Zhao D, Wang X. Agglomerative clustering via maximum incremental path integral. Pattern Recogn. 2013;46(11):3056–65.
 11.
Rokach L, Maimon O. Clustering methods. In: Maimon O, Rokach L, editors. Data mining and knowledge discovery handbook. US: Springer; 2005. p. 321–52.
 12.
Székely GJ, Rizzo ML. Hierarchical clustering via joint betweenwithin distances: extending Ward’s minimum variance method. J Classif. 2005;22:151–83.
 13.
Karypis G, Han EH, Kumar V. Chameleon: hierarchical clustering using dynamic modeling. Computer. 1999;32(8):68–75.
 14.
Murtagh F, Legendre P. Ward’s hierarchical agglomerative clustering method: which algorithms implement ward’s criterion? J Classif. 2014;31(3):274–95.
 15.
Bouguettaya A, Yu Q, Liu X, Zhou X, Song A. Efficient agglomerative hierarchical clustering. Expert Syst Appl Int J. 2015;42(5):2785–97.
 16.
Zhang X, Xu Z. Hesitant fuzzy agglomerative hierarchical clustering algorithms. Int J Syst Sci. 2015;46(3):562–76.
 17.
Murtagh F, Contreras P. Algorithms for hierarchical clustering: an overview. WIREs Data Min Knowl Discov. 2012;2(1):86–97.
 18.
Gómez D, Yáñez J, Guada C, Rodríguez JT, Montero J, Zarrazola E. Fuzzy image segmentation based upon hierarchical clustering. KnowlBased Syst. 2015;87:26–37.
 19.
Zhou GT, Hwang SJ, Schmidt M, Sigal L, Mori G. Hierarchical maximummargin clustering. Ithaca: Cornell University Library; 2015.
 20.
Wang J, Zhong J, Chen G, Li M, FX W, Pan Y. ClusterViz: a cytoscape APP for cluster analysis of biological network. IEEE/ACM Trans Comput Biol Bioinform. 2015;12(4):815–22.
 21.
Kaufman L, Rousseeuw PJ. Finding groups in data: an introduction to cluster analysis. 1st ed: WileyInterscience; Hoboken, New Jersey; 2005.
 22.
Duda RO, Hart PE, Stork DG. Pattern classification. 2nd ed: WileyInterscience; New York; 2000.
 23.
MacnaughtonSmith P, T Williams W B. Dale M G. Mockett L: Dissimilarity analysis: a new technique of hierarchical subdivision. Nature 1964, 202:10341035.
 24.
Roux M. A comparative study of divisive hierarchical clustering algorithms. Ithaca: Cornell University Library; 2015.
 25.
Roux M. Basic procedures in hierarchical cluster analysis. In: Devillers J, Karcher W, editors. Applied multivariate analysis in SAR and environmental studies. Dordrecht: Springer Netherlands; 1991. p. 115–35.
 26.
Guénoche A, Hansen P, Jaumard B. Efficient algorithms for divisive hierarchical clustering with the diameter criterion. J Classif. 1991;8(1):5–30.
 27.
Farrell S, Ludwig CJH. Bayesian and maximum likelihood estimation of hierarchical response time models. Psychon Bull Rev. 2008;15(6):1209–17.
 28.
Pinheiro J, Bates D. Mixedeffects models in S and SPLUS. 2000th ed. New York: Springer; 2009.
 29.
Raudenbush SW, Bryk AS. Hierarchical linear models: applications and data analysis methods. 2nd ed. Thousand Oaks: SAGE Publications; 2001.
 30.
Sibson R. SLINK: an optimally efficient algorithm for the singlelink cluster method. Comput J. 1973;16(1):30–4.
 31.
Sokal RR, Michener CD. A statistical method for evaluating systematic relationships. Univ Kansas Sci Bull. 1958;38(22):1409–38.
 32.
Legendre P, Legendre L. Numerical Ecology. 3rd ed: Elsevier; 2012.
 33.
Everitt BS, Landau S, Leese M, Stahl D. Cluster analysis. 5th ed: Wiley; 2011.
 34.
Ward JH. Hierarchical grouping to optimize an objective function. J Am Stat Assoc. 1963;58(301):236–44.
 35.
Podani J. Multivariate data analysis in ecology and systematics: SPB Academic Publishing; 1994.
 36.
McQuitty LL. Similarity analysis by reciprocal pairs for discrete and continuous data. Educ Psychol Meas. 1966;26(4):825–31.
 37.
McLachlan G, Krishnan T. The EM algorithm and extensions. 2nd ed: WileyInterscience; 2008.
 38.
Jordan MI, Jacobs RA. Hierarchical mixtures of experts and the EM algorithm. Neural Comput. 1994;6(2):181–214.
 39.
González M, Gutiérrez C, Martínez R. Expectationmaximization algorithm for determining natural selection of Ylinked genes through twosex branching processes. J Comput Biol. 2012;19(9):1015–26.
 40.
Do CB, Batzoglou S. What is the expectation maximization algorithm? Nat Biotechnol. 2008;26(8):897–9.
 41.
Holmes I. Using evolutionary expectation maximization to estimate indel rates. Bioinformatics. 2005;21(10):2294–300.
 42.
Sharma A, Paliwal KK. Fast principal component analysis using fixedpoint algorithm. Pattern Recogn Lett. 2007;28(10):1151–5.
 43.
Sharma A, Paliwal KK. A gradient linear discriminant analysis for small sample sized problem. Neural Process Lett. 2008;27(1):17–24.
 44.
Sharma A, Paliwal KK. Cancer classification by gradient LDA technique using microarray gene expression data. Data Knowl Eng. 2008;66(2):338–47.
 45.
Sharma A, Paliwal KK. Rotational linear discriminant analysis technique for dimensionality reduction. IEEE Trans Knowl Data Eng. 2008;20(10):1336–47.
 46.
Sharma A, Paliwal KK, Onwubolu GC. Classdependent PCA, MDC and LDA: a combined classifier for pattern classification. Pattern Recogn. 2006;39(7):1215–29.
 47.
Sharma A, Paliwal KK. A new perspective to null linear discriminant analysis method and its fast implementation using random matrix multiplication with scatter matrices. Pattern Recogn. 2012;45(6):2205–13.
 48.
Chang H, Yeung DY. Robust pathbased spectral clustering. Pattern Recogn. 2008;41(1):191–203.
 49.
Fu LM, Medico E. FLAME, a novel fuzzy clustering method for the analysis of DNA microarray data. BMC Bioinformatics. 2007;83. pp. 115, https://doi.org/10.1186/1471210583.
 50.
Gionis A, Mannila H, Tsaparas P. Clustering aggregation. ACM Trans Knowl Discov Data. 2007;1(1):1–30.
 51.
Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, et al. Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999;286(5439):531–7.
 52.
Armstrong SA, Staunton JE, Silverman LB, Pieters R, den Boer ML, Minden MD, Sallan SE, Lander ES, Golub TR, Korsmeyer SJ. MLL translocations specify a distinct gene expression profile that distinguishes a unique leukemia. Nat Genet. 2002;30(1):41–7.
Acknowledgements
We thank Maurice Roux for kindly providing us with his implementation of divisive hierarchical clustering algorithms. Highend computing resources were provided by the Advanced Center for Computing and Communication, RIKEN.
Funding
This study was funded by the CREST, JST Grant, Japan.
Availability of data and materials
The ALL and MLL datasets are publicly accessible and can be downloaded via the author’s webpage or the Kent Ridge Biomedical Data Set Repository. The mutation data can be retrieved from The Cancer Genome Atlas Project. The MATLAB package of DRAGON can be accessed at http://www.riken.jp/en/research/labs/ims/med_sci_math/ or http://www.alokailab.com
Declarations
The publication charges of this article were funded by JST CREST grant JPMJCR1412.
About this supplement
This article has been published as part of BMC Bioinformatics Volume 18 Supplement 16, 2017: 16th International Conference on Bioinformatics (InCoB 2017): Bioinformatics. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume18supplement16.
Author information
Affiliations
Contributions
AS developed the concept, conducted the experiments and wrote the first draft of the manuscript. YL carried out the experiments related to divisive hierarchical clustering and proofread the manuscript. TT supervised the entire project. All the authors read and approved the final manuscript.
Corresponding author
Correspondence to Tatsuhiko Tsunoda.
Ethics declarations
Ethics approval and consent to participate
None
Consent for publication
Not applicable
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional files
Additional file 1:
Divisive hierarchical maximum likelihood clustering. In this file an illustration of DRAGON method is given. Additionally, performance (in terms of Rand index) is given for synthetic data (Flame, Pathbased and Aggregation). (PDF 963 kb)
Additional file 2:
Computational consideration of DRAGON search. In this file derivation of computational complexity of DRAGON search is given. (PDF 84 kb)
Additional file 3:
Clustering accuracy using InfoGain feature selection method. In this file, InfoGain filtering method was used to perform feature selection. Thereafter, various clustering methods were used to evaluate the performance of DRAGON method. (PDF 68 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Sharma, A., López, Y. & Tsunoda, T. Divisive hierarchical maximum likelihood clustering. BMC Bioinformatics 18, 546 (2017). https://doi.org/10.1186/s1285901719655
Published:
Keywords
 Divisive approach
 Hierarchical clustering
 Maximum likelihood