# Fast dimension reduction and integrative clustering of multi-omics data using low-rank approximation: application to cancer molecular classification

- Dingming Wu
^{1}, - Dongfang Wang
^{1}, - Michael Q. Zhang
^{1, 2}Email author and - Jin Gu
^{1}Email author

**Received: **22 May 2015

**Accepted: **16 November 2015

**Published: **1 December 2015

## Abstract

### Background

One major goal of large-scale cancer omics study is to identify molecular subtypes for more accurate cancer diagnoses and treatments. To deal with high-dimensional cancer multi-omics data, a promising strategy is to find an effective low-dimensional subspace of the original data and then cluster cancer samples in the reduced subspace. However, due to data-type diversity and big data volume, few methods can integrative and efficiently find the principal low-dimensional manifold of the high-dimensional cancer multi-omics data.

### Results

In this study, we proposed a novel low-rank approximation based integrative probabilistic model to fast find the shared principal subspace across multiple data types: the convexity of the low-rank regularized likelihood function of the probabilistic model ensures efficient and stable model fitting. Candidate molecular subtypes can be identified by unsupervised clustering hundreds of cancer samples in the reduced low-dimensional subspace. On testing datasets, our method LRAcluster (low-rank approximation based multi-omics data clustering) runs much faster with better clustering performances than the existing method. Then, we applied LRAcluster on large-scale cancer multi-omics data from TCGA. The pan-cancer analysis results show that the cancers of different tissue origins are generally grouped as independent clusters, except squamous-like carcinomas. While the single cancer type analysis suggests that the omics data have different subtyping abilities for different cancer types.

### Conclusions

LRAcluster is a very useful method for fast dimension reduction and unsupervised clustering of large-scale multi-omics data. LRAcluster is implemented in R and freely available via http://bioinfo.au.tsinghua.edu.cn/software/lracluster/.

### Keywords

Mutli-omics Cancer Low-rank approximation Clustering Dimension reduction Algorithm## Background

Cancer is a large family of lethal diseases which are killing millions of lives each year [1, 2]. Highly genetic heterogeneity makes it hard to develop general and effective treatments against cancer [3, 4]. One of the major goal of cancer multi-omics study is to discover possible cancer subtypes using molecule-level signatures, which can be used for more accurate diagnoses and treatments [5–8]. Several international collaborated projects, such as TCGA [9], ICGC [10], and CCLE [11] generated tons of cancer multi-omics data. However, we still face several challenges for analyzing such large-scale cancer multi-omics data: 1) need to handle different data types of different platforms at the same time, such as count based data of sequencing, continuous data of microarray and binary data of genetic variations; 2) the data dimension (the number of the molecular features) is much higher than the sample number; and 3) the big data volumes require efficient and robust computational algorithms.

The molecules involved in the same biological processes are usually highly correlated. It is commonly believed that the high-dimensional cancer genomic data can be reduced to a low-dimensional subspace associated to a few major biological processes [12–15], such as sustainable proliferation, apoptosis resistance, activated invasion and immune avoidance [16, 17]. Several efforts have been made to do such integration analysis [18–22]. To find the shared low-dimensional subspace across multiple data types, Shen et al. proposed a latent model iCluster + based on probabilistic principal component analysis, which used generalized linear models to transform continuous, discretized and count variables as a sparse linear regression on a set of latent driving factors. Then, cancer subtyping can be done in the reduced subspace consisting of the latent driving factors [21, 22]. Lock et al. proposed another Bayesian latent model (Bayesian consensus clustering, BCC) to simultaneously find the latent low-dimension subspaces and assign samples into different clusters [23]. However, the low computational efficiency limits its applications on large-scale cancer omics dataset.

