Calinski-Harabasz index for estimating the number of clusters,
based on an observations/variables-matrix here. A distance based
version is available through cluster.stats.
Usage
calinhara(x,clustering,cn=max(clustering))
Value
Calinski-Harabasz statistic, which is
(n-cn)*sum(diag(B))/((cn-1)*sum(diag(W))). B being the
between-cluster means,
and W being the within-clusters covariance matrix.