code: def est_cluster_metric(self, cluster_opts=None):
          
          wcv, bcv, tv = self.est_cluster_covs(cluster_opts)
          return np.diag(bcv) / np.diag(tv)
docs: Returns an estimate of how much of the variance in the current posterior
          can be explained by a separation between *clusters*.