In recent years, low-rank approximation (LRA) is becoming one kind of promising dimension reduction methods [20, 24]. In most cases, LRA is convex and can be solved using fast algorithm [25–27]. A few studies show the advantages of LRA for single data type analysis, such as cancer copy number variations [20, 28]. In this study, we formulated a novel low-rank approximation based integrative probabilistic model, which can deal with different data types with high computational efficiency and stability. It assumes that a few major biological factors determine a set of high-dimensional but low-rank systems parameters and the observed cancer omics data are generated based on these parameters. Results show that our method LRAcluster can run much faster than iCluster + with stable model fitting, which makes it possible to analyze large-scale cancer multi-omics data on a small server or even a personal computer.

Then, LRAcluster is applied on a large-scale TCGA multi-omics dataset of 11 different cancer types with four different data types, which is hard to be processed by previous methods. The pan-cancer analysis results suggest that different cancer types (or different tissue origins) can be generally grouped into independent clusters except squamous-like carcinomas in the reduced low-dimensional subspace. While, the single cancer type analysis results show that the multi-omics data have different subtyping capabilities for different cancer types.

## Methods

### LRAcluster overview

### Probabilistic model

*k*-th type of omics data are denoted as

*X*

_{ ij }

^{(k)}(the row index represents the

*i*-th molecular feature and the column index represent the

*j*-th sample), while

*Θ*

^{(k)}denotes the size-matched parameter matrix of

*X*

^{(k)}. The probabilistic model specifies the probability density (mass) function of the observations given the parameters for each data type as below:

- a)
\( \Pr \left({X}_{ij}^{(k)}\Big|{\varTheta}_{ij}^{(k)}\right)\propto \exp \left(-\frac{1}{2}{\left({X}_{ij}^{(k)}-{\varTheta}_{ij}^{(k)}\right)}^2\right) \) for real-type data, Gaussian distribution (CNV and DNA methylation data in this study);

- b)
\( \Pr \left({X}_{ij}^{(k)}\Big|{\varTheta}_{ij}^{(k)}\right)=\frac{e^{\varTheta_{ij}^{(k)}}}{1+{e}^{\varTheta_{ij}^{(k)}}}I\left({X}_{ij}^{(k)}=1\right)+\frac{1}{1+{e}^{\varTheta_{ij}^{(k)}}}I\left({X}_{ij}^{(k)}=0\right) \) for binary data, Bernoulli distribution (somatic mutation data in this study);

- c)
\( \Pr \left({X}_{ij}^{(k)}\Big|{\varTheta}_{ij}^{(k)}\right)\propto {\left({\lambda}_{ij}^{(k)}\right)}^{X_{ij}^{(k)}}{e}^{\left(-{\lambda}_{ij}^{(k)}\right)},\kern0.5em {\lambda}_{ij}^{(k)}={e}^{\varTheta_{ij}^{(k)}} \) for count data, Poisson distribution (RNAseq normalized count data in this study).

Categorical data can be transformed using dummy code and thus can be treated as binary variables.

*X*

^{(k)}(

*k*= 1, 2, …,

*K*). Thus the overall parameter matrix

*Θ*stacks all the parameter matrices (

*Θ*

^{(k)}) used for each observed data matrix. The overall likelihood function is the sum of the likelihood functions of different data types:

*X*

_{ ij }are independently distributed conditional on the ultrahigh dimensional parameter matrix

*Θ*. The prior assumption of the model is that

*Θ*has low-rank structure. The low-rank assumption is used to penalize the freedom of the model and eventually leads to the following optimization problem:

where *μ* is a tuning parameter and |•|* denotes the nuclear norm of the matrix [25].

### Fast low-rank approximation

*f*is a convex function. The iterative solution framework can be briefly expressed as the following steps:

- 1)
initialize

*Θ*^{0}and iterate the following two steps until convergence - 2)$$ {\varTheta}^{2n+1}={\varTheta}^{2n}-{\delta}_n\nabla f $$
- 3)$$ {\varTheta}^{2n+2}={\mathrm{D}}_{\mu}\left({\varTheta}^{2n+1}\right) $$

