Title: Cluster Analysis Part I
1Cluster AnalysisPart I
2Learning Objectives
- What is Cluster Analysis?
- Types of Data in Cluster Analysis
- A Categorization of Major Clustering Methods
- Partitioning Methods
3What is Cluster Analysis?
- Cluster a collection of data objects
- Similar to one another within the same cluster
- Dissimilar to the objects in other clusters
- Cluster analysis
- Grouping a set of data objects into clusters
- Clustering is unsupervised classification no
predefined classes - Typical applications
- As a stand-alone tool to get insight into data
distribution - As a preprocessing step for other algorithms
4General Applications of Clustering
- Pattern Recognition
- Spatial Data Analysis
- create thematic maps in GIS by clustering feature
spaces - detect spatial clusters and explain them in
spatial data mining - Image Processing
- Economic Science (especially market research)
- WWW
- Document classification
- Cluster Weblog data to discover groups of similar
access patterns
5Examples of Clustering Applications
- Marketing Help marketers discover distinct
groups in their customer bases, and then use this
knowledge to develop targeted marketing programs - Land use Identification of areas of similar land
use in an earth observation database - Insurance Identifying groups of motor insurance
policy holders with a high average claim cost - City-planning Identifying groups of houses
according to their house type, value, and
geographical location - Earth-quake studies Observed earth quake
epicenters should be clustered along continent
faults
6What Is Good Clustering?
- A good clustering method will produce high
quality clusters with - high intra-class similarity
- low inter-class similarity
- The quality of a clustering result depends on
both the similarity measure used by the method
and its implementation. - The quality of a clustering method is also
measured by its ability to discover some or all
of the hidden patterns.
7Requirements of Clustering in Data Mining
- Scalability
- Ability to deal with different types of
attributes - Discovery of clusters with arbitrary shape
- Minimal requirements for domain knowledge to
determine input parameters - Able to deal with noise and outliers
- Insensitive to order of input records
- High dimensionality
- Incorporation of user-specified constraints
- Interpretability and usability
8- What is Cluster Analysis?
- Types of Data in Cluster Analysis
- A Categorization of Major Clustering Methods
- Partitioning Methods
- Hierarchical Methods
- Density-Based Methods
- Grid-Based Methods
- Model-Based Clustering Methods
- Outlier Analysis
- Summary
9Data Structures
- Data matrix
- (two modes)
- Dissimilarity matrix
- (one mode)
10Measure the Quality of Clustering
- Dissimilarity/Similarity metric Similarity is
expressed in terms of a distance function, which
is typically metric d(i, j) - There is a separate quality function that
measures the goodness of a cluster. - The definitions of distance functions are usually
very different for interval-scaled, boolean,
categorical, ordinal and ratio variables. - Weights should be associated with different
variables based on applications and data
semantics. - It is hard to define similar enough or good
enough - the answer is typically highly subjective.
11Type of data in clustering analysis
- Interval-scaled variables
- Binary variables
- Nominal, ordinal, and ratio variables
- Variables of mixed types
12Interval-valued variables
- Standardize data
- Calculate the mean absolute deviation
- where
- Calculate the standardized measurement (z-score)
- Using mean absolute deviation is more robust than
using standard deviation
13Similarity and Dissimilarity Between Objects
- Distances are normally used to measure the
similarity or dissimilarity between two data
objects - Some popular ones include Minkowski distance
- where i (xi1, xi2, , xip) and j (xj1, xj2,
, xjp) are two p-dimensional data objects, and q
is a positive integer - If q 1, d is Manhattan distance
14Similarity and Dissimilarity Between Objects
(Cont.)
- If q 2, d is Euclidean distance
- Properties
- d(i,j) ? 0
- d(i,i) 0
- d(i,j) d(j,i)
- d(i,j) ? d(i,k) d(k,j)
- Also one can use weighted distance, parametric
Pearson product moment correlation, or other
disimilarity measures.
15Binary Variables
- A contingency table for binary data
- Simple matching coefficient (invariant, if the
binary variable is symmetric) - Jaccard coefficient (noninvariant if the binary
variable is asymmetric)
Object j
Object i
16Dissimilarity between Binary Variables
- Example
- gender is a symmetric attribute
- the remaining attributes are asymmetric binary
- let the values Y and P be set to 1, and the value
N be set to 0
17Nominal Variables
- A generalization of the binary variable in that
it can take more than 2 states, e.g., red,
yellow, blue, green - Method 1 Simple matching
- m of matches, p total of variables
- Method 2 use a large number of binary variables
- creating a new binary variable for each of the M
nominal states
18Ordinal Variables
- An ordinal variable can be discrete or continuous
- order is important, e.g., rank
- Can be treated like interval-scaled
- replacing xif by their rank
- map the range of each variable onto 0, 1 by
replacing i-th object in the f-th variable by - compute the dissimilarity using methods for
interval-scaled variables
19Ratio-Scaled Variables
- Ratio-scaled variable a positive measurement on
a nonlinear scale, approximately at exponential
scale, such as AeBt or Ae-Bt - Methods
- treat them like interval-scaled variables not a
good choice! (why?) - apply logarithmic transformation
- yif log(xif)
- treat them as continuous ordinal data treat their
rank as interval-scaled.
20Variables of Mixed Types
- A database may contain all the six types of
variables - symmetric binary, asymmetric binary, nominal,
ordinal, interval and ratio. - One may use a weighted formula to combine their
effects. - f is binary or nominal
- dij(f) 0 if xif xjf , or dij(f) 1 o.w.
- f is interval-based use the normalized distance
- f is ordinal or ratio-scaled
- compute ranks rif and
- and treat zif as interval-scaled
21- What is Cluster Analysis?
- Types of Data in Cluster Analysis
- A Categorization of Major Clustering Methods
- Partitioning Methods
- Hierarchical Methods
- Density-Based Methods
- Grid-Based Methods
- Model-Based Clustering Methods
- Outlier Analysis
- Summary
22Major Clustering Approaches
- Partitioning algorithms Construct various
partitions and then evaluate them by some
criterion - Hierarchy algorithms Create a hierarchical
decomposition of the set of data (or objects)
using some criterion - Density-based based on connectivity and density
functions - Grid-based based on a multiple-level granularity
structure - Model-based A model is hypothesized for each of
the clusters and the idea is to find the best fit
of that model to each other
23- What is Cluster Analysis?
- Types of Data in Cluster Analysis
- A Categorization of Major Clustering Methods
- Partitioning Methods
- Hierarchical Methods
- Density-Based Methods
- Grid-Based Methods
- Model-Based Clustering Methods
- Outlier Analysis
- Summary
24Partitioning Algorithms Basic Concept
- Partitioning method Construct a partition of a
database D of n objects into a set of k clusters - Given a k, find a partition of k clusters that
optimizes the chosen partitioning criterion - Global optimal exhaustively enumerate all
partitions - Heuristic methods k-means and k-medoids
algorithms - k-means (MacQueen67) Each cluster is
represented by the center of the cluster - k-medoids or PAM (Partition around medoids)
(Kaufman Rousseeuw87) Each cluster is
represented by one of the objects in the cluster
25The K-Means Clustering Method
- Given k, the k-means algorithm is implemented in
4 steps - Partition objects into k nonempty subsets
- Compute seed points as the centroids of the
clusters of the current partition. The centroid
is the center (mean point) of the cluster. - Assign each object to the cluster with the
nearest seed point. - Go back to Step 2, stop when no more new
assignment.
26The K-Means Clustering Method
27Comments on the K-Means Method
- Strength
- Relatively efficient O(tkn), where n is
objects, k is clusters, and t is iterations.
Normally, k, t ltlt n. - Often terminates at a local optimum. The global
optimum may be found using techniques such as
deterministic annealing and genetic algorithms - Weakness
- Applicable only when mean is defined, then what
about categorical data? - Need to specify k, the number of clusters, in
advance - Unable to handle noisy data and outliers
- Not suitable to discover clusters with non-convex
shapes
28Variations of the K-Means Method
- A few variants of the k-means which differ in
- Selection of the initial k means
- Dissimilarity calculations
- Strategies to calculate cluster means
- Handling categorical data k-modes (Huang98)
- Replacing means of clusters with modes
- Using new dissimilarity measures to deal with
categorical objects - Using a frequency-based method to update modes of
clusters - A mixture of categorical and numerical data
k-prototype method
29The K-Medoids Clustering Method
- Find representative objects, called medoids, in
clusters - PAM (Partitioning Around Medoids, 1987)
- starts from an initial set of medoids and
iteratively replaces one of the medoids by one of
the non-medoids if it improves the total distance
of the resulting clustering - PAM works effectively for small data sets, but
does not scale well for large data sets - CLARA (Kaufmann Rousseeuw, 1990)
- CLARANS (Ng Han, 1994) Randomized sampling
- Focusing spatial data structure (Ester et al.,
1995)
30PAM (Partitioning Around Medoids) (1987)
- PAM (Kaufman and Rousseeuw, 1987), built in Splus
- Use real object to represent the cluster
- Select k representative objects arbitrarily
- For each pair of non-selected object h and
selected object i, calculate the total swapping
cost TCih - For each pair of i and h,
- If TCih lt 0, i is replaced by h
- Then assign each non-selected object to the most
similar representative object - repeat steps 2-3 until there is no change
31PAM Clustering Total swapping cost TCih?jCjih
32CLARA (Clustering Large Applications) (1990)
- CLARA (Kaufmann and Rousseeuw in 1990)
- Built in statistical analysis packages, such as
S - It draws multiple samples of the data set,
applies PAM on each sample, and gives the best
clustering as the output - Strength deals with larger data sets than PAM
- Weakness
- Efficiency depends on the sample size
- A good clustering based on samples will not
necessarily represent a good clustering of the
whole data set if the sample is biased
33CLARANS (Randomized CLARA) (1994)
- CLARANS (A Clustering Algorithm based on
Randomized Search) (Ng and Han94) - CLARANS draws sample of neighbors dynamically
- The clustering process can be presented as
searching a graph where every node is a potential
solution, that is, a set of k medoids - If the local optimum is found, CLARANS starts
with new randomly selected node in search for a
new local optimum - It is more efficient and scalable than both PAM
and CLARA - Focusing techniques and spatial access structures
may further improve its performance (Ester et
al.95)