Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

Open Access Highly Accessed Methodology article

Consensus embedding: theory, algorithms and application to segmentation and classification of biomedical data

Satish Viswanath and Anant Madabhushi

Author Affiliations

Dept. of Biomedical Engineering, Rutgers University, 599 Taylor Road, Piscataway, New Jersey 08854, USA

BMC Bioinformatics 2012, 13:26  doi:10.1186/1471-2105-13-26

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/13/26


Received:23 August 2011
Accepted:8 February 2012
Published:8 February 2012

© 2012 Viswanath and Madabhushi; licensee BioMed Central Ltd.

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

Background

Dimensionality reduction (DR) enables the construction of a lower dimensional space (embedding) from a higher dimensional feature space while preserving object-class discriminability. However several popular DR approaches suffer from sensitivity to choice of parameters and/or presence of noise in the data. In this paper, we present a novel DR technique known as consensus embedding that aims to overcome these problems by generating and combining multiple low-dimensional embeddings, hence exploiting the variance among them in a manner similar to ensemble classifier schemes such as Bagging. We demonstrate theoretical properties of consensus embedding which show that it will result in a single stable embedding solution that preserves information more accurately as compared to any individual embedding (generated via DR schemes such as Principal Component Analysis, Graph Embedding, or Locally Linear Embedding). Intelligent sub-sampling (via mean-shift) and code parallelization are utilized to provide for an efficient implementation of the scheme.

Results

Applications of consensus embedding are shown in the context of classification and clustering as applied to: (1) image partitioning of white matter and gray matter on 10 different synthetic brain MRI images corrupted with 18 different combinations of noise and bias field inhomogeneity, (2) classification of 4 high-dimensional gene-expression datasets, (3) cancer detection (at a pixel-level) on 16 image slices obtained from 2 different high-resolution prostate MRI datasets. In over 200 different experiments concerning classification and segmentation of biomedical data, consensus embedding was found to consistently outperform both linear and non-linear DR methods within all applications considered.

Conclusions

We have presented a novel framework termed consensus embedding which leverages ensemble classification theory within dimensionality reduction, allowing for application to a wide range of high-dimensional biomedical data classification and segmentation problems. Our generalizable framework allows for improved representation and classification in the context of both imaging and non-imaging data. The algorithm offers a promising solution to problems that currently plague DR methods, and may allow for extension to other areas of biomedical data analysis.

Background

The analysis and classification of high-dimensional biomedical data has been significantly facilitated via the use of dimensionality reduction techniques, which allow classifier schemes to overcome issues such as the curse of dimensionality. This is an issue where the number of variables (features) is disproportionately large compared to the number of training instances (objects) [1]. Dimensionality reduction (DR) involves the projection of data originally represented in a N-dimensional (N-D) space into a lower n-dimensional (n-D) space (known as an embedding) such that n << N. DR techniques are broadly categorized as linear or non-linear, based on the type of projection method used.

Linear DR techniques make use of simple linear projections and consequently linear cost functions. An example of a linear DR scheme is Principal Component Analysis [2] (PCA) which projects data objects onto the axes of maximum variance. However, maximizing the variance within the data best preserves class discrimination only when distinct separable clusters are present within the data, as shown in [3]. In contrast, non-linear DR involves a non-linear mapping of the data into a reduced dimensional space. Typically these methods attempt to project data so that relative local adjacencies between high dimensional data objects, rather than some global measure such as variance, are best preserved during data reduction from N- to n-D space [4]. This tends to better retain class-discriminatory information and may also account for any non-linear structures that exist in the data (such as manifolds), as illustrated in [5]. Examples of these techniques include locally linear embedding [5] (LLE), graph embedding [6] (GE), and isometric mapping [7] (ISOMAP). Recent work has shown that in several scenarios, classification accuracy may be improved via the use of non-linear DR schemes (rather than linear DR) for gene-expression data [4,8] as well as medical imagery [9,10].

However, typical DR techniques such as PCA, GE, or LLE may not guarantee an optimum result due to one or both of the following reasons:

• Noise in the original N-D space tends to adversely affect class discrimination, even if robust features are used (as shown in [11]). A single DR projection may also fail to account for such artifacts (demonstrated in [12,13]).

• Sensitivity to choice of parameters being specified during projection; e.g. in [14] it was shown that varying the neighborhood parameter in ISOMAP can lead to significantly different embeddings.

In this paper, we present a novel DR scheme known as consensus embedding which aims to overcome the problems of sensitivity to noise and choice of parameters that plague several popular DR schemes [12-14]. The spirit behind consensus embedding is to construct a single stable embedding by generating and combining multiple uncorrelated, independent embeddings; the hypothesis being that this single stable embedding will better preserve specific types of information in the data (such as class-based separation) as compared to any of the individual embeddings. Consensus embedding may be used in conjunction with either linear or non-linear DR methods and, as we will show, is intended to be easily generalizable to a large number of applications and problem domains. In this work, we will demonstrate the superiority of the consensus embedding representation for a variety of classification and clustering applications.

Figure 1 illustrates an application of consensus embedding in separating foreground (green) and background (red) regions via pixel-level classification. Figure 1(a) shows a simple RGB image to which Gaussian noise was added to the G and B color channels (see Figure 1(b)). We now consider each of the 3 color channels as features (i.e. N = 3) for all of the image objects (pixels). Classification via replicated k-means clustering [15] of all the objects (without considering class information) was first performed using the noisy RGB feature information (Figure 1(b)), in order to distinguish the foreground from background.

thumbnailFigure 1. Region partitioning of toy image data. (a) Original RGB image to which Gaussian noise was added to create (b) noisy RGB image. Image visualization of classes obtained by replicated k-means clustering [15] of all the pixels via (c) original noisy RGB space, and (d) graph embedding [6] of noisy RGB data. 2D plots of (e) R-G, (f) R-B, and (g)G-B planes are also shown where colors of objects plotted correspond to the region in (b) that they are derived from. The discriminatory 2D spaces ((e) and (f)) are combined via consensus embedding, and the visualized classification result is shown in (h). Note the significantly better image partitioning into foreground and background of (h) compared to (c) and (d).

The labels so obtained for each object (pixel) are then visualized in the image shown in Figure 1(c), where the color of the pixel corresponds to its cluster label. The 2 colors in Figure 1(c) hence correspond to the 2 classes (clusters) obtained. No discernible regions are observable in this figure. Application of DR (via GE) reduces the data to a n = 2-D space, where the graph embedding algorithm [6] non-linearly projects the data such that the object classes are maximally discriminable in the reduced dimensional space. However, as seen in Figure 1(d), clustering this reduced embedding space does not yield any obviously discernible image partitions either.

By plotting all the objects onto 2D plots using only the R-G (Figure 1(e)) and R-B (Figure 1(f)) color channels respectively, we can see that separation between the two classes exists only along the R axis. In contrast, the 2D G-B plot (Figure 1(g)) shows no apparent separation between the classes. Combining 1D embeddings obtained via applying graph embedding to Figures 1(e) and 1(f), followed by unsupervised clustering, yields the consensus embedding result shown in Figure 1(h). Consensus embedding clearly results in superior background/foreground partitioning compared to the results shown in Figures 1(c),(d).

Related Work and Significance

Classifier and clustering ensembles

Researchers have attempted to address problems of classifier sensitivity to noise and choice of parameters via the development of classifier ensemble schemes, such as Boosting [16] and Bagging [17]. These classifier ensembles guarantee a lower error rate as compared to any of the individual members (known as "weak" classifiers), assuming that the individual weak classifiers are all uncorrelated [18]. Similarly a consensus-based algorithm has been presented [15] to find a stable unsupervised clustering of data using unstable methods such as k-means [19]. Multiple "uncorrelated" clusterings of the data were generated and used to construct a co-association matrix based on cluster membership of all the points in each clustering. Naturally occurring partitions in the data were then identified. This idea was further extended in [20] where a combination of clusterings based on simple linear transformations of high-dimensional data was considered. Note that ensemble techniques thus (1) make use of uncorrelated, or relatively independent, analyses (such as classifications or projections) of the data, and (2) combine multiple analyses (such as classifications or projections) to enable a more stable result.

Improved DR schemes to overcome parameter sensitivity

As shown by [7], linear DR methods such as classical multi-dimensional scaling [21] are unable to account for non-linear proximities and structures when calculating an embedding that best preserves pairwise distances between data objects. This led to the development of non-linear DR methods such as LLE [5] and ISOMAP [7] which make use of local neighborhoods to better calculate such proximities. As previously mentioned, DR methods are known to suffer from certain shortcomings (sensitivity to noise and/or change in parameters). A number of techniques have recently been proposed to overcome these shortcomings.

In [22,23] methods were proposed to choose the optimal neighborhood parameter for ISOMAP and LLE respectively. This was done by first constructing multiple embeddings based on an intelligently selected subset of parameter values, and then choosing the embedding with the minimum residual variance. Attempts have been made to overcome problems due to noisy data by selecting data objects known to be most representative of their local neighborhood (landmarks) in ISOMAP [24], or estimating neighborhoods in LLE via selection of data objects that are unlikely to be outliers (noise) [13]. Similarly, graph embedding has also been explored with respect to issues such as the scale of analysis and determining accurate groups in the data [25]. However, all of these methods require an exhaustive search of the parameter space in order to best solve the specific problem being addressed. Alternatively, one may utilize class information within the supervised variants [26,27] of ISOMAP and LLE which attempt to construct weighted neighborhood graphs that explicitly preserve class information while embedding the data.

Learning in the context of dimensionality reduction