∇*f* is the gradient of the un-regularized likelihood function (2) and *δ*
_{n} is the step length. D_{
μ
} represents the “singular value shrinkage operator”: let us denote the singular value decomposition (SVD) of a matrix *Θ* as *Θ* = *UΣV*
^{
T
}, then D_{
μ
}(*Θ*) = *U*D_{
μ
}(*Σ*)*V*
^{T}. D_{
μ
}(*Σ*) is a diagonal matrix with the same size as *Σ* and each diagonal element is the shrinkage of the singular value of *Σ*. For a positive singular value *λ*, the shrinkage result is (*λ*–*μ*) when *λ* > *μ* and 0 when *λ* ≤ *μ*.

The objective function of LRAcluster is convex, so any initial value of the iteration will converge to the global minimum. LRAcluster simply initializes *Θ* as a zero matrix. The original framework needs a user defined constraint parameter *μ* which is hard to choose in practical use. Instead of *μ*, LRAcluster receives the rank *r* (also the target dimension) as the user defined constraint parameter. *μ* is automatically chosen as the rank *r* + 1 largest singular value in each iteration. The choice of *μ* is to guarantee that *Θ* has rank *r* and the shrinkage has minimal effect on *Θ*. For simple “matrix completion problem”, Cai et al. proves that when the step length *δ* is between 0.5 and 2, the algorithm converges definitely [26]. LRAcluster set *δ* as 0.5, which ensures convergence for real applications in this study.

The target rank (or dimension) *r* is the only user-defined parameter in dimension reduction step. The log likelihood − L(*θ*; *X*) corresponding to the optimized solution *θ** (denoted as ℒ
_{
r
}
^{*}
) is used for guiding the choice of parameter *r*: for the same dataset, larger *r* means weaker penalization of the model freedom and leads to better data fitting (larger likelihood ℒ
_{
r
}
^{*}
). Thus, ℒ
_{
r = 0}
^{*}
is the minimum and ℒ
_{
r = + ∞}
^{*}
is the maximum among all the ℒ
_{
r
}
^{*}
. The quantity ℒ
_{
r
}
^{*}
describes to what extend the model fits the data. As LRAcluser mainly deals with large dataset, ℒ
_{
r
}
^{*}
is usually a big value. So, instead of ℒ
_{
r
}
^{*}
, LRAcluster uses the normalized quantity \( \frac{{\mathrm{\mathcal{L}}}_{r=+\infty}^{*}-{\mathrm{\mathcal{L}}}_r^{*}}{{\mathrm{\mathcal{L}}}_{r=+\infty}^{*}-{\mathrm{\mathcal{L}}}_{r=0}^{*}} \) (between 0 and 1) as “explained variation” for choosing a desirable rank *r*. We will describe the basic principles for the choice of rank *r* in Results section.

### Dimension reduction and clustering

The dimension reduction is straightforward after getting the low-rank matrix *Θ*. As the rank of *Θ* is no more than *r*, the singular vector decomposition (SVD) of that matrix *Θ* = *UΣV*
^{T} has *Σ* with no more than *r* non-zero singular values. Thus the first *r* columns of *ΣV*
^{T} are just the dimension reduction result of the original data matrix *X* with the target dimension (rank) *r*.

LRAcluster uses *k*-means to identify the candidate molecular subtypes in the reduced low-dimensional subspace. Silhouette values [29] is used to determine the cluster number *k*. Any other unsupervised clustering algorithm can be used instead of *k*-means.

### Datasets

In this study, all the datasets were downloaded from publicly released TCGA level 3 data (processed data from UCSC Cancer Genome Browser [30]). No ethics approval is required for this study. The whole dataset consists of 11 types of cancer (BRCA, COAD, GBM, HNSC, KIRC, LGG, LUAD, LUSC, PRAD, STAD, and THCA) with somatic mutations, copy number variations, DNA methylations and gene expressions. For somatic mutation and copy number variation data, our preliminary studies indicate that the massive passenger variations of the complete datasets deteriorated the clustering stability. Thus, only the somatic mutations and copy number variations of the ~500 genes reported as “causally implicated in cancer” in COSMIC [31] were included in this study. For DNA methylation data using Illumina HumanMethylation450 BeadChip (450 k array), probes annotated as “promoter_associated” (based on the annotations of IlluminaHumanMethylation450k.db [32]) were selected (if a gene has multiple promoter associated probes, only one of them was chosen). Overall, ~8,000 probes were used. The normalized count-based data from RNA-Seq were all included with ~20,000 genes.

