Comparisons of Non-Gaussian Statistical Models in DNA Methylation Analysis
2014

Zhanyu Ma
Pattern Recognition and Intelligent System Lab
Beijing University of Posts and Telecommunications
No. 10 Xitucheng Road100876BeijingChina

Andrew E Teschendorff firstname.lastname@example.com
Shanghai Institute for Biological Sciences
Computational Systems Genomics, CAS-MPG Partner Institute for Computational Biology
Chinese Academy of Sciences
320 Yue Yang Road200031ShanghaiChina

Statistical Genomics Group
UCL Cancer Institute
Paul O'Gorman Building

University College London
72 Huntley StreetWC1E 6BTLondonUK

Hong Yu
Pattern Recognition and Intelligent System Lab
Beijing University of Posts and Telecommunications
No. 10 Xitucheng Road100876BeijingChina

Jalil Taghia firstname.lastname@example.com
Communication Theory Lab
KTH -Royal Institute of Technology
Osquldas väg 1010044StockholmSweden

Jun Guo email@example.com.
Pattern Recognition and Intelligent System Lab
Beijing University of Posts and Telecommunications
No. 10 Xitucheng Road100876BeijingChina

Comparisons of Non-Gaussian Statistical Models in DNA Methylation Analysis

Int. J. Mol. Sci
15201410.3390/ijms150610835Received: 24 March 2014; in revised form: 12 May 2014 / Accepted: 10 June 2014 /OPEN ACCESS International Journal of Molecular Sciences * Author to whom correspondence should be addressed; E-Mail: email@example.com;non-Gaussian statistical modelsdimension reductionunsupervised learningfeature selectionDNA methylation analysis
As a key regulatory mechanism of gene expression, DNA methylation patterns are widely altered in many complex genetic diseases, including cancer. DNA methylation is naturally quantified by bounded support data; therefore, it is non-Gaussian distributed. In order to capture such properties, we introduce some non-Gaussian statistical models to perform dimension reduction on DNA methylation data. Afterwards, non-Gaussian statistical model-based unsupervised clustering strategies are applied to cluster the data. Comparisons and analysis of different dimension reduction strategies and unsupervised clustering methods are presented. Experimental results show that the non-Gaussian statistical model-based methods are superior to the conventional Gaussian distribution-based method. They are meaningful tools for DNA methylation analysis. Moreover, among several non-Gaussian methods, the one that captures the bounded nature of DNA methylation data reveals the best clustering performance.

Keywords: non-Gaussian statistical models; dimension reduction; unsupervised learning; feature selection; DNA methylation analysis

Introduction

DNA methylation is a covalent modification of DNA, which can regulate the expression of genes [1]. Recently, DNA methylation has attracted considerable interest, due its role in the etiology of complex diseases, specially cancer [2], but also because it can be easily measured genome-wide from limited amounts of DNA, allowing measurements in clinical specimens [3].

Relatively little is known about the taxonomy of cancers at the DNA methylation level. Hence, there is currently a strong interest in performing unsupervised clustering of large-scale DNA methylation data sets in order to identify novel cancer subtypes. DNA methylation data is quantified naturally in terms of a beta-distribution. The DNA methylation beta-value, β, at a specific genomic locus provides an estimate of the fraction of cells that have that locus methylated. Although studies have considered using the logit-transform y = log 2 β/(1 − β) instead of β for subsequent statistical inferences [4], it was shown in Zhuang et al. [5] that the logit-basis can, under certain circumstances, lead to worse inference, as it can aggravate the effects of outliers (i.e., β values close to zero or one): from a biological perspective, an outlier at β = 0.999 is not more interesting than one at β = 0.9, yet on the logit scale, they would be widely separated. Moreover, analyzing DNA methylation data in terms of beta-values helps interpretability. Thus, there is considerable benefit in performing statistical inferences from the original beta-valued distributions. As a result of this, normalization, feature selection and clustering methods designed for beta-valued DNA methylation data have recently been investigated [6][7][8][9][10][11][12][13][14][15].

However, there still remains a significant shortage of methods, specially for the dimensional reduction of large DNA methylation data sets. For instance, blind source separation (BSS) [16,17], independent/principal component analysis (ICA/PCA) [18,19] and nonnegative matrix factorization (NMF) [20,21] techniques have been extensively studied in the gene expression field. Due to the bounded support property of the beta-valued data, the DNA methylation level cannot be efficiently described by these existing dimension reduction methods, which mainly assume the Gaussian distribution of the data. Therefore, the analysis results based on such mismatched models are not sufficiently promising.

Gaussian distribution is the ubiquitous probability distribution used in statistics [22][23][24]. It has an analytically tractable probability density function (pdf), and analysis based on it can be derived in an explicit form. In practice, not all of the data we need to model are Gaussian distributed [25,26]. Recent research showed that, when processing the non-Gaussian distributed data, applying suitable non-Gaussian distribution to model the data can lead to better performance than that obtained by a conventional Gaussian distribution. The advantages of applying a non-Gaussian distribution have been demonstrated in various real-life applications, including image processing [27,28], speech coding [29], document analysis [30], communication and compressive sensing systems [31], complex network analysis [32], decision-making in expert systems [33,34] and biomedical signal processing [35,36]. In particular, non-Gaussian statistical models have also been widely applied in bioinformatics [12,14,15], especially to the analysis of omics data. Indeed, omics data types are rarely Gaussianly distributed.

In this paper, we will introduce and compare some machine learning methods, which are based on non-Gaussian statistical models, for DNA methylation data analysis. The analysis of DNA methylation data includes two parts. (1) Dimensional reduction: DNA methylation array data is high-dimensional, typically involving on the order of 25 k up to 500 k dimensions (and even higher). As with other omics data, the number of samples is typically on the order of 100. However, typically, most of the salient variability in the data, e.g., variation distinguishing cancer from normal samples, or distinguishing different cancer phenotypes, is captured by a much lower-dimensional space. Hence, we need to perform some forms of dimension reduction; (2) Unsupervised clustering: Cancers especially are known to be highly heterogeneous [37]. Hence, we expect an effective and accurate unsupervised learning method to reveal such heterogeneity. Thus, although normal and cancer should be well discriminated, cancer samples may form multiple distinct clusters.

Results and Discussion

Data Description and Preprocessing

The DNA methylation data is obtained from Gene Expression Omnibus (GEO) website [38]. GEO is a public functional genomics data repository supporting MIAME-compliant data submissions. Dataset GSE32393 [39] is used for the evaluation of the above-mentioned dimension reduction and unsupervised clustering methods. DNA methylation profiles were obtained across approximately 27, 578 CpGs in breast tissues from women with and without breast cancer. Breast tissue samples were drawn from 113 breast cancers and 23 non-neoplastic breast tissues.

We considered a DNA methylation data matrix over 5000 dimensions (specifically, CpG dinucleotides). The 5000 CpGs were selected as those with the highest variance across the 136 samples. A data matrix X with a size of 5000 × 136 was obtained. Figure 1 illustrates the effect of dimension reduction. The data points in Figure 1b are separated better than those in Figure 1a, which indicates that dimension reduction can potentially improve the clustering performance. From here and thereafter, the stochastic neighbor embedding (t-SNE) method [40] (the t-SNE method can only approximately illustrate the high-dimensional data for visualization convenience; it is not the exact representation of the relations between data points) is applied to visualize the high-dimensional data in two-dimensional space.