The application of classification theory to DR has begun to be explored recently. Athitsos et al presented a nearest neighbor retrieval method known as BoostMap [28], in which distances from different reference objects are combined via boosting. The problem of selecting and weighting the most relevant distances to reference objects was posed in terms of classification in order to utilize the Adaboost algorithm [16], and BoostMap was shown to improve the accuracy and speed of overall nearest neighbor discovery compared to traditional methods. DR has also previously been formulated in terms of maximizing the entropy [29] or via a simultaneous dimensionality reduction and regression methodology involving Bayesian mixture modeling [30]. The goal in such methods is to probabilistically estimate the relationships between points based on objective functions that are dependent on the data labels [29]. These methods have been demonstrated in the context of application of PCA to non-linear datasets [30]. More recently, multi-view learning algorithms [31] have attempted to address the problem of improving the learning ability of a system by considering several disjoint subsets of features (views) of the data. The work most closely related to our own is that of [32] in the context of web data mining via multi-view learning. Given that a hidden pattern exists in a dataset, different views of this data are each embedded and transformed such that known domain information (encoded via pairwise link constraints) is preserved within a common frame of reference. The authors then solve for a consensus pattern which is considered the best approximation of the underlying hidden pattern being solved for. A similar idea was examined in [33,34] where 1D projections of image data were co-registered in order to better perform operations such as image-based breathing gating as well as multi-modal registration. Unlike consensus embedding, these algorithms involve explicit transformations of embedding data to a target frame of reference, as well as being semi-supervised in encoding specific link constraints in the data.

Intuition and significance of consensus embedding

In this paper we present a novel DR scheme (consensus embedding) that involves first generating and then combining multiple uncorrelated, independent (or base) n-D embeddings. These base embeddings may be obtained via either linear or non-linear DR techniques being applied to a large N-D feature space. Note that we use the terms "uncorrelated, independent" with reference to the method of constructing base embeddings; similar to their usage in ensemble classification literature [18]. Indeed, techniques to generate multiple base embeddings may be seen to be analogous to those for constructing classifier ensembles. In the latter, base classifiers with significant variance can be generated by varying the parameter associated with the classification method (k in kNN classifiers [35]) or by varying the training data (combining decision trees via Bagging [17]). Previously, a consensus method for LLE was examined in [36] with the underlying hypothesis that varying the neighborhood parameter (κ) will effectively generate multiple uncorrelated, independent embeddings for the purposes of constructing a consensus embedding. The combination of such base embeddings for magnetic resonance spectroscopy data was found to result in a low-dimensional data representation which enabled improved discrimination of cancerous and benign spectra compared to using any single application of LLE. In this work we shall consider an approach inspired by random forests [37] (which in turn is a modification of the Bagging algorithm [17]), where variations within the feature data are used to generate multiple embeddings which are then combined via our consensus embedding scheme. Additionally, unlike most current DR approaches which require tuning of associated parameters for optimal performance in different datasets, consensus embedding offers a methodology that is not significantly sensitive to parameter choice or dataset type.

The major contributions of our work are:

■ A novel DR approach which generates and combines embeddings.

■ A largely parameter invariant scheme for dimensionality reduction.

■ A DR scheme easily applicable to a wide variety of pattern recognition problems including image partitioning, data mining, and high dimensional data classification.

The organization of the rest of this paper is as follows. In Section 2 we will examine the theoretical grounding and properties of consensus embedding, followed by algorithms to efficiently implement the consensus embedding scheme. In Section 3 we show the application of consensus embedding in the context of (1) partitioning of synthetic as well as clinical images, and (2) classification of gene-expression studies. Quantitative and qualitative results of this evaluation, as well as discussion of the results and concluding remarks, are presented in Section 4.

Methods

Theory of Consensus Embedding

The spirit of consensus embedding lies in the generation and combination of multiple embeddings in order to construct a more stable, stronger result. Thus we will first define various terms associated with embedding construction. Based on these, we can mathematically formalize the concept of generating and combining multiple base embeddings, which will in turn allow us to derive necessary and sufficient conditions that must be satisfied when constructing a consensus embedding. Based on these conditions we will describe the specific algorithmic steps in more detail. Notation that is used in this section is summarized in Table 1.

Table 1. Notation and symbols

Preliminaries

An object shall be referred to by its label c and is defined as a point in an N-dimensional space ℝN. It is represented by an N-tuple F(c) comprising its unique N-dimensional co-ordinates. In a sub-space ℝn ⊂ ℝN such that n << N, this object c in a set C is represented by an n-tuple of its unique n-dimensional coordinates X(c). ℝn is also known as the embedding of objects c C and is always calculated via some projection of ℝN. For example in the case of ℝ3, we can define F(c) = {f1, f2, f3} based on the co-ordinate locations (f1, f2, f3) on each of the 3 axes for object c C. The corresponding embedding vector of c C in ℝ2 will be X(c) = {e1, e2} with co-ordinate axes locations (e1, e2). Note that in general, determining the target dimensionality (n) for any ℝN may be done by a number of algorithms such as the one used in this work [38].

The notation Λcd, henceforth referred to as the pairwise relationship, will represent the relationship between two objects c, d C with corresponding vectors F(c), F(d) ∈ ℝN. Similarly, the notation δcd will be used to represent the pairwise relationship between two objects c, d C with embedding vectors X(c), X(d) ∈ ℝn. We assume that this relationship satisfies the three properties of a metric (e.g. Euclidean distance). Finally, a triplet of objects c, d, e C is referred to as an unique triplet if c d, d e, and c e. Unique triplets will be denoted simply as (c, d, e).

Definitions

Definition 1 The function Δ defined on a unique triplet (c, d, e) is called a triangle relationship, Δ(c, d, e), if when Λcd < Λce and Λcd < Λde, then δcd < δce and δcd < δde.

For objects c, d, e C whose relative pairwise relationships in ℝN are preserved in ℝn, the triangle relationship Δ(c, d, e) = 1. For ease of notation, the triangle relationship Δ(c, d, e) will be referred to as Δ where appropriate. Note that for a set of R unique objects (R = |C|, |.| is cardinality of a set), <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M1">View MathML</a> unique triplets may be formed.

Definition 2 Given Z unique triplets (c, d, e) ∈ C and an embedding n of all objects c, d, e C, the associated embedding strength <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M2">View MathML</a>.

The embedding strength (ES) of an embedding ℝn, denoted ψES(ℝn), is hence the fraction of unique triplets (c, d, e) ∈ C for which Δ(c, d, e) = 1.

Definition 3 A true embedding, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M3">View MathML</a>, is an embedding for which <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M4">View MathML</a>.

A true embedding <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M3">View MathML</a> is one for which the triangle relationship is satisfied for all unique triplets (c, d, e) ∈ C, hence perfectly preserving all pairwise relationships from ℝN to <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M3">View MathML</a>. Additionally, for all objects c, d C in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M3">View MathML</a>, the pairwise relationship is denoted as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M5">View MathML</a>.

Note that according to Definition 3, the most optimal true embedding may be considered to be the original ℝN itself, i.e. <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M6">View MathML</a>. However, as ℝN may not be optimal for classification (due to the curse of dimensionality), we are attempting to approximate a true embedding as best possible in n-D space. Note that multiple true embeddings in n-D space may be calculated from a single ℝN; any one of these may be chosen to calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M5">View MathML</a>.

Practically speaking, any ℝn will be associated with some degree of error compared to the original ℝN. This is almost a given since some loss of information and concomitant error can be expected to occur in going from a high- to a low-dimensional space. We can calculate the probability of pairwise relationships being accurately preserved from ℝN to ℝn i.e. the probability that Δ(c, d, e) = 1 for any unique triplet (c, d, e) ∈ C in any ℝn as,

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M7">View MathML</a>

(1)

More details on this formulation may be found in the Appendix. Note that the probability in Equation 1 is binomial as the complementary probability to p(Δ|c, d, e, ℝn) (i.e. the probability that Δ(c, d, e) ≠ 1 for any unique triplet (c, d, e) ∈ C in any ℝn) is given by 1 - p(Δ|c, d, e, ℝn) (in the case of binomial probabilities, event outcomes can be broken down into two probabilities which are complementary, i.e. they sum to 1).

Definition 4 A strong embedding, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M8">View MathML</a>, is an embedding for which <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M9">View MathML</a>.

In other words, a strong embedding is defined as one which accurately preserves the triangle relationship for more than some significant fraction (θ) of the unique triplets of objects c, d, e C that exist. An embedding ℝn which is not a strong embedding is referred to as a weak embedding, denoted as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M10">View MathML</a>.

We can calculate multiple uncorrelated (i.e. independent) embeddings from a single ℝN which may be denoted as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M11">View MathML</a>, where M is total number of possible uncorrelated embeddings. Note that both strong and weak embeddings will be present among all of the M possible embeddings. All objects c, d C can then be characterized by corresponding embedding vectors <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M12">View MathML</a> with corresponding pairwise relationship <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M13">View MathML</a>. Given multiple <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M13">View MathML</a>, we can form a distribution <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M14">View MathML</a>, over all M embeddings. Our hypothesis is that the maximum likelihood estimate (MLE) of <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M14">View MathML</a>, denoted as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M15">View MathML</a>, will approximate the true pairwise relationship <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M5">View MathML</a> for objects c, d C.

Definition 5 An embedding n is called a consensus embedding, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a>, if for all objects <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M17">View MathML</a>.

We denote the consensus embedding vectors for all objects c C by <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M18">View MathML</a>. Additionally, from Equation 1, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M19">View MathML</a> represents the probability that Δ(c, d, e) = 1 for any (c, d, e) ∈ C in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a>.

Necessary and sufficient conditions for consensus embedding