The three cancer-type testing dataset consists of BRCA, COAD, LUAD cancer types with RNA-Seq and and 450 k DNA methylation data. The other datasets consists of all the four data types described as above.

## Results

LRAcluster is a computational-efficient method for fast dimension reduction and integrative clustering of large-scale cancer multi-omics data. We first show the performances and parameter tuning of LRAcluster on a three cancer-type testing dataset and a breast cancer dataset labeled with ER+/ER- subtypes. Then, it was applied on the large-scale TCGA pan-cancer dataset.

### The computational performances of LRAcluster

A three cancer-type dataset was used to compare the clustering performances and time consumption between LRAcluster and iCluster+. The molecular features (genes for expression data and probes for DNA methylation data) with largest variances across all samples are selected to construct datasets of different sizes. The smallest dataset containing top 100 molecular features of each data type is used to test LRAcluster and iCluster+’s clustering performances with different target dimension (from 2 to 10). Time consumption of the two methods was recorded for datasets with different feature sizes (from 100 to 5000 features). iCluster + runs under both normal model (random initialization of penalty parameter for better model) and simple model (fixed penalty parameter).

*μ*demonstrated that LRAcluster can quickly converge within only a few iterations (Fig. 2c & d). There are two important parameters in LRAcluster: the rank (or dimension) of the reduced subspace

*r*and the cluster number

*c*. To illustrate how to choose these parameters empirically, we used the BRCA dataset with known ER+/ER- subtypes as an example: the dimension

*r*can be chosed according to the curve of “explained variance” (Fig. 3a) and the cluster number

*c*can be chosen according to the curve of silhouette value (s-value) (Fig. 3b). For the BRCA dataset, dimension

*r*should be chosen as 2, because there was a turning point at 2 on the curve of the “explained variance” (Fig. 3a). This empirical rule is based on the principle that the increase of model fitness is much slower after the changing point. The choice of cluster number

*c*is straightforward: larger s-value indicates better clustering performance. For the BRCA dataset, the largest s-value was achieved when

*c*= 2 (Fig. 3b). Results show that LRAcluster can find two subtypes in the reduced 2-dimensional subspace and the identified subtypes are highly consistent with known ER+/ER- subtypes (accuracy 92.1 %) (Fig. 3c).

### Application on the large-scale TCGA pan-cancer dataset

The unsupervised clustering results of pan-cancer analysis

BRCA | COAD | GBM | HNSC | KIRC | LGG | LUAD | LUSC | PRAD | STAD | THCA | Total | |
---|---|---|---|---|---|---|---|---|---|---|---|---|

C1 | 1 | 0 | 0 | 286 | 0 | 0 | 0 | 6 | 0 | 0 | 0 | 293 |

C2 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 411 | 412 |

C3 | 0 | 0 | 41 | 0 | 0 | 451 | 0 | 0 | 0 | 0 | 0 | 492 |

C4 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 231 | 0 | 231 |

C5 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 293 | 0 | 0 | 293 |

C6 | 0 | 190 | 0 | 1 | 0 | 0 | 2 | 0 | 1 | 0 | 0 | 194 |

C7 | 3 | 17 | 0 | 0 | 1 | 0 | 406 | 7 | 0 | 0 | 3 | 437 |

C8 | 0 | 0 | 0 | 0 | 240 | 0 | 0 | 0 | 0 | 0 | 0 | 240 |

C9 | 448 | 0 | 1 | 2 | 1 | 0 | 4 | 1 | 0 | 0 | 0 | 457 |