In each dimension reduction method, we need to specify the number of dimensions for which to search. The random matrix theory (RMT) [41] was used for determining the underlying number of dimensions. Although data is distinctly non-Gaussian (even after mean-centering each CpG), RMT provides a reasonable approximation of the dimensionality, as shown by us previously in [16]. For our data matrix of 5000 × 136, we estimated a total of 14 dimensions. Hence, we set K = 14 in all of the following experiments.

Unsupervised Clustering

The level of DNA methylation can discriminate normal and cancer samples [42,43]. Cancers especially are also known to be highly heterogeneous [37], which means that the cancer samples should be grouped into more than one cluster. The above obtained matrix X still has more features than samples. Therefore, we need to further reduce the dimension and then apply an advanced clustering method to group samples appropriately.

PCA (Principal Component Analysis) + VBGMM (Variational Bayesian Gaussian Mixture Model)

When implementing dimension reduction with Gaussian assumptions, PCA [22] is the widely used method. After taking eigenvalue decomposition on the covariance matrix of the observed data, PCA keeps K eigenvectors that correspond to the K largest eigenvalues. With PCA, a matrix of size 14 × 136 was obtained. For this reduced feature matrix, the variational Bayesian Gaussian mixture model (VBGMM) [22] was applied to estimate the number of clusters, as well as to cluster the data. The VBGMM estimated 10 clusters. Figure 2a illustrates the reduced features' distribution, and the clustering results are shown in Figure 2b. The PCA + VBGMM misclustered nine samples, out of which nine cancer samples are clustered as normal ones and no normal sample is recognized as a cancer one. We applied the BG-NMF method (this work was presented in; it is used as the benchmark for non-Gaussian methods in our paper) to the above mentioned 5000 × 136 matrix. We set the number of basis vectors equal to 14 when applying the BG-NMF method.

Setting the number of basis vectors equal to 14 and applying BG-NMF to X T resulted in a 136 × 14 pseudo-basis matrix and a 14 × 5000 excitation matrix. The hypothesis is that the dimensionally-reduced basis matrix, whose element remains bounded, supported and is assumed to be beta distributed, captures the salient patterns of variation.

The benchmarked RPBMM algorithm was applied to estimate the final clusters of the reduced 136 × 14 matrix, which is illustrated in Figure 3a. The clustering is carried out on a 14-dimensional space. Table 1 lists the final results. RPBMM inferred a total of nine clusters. In summary, five samples (out of 136 samples) were misclustered. Four cancer samples were identified as normal samples, and only one normal sample is identified as cancer. This result is the same as the one obtained by applying RPBMM directly on the original 136 × 5000 data matrix. However, the overall time consumption is reduced from 1275 s to 139 s. Hence, dimension reduction via BG-NMF cannot only yield a convincing clustering result, but also facilitate the calculation.

BG-NMF + VBBMM (Variational Bayesian Estimation Framework for BMM)

One disadvantage of the RPBMM method is that it will estimate the number of mixture components in a recursive manner. It employs the wtdBIC to decide whether to further split one mixture component into two or not. The VBBMM can potentially estimate the model complexity automatically. After convergence, each mixture component corresponds to one cluster. Hence, we applied the VBBMM method on the 136 × 14 pseudo-basis matrix to cluster the samples. The initial settings of the number of mixture component is 15. Eventually, the VBBMM method estimated nine clusters (by removing components whose mixture weights are smaller than 0.01), which is the same as that inferred by RPBMM. Five samples are misclustered, out of which four cancer samples are classified as normal ones, and only one normal sample is recognized as a cancer one. The overall computational time in this scenario is about 124 seconds, which is faster than BG-NMF + RPBMM. This is the main advantage of applying the BG-NMF + RPBMM method. The clustering results are illustrated in Figure 3.

SC (Spectral Clustering) + VBvMM (Variational Inference Framework-based Bayesian Analysis of the vMF Mixture Model)

An alternative way for dimension reduction is to apply the SC method. As the RMT estimated K = 14, we set the target dimension K in the SC method to be 14. With the method described in Algorithm 1, a feature matrix of a size of 14 × 136 was obtained. This matrix is visualized in Figure 4a. The parameter σ was empirically set as the squared root of the data's variance in matrix B, where
B ij = ∥x i − x j ∥ 2 .
This choice is due to the fact that the affinity matrix A in Algorithm 1 takes the form of a Gaussian kernel function. Therefore, it is natural to consider the variance of data as σ 2 . The choice of the optimal σ in the SC method is an open topic for future research. The reduced feature from SC has a unit length, which means the l 2 norm of each feature equals one. To model such property efficiently, the VBvMM is used to realize the unsupervised clustering. After convergence, VBvMM determined 14 clusters. Seven normal samples are clustered as cancer ones, while zero cancer samples are clustered as normal ones. The overall number of misclustered samples is seven. Figure 4b illustrates the clustering result.

SC + VBWMM (Variational Bayesian Estimation of WMM)

Assuming what we observed from the SC resulting feature is only one of the mirrored pairs (with respect to the origin), the VBWMM can be applied to model the axially symmetric data. With a similar approach as the SC + VBvMM, we have six cancer samples clustered as normal ones and none of the normal samples inferred as the cancer one. The overall misclustered samples are six. The VBWMM inferred seven clusters in the end. The clustering results are shown in Figure 4c.

Discussion

The comparisons of the above-mentioned four methods are listed in Table 2. All of the above-mentioned algorithms can yield appropriate unsupervised clustering results. In general, the non-Gaussian distribution-based methods are superior to the Gaussian distribution-based method. This is due to the fact that the Gaussian distribution cannot describe the bounded/unit length property of the features properly. Among all four non-Gaussian distribution-based methods, the BG-NMF + RPBMM/BG-NMF + VBBMM methods outperform the SC + VBvMM/SC + VBWMM methods. This is mainly due to the fact that bounded support is an important property of DNA methylation data [5]. The BG-NMF methods yield reduced features that retain the bounded property of the data, while the SC-related methods produced a feature (which has a unit length) that does not directly reflect the bounded support property. The BG-NMF + VBBMM performs the best among all three proposed methods. It also outperforms the benchmarked BG-NMF + RPBMM method in terms of computational cost, which is because of the advantage that the variational Bayesian method can automatically determine the complexity. When looking at the misclustered samples, all of the BG-NMF related clustering methods miscluster four or five cancer samples to normal and miscluster one normal sample to cancer, while the SC-related method estimated six or seven cancer samples to normal, but no normal sample to cancer. Misclustering happens since the data are highly-dimensionally correlated. Although we have reduced the dimensions to remove redundant features, it is still difficult to separate one type of data from the other. The SC-related methods, however, do not miscluster any normal sample to cancer. We speculate that this is because the SC method embedded the data in a tight manner, so that a relatively "clearer" positive/negative boundary can be obtained than the BG-NMF method. On the one hand, BGNMF-related methods have overall better clustering performance than the SC-related methods, but misclustered data in both ways. On the other hand, SC-related methods do not cluster any normal data to cancer, but have relatively worse overall accuracy. These observations motivate us to improve the unsupervised clustering method so that better clustering results can be obtained.