While <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> is expected to approximate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M3">View MathML</a> as best possible, it cannot be guaranteed that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M20">View MathML</a> as this is dependent on how well <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M15">View MathML</a> approximates <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M5">View MathML</a>, for all objects c, d C. <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M15">View MathML</a> may be calculated inaccurately as a result of considering pairwise relationships derived from weak embeddings, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M10">View MathML</a>, present amongst the M embeddings that are generated. As Proposition 1 and Lemma 1 below demonstrate, in order to ensure that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M21">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> must be constructed from a combination of multiple strong embeddings <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M8">View MathML</a> alone, so as to avoid including weak embeddings.

Proposition 1 If K M independent, strong embeddings <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M22">View MathML</a>, with a constant <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M23">View MathML</a>that Δ(c, d, e) = 1 for all (c, d, e) ∈ C, are used to calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M21">View MathML</a>as K → ∞.

Proof. If K M independent, strong embeddings alone are utilized in the construction of <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a>, then the number of weak embeddings is (M - K). As Equation 1 represents a binomial probability, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M19">View MathML</a> can be approximated via the binomial formulation of Equation 1 as,

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M24">View MathML</a>

(2)

where <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M25">View MathML</a> (Equation 1) is considered to be constant. Based on Equation 2, as K → ∞, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M26">View MathML</a>, which in turn implies that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M21">View MathML</a>; therefore <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> approaches <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M3">View MathML</a>.   □

Proposition 1 demonstrates that for a consensus embedding to be strong, it is sufficient that strong embeddings be used to construct it. Note that as K M, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M27">View MathML</a>. In other words, if <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M28">View MathML</a>. Based on Equation 1 and Definitions 2, 4, this implies that as K M, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M29">View MathML</a>. Lemma 1 below demonstrates the necessary nature of this condition i.e. if weak embeddings are considered when constructing <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M30">View MathML</a> (it will be a weak embedding).

Lemma 1 If K M independent, weak embeddings <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M31">View MathML</a>, with <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M32">View MathML</a>, are used to calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a>, then <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M30">View MathML</a>.

Proof. From Equation 1 and Definitions 2, 4, if <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M32">View MathML</a>, then <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M33">View MathML</a>. Substituting <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M34">View MathML</a> in Equation 2, will result in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M35">View MathML</a>. Thus <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M30">View MathML</a>, and <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> will be weak.   □

Proposition 1 and Lemma 1 together demonstrate the necessary and sufficient nature of the conditions required to construct a consensus embedding: that if a total of M base embeddings are calculated from a single ℝN, some minimum number of strong embeddings (K M) must be considered to construct a <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> that is a strong embedding. Further, a <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> so constructed will have an embedding strength <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M36">View MathML</a> that will increase significantly as we include more strong embeddings in its computation. Appendix B demonstrates an additional property of <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> showing that it preserves information from ℝN with less inherent error than any ℝn used in its construction.

Algorithms and Implementation

Based on Proposition 1, 3 distinct steps are typically required for calculating a consensus embedding. First, we must generate a number of base embeddings (M), the steps for which are described in CreateEmbed. We then select for strong embeddings from amongst M base embeddings generated, described in SelEmbed. We will also discuss criteria for selecting strong embeddings. Finally, selected embeddings are combined to result in the final consensus embedding representation as explained in CalcConsEmbed. We also discuss some of the computational considerations of our implementation.

Creating n-dimensional data embeddings

One of the requirements for consensus embedding is the calculation of multiple uncorrelated, independent embeddings ℝn from a single ℝN. This is also true of ensemble classification systems such as Boosting [16] and Bagging [17] which require multiple uncorrelated, independent classifications of the data to be generated prior to combination. As discussed previously, the terms "uncorrelated, independent" are used by us with reference to the method of constructing embeddings, as borrowed from ensemble classification literature [18]. Similar to random forests [37], we make use of a feature space perturbation technique to generate uncorrelated (base) embeddings. This is implemented by first creating M bootstrapped feature subsets of V features each (every subset ηm, m ∈ {1, . . . , M} containing <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M37">View MathML</a> features, no DR involved). Note, that the number of samples in each V-dimensional subset is the same as in the original N-dimensional space. Each V-dimensional ηm is then embedded in n-D space via DR (i.e. projecting from ℝV to ℝn). M is chosen such that each of N dimensions appears in at least one ηm.

Algorithm CreateEmbed

Input: F(c) ∈ ℝN for all objects c C, n

Output: <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M38">View MathML</a>

Data Structures: Feature subsets ηm, total number of subsets M, number of features in each subset V, DR method Φ

begin

   0. for m = 1 to M do

   1.    Select V < N features from ℝN, forming subset ηm;

   2.    Calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M39">View MathML</a>, for all c C using ηm and method Φ;

   3. endfor

end

As discussed in the introduction, multiple methods exist to generate base embeddings, such as varying a parameter associated with a method (e.g. neighborhood parameter in LLE, as shown in [36]) as well as the method explored in this paper (feature space perturbation). These methods are analogous to methods in the literature for generating base classifiers in a classifier ensemble [18], such as varying k in kNN classifiers (changing associated parameter) [39], or varying the training set for decision trees (perturbing the feature space) [37].

Selection of strong embeddings

Having generated M base embeddings, we first calculate their embedding strengths <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M40">View MathML</a> for all <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M41">View MathML</a>. The calculation of ψES can be done via performance evaluation measures such as those described below, based on the application and prior domain knowledge. Embeddings for which <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M42">View MathML</a> are then selected as strong embeddings, where θ is a pre-specified threshold.

Algorithm SelEmbed

Input: <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M39">View MathML</a> for all objects c C, m ∈ {1, . . . , M}

Output: <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M43">View MathML</a>

Data Structures: A list Q, embedding strength function ψES, embedding strength threshold θ

begin

   0. for m = 1 to M do

   1.    Calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M40">View MathML</a>;

   2.    if <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M42">View MathML</a>

   3.       Put m in Q;

   4.    endif

   5. endfor

   6. For each element k of Q, store <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M44">View MathML</a> for all objects c C;

end

Note that while θ may be considered to be a parameter which needs to be specified to construct the consensus embedding, we have found in our experiments that the results are relatively robust to variations in θ. In general, θ may be defined based on the manner of evaluating the embedding strength, as discussed in the next section.

Evaluation of embedding strength

We present two performance measures in order to evaluate embedding strength: one measure being supervised and relying on label information; the other being unsupervised and driven by the separability of distinct clusters in the reduced dimensional embedding space. In Experiment 4 we compare the two performance measures against each other to determine their relative effectiveness in constructing a strong consensus embedding.

Supervised evaluation of embedding strength

We have demonstrated that embedding strength increases as a function of classification accuracy (Theorem 1, Appendix), implying that strong embeddings will have high classification accuracies. Intuitively, this can be explained as strong embeddings showing greater class separation compared to weak embeddings. Given a binary labeled set of samples C, we denote the sets of objects corresponding to the two classes as S+ and S-, such that C = S+ S- and S+ S- = ∅. When using a classification algorithm that does not consider class labels, we can evaluate classification accuracy as follows:

1. Apply classification algorithm to C (embedded in ℝn) to find T clusters (unordered, labeled set of objects), denoted via <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M45">View MathML</a>.

2. For each <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M46">View MathML</a>

(a) Calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M47">View MathML</a>.

(b) Calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M48">View MathML</a>.

(c) Calculate classification accuracy for <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M46">View MathML</a>, as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M49">View MathML</a>.

3. Calculate classification accuracy of ℝn as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M50">View MathML</a>.

As classification has been done without considering label information, we must evaluate which of the clusters so obtained shows the greatest overlap with S+ (the class of interest). We therefore consider the classification accuracy of the cluster showing the most overlap with S+ as an approximation of the embedding strength of ℝn, i.e. ψES(ℝn) ≈ ϕAcc(ℝn).

Unsupervised evaluation of embedding strength

We utilize a measure known as the R-squared index (RSI), based off cluster validity measures [40], which can be calculated as follows:

1. Apply classification algorithm to C (embedded in ℝn) to find T clusters (unordered, labeled set of objects), denoted via <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M51">View MathML</a>.

2. Calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M52">View MathML</a> (where X(cj) is the mean of data values in the jth dimension).

3. Calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M53">View MathML</a>.

4. Calculate R-squared index of ℝn as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M54">View MathML</a>.

RSI may be considered both a measure of the degree of difference between clusters found in a dataset as well as measurement of the degree of homogeneity between them. The value of ϕRS ranges between 0 and 1, where if ϕRS = 0, no difference exists among clusters. Conversely, a value close to ϕRS = 1 suggests well-defined, separable clusters in the embedding space. Note that when using RSI to evaluate embedding strength, it will be difficult to ensure that all selected embeddings are strong without utilizing a priori information. In such a case we can attempt to ensure that a significant majority of the embeddings selected are strong, which will also ensure that the consensus embedding <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> is strong (based off Proposition 1).

Constructing the consensus embedding

Given K selected embeddings <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M55">View MathML</a>, we quantify pairwise relationships between all the objects in each <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M56">View MathML</a> via Euclidean pairwise distances. Euclidean distances were chosen for our implementation as they are well understood, satisfy the metric assumption of the pairwise relationship, as well as being directly usable within the other methods used in this work. Ωk denotes the ML estimator used for calculating <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M15">View MathML</a> from K observations <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M57">View MathML</a> for all objects c, d C.

Algorithm CalcConsEmbed

Input: <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M58">View MathML</a> for all objects c C, k ∈ {1, . . . , K}

Output: <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M59">View MathML</a>

Data Structures: Confusion matrix W, ML estimator Ω, projection method γ

begin

   0. for k = 1 to K do

   1.    Calculate Wk(i, j) = ||Xk(c) - Xk(d)||2 for all objects c, d C with indices i, j;

   2. endfor

   3. Apply normalization to all Wk, k ∈ {1, . . . , K};

   4. Obtain <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M60">View MathML</a>;

   5. Apply projection method γ to <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M61">View MathML</a> to obtain final consensus embedding <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a>;

end

Corresponding entries across all Wk (after any necessary normalization) are used to estimate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M15">View MathML</a> (and stored in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M61">View MathML</a>). In our implementation, we have used the median as the ML estimator as (1) the median is less corruptible to outliers, and (2) the median and the expectation are interchangeable if one assumes a normal distribution [41]. In Section 3 we compare classification results using both the mean and median individually as the ML estimator. We apply a projection method γ, such as multi-dimensional scaling (MDS) [21], to the resulting <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M61">View MathML</a> to embed the objects in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> while preserving the pairwise distances between all objects c C. The underlying intuition for this final step is based on a similar approach adopted in [15] where MDS was applied to the co-association matrix (obtained by accumulating multiple weak clusterings of the data) in order to visualize the clustering results. As <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M61">View MathML</a> is analogous to the co-association matrix, the projection method γ will allow us to construct the consensus embedding space <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a>.