C10 | 8 | 1 | 0 | 195 | 0 | 0 | 6 | 60 | 0 | 0 | 0 | 270 |

Total | 460 | 208 | 42 | 484 | 242 | 452 | 418 | 74 | 294 | 231 | 414 | 3319 |

Results show that most samples from the same cancer types are grouped as independent clusters. These results are similar with a recent pan-cancer study [8]. The two brain cancers (LGG and GBM) are grouped together (Cluster C3). Only HNSC are separated into two major clusters (Cluster C1 & C10) and the samples (40.3 % of HNSC) in Cluster C10 are clustered together with LUSC samples (81.1 % of LUSC), which indicates that the squamous carcinomas of different tissue origins may share some common molecular mechanisms. A recent work also reported an integrative network-based stratification (jNBS) pan-cancer clustering analysis on TCGA dataset, which incorporated multi-omics data with the information of a pre-given gene network [33]. Generally speaking, it reported similar results with LRAcluster: most of cancer types are separately clustered according to their tissue origin, and two types of squamous carcinomas, head/neck squamous carcinoma and lung squamous carcinoma are cluster together. But it found more cross tissue type clusters. Because the jNBS analysis only used genetic (mutation & CNV) and epigenetic (DNA methylation) data, the results are hard to be directly compared. The molecular signatures associated with the pan-cancer clusters were shown in Additional file 1: Figure S3.

The results of single-cancer analysis

Cancer | Dimension | #Cluster | Silhouette values |
---|---|---|---|

BRCA | 2 | 2 | 0.55 |

COAD | 4 | 4 | 0.40 |

GBM | 8 | 2 | 0.35 |

HNSC | 7 | 3 | 0.26 |

KIRC | 6 | 2 | 0.36 |

LGG | 2 | 3 | 0.44 |

LUAD | 5 | 2 | 0.34 |

LUSC | 5 | 4 | 0.32 |

PRAD | 2 | 4 | 0.41 |

STAD | 4 | 3 | 0.37 |

THCA | 2 | 2 | 0.61 |

## Conclusion

LRAcluster probabilistically models the observed data conditional on the size-matched parameters. The low-rank constraint is the key to get the low-dimensional representation of the original data. And the convexity of the regularized likelihood function provides efficient gradient-descent algorithm for model fitting. Results show that LRAcluster runs fast with high classification accuracy and it is suitable for large-scale cancer multi-omics analysis.

## Discussions

In LRAcluster probabilistic model, the real-type data are modeled as Gaussian-distributed random variables with variance 1. Though the assumption of all features having the same variance seems unnatural for any dataset as the different features should have different variance, it is consistent with the simple method of principle component analysis. Minus log likelihood function of the real-type data is \( \frac{1}{2}{\left({X}_{ij}-{\varTheta}_{ij}\right)}^2 \) which is the same as the loss function of principle component analysis (PCA). So, if there are only real-type data as input, the LRAcluster solution is in principle the same as the PCA. The only difference is the scale of each principle component because the LRAcluster considers the L_{1} norm but PCA considers the L_{0} norm.

LRAcluster receives the rank *r* of the matrix *Θ* as the user-defined parameter instead of the original parameter *μ*. This setting makes the dimension reduction more straightforward: *r* is just the target dimension of the reduced subspace. From computational view, *μ* and *r* have the same function as they are both used to penalize the complexity of the probabilistic model.

LRAcluster does not penalize the association between molecular features and the reduced subspace (latent factors) via sparsity assumption. It is a better strategy to find the molecular features associated the identified clusters or subtypes by molecular signature analysis: find the significantly differential features between the samples in that cluster and all the other samples (please see the heatmap of the selected molecular features of TCGA pan-cancer analysis in Additional file 1: Figure S3). Besides, LRAcluster will prefer the inter-omics features with large co-variances implied by the low-rank assumption (for example, the significantly correlated CNVs and mRNA expressions). The inter-omics regulatory information can be modeled as a separate pre-processing step to find the cancer driving factors and then only the molecular features significantly associated with these drivers are used as the input of LRAcluster.

