clustering multidimensional data - computer science- uc...
Post on 23-Apr-2019
227 Views
Preview:
TRANSCRIPT
ECS 234
What is Clustering?
Given n objects, assign them to groups (clusters) based on their similarity
• Unsupervised Machine Learning• Class Discovery• Difficult, and maybe ill-posed problem!
ECS 234
Impossibility of Clustering
• Scale-invariance: meters vs inches• Richness: all partitions as possible solutions• Consistency: increasing distances between
clusters and decreasing distances within clusters should yield the same solution
No function exists that satisfies all three.No function exists that satisfies all three.Kleinberg, NIPS 2002
ECS 234
Clustering Microarray DataClustering reveals similar expression patterns, in particular in time-series expression data
Guilt-by-association: a gene of unknown function has the same function as a similarly expressed gene of known function
Genes of similar expression might be similarly regulated
ECS 234
Clustering ApproachesNon-Parametric Parametric
GenerativeReconstructive(Hierarchical)
Agglomerative Divisive
Gaussian MixtureModels
Fuzzy C-Means
K-Means
K-Medoids (PAM)Single Link
Average Link
Complete Link
Ward Method
Divisive Set Partitioning
SOM
Graph Models
CorruptedClique
Bayesian Models
Hard ClusteringSoft Clustering
Multi-feature
BiclusteringPlaid Models
ECS 234
How To Choose the Right Clustering?
• Data Type– Independent Experiments (e.g. knockouts)– Dependent experiments (e.g. time series)
• Parametric vs. non-parametric clustering• Quality of Clustering• Software Availability• Features of the Methods
– Computing averages (sometimes impossible or too slow)– Stability analysis– Properties of the clusters– Speed– Memory
ECS 234
Clustering Meta-Procedure
1. Compare the similarity of all pairs of objects
2. Group the most similar ones together into clusters
3. Reason about the resulting groups of clusters
ECS 234
Certain properties are expected from distance measures
1. d(x,y)=0– d(x,y)>0, x≠y3. d(x,y)=d(y,x)– d(x,y)≤d(x,z)+d(z,y) the triangle inequality
If properties 1-4 are satisfied, the distance measure is a metric
Distance Measures, d(x,y)
ECS 234
The Lp norm
p = 2, Euclidean Dist.p = ∞, Manhattan Dist.(downtown Davis distance)
p pnn
p yxyxyxd ||||),( 11 −++−=
p=1 p=2 p=3 p=4 p=20
Equidistant points from a center, for different norms
ECS 234
Pearson Correlation Coefficient
))(
)()(
(
),(2
2
2
2
n
yy
n
xx
n
yxyx
yxr
kk
kk
kk
kk
k kkk
kkk
∑∑
∑∑
∑ ∑∑
−−
−=
(Normalized vector dot product)
Good for comparing expression profiles because it is insensitive to scaling (but data should be normally distributed, e.g. log expression)!
Not a metric!
ECS 234
Hierarchical Clustering• Input: Data Points, x1,x2,…,xn
• Output:Tree – the data points are leaves– Branching points indicate similarity between sub-trees– Horizontal cut in the tree produces data clusters
1 2
5
3 7
4
61 2 57 4 6 C
lust
er M
ergi
ng C
ost
3
ECS 234
1 2 57 4 6 Clu
ster
Mer
ging
Cos
t
3
Maximum iterations: n-1
General Algorithm• Place each element in its own cluster, Ci={xi}• Compute (update) the merging cost between every pair
of elements in the set of clusters to find the two cheapest to merge clusters Ci, Cj,
• Merge Ci and Cj in a new cluster Cij which will be the parent of Ci and Cj in the result tree.
4. Go to (2) until there is only one set remaining
ECS 234
Different Types of Algorithms Based on The Merging Cost
• Single Link,
• Average Link,
• Complete Link,
• Others (Ward method-least squares)
),(min,
yxdji CyCx ∈∈
),(max,
yxdji CyCx ∈∈
),(||||
1 yxdCjCi
i jCx Cy∑ ∑
∈ ∈
ECS 234
Characteristics of Hierarchical Clustering
• Greedy Algorithms – suffer from local optima, and build a few big clusters
• A lot of guesswork involved:– Number of clusters– Cutoff coefficient– Size of clusters
• Average Link is fast and not too bad: biologically meaningful clusters are retrieved
ECS 234
K-MeansInput: Data Points, Number of Clusters (K)Output: K clustersAlgorithm: Starting from k-centroids assign data points to them based on proximity, updating the centroids iteratively
• Select K initial cluster centroids, c1, c2, c3, ..., ck
• Assign each element x to nearest centroid7. For each cluster, re-compute its centroid by averaging the
data points in it 8. Go to (2) until convergence is achieved
ECS 234
K-Means Properties
• Must know the number of clusters before hand
• Sensitive to perturbations• Clusters formed ad hoc with no indication
of relationships among them• Results depend on initial choice for centers• In general, betters average link clustering
ECS 234
Self Organizing Maps ClusteringInput: Data Points, SOM Topology (K nodes and a distance function)Output: K clusters, (near clusters are similar) Algorithm: Starting with a simple topology (connected nodes) iteratively move the nodes “closer” to the data
1. Select initial topology 2. Select a random data point P3. Move all the nodes towards P by varying amounts4. Go to (2) until convergence is achieved.
ECS 234
))()(),,(()()(1 NfPiNNdNfNf ipii −+=+ τ
iiteration at N node ofPosition )(
N and Nbetween Distance),(
P closest to Node P;point Random Node;
p
=
=
===
NfNNd
NPN
i
p
p
τ is the learning rate (decreases with d and I)
Initial Topology
ECS 234
SOM Properties
• Neighbouring clusters are similar• Element on the borders belong to both
clusters• Very robust• Works for short profile data too
ECS 234
What if the number of clusters is not known? • Elbow criterion: look for a clustering that explains most of
the variance or stability in data with the fewest clusters• Information theoretic: maximize (or minimize) some
Information Criterion (like BIC or AIC or MDL)• Within/between cluster distance/separation: silhouettes
ECS 234
Note on Missing Values• Microarray experiments often have missing values,
as a result of experimental error, values out of bound, spot reading error, batch errors, etc.
• Many clustering algorithms (all of the ones presented here) are sensitive to missing data
• Filling in the holes:– All 0s– Average– Better: weighted K-nearest neighbor, or SVD based
methods (SVDimpute, KNNimpute) Troyanskaya et al. 2000 (AVAILABLE FOR DOWNLOAD)
• Robust• Do better than average
ECS 234
Cluster Visualization• How to “see” the clusters effectively?• Present gene expressions in different colors• Plot similar genes close to each other• R• GeneXPress• Expander• CytoScape
ECS 234
Algorithm Comparison and Cluster Validation
• Paper: Chen et al. 2001
• Data: embryonic stem cells expression data• Results: evaluated advantages and
weaknesses of algorithms w/respect to both internal and external quality measures
• Used known and developed novel indices to measure clustering efficacy
ECS 234
Algorithms Compared
• Average Link Hierarchical Clustering, • K-Means and PAM , and • SOM, two different neighborhood radii
– R=0 (theoretically approaches K-Means)– R=1
• Compared them for different numbers of clusters
ECS 234
Clustering Quality Indices• Homogeneity and Separation
– Homogeneity is calculated as the average distance between each gene expression profile and the center of the cluster it belongs to
– Separation is calculated as the weighted average distance between cluster centers
– H reflects the compactness of the clusters while S reflects the overall distance between clusters
– Decreasing H or increasing S suggest an improvement in the clustering results
ECS 234
Results:
•K-Means and PAM scored identically
•SOM_r0 very close to both above
•All three beat ALHC
•SOM_r1 worst
ECS 234
• Silhouette Width– A composite index reflecting the compactness and
separation of the clusters, and can be applied to different distance metrics
– A larger value indicates a better overall quality of the clustersResults:•All had low scores indicating underlying “blurriness” of the data
•K-Means, PAM, SOM_r0 very close
•All three slightly better than ALHC
•SOM_r1 had the lowest score
ECS 234
• Redundant Scores (external validation)– Almost every microarray data set has a small portion of duplicates,
i.e. redundant genes (check genes)– A good clustering algorithm should cluster the redundant genes’
expressions in the same clusters with high probability– DRRS (difference of redundant separation scores) between control
and redundant genes was used as a measure of cluster quality– High DRRS suggests the redundant genes are more likely to be
clustered together than randomly chosen genes
Results:
- K-means consistently better than ALHC
- PAM and SOM_r0 close to the above
- SOM_r1 was consistently the worst
ECS 234
• WADP – Measure of Robustness– If the input data deviate slightly from their current
value, will we get the same clustering?– Important in Microarray expression data analysis
because of constant noise– Experiment:
• each gene expression profile was perturbed by adding to it a random vector of the same dimension
• values for the random vector generated from a Gaussian distr. (mean zero, and stand. dev.=0.01)
• data was renormalized and clustered• WADP Cluster discrepancy: measure of inconsistent
clusterings after noise. WADP=0 is perfect.
ECS 234
Results:
•SOM_r1 clusters are the most robust of all
•K-means and ALHC were high through all cluster numbers
•PAM and SOM_r1 were better for small number of clusters
ECS 234
Comparison of Cluster Content• How similar are two clusterings in all the methods?
– WADP
• Other measures of similarity based on co-clusteredness of elements– Rand index– Adjusted Rand– Jaccard
ECS 234
Clustering: Conclusions• K-means outperforms ALHC• SOM_r0 is almost K-means and PAM• Tradeoff between robustness and cluster quality:
SOM_r1 vs SOM_r0, based on the topological neighborhood
• Whan should we use which? Depends on what we know about the data– Hierarchical data – ALHC– Cannot compute mean – PAM– General quantitative data - K-Means– Need for robustness – SOM_r1– Soft clustering: Fuzzy C-Means– Clustering genes and experiments - Biclustering
ECS 234
Biclustering• Problem with clustering:
– Clustering the same genes under different subsets of conditions can result in very different clusterings
• Additional Motivation– sometimes only subset of genes are interesting and one
wants to cluster those– Genes expressed differentially in different conditions
and pathways• Proposed solutions: cluster simultaneously the
genes and the conditions
ECS 234
Clustering conditions
Clustering Genes Biclustering
The biclustering methods look for submatrices in the expression matrix which show coordinated differential expression of subsets of genes in subsets of conditions. The biclusters are also statistically significant.
Clustering is a global similarity method, while biclustering is a local one.
ECS 234
Biclustering Methods
• Biclustering• Coupled Two-way Clustering• Iterative Signature Algorithm• SAMBA• Spectral Biclustering• Plaid Models
ECS 234
Other Dimension Reduction Techniques
• All (including clustering) are based on the premise that not all genes (or experiments) show different behavior, so groups of similar genes (experiments) are sought
• Principal Component Analysis– Identifies the underlying classes or “base” genes of the data
representing most variability (best separating the genes)– All other genes expressions are linear combination of those– Classes are built around a few top “base” genes– Typically used for 2D or 3D data visualization and seeding k-means
• Independent Component Analysis– Similar as PCA but here the “base” components are required to be
statistically independent• Non-zero Matrix Factorization
top related