One can hypothesize that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M61">View MathML</a> is an approximation of distances calculated in the original feature space. Distances in the original feature space can be denoted as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M62">View MathML</a> with indices i, j. An alternative approach could therefore be to calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M63">View MathML</a> in the original feature space and apply γ to it instead. However, noise artifacts in the original feature space may prevent it from being truly optimal for analysis [11]. As we will demonstrate in Section 3, simple DR, as well as consensus DR, provide superior representations of the data (by accounting for noise artifacts) as compared to using the original feature space directly.

Computational efficiency of Consensus Embedding

The most computationally expensive operations in consensus embedding are (1) calculation of multiple uncorrelated embeddings (solved as an eigenvalue problem in O(n3) time for n objects), and (2) computation of pairwise distances between all the objects in each strong embedding space (computed in time O(n2) for n objects). A slight reduction in both time and memory complexity can be achieved based on the fact that distance matrices will be symmetric (hence only the upper triangular need be calculated). Additionally, multiple embeddings and distance matrices can be computed via code parallelization. However these operations still scale polynomially based on the number of objects n.

To further reduce the computational burden we embed the consensus embedding paradigm within an intelligent sub-sampling framework. We make use of a fast implementation [42] of the popular mean shift algorithm [43] (MS) to iteratively represent data objects via their most representative cluster centers. As a result, the space retains its original dimensionality, but now comprises only some fractional number (n/t) of the original objects. These n/t objects are used in the calculations of consensus embedding as well as for any additional analysis. A mapping (Map) is retained from all n original objects to the final n/t representative objects. We can therefore map back results and analyses from the lowest resolution (n/t objects) to the highest resolution (n objects) easily. The fewer number of objects (n/t << n) ensures that consensus embedding is computationally feasible. In our implementation, t was determined automatically based on the number of stable cluster centers detected by MS.

Algorithm ConsEmbedMS

Input: F(c) ∈ ℝN for all objects c C, n

Output: <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M64">View MathML</a>

Data Structures: Reduced set of objects <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M65">View MathML</a>

begin

   0. Apply MS [42] to ℝN resulting in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M66">View MathML</a> for sub-sampled set of objects <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M65">View MathML</a>;

   1. Save Map from sub-sampled set of objects <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M65">View MathML</a> to original set of objects c C;

   2. <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M67">View MathML</a>;

   3. <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M68">View MathML</a>;

   4. <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M69">View MathML</a>;

   5. Use MS and Map to calculate <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M64">View MathML</a> from <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M70">View MathML</a> for all objects c C;

end

For an MRI image comprising 5589 pixels (objects) for analysis, the individual algorithms CreateEmbed, SelEmbed and CalcConsEmbed took 121.33, 12.22, and 35.75 seconds respectively to complete (on average). By implementing our mean-shift optimization it took only 119 seconds (on average) for ConsEmbedMS to complete analysis of an MRI image comprising between 15,000 and 40,000 pixels (objects); a calculation that would have been computationally intractable otherwise. All experiments were conducted using MATLAB 7.10 (Mathworks, Inc.) on a 72 GB RAM, 2 quad core 2.33 GHz 64-bit Intel Core 2 processor machine.

Experimental Design for Evaluating Consensus Embedding

Dataset description

The different datasets used in this work included: (1) synthetic brain image data, (2) clinical prostate image data, and (3) gene-expression data (comprehensively summarized in Table 2). The overarching goal in each experiment described was to determine the degree of improvement in class-based separation via the consensus embedding representation as compared to alternative representations (quantified in terms of classification accuracy). Note that in the case of prostate images as well as gene-expression data we have tested the robustness of the consensus embedding framework via the use of independent training and testing sets.

Table 2. Datasets

In the case of image data (brain, prostate MRI), we have derived texture features [44] on a per-pixel basis from each image. These features are based on calculating statistics from a gray level intensity co-occurrence matrix constructed from the image, and were chosen due to previously demonstrated discriminability between cancerous and non-cancerous regions in the prostate [45] and different types of brain matter [46] for MRI data. Following feature extraction, each pixel c in the MR image is associated with a N dimensional feature vector F(c) = [fu(c)|u ∈ {1, . . . , N}] ∈ ℝN, where fu(c) is the response to a feature operator for pixel c. In the case of gene-expression data, every sample c is considered to be associated with a high-dimensional gene-expression vector, also denoted F(c) ∈ ℝN.

DR methods utilized to reduce ℝN to ℝn were graph embedding (GE) [6] and PCA [2]. These methods were chosen in order to demonstrate instantiations of consensus embedding using representative linear and non-linear DR schemes. Additionally, these methods have been leveraged both for segmentation as well as classification of similar biomedical image and bioinformatics datasets in previous work [47,48]. The dimensionality of the embedding space, n, is calculated as the intrinsic dimensionality of ℝN via the method of [38]. To remain consistent with notation defined previously, the result of DR on F(c) ∈ ℝN is denoted XΦ(c) ∈ ℝn, while the result of consensus DR will be denoted <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M71">View MathML</a>. The subscript Φ corresponds to the DR method used, Φ ∈ {GE, PCA}. For ease of description, the corresponding classification results are denoted Ψ(F), Ψ(XΦ), <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M72">View MathML</a>, respectively.

Experiment 1: Synthetic MNI brain data

Synthetic brain data [49] was acquired from BrainWeb1, consisting of simulated proton density (PD) MRI brain volumes at various noise and bias field inhomogeneity levels. Gaussian noise artifacts have been added to each pixel in the image, while inhomogeneity artifacts were added via pixel-wise multiplication of the image with an intensity non-uniformity field. Corresponding labels for each of the separate regions within the brain, including white matter (WM) and grey matter (GM), were also available. Images comprising WM and GM alone were obtained from 10 sample slices (ignoring other brain tissue classes). The objective was to successfully partition GM and WM regions on these images across all 18 combinations of noise and inhomogeneity, via pixel-level classification (an application similar to Figure 1). Classification is done for all pixels c C based on each of,

(i) the high-dimensional feature space F(c) ∈ ℝN, N = 14,

(ii) simple GE on F(c), denoted XGE(c) ∈ ℝn, n = 3,

(iii) multi-dimensional scaling (MDS) on distances calculated directly in ℝN, denoted as XMDS(c) ∈ ℝn, n = 3 (alternative to consensus embedding, explained in Section 2),

(iv) consensus embedding, denoted <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M73">View MathML</a>, n = 3.

The final slice classification results obtained for each of these spaces are denoted as Ψ(F), Ψ(XGE), Ψ(XMDS), <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a>, respectively.

Experiment 2: Comparison of ML estimators in consensus embedding

For the synthetic brain data [49], over all 18 combinations of noise and inhomogeneity and over all 10 images, we compare the use of mean and median as ML estimators in CalcConsEmbed. This is done by preserving outputs from SelEmbed and only changing the ML estimator in the CalcConsEmbed. We then compare classification accuracies for detection of white matter in each of the resulting consensus embedding representations, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M75">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M76">View MathML</a> (superscript denotes choice of ML estimator).

Experiment 3: Clinical prostate MRI data

Two different prostates were imaged ex vivo using a 4 Tesla MRI scanner following surgical resection. The excised glands were then sectioned into 2D histological slices which were digitized using a whole slide scanner. Regions of cancer were determined via Haemotoxylin and Eosin (H&E) staining of the histology sections. The cancer areas were then mapped onto corresponding MRI sections via a deformable registration scheme [50]. Additional details of data acquisition are described in [45].

For this experiment, a total of 16 4 Tesla ex vivo T2-weighted MRI and corresponding digitized histology images were considered. The purpose of this experiment was to accurately identify cancerous regions on prostate MRI data via pixel-level classification, based on exploiting textural differences between diseased and normal regions on T2-weighted MRI [45]. For each MRI image, M embeddings, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M77">View MathML</a>, were first computed (via CreateEmbed) along with their corresponding embedding strengths <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M78">View MathML</a> (based on clustering classification accuracy). Construction of the consensus embedding was performed via a supervised cross-validation framework, which utilized independent training and testing sets for selection of strong embeddings (SelEmbed). The algorithm proceeds as follows,

(a) Training (Str) and testing (Ste) sets of the data (MRI images) were created.

(b) For each element (image) of Str, strong embeddings were identified based on <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M79">View MathML</a>.

(c) Those embeddings voted as being strong across all the elements (images) in Str were then identified and selected.

(d) For the data (images) in Ste, corresponding embeddings were then combined (via CalcConsEmbed) to yield the final consensus embedding result.