In summary, for DNA methylation analysis, the bounded nature of the data plays an important role. Thus, such a property should be retained in both the dimension reduction and clustering methods. Furthermore, an appropriate unsupervised learning method is required for revealing the heterogeneity more accurately.

Experimental Section

Non-Gaussian Statistical Distributions

The Gaussian distribution (both univariate and multivariate) has a symmetrical "bell" shape, and the variable's definition is on the interval (−∞, ∞). Non-Gaussian statistical distributions refer to a set of distributions that have special properties that the Gaussian distribution cannot characterize. For example, the beta distribution is defined on the interval [0, 1] (in a general form, the beta distribution could have definition on any interval [a,b]; after linear scaling, it can be represented with the standard beta distribution [44]) and can have a symmetric or asymmetric shape [27]. The Dirichlet distribution, which is a multivariate generalization of the beta distribution, has a pdf with respect to the Lebesgue measure on the Euclidean space [45]. The gamma distribution is defined on the interval (0, ∞), and the shape cannot be symmetric [46]. To model data whose l 2 norm equals one, the von Mises-Fisher (vMF) distribution [47] and the Watson distribution [48] are usually applied. These distributions show characteristics that are significantly different from a Gaussian distribution.

In the remaining part of this section, we will introduce some typical non-Gaussian distributions that can be applied in DNA methylation analysis.

Beta Distribution

The beta distribution is characterized by two positive shape parameters u and v. The pdf of the beta distribution is:
Beta(x; u, v) = Γ(u + v) Γ(u)Γ(v) x u−1 (1 − x) v−1(1)
where Γ(·) is the gamma function. The beta distribution has a flexible shape, which is shown in Figure 5.

In real applications, the beta distribution can be applied to model the distribution of a gray image pixel [49], to describe the probability of human immunodeficiency virus (HIV) transmission [50] and to capture the bounded property of the DNA methylation level [8,15].

vMF (von Mises-Fisher) Distribution

The vMF distribution is considered a popular distribution in the family of directional distributions [47,51]. The data following a vMF distribution are located on a unit hypersphere. Hence, the vMF variable's l 2 norm equals one, i.e., ∥x∥ 2 = 1. The vMF distribution contains two parameters, namely the mean direction µ and the concentration parameter λ. The pdf of a K-dimensional vMF distribution can be expressed as:
F(x | µ, λ) = c K (λ) e λµ T x(2)

c K (λ) = λ K−2 2 (2π) K 2 I K−2 2 (λ)(3)
where I ν (·) represents the modified Bessel function of the first kind of order ν [52]. The pdf of the vMF distribution is illustrated in Figure 6. In information retrieval applications, the vMF distribution can be applied to model the cosine similarity for the clustering of text documents [53,54]. It can also be applied in modeling the gene expression data, which has been shown to have directional characteristics [55]. Samples generated from F(µ, 400) (shown by green colors) are highly concentrated around the mean direction, while for samples generated from F(µ, 4) (shown by blue colors), the distribution of samples on the sphere is more uniform around the mean direction.

Watson Distribution

Observations on the sphere might have an additional structure, such that the unit vectors x and −x are equivalent. In other words, that is ±x that are observed. Here, we need probability density functions for x on S p−1 which are axially symmetric, that is f (−x) = f (x). In such cases, the p-dimensional observation ±x can be regarded as being on the projective space P p−1 , which is obtained by identifying opposite points on the sphere S p−1 .

One of the simplest distributions for axial data, with a rotational symmetry property, is the (Dimroth-Scheidegger-) Watson distribution. The Watson distribution is a special case of the Bingham distribution [56], which is developed for axial data with no rotational symmetry property.

A random vector x ∈ P p−1 , or equivalently ±x ∈ S p−1 , has the (p − 1)-dimensional Watson distribution W p (µ, κ), with the mean direction µ and the concentration parameter κ, if its probability density function is:
W p (x | µ, κ) = Γ(p) (2π) p 1 F 1 (r, p, κ) (κ)e κ∥µ † x∥ 2(4)

1 F 1 (r, p, κ) = ∑ j≥0 r (j) p (j) κ j j!(5)
where r (j) ≡ Γ(r+j) Γ(r) is the raising factorial. Similar to the case of vMF distributions, for κ > 0, as κ → 0, W p (µ, κ) reduces to the uniform density, and as κ → ∞, W p (µ, κ) tends to a point density. For κ < 0, as κ → −∞, the density concentrates around the great circle orthogonal to the mean direction ([59] (Chapter 9.4)). The samples generated from Watson distribution are shown in Figure 7.

Non-Gaussian Dimension Reduction Methods

When analyzing DNA methylation data, the high-dimensional property presents mathematical challenges, as well as opportunities. The main purpose of applying dimension reduction methods on microarray data is to extract the core features driving interesting biological variability [61]. Such methods include principal component analysis (PCA) [22,62], nonnegative matrix factorization (NMF) [5,63] and singular value decomposition (SVD) [64].

Nonnegative Matrix Factorization for Bounded Support Data

Unlike PCA or ICA, NMF reveals the data's nonnegativity during dimension reduction. Traditional NMF decomposes the data matrix into a product of two nonnegative matrices as:
X P ×T ≈ W P ×K V K×T(6)
where X P ×T , W P ×K and V K×T contain nonnegative values X pt , W pk and V kt , respectively, and p = 1, . . . , P , t = 1, . . . , T , k = 1, . . . , K, K << T . The DNA methylation data are naturally bounded on interval [0, 1]. Conventional NMF strategies do not take such a nature into account. In order to capture such a bounded feature explicitly, we proposed an NMF for bounded support data [65]. Each bounded support element X pt is assumed to be generated from a beta distribution with parameters a pt and b pt . With an observation matrix X P ×T , two parameter matrices a and b of size P × T are obtained, respectively. Each parameter matrix, rather than the observation matrix, is decomposed into a product of a basis matrix and an excitation matrix as:
a P ×T ≈ A P ×K H K×T b P ×T ≈ B P ×K H K×T(7)
With the above description, we assume that the matrix X (with element X pt ∈ [0, 1]) is drawn according to the following generative model:
A pk ∼ Gamma(A pk ; µ 0 , α 0 ) B pk ∼ Gamma(B pk ; ν 0 , β 0 ) H kt ∼ Gamma(H kt ; ρ 0 , ζ 0 ) X pt ∼ Beta(X pt ; ∑ k A pk H kt , ∑ k B pk H kt )(8)
where Gamma(x; k, θ) is the gamma density with parameters k, θ defined as:
Gamma(x; k, θ) = θ k Γ(k) x k−1 e −θx , k, θ > 0(9)
As the data is assumed to be beta distributed and the parameters of the beta distribution are assumed to be gamma distributed, this model is named BG-NMF.

