方法一：

K平均算法（K-means聚類分析）
在下面的誤差平方和圖中，拐點（bend or elbow）的位置對應的x軸即k-means聚類給出的合適的類的個數(shù)。

> n = 100
> g=6
> set.seed(g)
> d <- data.frame(x = unlist(lapply(1:g, function(i) rnorm(n/g, runif(1)*i^2))), y = unlist(lapply(1:g, function(i) rnorm(n/g, runif(1)*i^2))))
> plot(d)
>

> mydata <- d
>
> wss <- (nrow(mydata)-1)*sum(apply(mydata,2,var))
> for (i in 2:15)
+ wss[i] <- sum(kmeans(mydata,centers=i)$withinss)
> ###這里的wss(within-cluster sum of squares)是組內(nèi)平方和
> plot(1:15, wss, type="b", xlab="Number of Clusters",ylab="Within groups sum of squares")
>

由上圖可以看出，該方法給出合理的類別個數(shù)是4個。
方法二：

K中心聚類算法（K-mediods）
使用fpc包里的pamk函數(shù)來估計類的個數(shù)：

> library(cluster)
Warning message:
程輯包‘cluster’是用R版本3.2.3 來建造的
> library(fpc)
> pamk.best <- pamk(d)
> cat("number of clusters estimated by optimum average silhouette width:", pamk.best$nc, "\n")
number of clusters estimated by optimum average silhouette width: 4
> plot(pam(d, pamk.best$nc))

sihouette值是用來表示某一個對象和它所屬類的凝合力強度以及和其他類分離強度的，值范圍為-1到1，值越大表示該對象越匹配所屬類以及和鄰近類有多不匹配。
所以從上圖sihouette plot中可以看出，該方法給出的合理類的個數(shù)為4個。
方法三：
基于Calinsky Criterion
> require(vegan)
載入需要的程輯包：vegan
載入需要的程輯包：permute
載入需要的程輯包：lattice
This is vegan 2.4-0
Warning messages:
1: 程輯包‘vegan’是用R版本3.2.5 來建造的
2: 程輯包‘permute’是用R版本3.2.5 來建造的
3: 程輯包‘lattice’是用R版本3.2.3 來建造的
> fit <- cascadeKM(scale(d, center = TRUE,  scale = TRUE), 1, 10, iter = 1000)
> plot(fit, sortg = TRUE, grpmts.plot = TRUE)
> calinski.best <- as.numeric(which.max(fit$results[2,]))
> cat("Calinski criterion optimal number of clusters:", calinski.best, "\n")
Calinski criterion optimal number of clusters: 5
>

由上圖我們可以看到，根據(jù)Calinsky標準，得到類的個數(shù)是5個。
方法四：

基于模型假設的聚類，利用的是mclust包：

> library(mclust)
    __  ___________    __  _____________
   /  |/  / ____/ /   / / / / ___/_  __/
  / /|_/ / /   / /   / / / /\__ \ / /   
 / /  / / /___/ /___/ /_/ /___/ // /    
/_/  /_/\____/_____/\____//____//_/    version 5.1
Type 'citation("mclust")' for citing this R package in publications.
Warning message:
程輯包‘mclust’是用R版本3.2.4 來建造的
> d_clust <- Mclust(as.matrix(d), G=1:20)
> m.best <- dim(d_clust$z)[2]
> cat("model-based optimal number of clusters:", m.best, "\n")
model-based optimal number of clusters: 4
> plot(d_clust)
Model-based clustering plots:

1: BIC
2: classification
3: uncertainty
4: density

方法五：

基于AP算法的聚類

> library(apcluster)

載入程輯包：‘a(chǎn)pcluster’

The following object is masked from ‘package:stats’:

    heatmap

Warning message:
程輯包‘a(chǎn)pcluster’是用R版本3.2.5 來建造的
> d.apclus <- apcluster(negDistMat(r=2), d)
> cat("affinity propogation optimal number of clusters:", length(d.apclus@clusters), "\n")
affinity propogation optimal number of clusters: 4
> #4 得出的分類個數(shù)
> heatmap(d.apclus)
> plot(d.apclus, d)
>