A leave-one-out cross-validation strategy was employed in this experiment. A comparison is made between the pixel-level classifications for (1) simple GE denoted as Ψ(XGE), and (2) consensus GE denoted as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a>.

Experiment 4: Gene-expression data

Four publicly available binary class gene-expression datasets were obtained2 with corresponding class labels for each sample [4]; the purpose of the experiment being to differentiate the two classes in each dataset. This data comprises the gene-expression vectorial data profiles of normal and cancerous samples for each disease listed in Table 2, where the total number of samples range from 72 to 181 patients and the number of corresponding features range from 7130 to 24,481 genes or peptides. All 4 data sets comprise independent training (Str) and testing (Ste) subsets, and these were utilized within a supervised framework for constructing and evaluating the consensus embedding representation.

Prior to analysis, each dataset was first pruned to the 300 most class-informative features based on t-statistics as described in [51]. The supervised cross-validation methodology for constructing the consensus embedding using independent training and testing sets is as follows,

(a) First, CreateEmbed is run concurrently on data in Str and Ste, such that the same subsets of features are utilized when generating base embeddings for each of Str and Ste.

(b) SelEmbed is then executed on base embeddings generated from Str alone, thus selecting strong embeddings from amongst those generated. Strong embeddings were defined based on <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M79">View MathML</a>.

(c) Corresponding (selected) embeddings for data in Ste are then combined within CalcConsEmbed to obtain the final consensus embedding vectors denoted as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M80">View MathML</a>, Φ ∈ {GE, PCA}, n = 4.

For this dataset, both supervised (via clustering classification accuracy, superscript S) and unsupervised (via RSI, superscript US) measures of embedding strength were evaluated in terms of the classification accuracy of the corresponding consensus embedding representations.

In lieu of comparative DR strategies, a semi-supervised variant of GE [52] (termed SSAGE) was implemented, which utilizes label information when constructing the embedding. Within this scheme, higher weights are given to within-class points and lower weights to points from different classes. When running SSAGE, both Str and Ste were combined into a single cohort of data, and labels corresponding to Str alone were revealed to the SSAGE algorithm.

An additional comparison was conducted against a supervised random forest-based kNN classifier operating in the original feature space to determine whether DR provided any advantages in the context of high-dimensional biomedical data. This was implemented by training a kNN classifier on each of the feature subsets for Str (that were utilized in Create Embed), but without performing DR on the data. Each such kNN classifier was then used to classify corresponding data in Ste. The final classification result for each sample in Ste is based on ensemble averaging to calculate the probability of a sample belonging to the target class. Classifications compared in this experiment were Ψ(F), Ψ(XSSGE), <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M81">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M82">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M83">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M84">View MathML</a>, respectively.

Classification

For image data (brain, prostate MRI), classification was done via replicated k-means clustering [15], while for gene-expression data, classification was done via hierarchical clustering [53]. The choice of clustering algorithm was made based on the type of data being considered in each of the different experiments, as well as previous work in the field. Note that both these clustering techniques do not consider class label information while classifying the data, and have been demonstrated as being deterministic in nature (hence ensuring reproducible results). The motivation in using such techniques for classification was to ensure that no classifier bias or fitting optimization was introduced during evaluation. As our experimental intent was purely to examine improvements in class separation offered by the different data representations, all improvements in corresponding classification accuracies may be directly attributed to improved class discriminability in the corresponding space being evaluated (without being dependent on optimizing the technique used for classification).

Evaluating and visualizing results

To visualize classification results as region partitions on the images (brain, prostate MRI), all the pixels were plotted back onto the image and assigned colors based on their classification label membership. Similar to the partitioning results shown in Figure 1, pixels of the same color were considered to form specific regions. For example, in Figure 1(h), pixels colored green were considered to form the foreground region, while pixels colored red were considered to form the background.

Classification accuracy of clustering results for images as well as gene-expression data can be quantitatively evaluated as described previously (Section 2). Image region partitioning results as well as corresponding classification accuracies of the different methods (GE, PCA, consensus embedding) were used to determine what improvements are offered by consensus embedding.

Results and Discussion

Experiment 1: Synthetic MNI Brain data

Figure 2 shows qualitative pixel-level WM detection results on MNI brain data for comparisons to be made across 3 different noise and inhomogeneity combinations (out of 18 possible combinations). The original PD MRI image for selected combinations of noise and inhomogeneity with the ground truth for WM superposed as a red contour is shown in Figures 2(a), (f), (k). Note that this is a 2 class problem, and GM (red) and WM (green) region partitions are visualized together in all the result images, as explained previously. Other brain tissue classes were ignored. Comparing the different methods used, when only noise (1%) is added to the data, all three of Ψ(F) (Figure 2(b)), Ψ(XMDS) (Figure 2(c)), and Ψ(XGE) (Figure 2(d)) are only able to identify the outer boundary of the WM region. However, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> (Figure 2(e)) shows more accurate detail of the WM region in the image (compare with the ground truth WM region outlined in red in Figure 2(a)). When RF inhomogeneity (20%) is added to the data for intermediate levels of noise (3%), note the poor WM detection results for Ψ(F) (Figure 2(g)), Ψ(XMDS) (Figure 2(h)), and Ψ(XGE) (Figure 2(i)). <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> (Figure 2(j)), however, yields a more accurate WM detection result (compared to the ground truth WM region in Figure 2(f)). Increasing the levels of noise (7%) and inhomogeneity (40%) results in further degradation of WM detection performance for Ψ(F) (Figure 2(l)), Ψ(XMDS) (Figure 2(m)), and Ψ(XGE) (Figure 2(n)). Note from Figure 2(o) that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> appears to fare far better than Ψ(F), Ψ(XMDS), and Ψ(XGE).

thumbnailFigure 2. WM detection results on synthetic BrainWeb data. Figure 2: Pixel-level WM detection results visualized for one image from the MNI brain MRI dataset, each row corresponding to a different combination of noise and inhomogeneity: (a)-(e) 1% noise, 0% inhomogeneity, (f)-(j) 3% noise, 20% inhomogeneity, (k)-(o) 7% noise, 40% inhomogeneity. The first column shows the original PD MRI image with the ground truth for WM outlined in red, while the second, third, fourth, and fifth columns show the pixel-level WM classification results for Ψ(F), Ψ(XMDS), Ψ(XGE), and <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a>, respectively. The red and green colors in (b)-(e), (g)-(j), (l)-(o) denote the GM and WM regions identified in each result image.

For each of the 18 combinations of noise and inhomogeneity, we averaged the WM detection accuracies ϕAcc(F), ϕAcc(XMDS), ϕAcc(XGE), <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M85">View MathML</a>(calculated as described in Section 2) over all 10 images considered (a total of 180 experiments). These results are summarized in Table 3 (corresponding trend visualization in Figure 3) with accompanying standard deviations in accuracy. Note that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M85">View MathML</a> shows a consistently better performance than the remaining methods (ϕAcc(F), ϕAcc(XMDS), ϕAcc(XGE)) in 17 out of 18 combinations of noise and inhomogeneity. This trend is also visible in Figure 3.

Table 3. WM detection results for synthetic BrainWeb data

thumbnailFigure 3. Trends in WM detection accuracy across Experiments 1 and 2. Visualization of classification accuracy trends (Tables 3 and 4). <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> (consensus embedding) performs significantly better than comparative strategies (original feature space, GE, MDS); using median as ML estimator (purple) may be marginally more consistent than using mean as ML estimator (orange). Ψ(F) (blue) and Ψ(XMDS) (red) perform similarly (corresponding trends directly superposed on one another).

For each combination of noise and inhomogeneity, a paired Students' t-test was conducted between <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M85">View MathML</a> and each of ϕAcc(F), ϕAcc(XMDS), and ϕAcc(XGE), with the null hypothesis being that there was no improvement via <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> over all 10 brain images considered. <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> was found to perform significantly better (p < 0.05) than all of Ψ(F), Ψ(XMDS), and Ψ(XGE) in 16 out of 18 combinations of noise and inhomogeneity.

Comparing ϕAcc(F), ϕAcc(XMDS), and ϕAcc(XGE), it can be observed that Ψ(F) and Ψ(XMDS) perform similarly for all combinations of noise and inhomogeneity (note that the corresponding red and blue trend-lines completely overlap in Figure 3). In contrast, Ψ(XGE) shows improved performance at every combination of noise and inhomogeneity as compared to either of Ψ(F) and Ψ(XMDS). <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> was seen to significantly improve over all of Ψ(F), Ψ(XMDS), and Ψ(XGE), reflecting the advantages of consensus embedding.

Experiment 2: Comparison of ML estimators

WM pixel-level detection accuracy results for consensus embedding using two different ML estimators (median and mean) were averaged over all 10 MNI brain images considered and summarized in Table 4, for each of the 18 combinations of noise and inhomogeneity (total of 180 experiments). We see that the accuracy values are generally consistent across all the experiments conducted. No statistically significant difference in classifier performance was observed when using <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M86">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M87">View MathML</a>. It would appear that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M86">View MathML</a> is less susceptible to higher noise and bias field levels compared to <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M87">View MathML</a> (trends in Figure 3).

Table 4. Comparing the mean and median as ML estimators within CalcConsEmbed

Experiment 3: Clinical Prostate MRI data

Figure 4 shows qualitative results of the ConsEmbedMS algorithm in detecting prostate cancer (CaP) on T2-weighted MRI, each row corresponding to a different 2D MRI image. Comparing the pixel-level CaP detection results (visualized in green) in Figures 4(c) and 4(g) to the green CaP masks in Figures 4(b) and 4(f), obtained by registering the MRI images with corresponding histology images [50] (not shown), reveals that Ψ(XGE) results in a large false positive error. In contrast, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> (Figures 4(d) and 4(h)) appears to better identify the CaP region when compared to the ground truth for CaP extent in Figures 4(b) and 4(f). Figure 5 illustrates the relative pixel-level prostate cancer detection accuracies averaged across 16 MRI slices for the 2 methods compared. <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> was found to significantly (p < 0.05) outperform Ψ(XGE) in terms of accuracy and specificity of CaP segmentation over all 16 slices considered.