For BG-NMF, the variational inference (VI) method [66] is applied to estimate the posterior distributions. The expected value of X pt is X pt = apt apt+bpt . If we take the point estimate to A pk , B pk and H kt , then the expected value of X pt can be approximated as:
X pt ≈ ∑ k A pk H kt ∑ k A pk H kt + ∑ k B pk H kt(10)
which can be expressed in matrix form as:
X ≈ (A H) ⊘ (A H + B H)(11)
where ⊘ means element-wise division. When placing sparsity constraints on the columns in H, the reconstruction in Equation (11) could be approximated as:
X ≈ [ A ⊘ (A + B) ] H = W H(12)
Hence, the resulting pseudo-basis matrix W is low-dimensional while retaining the bounded support constraint.

Spectral Clustering for Non-Gaussian Reduced Features

Recently, spectral clustering (SC) has become one of the most popular clustering algorithms [67]. It is an alternative method for the K-means algorithm. When the natural clusters in R L are not corresponding to the convex region, the K-means algorithm cannot provide satisfactory clustering results. However, when mapping the data points to R K space via SC, they may form tight clusters [68]. SC analyzes the affinity matrix of data. Assuming that the data are likely to be clustered in K-dimensional space, the reduced features, each of which is K-dimensional, are extracted by taking eigenvalue analysis of an intermediate matrix M. The reduced features will then be used for clustering, with conventional methods like K-means. The feature extraction procedure via the SC method [68] is summarized in Algorithm 1. With the above extracted features, the task of data clustering can be carried out in the reduced R K space.