Joint non-negative matrix factorization (jNMF) is another strategy to find the shared principal subspace across multiple omics datasets [34, 35]. Theoretically, NMF can be treated as a matrix version of latent factor analysis. jNMF will also encounter the optimization difficulty of non-convey loss function. But the advantage of jNMF is that the model can also get the molecular features (or called as modules) significantly associated each dimension.

## Declarations

### Acknowledgements

We thank Songpeng Zu, Zijian Ding and Qiuyu Lian for their kind discussions and method testing. This work is supported by National Basic Research Program of China [2012CB316503], National Natural Science Foundation of China [61370035 and 31361163004] and Tsinghua University Initiative Scientific Research Program.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359–386.View ArticlePubMedGoogle Scholar
- Siegel R, Ma J, Zou Z, Jemal A. Cancer statistics, 2014. CA Cancer J Clin. 2014;64(1):9–29.View ArticlePubMedGoogle Scholar
- Bedard PL, Hansen AR, Ratain MJ, Siu LL. Tumour heterogeneity in the clinic. Nature. 2013;501:355–64.View ArticlePubMedGoogle Scholar
- Burrell RA, McGranahan N, Bartek J, Swanton C. The causes and consequences of genetic heterogeneity in cancer evolution. Nature. 2013;501:338–45.View ArticlePubMedGoogle Scholar
- Hayhoe FG. Classification of acute leukaemias. Blood Rev. 1988;2:186–93.View ArticlePubMedGoogle Scholar
- Yan H, Peng Z-G, Wu Y-L, Jiang Y, Yu Y, Huang Y, et al. Hypoxia-simulating agents and selective stimulation of arsenic trioxide-induced growth arrest and cell differentiation in acute promyelocytic leukemic cells. Haematologica. 2005;90:1607–16.PubMedGoogle Scholar
- Yersal O, Barutca S. Biological subtypes of breast cancer: Prognostic and therapeutic implications. World J Clin Oncol. 2014;5:412–24.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoadley KA, Yau C, Wolf DM, Cherniack AD, Tamborero D, Ng S, et al. Multiplatform analysis of 12 cancer types reveals molecular classification within and across tissues of origin. Cell. 2014;158:929–44.PubMed CentralView ArticlePubMedGoogle Scholar
- The Cancer Genome Atlas [http://cancergenome.nih.gov/]
- Hudson TJ, Anderson W, Artez A, Barker AD, Bell C, Bernabe RR, et al. International network of cancer genome projects. Nature. 2010;464:993–8.View ArticlePubMedGoogle Scholar
- Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483:603–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang E, Ishida S, Pittman J, Dressman H, Bild A, Kloos M, et al. Gene expression phenotypic models that predict the activity of oncogenic pathways. Nat Genet. 2003;34:226–30.View ArticlePubMedGoogle Scholar
- Li L, Li H. Dimension reduction methods for microarrays with application to censored survival data. Bioinforma Oxf Engl. 2004;20:3406–12.View ArticleGoogle Scholar
- Li H, Gui J. Partial Cox regression analysis for high-dimensional microarray gene expression data. Bioinforma Oxf Engl. 2004;20 Suppl 1:i208–215.View ArticleGoogle Scholar
- Jia P, Pao W, Zhao Z. Patterns and processes of somatic mutations in nine major cancers. BMC Med Genomics. 2014;7:11.PubMed CentralView ArticlePubMedGoogle Scholar
- Hanahan D, Weinberg RA. The hallmarks of cancer. Cell. 2000;100(1):57–70.View ArticlePubMedGoogle Scholar
- Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646–74.View ArticlePubMedGoogle Scholar
- Alter O, Golub GH. Integrative analysis of genome-scale data by using pseudoinverse projection predicts novel correlation between DNA replication and RNA transcription. Proc Natl Acad Sci U S A. 2004;101:16577–82.PubMed CentralView ArticlePubMedGoogle Scholar
- Yuan Y, Savage RS, Markowetz F. Patient-specific data fusion defines prognostic cancer subtypes. PLoS Comput Biol. 2011;7:e1002227.PubMed CentralView ArticlePubMedGoogle Scholar
- Lock EF, Hoadley KA, Marron JS, Nobel AB. Joint and individual variation explained (JIVE) for integrated analysis of multiple data types. Ann Appl Stat. 2013;7(1):523–42.PubMed CentralView ArticlePubMedGoogle Scholar
- Mo Q, Wang S, Seshan VE, Olshen AB, Schultz N, Sander C, et al. Pattern discovery and cancer gene identification in integrated cancer genomic data. Proc Natl Acad Sci U S A. 2013;110:4245–50.PubMed CentralView ArticlePubMedGoogle Scholar
- Shen R, Olshen AB, Ladanyi M. Integrative clustering of multiple genomic data types using a joint latent variable model with application to breast and lung cancer subtype analysis. Bioinforma Oxf Engl. 2009;25:2906–12.View ArticleGoogle Scholar
- Lock EF, Dunson DB. Bayesian consensus clustering. Bioinforma Oxf Engl. 2013;29:2610–6.View ArticleGoogle Scholar
- Yuan M, Ekici A, Lu Z, Monteiro R. Dimension reduction and coefficient estimation in multivariate linear regression. J R Stat Soc Ser B-Stat Methodol. 2007;69:329–46.View ArticleGoogle Scholar
- Candes EJ, Recht B. Exact Matrix Completion via Convex Optimization. Found Comput Math. 2009;9:717–72.View ArticleGoogle Scholar
- Cai JF, Candès EJ, Shen Z. A singular value thresholding algorithm for matrix completion. SIAM J Optim. 2010;20:1956–82.View ArticleGoogle Scholar
- Hsieh CJ, Olsen PA. Nuclear Norm Minimization via Active Subspace Selection. Proc 31st Int Conf Mach Learn. 2014.Google Scholar
- Zhou X, Liu J, Wan X, Yu W. Piecewise-constant and low-rank approximation for identification of recurrent copy number variations. Bioinforma Oxf Engl. 2014;30:1943–9.View ArticleGoogle Scholar
- Rousseeuw P. silhouettes - A graphical aid to the integration of cluster-analysis. J Comput Appl Math. 1987;20:53–65.View ArticleGoogle Scholar
- Goldman M, Craft B, Swatloski T, Cline M, Morozova O, Diekhans M, et al. The UCSC Cancer Genomics Browser: update 2015. Nucleic Acids Res. 2015;43(Database issue):D812–817.PubMed CentralView ArticlePubMedGoogle Scholar
- Forbes SA, Beare D, Gunasekaran P, Leung K, Bindal N, Boutselakis H, et al. COSMIC: exploring the world's knowledge of somatic mutations in human cancer. Nucleic Acids Res. 2015;43(Database issue):D805–811.PubMed CentralView ArticlePubMedGoogle Scholar
- Triche T, Jr. IlluminaHumanMethylation450k.db: Illumina Human Methylation 450k annotation data.Google Scholar
- Liu Z, Zhang S. Tumor characterization and stratification by integrated molecular profiles reveals essential pan-cancer features. BMC Genomics. 2015;16:503.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang S, Li Q, Liu J, Zhou XJ. A novel computational framework for simultaneous integration of multiple types of genomic data to identify microRNA-gene regulatory modules. Bioinforma Oxf Engl. 2011;27:i401–409.View ArticleGoogle Scholar
- Zhang S, Liu C-C, Li W, Shen H, Laird PW, Zhou XJ. Discovery of multi-dimensional modules by integrative analysis of cancer genomic data. Nucleic Acids Res. 2012;40:9379–91.PubMed CentralView ArticlePubMedGoogle Scholar