thumbnailFigure 4. Prostate cancer detection results on ex vivo MRI data. (a), (e) 2D sections from 3D prostate MRI data, and (b), (f) corresponding CaP masks superposed on the MRI, obtained via deformable registration with the corresponding histology slice (not shown) [50]. Corresponding CaP detection results via (c), (g) Ψ(XGE) (graph embedding), and (d), (h) <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M112','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M112">View MathML</a> (consensus embedding) are superposed back onto the original MRI sections ((a), (e)). In each of (b)-(d) and (f)-(h), green denotes the CaP segmentation region. Note the significantly fewer false positives in (d) and (h) compared to (c) and (g) respectively.

thumbnailFigure 5. Prostate cancer detection accuracy on ex vivo MRI data. Pixel-level classification accuracy in identifying prostate cancer on T2-weighted MRI, averaged over 16 2D MRI slices for Ψ(XGE) (blue) and <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M74">View MathML</a> (red).

Experiment 4: Classification of Gene-Expression Data

Table 5 summarizes classification accuracies for each of the strategies compared: supervised consensus-PCA and consensus-GE (<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M82">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M81">View MathML</a>, respectively), unsupervised consensus-PCA and consensus-GE (<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M84">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M83">View MathML</a>, respectively), SSAGE (Ψ(XSSGE)), as well as supervised classification of the original feature space (Ψ(F)). These results suggest that consensus embedding yields a superior classification accuracy compared to alternative strategies. We posit that this improved performance is due to the more accurate representation of the data obtained via consensus embedding.

Table 5. Classification accuracies for different representation strategies for gene-expression data

The presence of a large noisy, high-dimensional space was seen to adversely affect supervised classification performance of F, which yielded a worse classification accuracy than unsupervised classification (of consensus-GE and consensus-PCA) in 3 out of the 4 datasets. Moreover, semi-supervised DR, which utilized label information to construct XSSGE, was also seen to perform worse than consensus embedding (both supervised and unsupervised variants). We posit that this is because SSAGE does not explicitly account for noise, but only modifies the pairwise relationships between points based on label information (possibly exacerbating the effects of noise). By contrast, any label information used by consensus embedding is used to account for noisy samples when approximating the "true" pairwise relationships between points. The difference in the final embedding representations can be visualized in 3D in Figure 6, obtained by plotting all the samples in the lung cancer gene-expression dataset in 3D Eigen space. Note that consensus DR (Figures 6(b)-(e)) shows significantly better separation between the classes with more distinct, tighter clusters as well as fewer false positives compared to SSAGE (Figure 6(a)).

thumbnailFigure 6. Visualization of 3D embeddings for gene-expression data (breast cancer). 3D visualization of embedding results for lung cancer gene-expression data: (a) XSSGE, (b) <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M113','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M113">View MathML</a>, (c) <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M114','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M114">View MathML</a>, (d) <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M115','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M115">View MathML</a>, (e) <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M116','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M116">View MathML</a>. The 3 axes correspond to the primary 3 eigenvalues obtained via different DR methods (SSAGE, consensus-GE and consensus-PCA), while the colors of the objects (red and blue) are based on known class information (cancer and non-cancer, respectively). Note the relatively poor performance of (a) semi-supervised DR compared to (b)-(e) consensus DR. Both supervised ((b) and (d)) and unsupervised ((c) and (e)) consensus DR show relatively consistent separation between the classes with distinct, tight clusters. The best clustering accuracy for this dataset was achieved by (c) unsupervised consensus GE <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M117','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M117">View MathML</a>.

Further, comparing the performance of supervised <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M88">View MathML</a> and unsupervised <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M89">View MathML</a> variants of consensus embedding demonstrates comparable performance between them, though a supervised measure of embedding strength shows a trend towards being more consistent. The relatively high performance of <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M84">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M83">View MathML</a> demonstrate the feasibility of a completely unsupervised framework for consensus embedding.

For both consensus PCA and consensus GE we tested the parameter sensitivity of our scheme by varying the number of feature subsets generated (M ∈ {200, 500, 1000}) in the CreateEmbed algorithm (Tables 6 & 7). The relatively low variance in classification accuracy as a function of M reflects the invariance to parameters of consensus embedding. No consistent trend was seen in terms of either of consensus-GE or consensus-PCA outperforming the other.

Table 6. Variation in classification accuracy as a function of parameters for consensus-PCA on gene-expression data

Table 7. Variation in classification accuracy as a function of parameters for consensus-GE on gene-expression data

Conclusions

We have presented a novel dimensionality reduction scheme called consensus embedding which can be used in conjunction with a variety of DR methods for a wide range of high-dimensional biomedical data classification and segmentation problems. Consensus embedding exploits the variance within multiple base embeddings and combines them to produce a single stable solution that is superior to any of the individual embeddings, from a classification perspective. Specifically, consensus embedding is able to preserve pairwise object-class relationships from the high- to the low-dimensional space more accurately compared to any single embedding technique. Using an intelligent sub-sampling approach (via mean-shift) and code parallelization, computational feasibility and practicability of our method is ensured.

Results of quantitative and qualitative evaluation in over 200 experiments on toy, synthetic, and clinical images in terms of detection and classification accuracy demonstrated that consensus embedding shows significant improvements compared to traditional DR methods such as PCA. We also compared consensus embedding to using the feature space directly, as well as to using an embedding based on distance preservation directly from the feature space (via MDS [21]), and found significant performance improvements when using consensus embedding. Even though the features and classifier used in these experiments were not optimized for image segmentation purposes, consensus embedding outperforms state-of-the-art segmentation schemes (graph embedding, also known as normalized cuts [6]), differences being statistically significant in all cases. Incorporating spatial constraints via algorithms such as Markov Random Fields [54] could be used to further bolster the image segmentation results via consensus embedding.

In experiments for high-dimensional biomedical data analysis using gene-expression signatures, consensus embedding also demonstrated improved results compared to semi-supervised DR methods (SSAGE [52]). Evaluating these results further illustrates properties of consensus embedding: (1) the consensus of multiple projections improves upon any single projection (via either linear PCA or non-linear GE), (2) the error rate for consensus embedding is not significantly affected by parameters associated with the method, as compared to traditional DR. Finally, the lower performance of a supervised classifier using the original noisy feature space as compared to using the consensus embedding representation demonstrates the utility of DR to obtain improved representations of the data for classification.

It is however worth noting that in certain scenarios, consensus embedding may not yield optimal results. For instance, if very few embeddings are selected for consensus, the improvement in performance via consensus embedding over simple DR techniques may not be as significant. This translates to having a sparsely populated distribution for the estimation of the consensus pairwise relationship, resulting in a lower confidence being associated with it. Such a scenario may arise due to incorrectly specified selection criteria for embeddings; however, it is relatively simple to implement self-tuning for the selection parameter (θ). Note we have reported results for a fixed value of θ in all our experiments/applications, further demonstrating the robustness of our methodology to choice of parameters.

In this work, consensus embedding has primarily been presented within a supervised framework (using class label information to evaluate embedding strength). Preliminary results in developing an unsupervised evaluation measure using the R-squared cluster validity index [40] are extremely promising. However, additional tuning and testing of the measure is required to ensure robustness.

Another area of future work is developing algorithms for the generating uncorrelated, independent embeddings. This is of great importance as generating truly uncorrelated, independent embeddings will allow us to capture the information from the data better, hence ensuring in an improved consensus embedding result. As mentioned previously, methods to achieve this could include varying the parameter associated with the DR method (e.g. neighborhood parameter in LLE [5]) as well as the feature space perturbation method explored in this paper. These approaches are analogous to methods of generating weak classifiers within a classifier ensemble [18], such as varying the k parameter in kNN classifiers [39] or varying the training set for decision trees [37]. Note our feature space perturbation method to generate multiple, uncorrelated independent embeddings is closely related to the method used in random forests [37] to generate multiple weak, uncorrelated classifiers. Thus the embeddings we generate, as with the multiple classifiers generated in ensemble classifier schemes, are not intended to be independent in terms of information content, but rather in their method of construction.

The overarching goal of consensus embedding is to optimally preserve pairwise relationships when projecting from high- to low-dimensional space. In this work, pairwise relationships were quantified by us using the popular Euclidean distance metric. This was chosen as it is well understood in the context of these methods used within our algorithm (e.g. the use of MDS). Alternative pairwise relationship measures could include the geodesic distance [7] or the symmetrized Kullback-Leibler divergence [55]. It is important to note that such measures will need to satisfy the properties of a metric to ensure that they correctly quantify both triangle as well as pairwise relationships. We currently use MDS [21] to calculate the final consensus embedding (based on the consensus pairwise distance matrix). We have chosen to use MDS due to ease of computational complexity, but this method could be replaced by a non-linear variant instead. Finally, our intelligent sub-sampling approach to ensure computational feasibility comes with a caveat of the out-of-sample extension problem [56]. We currently handle this using a mapping of results from high to low-resolutions, but are currently identifying more sophisticated solutions. We intend to study these areas in greater detail to further validate the generalizability of consensus embedding.

Authors' contributions

AM and SV co-conceived the core algorithm and theoretical justifications of consensus embedding. SV further developed, evaluated, and refined the implementation and experiments. AM directed the research and the development of the manuscript. Both authors contributed to writing and editing, and have read and approved the final manuscript.

Appendix A: Correspondence between Equation 1 and Definition 2

