calinhara function

Calinski-Harabasz index