Create the affinity matrix
A N ×N , where A ij = { e −∥x i −x j ∥ 2 2σ 2 i ̸ = j 0 i = j ; 2. Construct the intermediate matrix M = D − 1 2 AD − 1 2 ,
where D is a diagonal matrix whose (i, i)th element is the summation of the i-th row of A;

3. Apply eigenvalue analysis on M and create a matrix Y K×N , which contains K eigenvectors corresponding to the largest K eigenvalues; 4. Form a matrix Z from Y by normalizing each column of Y. Each column of Z has a unit length.

Output: Reduced K-dimensional feature z n for each data point x n .

Non-Gaussian Statistical Models for Unsupervised Clustering

The ultimate goal of dimension reduction is to benefit the clustering of DNA methylation data. The dimension reduction methods introduced above yield features with special properties. The BG-NMF method provides a basis matrix (see Equation (12)), every element of which is on the interval [0, 1]. The SC method constructs feature vectors with a unit length in the reduced K-dimensional space. These data are obviously non-Gaussian distributed. It has been shown in several studies that applying a suitable non-Gaussian distribution to model the non-Gaussian distributed data can improve the corresponding performance [25,26,31]. To this end, we apply the non-Gaussian distribution that can nicely describe the data distribution to realize the clustering task in the reduced feature space.

Recursive Partitioning Beta Mixture Model

For the beta-valued DNA methylation data, it is natural to consider the beta distribution as a candidate to model the underlying distribution. Since the DNA methylation data that show an obvious normal/cancer status are multi-modal, a beta mixture model (BMM) can be applied for modeling. One mixture component represents one cluster. In unsupervised clustering, selecting the optimal number of clusters is a big challenge. One popular method designed for such purpose is the recursive partitioning beta mixture model (RPBMM) [12]. The RPBMM navigates clusters in a BMM. It treats each mixture component as a node and takes the weighted version of Bayesian information criterion (wtdBIC) to decide whether to stop further recursive partitioning of the current node. In the end, the optimal number of clusters is determined. This is a recursive unsupervised learning method. In this paper, the clustering results based on the RPBMM method are taken as the benchmark.

Variational Beta Mixture Model

Another way of carrying out model selection is to employ the variational Bayesian estimation framework for BMM (VBBMM). Under this circumstance, the joint posterior distribution of the weighting factors is modeled by a sparse Dirichlet distribution, so that the component with a very small weight will be pruned from the mixture model. In [27], the extended variational inference strategy is applied to derive an analytically tractable solution for estimating the parameters in a BMM. The proposed algorithm can remove the redundant component. This closed-form solution avoids the cross-validation in the methods using BIC as the criterion (e.g., [12,44]). Hence, it is computationally efficient, and the unsupervised clustering is facilitated.

Variational von Mises-Fisher Mixture Model

The data with its l 2 norm equaling one has a directional property. The von Mises-Fisher (vMF) distribution is suitable for such a type of data [47,69]. In order to decide the model complexity (in terms of free parameters) automatically based on the data, we proposed a variational inference framework-based Bayesian analysis of the vMF mixture model (VBvMM) in [47]. This method can potentially determine the model complexity and avoid the over-fitting problem associated with conventional approaches based on the expectation maximization. This variational vMM is a suitable model for the unsupervised clustering of directional data.

Variational Watson Mixture Model

The Watson distribution is a simple distribution for modeling axially symmetric data on the unit hypersphere ([59] (Chapter 9.4)). By assuming that any data point has its axial mirror, it is natural to model the distribution of the data with the unit length (i.e., l 2 norm equals one) and its axial mirror by the Watson distribution. Similarly, when such data are multi-modally distributed, a Watson mixture model (WMM) can be applied. With a variational inference framework, Taghia et al. [60] proposed the variational Bayesian estimation of WMM (VBWMM), where the model complexity can be determined by pruning the irrelevant components. This variational WMM can also be applied for the purpose of unsupervised clustering.

Conclusions

Cancer is characterized by alterations at the DNA methylation level. A Gaussian distribution, in general, cannot describe the DNA methylation data appropriately.

Hence, the Gaussian distribution-based unsupervised clustering does not provide convincing performance.

For the purpose of efficiently clustering DNA methylation data, we proposed several dimension reduction methods and consequent unsupervised learning methods, which are all based on non-Gaussian distributions. They all perform better than the Gaussian distribution-based method. In the dimension reduction step, both the BG-NMF and the SC methods can remove the redundant dimensions efficiently. In unsupervised clustering, the VBBMM method, the VBvMM method and the VBWMM method can all reveal the heterogeneity of the DNA methylation data appropriately. Clustering performance demonstrates that the proposed non-Gaussian distribution-based methods are meaningful tools for analyzing DNA methylation data. Experimental results also show that the BG-NMF + VBBMM method performs the best among all of the proposed methods and is faster than the benchmarked BG-NMF + RPBMM method. Furthermore, for the reduced features inferred from both the BG-NMF method and the SC method, the consequent unsupervised clustering method needs to be improved, so that better clustering accuracy can be obtained.

Moreover, the methodology introduced in this paper can be easily extended to analyze other DNA methylation data sets. Some other non-Gaussian statistical models can also be applied for such purposes.

Figure 1 .
1Comparisons of DNA methylation data before and after dimension reduction. (a) Visualization of the original 27, 578 dimensional data; (b) visualization of the reduced 5000 dimensional data.

Figure 2 .
2Illustration of the clustering result via PCA (principal component analysis) + the variational Bayesian Gaussian mixture model (VBGMM). The clusters are color coded. The normal data are marked with dots, and the cancer data are marked with crosses. Samples in a larger size are those misclustered. 2. BG-NMF (Beta-Gamma-Nonnegative Matrix Factorization) + RPBMM (Recursive Partitioning Beta Mixture Model)

Figure 3 .
3Illustration of the clustering result via beta-gamma (BG)-nonnegative matrix factorization (NMF) + VBBMM (variational Bayesian estimation framework for BMM). The clusters are color coded. The normal data are marked with dots, and the cancer data are marked with crosses. Samples in the larger size are those misclustered. (a) Data visualization; (b) clustering result via BG-NMF + VBBMM.

Figure 4 .
4Illustration of the clustering result via spectral clustering (SC) + variational inference framework-based Bayesian analysis of the vMF mixture model (VBvMM) and SC + VBWMM (Watson mixture model (WMM)). The clusters are color coded. The normal data are marked with dots, and the cancer data are marked with crosses. Samples in the larger size are those misclustered. (a) Data visualization; (b) clustering result via SC + VBvMM; (c) clustering result via SC + VBWMM.

Figure 5 .
5Beta distributions for different pairs of parameters. (a) u = 5,

Figure 6 .
6Scatter plot of samples from a single von Mises-Fisher (vMF) distribution on the sphere for different concentration parameters, λ = {4, 40, 400}, and around the same mean direction µ = [0, 0, 1] T .

Figure 7 .
7Scatter plot of samples from a single distribution, W p (µ, κ), on the sphere for positive and negative concentration parameters κ, around the same mean direction, µ = [0, 1] T[60]. For larger concentration parameters, i.e., κ = 40 or κ = −40, samples are more concentrated around the mean direction (shown by the red color). For smaller concentration parameters, i.e., κ = 4 or κ = −4, samples are more uniformly distributed around the mean direction (shown by the blue color). (a) κ > 0, κ ∈ {+4, +40}; (b) κ < 0, κ ∈ {−4, −40}.

Table 1 .
1Recursive partitioning beta mixture model (RPBMM) clustering details for all of the 136 samples over 14 BG-NMF pseudo-basis vectors.Cluster Normal Cancer

rLLLL
16
4
rLLLR
6
0
rLLR
1
32
rLR
0
31
rRLLL
0
4
rRLLR
0
3
rRLRL
0
10
rRLRR
0
5
rRR
0
34

Table 2 .
2Comparisons of the clustering performance of different methods.Method
Error Rate Cancer→Normal Normal→Cancer

PCA + VBGMM
6.62%
9
0
BGNMF + RPBMM
3.68%
4
1
BGNMF + VBBMM
3.68%
4
1
SC + VBvMM
5.15%
7
0
SC + VBWMM
4.41%
6
0

AcknowledgmentsThe authors would like to thank the editor for organizing the review process and thank the anonymous reviewers for their efforts in reviewing this manuscript and providing fruitful suggestions.This work is partly supported by the "Fundamental Research Funds for the Central UniversitiesConflicts of InterestThe authors declare no conflict of interest.
CpG islands and the regulation of transcription. A M Deaton, A Bird, Genes Dev. 25Deaton, A.M.; Bird, A. CpG islands and the regulation of transcription. Genes Dev. 2011, 25, 1010-1022.

Epigenetics as a unifying principle in the aetiology of complex traits and diseases. A Petronis, Nature. 465Petronis, A. Epigenetics as a unifying principle in the aetiology of complex traits and diseases. Nature 2010, 465, 721-727.

Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. J Sandoval, H Heyn, S Moran, J Serra-Musach, M A Pujana, M Bibikova, M Esteller, Epigenetics. 6Sandoval, J.; Heyn, H.; Moran, S.; Serra-Musach, J.; Pujana, M.A.; Bibikova, M.; Esteller, M. Validation of a DNA methylation microarray for 450,000 CpG sites in the human genome. Epigenetics 2011, 6, 692-702.

Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. P Du, X Zhang, C C Huang, N Jafari, W A Kibbe, L Hou, S M Lin, 10.1186/1471-2105-11-587BMC Bioinform. 11Du, P.; Zhang, X.; Huang, C.C.; Jafari, N.; Kibbe, W.A.; Hou, L.; Lin, S.M. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinform. 2010, 11, doi:10.1186/1471-2105-11-587.

A comparison of feature selection and classification methods in DNA methylation studies using the Illumina 27k platform. J Zhuang, M Widschwendter, A E Teschendorff, 10.1186/1471-2105-13-59BMC Bioinform. 13Zhuang, J.; Widschwendter, M.; Teschendorff, A.E. A comparison of feature selection and classification methods in DNA methylation studies using the Illumina 27k platform. BMC Bioinform. 2012, 13, doi:10.1186/1471-2105-13-59.

CpGassoc: An R function for analysis of DNA methylation microarray data. R T Barfield, V Kilaru, A K Smith, K N Conneely, Bioinformatics. 28Barfield, R.T.; Kilaru, V.; Smith, A.K.; Conneely, K.N. CpGassoc: An R function for analysis of DNA methylation microarray data. Bioinformatics 2012, 28, 1280-1281.

A graphical user interface package for the analysis of array-based DNA methylation data. V Kilaru, R T Barfield, J W Schroeder, A K Smith, K N Conneely, Methlab, Epigenetics. 7Kilaru, V.; Barfield, R.T.; Schroeder, J.W.; Smith, A.K.; Conneely, K.N. MethLAB: A graphical user interface package for the analysis of array-based DNA methylation data. Epigenetics 2012, 7, 225-229.

A beta-mixture model for dimensionality reduction, sample classification and analysis. K Laurila, B Oster, C Andersen, P Lamy, T Orntoft, O Yli-Harja, C Wiuf, 10.1186/1471-2105-12-215BMC Bioinform. 12Laurila, K.; Oster, B.; Andersen, C.; Lamy, P.; Orntoft, T.; Yli-Harja, O.; Wiuf, C. A beta-mixture model for dimensionality reduction, sample classification and analysis. BMC Bioinform. 2011, 12, doi:10.1186/1471-2105-12-215.

Semi-supervised recursively partitioned mixture models for identifying cancer subtypes. D C Koestler, C J Marsit, B C Christensen, M R Karagas, R Bueno, D J Sugarbaker, K T Kelsey, E A Houseman, Bioinformatics. 26Koestler, D.C.; Marsit, C.J.; Christensen, B.C.; Karagas, M.R.; Bueno, R.; Sugarbaker, D.J.; Kelsey, K.T.; Houseman, E.A. Semi-supervised recursively partitioned mixture models for identifying cancer subtypes. Bioinformatics 2010, 26, 2578-2585.

A statistical framework for Illumina DNA methylation arrays. P F Kuan, S Wang, X Zhou, H Chu, Bioinformatics. 26Kuan, P.F.; Wang, S.; Zhou, X.; Chu, H. A statistical framework for Illumina DNA methylation arrays. Bioinformatics 2010, 26, 2849-2855.

Copy number variation has little impact on bead-array-based measures of DNA methylation. E A Houseman, B C Christensen, M R Karagas, M R Wrensch, H H Nelson, J L Wiemels, S Zheng, J K Wiencke, K T Kelsey, C J Marsit, Bioinformatics. 25Houseman, E.A.; Christensen, B.C.; Karagas, M.R.; Wrensch, M.R.; Nelson, H.H.; Wiemels, J.L.; Zheng, S.; Wiencke, J.K.; Kelsey, K.T.; Marsit, C.J. Copy number variation has little impact on bead-array-based measures of DNA methylation. Bioinformatics 2009, 25, 1999-2005.

Model-based clustering of DNA methylation array data: A recursive-partitioning algorithm for high-dimensional data arising as a mixture of beta distributions. E A Houseman, B C Christensen, R F Yeh, C J Marsit, M R Karagas, M Wrensch, H H Nelson, J Wiemels, S Zheng, J K Wiencke, 10.1186/1471-2105-9-365BMC Bioinform. 9Houseman, E.A.; Christensen, B.C.; Yeh, R.F.; Marsit, C.J.; Karagas, M.R.; Wrensch, M.; Nelson, H.H.; Wiemels, J.; Zheng, S.; Wiencke, J.K.; et al. Model-based clustering of DNA methylation array data: A recursive-partitioning algorithm for high-dimensional data arising as a mixture of beta distributions. BMC Bioinform. 2008, 9, doi:10.1186/1471-2105-9-365.

Applications of beta-mixture models in bioinformatics. Y Ji, C Wu, P Liu, J Wang, K R Coombes, Bioinformatics. 21Ji, Y.; Wu, C.; Liu, P.; Wang, J.; Coombes, K.R. Applications of beta-mixture models in bioinformatics. Bioinformatics 2005, 21, 2118-2122.

A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. A E Teschendorff, F Marabita, M Lechner, T Bartlett, J Tegner, D Gomez-Cabrero, S Beck, Bioinformatics. 29Teschendorff, A.E.; Marabita, F.; Lechner, M.; Bartlett, T.; Tegner, J.; Gomez-Cabrero, D.; Beck, S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics 2013, 29, 189-196.

A variational Bayes beta Mixture Model for Feature Selection in DNA methylation Studies. Z Ma, A E Teschendorff, 10.1142/S0219720013500054J. Bioinform. Comput. Biol. 11Ma, Z.; Teschendorff, A.E. A variational Bayes beta Mixture Model for Feature Selection in DNA methylation Studies. J. Bioinform. Comput. Biol. 2013, 11. doi:10.1142/S0219720013500054.

Independent surrogate variable analysis to deconvolve confounding factors in large-scale microarray profiling studies. A E Teschendorff, J Zhuang, M Widschwendter, Bioinformatics. 27Teschendorff, A.E.; Zhuang, J.; Widschwendter, M. Independent surrogate variable analysis to deconvolve confounding factors in large-scale microarray profiling studies. Bioinformatics 2011, 27, 1496-1505.

Blind source separation methods for deconvolution of complex signals in cancer biology. A Zinovyev, U Kairov, T Karpenyuk, E Ramanculov, Biochem. Biophys. Res. Commun. 430Zinovyev, A.; Kairov, U.; Karpenyuk, T.; Ramanculov, E. Blind source separation methods for deconvolution of complex signals in cancer biology. Biochem. Biophys. Res. Commun. 2013, 430, 1182-1187.

Elucidating the altered transcriptional programs in breast cancer using independent component analysis. A E Teschendorff, M Journée, P A Absil, R Sepulchre, C Caldas, PLoS Comput. Biol. 3161Teschendorff, A.E.; Journée, M.; Absil, P.A.; Sepulchre, R.; Caldas, C. Elucidating the altered transcriptional programs in breast cancer using independent component analysis. PLoS Comput. Biol. 2007, 3, e161.

Direct detection of DNA methylation during single-molecule, real-time sequencing. B A Flusberg, D R Webster, J H Lee, K J Travers, E C Olivares, T A Clark, J Korlach, S W Turner, Nat. Methods. 7Flusberg, B.A.; Webster, D.R.; Lee, J.H.; Travers, K.J.; Olivares, E.C.; Clark, T.A.; Korlach, J.; Turner, S.W. Direct detection of DNA methylation during single-molecule, real-time sequencing. Nat. Methods 2010, 7, 461-465.

A modified non-negative matrix factorization algorithm utilizing uncertainty estimates. G Wang, A V Kossenkov, M F Ochs, Ls-Nmf, BMC Bioinform. 7175Wang, G.; Kossenkov, A.V.; Ochs, M.F. LS-NMF: A modified non-negative matrix factorization algorithm utilizing uncertainty estimates. BMC Bioinform. 2006, 7, 175.

A bioinformatic assay for pluripotency in human cells. F J Muller, B M Schuldt, R Williams, D Mason, G Altun, E P Papapetrou, S Danner, J E Goldmann, A Herbst, N O Schmidt, Nat. Methods. 8Muller, F.J.; Schuldt, B.M.; Williams, R.; Mason, D.; Altun, G.; Papapetrou, E.P.; Danner, S.; Goldmann, J.E.; Herbst, A.; Schmidt, N.O.; et al. A bioinformatic assay for pluripotency in human cells. Nat. Methods 2011, 8, 315-317.

. C M Bishop, Pattern Recognition and Machine Learning. SpringerBishop, C.M. Pattern Recognition and Machine Learning; Springer: Berlin/Heidelberg, Germany, 2006.

Gaussian Processes for Machine Learning; Adaptive computation and machine learning. C E Rasmussen, C K Williams, MIT PressCambridge, MA, USARasmussen, C.E.; Williams, C.K.I. Gaussian Processes for Machine Learning; Adaptive computation and machine learning; MIT Press: Cambridge, MA, USA, 2006.

Probability for Statistics and Machine Learning: Fundamentals and Advanced Topics; Springer Texts in Statistics. A Dasgupta, SpringerBerlin/Heidelberg, GermanyDasGupta, A. Probability for Statistics and Machine Learning: Fundamentals and Advanced Topics; Springer Texts in Statistics; Springer: Berlin/Heidelberg, Germany, 2011.

Topics in Non-Gaussian Signal Processing. E Wegman, S Schwartz, J Thomas, SpringerLondon; London, UKWegman, E.; Schwartz, S.; Thomas, J. Topics in Non-Gaussian Signal Processing; Springer London, Limited: London, UK, 2011.

Non-Gaussian Statistical Models and Their Applications. Z Ma, Stockholm, SwedenKTH -Royal Institute of TechnologyPh.D. ThesisMa, Z. Non-Gaussian Statistical Models and Their Applications. Ph.D. Thesis, KTH -Royal Institute of Technology, Stockholm, Sweden, 2011.

Bayesian estimation of beta mixture models with variational inference. Z Ma, A Leijon, IEEE Trans. Pattern Anal. Mach. Intell. 33Ma, Z.; Leijon, A. Bayesian estimation of beta mixture models with variational inference. IEEE Trans. Pattern Anal. Mach. Intell. 2011, 33, 2160-2173.

High-dimensional unsupervised selection and estimation of a finite generalized dirichlet mixture model based on minimum message length. N Bouguila, D Ziou, IEEE Trans. Pattern Anal. Mach. Intell. 29Bouguila, N.; Ziou, D. High-dimensional unsupervised selection and estimation of a finite generalized dirichlet mixture model based on minimum message length. IEEE Trans. Pattern Anal. Mach. Intell. 2007, 29, 1716-1731.

Vector quantization of LSF parameters with a mixture of dirichlet distributions. Z Ma, A Leijon, W B Kleijn, IEEE Trans. Audio, Speech, Lang. Process. 21Ma, Z.; Leijon, A.; Kleijn, W.B. Vector quantization of LSF parameters with a mixture of dirichlet distributions. IEEE Trans. Audio, Speech, Lang. Process. 2013, 21, 1777-1790.

Latent dirichlet allocation. D M.; Y Blei, N A Jordan, M I , J. Mach. Learn. Res. 3Blei, D.M.; Y., N.A.; Jordan, M.I. Latent dirichlet allocation. J. Mach. Learn. Res. 2003, 3, 993-1022.

Performance of Non-Gaussian Distribution Based Communication and Compressed Sensing Systems. H J Kwon, San Diego, CA, USAUniversity of CaliforniaPh.D. ThesisKwon, H.J. Performance of Non-Gaussian Distribution Based Communication and Compressed Sensing Systems. Ph.D. Thesis, University of California, San Diego, CA, USA, 2013.

An activation force-based affinity measure for analyzing complex networks. J Guo, H Guo, Z Wang, 10.1038/srep00113Sci. Rep. Guo, J.; Guo, H.; Wang, Z. An activation force-based affinity measure for analyzing complex networks. Sci. Rep. 2011, doi:10.1038/srep00113.

Toward reducing failure risk in an integrated vehicle health maintenance system: A fuzzy multi-sensor data fusion Kalman filter approach for IVHMS. J A Rodger, Expert Syst. Appl. 39Rodger, J.A. Toward reducing failure risk in an integrated vehicle health maintenance system: A fuzzy multi-sensor data fusion Kalman filter approach for IVHMS. Expert Syst. Appl. 2012, 39, 9821-9836.

A fuzzy nearest neighbor neural network statistical model for predicting demand for natural gas and energy cost savings in public buildings. J A Rodger, Expert Syst. Appl. 41Rodger, J.A. A fuzzy nearest neighbor neural network statistical model for predicting demand for natural gas and energy cost savings in public buildings. Expert Syst. Appl. 2014, 41, 1813-1829.

EEG signal classification with super-Dirichlet mixture model. Z Ma, Z H Tan, S Prasad, Proceedings of the IEEE Statistical Signal Processing Workshop (SSP). the IEEE Statistical Signal Processing Workshop (SSP)Ann Arbor, MI, USA, 5-8 AuguestMa, Z.; Tan, Z.H.; Prasad, S. EEG signal classification with super-Dirichlet mixture model. In Proceedings of the IEEE Statistical Signal Processing Workshop (SSP), Ann Arbor, MI, USA, 5-8 Auguest 2012; pp. 440-443.

Predictive distribution of the dirichlet mixture model by local variational inference. Z Ma, A Leijon, Z H Tan, S Gao, J. Signal Process. Syst. 74Ma, Z.; Leijon, A.; Tan, Z.H.; Gao, S. Predictive distribution of the dirichlet mixture model by local variational inference. J. Signal Process. Syst. 2014, 74, 1-16.

The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. C Curtis, S P Shah, S F Chin, G Turashvili, O M Rueda, M J Dunning, D Speed, A G Lynch, S Samarajiwa, Y Yuan, Nature. 486Curtis, C.; Shah, S.P.; Chin, S.F.; Turashvili, G.; Rueda, O.M.; Dunning, M.J.; Speed, D.; Lynch, A.G.; Samarajiwa, S.; Yuan, Y.; et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature 2012, 486, 346-352.

. Gene Expression Omnibus. Available online. 12GEO. Gene Expression Omnibus. Available online: http://www.ncbi.nlm.nih.gov/geo/ (accessed on 12 June 2014).

The dynamics and prognostic potential of DNA methylation changes at stem cell gene loci in women's cancer. J Zhuang, A Jones, S H Lee, E Ng, H Fiegl, M Zikan, D Cibula, A Sargent, H B Salvesen, I J Jacobs, PLoS Genet. 81002517Zhuang, J.; Jones, A.; Lee, S.H.; Ng, E.; Fiegl, H.; Zikan, M.; Cibula, D.; Sargent, A.; Salvesen, H.B.; Jacobs, I.J.; et al. The dynamics and prognostic potential of DNA methylation changes at stem cell gene loci in women's cancer. PLoS Genet. 2012, 8, e1002517.

Visualizing Data using t-SNE. L Van Der Maaten, G Hinton, J. Mach. Learn. Res. 9Van der Maaten, L.; Hinton, G. Visualizing Data using t-SNE. J. Mach. Learn. Res. 2008, 9, 2579-2605.

Random matrix approach to cross correlations in financial data. V Plerou, P Gopikrishnan, B Rosenow, L A Amaral, T Guhr, H E Stanley, Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 66126Plerou, V.; Gopikrishnan, P.; Rosenow, B.; Amaral, L.A.; Guhr, T.; Stanley, H.E. Random matrix approach to cross correlations in financial data. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2002, 65, 066126.

The epigenomics of cancer. P A Jones, S B Baylin, 128Jones, P.A.; Baylin, S.B. The epigenomics of cancer. Cell 2007, 128, 683-692.

Epigenetic stem cell signature in cancer. M Widschwendter, H Fiegl, D Egle, E Mueller-Holzner, G Spizzo, C Marth, D J Weisenberger, M Campan, J Young, I Jacobs, Nat. Genet. 39Widschwendter, M.; Fiegl, H.; Egle, D.; Mueller-Holzner, E.; Spizzo, G.; Marth, C.; Weisenberger, D.J.; Campan, M.; Young, J.; Jacobs, I.; et al. Epigenetic stem cell signature in cancer. Nat. Genet. 2007, 39, 157-158.

Practical Bayesian estimation of a finite Beta mixture through gibbs sampling and its applications. N Bouguila, D Ziou, E Monga, Stat. Comput. 16Bouguila, N.; Ziou, D.; Monga, E. Practical Bayesian estimation of a finite Beta mixture through gibbs sampling and its applications. Stat. Comput. 2006, 16, 215-225.

. S Kotz, N Balakrishnan, N L Johnson, Continuous Multivariate Distributions. 1WileyModels and ApplicationsKotz, S.; Balakrishnan, N.; Johnson, N.L. Continuous Multivariate Distributions. Volume 1: Models and Applications; Wiley: New York, NY, USA, 2004.

Handbook of Statistical Distributions with Applications. K Krishnamoorthy, Statistics: A Series of Textbooks and Monographs. Krishnamoorthy, K. Handbook of Statistical Distributions with Applications; Statistics: A Series of Textbooks and Monographs;

. & Taylor, Francis, Boca Ration, FL, USATaylor & Francis: Boca Ration, FL, USA, 2010.

Bayesian estimation of the von-Mises Fisher mixture model with variational inference. J Taghia, Z Ma, A Leijon, 10.1109/TPAMI.2014.2306426IEEE Trans. Pattern Anal. Mach. Intell. Taghia, J.; Ma, Z.; Leijon, A. Bayesian estimation of the von-Mises Fisher mixture model with variational inference. IEEE Trans. Pattern Anal. Mach. Intell. 2012, doi:10.1109/TPAMI.2014.2306426.

The multivariate Watson distribution: Maximum-likelihood estimation and other aspects. S Sra, D Karp, J. Multivar. Anal. 114Sra, S.; Karp, D. The multivariate Watson distribution: Maximum-likelihood estimation and other aspects. J. Multivar. Anal. 2013, 114, 256-269.

Beta mixture models and the application to image classification. Z Ma, A Leijon, Proceedings of the IEEE International Conference on Image Processing. the IEEE International Conference on Image ProcessingCairo, EgyptMa, Z.; Leijon, A. Beta mixture models and the application to image classification. In Proceedings of the IEEE International Conference on Image Processing, Cairo, Egypt, 7-10 November 2009, pp. 2045-2048.

The role of exposure history on HIV acquisition: Insights from repeated low-dose challenge studies. R R Regoes, 10.1371/journal.pcbi.1002767PLoS Comput. Biol. 8Regoes, R.R. The role of exposure history on HIV acquisition: Insights from repeated low-dose challenge studies. PLoS Comput. Biol. 2012, 8, doi:10.1371/journal.pcbi.1002767.

. K V Mardia, P E Jupp, Directional Statistics. Mardia, K.V.; Jupp, P.E. Directional Statistics;

. John Wiley, Sons, Hoboken, NJ, USAJohn Wiley and Sons: Hoboken, NJ, USA, 2009.

Handbook of Mathematical Functions. M Abramowitz, I A Stegun, Dover PublicationsNew York, NY, USAAbramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions; Dover Publications: New York, NY, USA, 1965.

Clustering on the unit hypersphere using von Mises-Fisher distributions. A Banerjee, I S Dhillon, J Ghosh, S Sra, J. Mach. Learn. Res. 6Banerjee, A.; Dhillon, I.S.; Ghosh, J.; Sra, S. Clustering on the unit hypersphere using von Mises-Fisher distributions. J. Mach. Learn. Res. 2005, 6, 1345-1382.

The von Mises Naive Bayes classifier for angular data. P L López-Cruz, C Bielza, P Larrañaga, Advances. López-Cruz, P.L.; Bielza, C.; Larrañaga, P. The von Mises Naive Bayes classifier for angular data. In Advances in Artificial Intelligence;

. J A Lozano, J A Gámez, J A Moreno, Springer7023Berlin/Heidelberg, GermanyLozano, J.A., Gámez, J.A., Moreno, J.A., Eds.; Springer: Berlin/Heidelberg, Germany, 2011; Volume 7023, pp. 145-154.

Diametrical clustering for identifying anti-correlated gene clusters. I S Dhillon, E M Marcotte, U Roshan, Bioinformatics. 19Dhillon, I.S.; Marcotte, E.M.; Roshan, U. Diametrical clustering for identifying anti-correlated gene clusters. Bioinformatics 2003, 19, 1612-1619.

An antipodally symmetric distribution on the sphere. C Bingham, Ann. Stat. 2Bingham, C. An antipodally symmetric distribution on the sphere. Ann. Stat. 1974, 2, 1201-1225.

Special Functions. G Andrews, R Askey, R Roy, Cambridge University PressCambridge, UKAndrews, G.; Askey, R.; Roy, R. Special Functions; Cambridge University Press: Cambridge, UK, 1999.

The NIST Handbook of Mathematical Functions. A B O Daalhuis, F W Olver, D W Lozier, R F Boisvert, C W Clark, Cambridge University PressCambridge, UKDaalhuis, A.B.O. The NIST Handbook of Mathematical Functions.; Olver, F.W., Lozier, D.W., Boisvert, R.F., Clark, C.W., Eds.; Cambridge University Press: Cambridge, UK, 2010; pp. 321-349.

. K V Mardia, P E Jupp, Directional Statistics. Mardia, K.V.; Jupp, P.E. Directional Statistics;

. John Wiley, Sons, Hoboken, NJ, USAJohn Wiley and Sons: Hoboken, NJ, USA, 2000.

Variational inference for Watson mixture model. J Taghia, A Leijon, IEEE Trans. Pattern Recognit. Mach. Intell. in pressTaghia, J.; Leijon, A. Variational inference for Watson mixture model. IEEE Trans. Pattern Recognit. Mach. Intell. 2014, in press.

Microarrays in cancer: Research and applications. M F Ochs, A K Godwin, Biotechniques. 34Suppl. POchs, M.F.; Godwin, A.K. Microarrays in cancer: Research and applications. Biotechniques 2003, 34 (Suppl. P), S4-S15.

A study of the influence of sex on genome wide methylation. J Liu, M Morgan, K Hutchison, V D Calhoun, PLoS One. 5Liu, J.; Morgan, M.; Hutchison, K.; Calhoun, V.D. A study of the influence of sex on genome wide methylation. PLoS One 2010, 5, e10028.

LinkNMF: Identification of histone modification modules in the human genome using nonnegative matrix factorization. I Jung, D Kim, Gene. 518Jung, I.; Kim, D. LinkNMF: Identification of histone modification modules in the human genome using nonnegative matrix factorization. Gene 2013, 518, 215-221.

Genome-wide DNA methylation analysis for diabetic nephropathy in type 1 diabetes mellitus. C Bell, A E Teschendorff, V K Rakyan, A P Maxwell, S Beck, D A Savage, 10.1186/1755-8794-3-33Med. Genomics. 3Bell, C.; Teschendorff, A.E.; Rakyan, V.K.; Maxwell, A.P.; Beck, S.; Savage, D.A. Genome-wide DNA methylation analysis for diabetic nephropathy in type 1 diabetes mellitus. Med. Genomics 2010, 3, doi:10.1186/1755-8794-3-33.

Variational bayesian matrix factorization for bounded support data. Z Ma, A E Teschendorff, A Leijon, IEEE Trans. Pattern Anal. Mach. Intell. in pressMa, Z.; Teschendorff, A.E.; Leijon, A. Variational bayesian matrix factorization for bounded support data. IEEE Trans. Pattern Anal. Mach. Intell. 2014; in press.

A variational bayesian framework for graphical models. H Attias, In Advances in Neural Information Processing Systems. 12MIT PressAttias, H. A variational bayesian framework for graphical models. In Advances in Neural Information Processing Systems 12; MIT Press: Cambridge, MA, USA, 2000; pp. 209-215.

A tutorial on spectral clustering. U Luxburg, Stat. Comput. 17Luxburg, U. A tutorial on spectral clustering. Stat. Comput. 2007, 17, 395-416.

On spectral clustering: Analysis and an algorithm. A Y Ng, M I Jordan, Y Weiss, Advances in Neural Information Processing Systems. Cambridge, MA, USAMIT PressNg, A.Y.; Jordan, M.I.; Weiss, Y. On spectral clustering: Analysis and an algorithm. In Advances in Neural Information Processing Systems; MIT Press: Cambridge, MA, USA, 2001; pp. 849-856.

A short note on parameter approximation for von Mises-Fisher distributions and a fast implementation of I s (x). S Sra, Comput. Stat. 27Sra, S. A short note on parameter approximation for von Mises-Fisher distributions and a fast implementation of I s (x). Comput. Stat. 2012, 27, 177-190.