In order to calculate the probability p(Δ|c, d, e, ℝn), we utilize the traditional formulation of a prior probability,

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M90">View MathML</a>

(3)

With reference to Equation 1, "instances" are triplets. Therefore Equation 3 becomes,

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M91','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M91">View MathML</a>

(4)

Independent of the above, we intuitively arrived at a mathematical formulation for embedding strength.

The strength of any embedding ℝn will depend on how well pairwise relationships are preserved from ℝN.

This in turn can written in terms of the triplet relationship as well,

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M92','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M92">View MathML</a>

(5)

Appendix B: Properties of consensus embedding

The following proposition will demonstrate that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> will have a lower inherent error in its pairwise relationships compared to the strong embeddings <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M93','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M93">View MathML</a>, used in its construction. Note that relevant notation and definitions have been carried over from Section 2.

We first define the mean squared error (MSE) in the pairwise relationship between every pair of objects c, d C in any embedding ℝn with respect to the true pairwise relationships in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M3">View MathML</a> as,

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M94">View MathML</a>

(6)

where Ecd is the expectation of the squared error in the pairwise relationships in ℝn calculated over all pairs of objects c, d C. We can hence calculate the expected MSE over all K base embeddings specified above as,

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M95">View MathML</a>

(7)

Given K observations <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M96">View MathML</a> (derived from selected base embeddings <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M97','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M97">View MathML</a>), we define the pairwise relationship in the consensus embedding <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> as <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M98','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M98">View MathML</a>, where EK is the expectation of <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M99','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M99">View MathML</a> over K observations. The MSE in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M15">View MathML</a> with respect to the true pairwise relationships in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M3">View MathML</a> may be defined as (similar to Equation 6),

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M100','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M100">View MathML</a>

(8)

where Ecd is the expectation of the squared error in the pairwise relationships in <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> calculated over over all pairs of objects c, d C. It is clear that if for all c, d C that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M101','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M101">View MathML</a>, then <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> is also a true embedding.

Proposition 2 Given K independent, strong embeddings, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M102">View MathML</a>, which are used to construct <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M103','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M103">View MathML</a>.

Proof.

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M104','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M104">View MathML</a>

Proposition 2 implies that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M16">View MathML</a> will never have a higher error than the maximum error associated with any individual strong embedding <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M102">View MathML</a>, involved in its construction. However if εK,X is low, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M105">View MathML</a> may not significantly improve on it. Similar to Bagging [17] where correlated errors across weak classifiers are preserved in the ensemble result, if the pairwise relationship <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M99','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M99">View MathML</a> is incorrect across all K embeddings, <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M15">View MathML</a> will be incorrect as well. However Proposition 2 guarantees that <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M105">View MathML</a> will never be worse than εK,X.

Appendix C: Practical implementation of embedding strength

While embedding strength may be seen as a generalized concept for evaluating embeddings, in this work we have examined applications of DR and consensus embedding to classifying biomedical data (Section 3). We now derive a direct relationship between embedding strength and classification accuracy, presented in Theorem 1 below.

For the purposes of the following discussion, all objects c, d, e C are considered to be associated with class labels l(c), l(d), l(e) ∈ {ω1, ω2}, respectively, such that if l(c) = l(d) = ω1 and l(e) = ω2 then Λcd < Λce and Λcd < Λde. Note that ω1, ω2 are binary class labels that can be assigned to all objects c C.

Definition 6 An unique triplet (c, d, e) ∈ C with l(c), l(d), l(e) ∈ {ω1, ω2} is called a class triplet, (c, d, e)l, if either l(c) ≠ l(d), or l(d) ≠ l(e), or l(c) ≠ l(e).

Thus, in a class triplet of objects, two out of three objects have the same class label but the third has a different class label, e.g. if l(c) = l(d) = ω1 and l(e) = ω2. Further, in the specific case that Δ(c, d, e) = 1 for a class triplet (c, d, e)l, it will be denoted as Δl(c, d, e). For the above example of a class triplet, we know that Λcd < Λce and Λcd < Λde (see above). If Δ(c, d, e) = 1, δcd < δce and δcd < δde. This implies that even after projection from ℝN to ℝn, the class-based pairwise relationships within the data are accurately preserved (a classifier can be constructed which will correctly classify objects c, d, e in ℝn).

Consider that if <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M106">View MathML</a> objects have class label ω1, then <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M107">View MathML</a> objects will have class label ω2. Based on the total number of unique triplets Z, the total number of triplets which are not class triplets is,

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M108','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M108">View MathML</a>

(9)

Y will be a constant for a given set of objects C, and is based on forming unique triplets (c, d, e) where l(c) = l(d) = l(e) (triplets which are not class triplets). U = (Z - Y) will correspond to the number of class triplets that may be formed for set C. If all U class triplets have Δl(c, d, e) = 1, then it is possible to construct U classifiers which correctly classify the corresponding objects in these class triplets.

Definition 7 Given U unique class triplets (c, d, e)l C and an embedding n of all objects c, d, e C, the associated classification accuracy <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M109','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M109">View MathML</a>

As illustrated previously, class triplets (c, d, e)l for which Δl(c, d, e) = 1 will correspond to those objects which will be classified correctly in ℝn. Therefore, the classification accuracy ϕAcc(ℝn) may simply be defined as the fraction of class triplets (c, d, e)l C for which Δl(c, d, e) = 1.

Theorem 1 For any n, the corresponding ψES(ℝn) increases monotonically as a function of ϕAcc(ℝn).

Proof.

<a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M110','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M110">View MathML</a>

As <a onClick="popup('http://www.biomedcentral.com/1471-2105/13/26/mathml/M111','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/13/26/mathml/M111">View MathML</a> is a constant, ψES(ℝn) increases monotonically with ϕAcc(ℝn).   □

Thus an embedding ℝn with a high embedding strength will have a high classification accuracy. Practically, this implies that ψES(ℝn) may be estimated via any measure of object-class discrimination such as classification accuracy or cluster-validity measures. We have exploited this relationship in our algorithmic implementation (Section 2).

Endnotes

1http://www.bic.mni.mcgill.ca/brainweb/ webcite

2These datasets were downloaded from the Biomedical Kent-Ridge Repositories at http://datam.i2r.a-star.edu.sg/datasets/krbd/ webcite

Acknowledgements

The authors would like to acknowledge Drs. Mark Rosen, John Tomasezewski, and Michael Feldman from the Hospital of the University of Pennsylvania for the use of ex vivo prostate MRI and histology data. They would also thank Andrew Janowczyk, Dr. Jonathan Chappelow, and Dr. James Monaco for results, discussions, and implementations used in this paper. This work was made possible via grants from the Department of Defense Prostate Cancer Research Program (W81XWH-08-1-0072), Wallace H. Coulter Foundation, National Cancer Institute (Grant Nos. R01CA136535-01, R01CA140772-01, and R03CA143991-01), The Cancer Institute of New Jersey, and the Society for Imaging Informatics in Medicine.

