estimate_k function

estimate_k: estimate K using the eigengap or multimodality gap heuristics