References

  1. Bellman R: Adaptive control processes: a guided tour. Princeton University Press; 1961. OpenURL

  2. Jolliffe I: Principal Component Analysis. Springer; 2002. OpenURL

  3. Lin T, Zha H: Riemannian Manifold Learning.

    IEEE Transactions on Pattern Analysis and Machine Intelligence 2008, 30(5):796-809. PubMed Abstract | Publisher Full Text OpenURL

  4. Lee G, Rodriguez C, Madabhushi A: Investigating the Efficacy of Nonlinear Dimensionality Reduction Schemes in Classifying Gene- and Protein-Expression Studies.

    IEEE/ACM Transactions on Computational Biology and Bioinformatics 2008, 5(3):1-17. PubMed Abstract | Publisher Full Text OpenURL

  5. Saul L, Roweis S: Think globally, fit locally: unsupervised learning of low dimensional manifolds.

    Journal of Machine Learning Research 2003, 4:119-155. OpenURL

  6. Shi J, Malik J: Normalized Cuts and Image Segmentation.

    IEEE Transactions on Pattern Analysis and Machine Intelligence 2000, 22(8):888-905. Publisher Full Text OpenURL

  7. Tenenbaum J, Silva V, Langford J: A Global Geometric Framework for Nonlinear Dimensionality Reduction.

    Science 2000, 290(5500):2319-2323. PubMed Abstract | Publisher Full Text OpenURL

  8. Dawson K, Rodriguez R, Malyj W: Sample phenotype clusters in high-density oligonucleotide microarray data sets are revealed using Isomap, a nonlinear algorithm.

    BMC Bioinformatics 2005, 6(1):195. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  9. Madabhushi A, Shi J, Rosen M, Tomaszeweski JE, Feldman MD: Graph embedding to improve supervised classification and novel class detection: application to prostate cancer.

    Proc 8th Int'l Conf Medical Image Computing and Computer-Assisted Intervention (MICCAI) 2005, 729-37. PubMed Abstract | Publisher Full Text OpenURL

  10. Varini C, Degenhard A, Nattkemper T: Visual exploratory analysis of DCE-MRI data in breast cancer by dimensional data reduction: A comparative study.

    Biomedical Signal Processing and Control 2006, 1(1):56-63. Publisher Full Text OpenURL

  11. Quinlan J: The effect of noise on concept learning. Morgan Kaufmann; 1986. OpenURL

  12. Balasubramanian M, Schwartz EL, Tenenbaum JB, de Silva V, Langford JC: The Isomap Algorithm and Topological Stability.

    Science 2002, 295(5552):7a. Publisher Full Text OpenURL

  13. Chang H, Yeung D: Robust locally linear embedding.

    Pattern Recognition 2006, 39(6):1053-1065. Publisher Full Text OpenURL

  14. Shao C, Huang H, Zhao L: P-ISOMAP: A New ISOMAP-Based Data Visualization Algorithm with Less Sensitivity to the Neighborhood Size.

    Dianzi Xuebao(Acta Electronica Sinica) 2006, 34(8):1497-1501. OpenURL

  15. Fred A, Jain A: Combining Multiple Clusterings Using Evidence Accumulation.

    IEEE Transactions on Pattern Analysis and Machine Intelligence 2005, 27(6):835-850. PubMed Abstract | Publisher Full Text OpenURL

  16. Freund Y, Schapire R: A decision-theoretic generalization of on-line learning and an application to boosting.

    Proc 2nd European Conf Computational Learning Theory 1995, 23-37. OpenURL

  17. Breiman L: Bagging predictors.

    Machine Learning 1996, 24(2):123-140. OpenURL

  18. Dietterich T: Ensemble Methods in Machine Learning.

    Proc 1st Int'l Workshop on Multiple Classifier Systems 2000, 1-15. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. MacQueen J: Some Methods for Classification and Analysis of MultiVariate Observations.

    Proc Fifth Berkeley Symposium on Mathematical Statistics and Probability 1967, 281-297. OpenURL

  20. Fern X, Brodley C: Random Projection for High Dimensional Data Clustering: A Cluster Ensemble Approach.

    Proc 20th Int'l Conf Machine Learning 2003, 186-193. PubMed Abstract OpenURL

  21. Venna J, Kaski S: Local multidimensional scaling.

    Neural Networks 2006, 19(6):889-899. PubMed Abstract | Publisher Full Text OpenURL

  22. Samko O, Marshall A, Rosin P: Selection of the optimal parameter value for the Isomap algorithm.

    Pattern Recognition Letters 2006, 27(9):968-979. Publisher Full Text OpenURL

  23. Kouropteva O, Okun O, Pietikainen M: Selection of the Optimal Parameter Value for the Locally Linear Embedding Algorithm.

    Proc 1st Int'l Conf Fuzzy Systems and Knowledge Discovery 2002, 359-363. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. de Silva V, Tenenbaum J: Global Versus Local Methods in Nonlinear Dimensionality Reduction.

    Proc 15th Conf Adv Neural Information Processing Systems (NIPS) 2003, 705-712. OpenURL

  25. Zelnik-Manor L, Perona P: Self-tuning spectral clustering.

    Proc 17th Conf Adv Neural Information Processing Systems (NIPS) 2004, 1601-1608. OpenURL

  26. Geng X, Zhan DC, Zhou ZH: Supervised nonlinear dimensionality reduction for visualization and classification.

    IEEE Transactions on Systems, Man, and Cybernetics: Part B, Cybernetics 2005, 35(6):1098-107. Publisher Full Text OpenURL

  27. de Ridder D, Kouropteva O, Okun O, Pietikainen M, Duin R: Supervised Locally Linear Embedding.

    Proc Artificial Neural Networks and Neural Information Processing 2003, 333-341. OpenURL

  28. Athitsos V, Alon J, Sclaroff S, Kollios G: BoostMap: An Embedding Method for Efficient Nearest Neighbor Retrieval.

    IEEE Transactions on Pattern Analysis and Machine Intelligence 2008, 30(1):89-104. PubMed Abstract | Publisher Full Text OpenURL

  29. Lawrence N: Spectral Dimensionality Reduction via Maximum Entropy.

    Proc 14th Intn'l Conf Artificial Intelligence and Statistics (AISTATS) 2011, 51-59. OpenURL

  30. Mao K, Liang F, Mukherjee S: Supervised Dimension Reduction Using Bayesian Mixture Modeling.

    Proc 13th Intn'l Conf Artificial Intelligence and Statistics (AISTATS) 2010, 501-508. OpenURL

  31. Blum A, Mitchell T: Combining labeled and unlabeled data with co-training.

    Proc 11th Annual Conf Computational Learning Theory 1998, 92-100. OpenURL

  32. Hou C, Zhang C, Wu Y, Nie F: Multiple view semi-supervised dimensionality reduction.

    Pattern Recognition 2009, 43(3):720-730. OpenURL

  33. Wachinger C, Yigitsoy M, Navab N: Manifold Learning for Image-Based Breathing Gating with Application to 4D Ultrasound.

    Proc 13th Int'l Conf Medical Image Computing and Computer-Assisted Intervention (MICCAI) 2010, 26-33. PubMed Abstract | Publisher Full Text OpenURL

  34. Wachinger C, Navab N: Manifold Learning for Multi-Modal Image Registration.

    Proc 11th British Machine Vision Conference (BMVC) 2010, 82.1-82.12. OpenURL

  35. Cover T, Hart P: Nearest neighbor pattern classification.

    IEEE Transactions on Information Theory 1967, 13(1):21-27. OpenURL

  36. Tiwari P, Rosen M, Madabhushi A: Consensus-locally linear embedding (C-LLE): application to prostate cancer detection on magnetic resonance spectroscopy.

    Proc 11th Int'l Conf Medical Image Computing and Computer-Assisted Intervention (MICCAI) 2008, 330-8. PubMed Abstract OpenURL

  37. Ho TK: The random subspace method for constructing decision forests.

    IEEE Transactions on Pattern Analysis and Machine Intelligence 1998, 20(8):832-844. Publisher Full Text OpenURL

  38. Levina E, Bickel P: Maximum Likelihood Estimation of Intrinsic Dimension.

    Proc 17th Conf Adv Neural Information Processing Systems (NIPS) 2005, 777-784. OpenURL

  39. Kuncheva L: Combining pattern classifiers: methods and algorithms. Wiley-Interscience; 2004. OpenURL

  40. Halkidi M, Batistakis Y, Vazirgiannis M: On Clustering Validation Techniques.

    Journal of Intelligent Information Systems 2001, 17(2):107-145. Publisher Full Text OpenURL

  41. Patel JK, Read CB: Handbook of the normal distribution. Marcel Dekker; 1996. OpenURL

  42. Yang C, Duraiswami R, Gumerov NA, Davis L: Improved fast gauss transform and efficient kernel density estimation.

    Proc 9th IEEE Int'l Conf Computer Vision (ICCV) 2003, 664-671. OpenURL

  43. Comaniciu D, Meer P: Mean shift: a robust approach toward feature space analysis.

    IEEE Transactions on Pattern Analysis and Machine Intelligence 2002, 24(5):603-619. Publisher Full Text OpenURL

  44. Haralick RM, Shanmugam K, Dinstein I: Textural Features for Image Classification.

    IEEE Transactions on Systems, Man and Cybernetics 1973, 3(6):610-621. OpenURL

  45. Madabhushi A, Feldman M, Metaxas D, Tomaszeweski J, Chute D: Automated Detection of Prostatic Adenocarcinoma from High-Resolution Ex Vivo MRI.

    IEEE Transactions on Medical Imaging 2005, 24(12):1611-1625. PubMed Abstract | Publisher Full Text OpenURL

  46. Herlidou-Meme S, Constans JM, Carsin B, Olivie D, Eliat PA, Nadal-Desbarats L, Gondry C, Le Rumeur E, Idy-Peretti I, de Certaines JD: MRI texture analysis on texture test objects, normal brain and intracranial tumors.

    Magnetic Resononance Imaging 2003, 21(9):989-993. Publisher Full Text OpenURL

  47. Dai JJ, Lieu L, Rocke D: Dimension reduction for classification with gene expression microarray data.

    Statistical Applications in Genetics and Molecular Biology 2006, 5(1):Article6. PubMed Abstract | Publisher Full Text OpenURL

  48. Carballido-Gamio J, Belongie S, Majumdar S: Normalized cuts in 3-D for spinal MRI segmentation.

    IEEE Transactions on Medical Imaging 2004, 23(1):36-44. PubMed Abstract | Publisher Full Text OpenURL

  49. Kwan R, Evans A, Pike G: MRI simulation-based evaluation of image-processing and classification methods.

    IEEE Transactions on Medical Imaging 1999, 18(11):1085-1097. PubMed Abstract | Publisher Full Text OpenURL

  50. Chappelow J, Bloch BN, Rofsky N, Genega E, Lenkinski R, DeWolf W, Madabhushi A: Elastic registration of multimodal prostate MRI and histology via multiattribute combined mutual information.

    Medical Physics 2011, 38(4):2005-2018. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Liu H, Li J, Wong L: A comparative study on feature selection and classification methods using gene expression profiles and proteomic patterns.

    Genome Informatics 2002, 13:51-60. PubMed Abstract OpenURL

  52. Zhao H: Combining labeled and unlabeled data with graph embedding.

    Neurocomputing 2006, 69(16):2385-2389. Publisher Full Text OpenURL

  53. Eisen M, Spellman P, Brown P, Botstein D: Cluster analysis and display of genome-wide expression patterns.

    Proceedings of the National Academy of Sciences of the United States of America 1998, 95(25):14863-14868. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Monaco JP, Tomaszewski JE, Feldman MD, Hagemann I, Moradi M, Mousavi P, Boag A, Davidson C, Abolmaesumi P, Madabhushi A: High-throughput detection of prostate cancer in histological sections using probabilistic pairwise Markov models.

    Medical Image Analysis 2010, 14(4):617-629. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Moakher M, Batchelor PG:

    Symmetric Positive-Definite Matrices: From Geometry to Applications and Visualization. 2006. OpenURL

  56. Bengio Y, Paiement J, Vincent P, Delalleau O, Le Roux N, Ouimet M: Out-of-Sample Extensions for LLE, Isomap, MDS, Eigenmaps, and Spectral Clustering.

    Proc 16th Conf Adv Neural Information Processing Systems (NIPS) 2004, 177-184. OpenURL