Author response
We thank you for your review. We appreciate that you have acknowledged the novelty of the problem and our contributions to interpretable approximations of kernel functions. We hope to resolve your remaining concerns in our response.
**Kernel $k$-means is not used as frequently in practice, spectral clustering or DBSCAN are preferred.**
Spectral clustering is equivalent to weighted kernel $k$-means [1]. Moreover, kernel $k$-means is closely related to clustering of distributions in the RKHS [2], which connects it to DBSCAN. Understanding its explainability is therefore not only important in its own right, but also a first step to extend the existing literature on interpretable clustering to more flexible methods.
[1] Dhillon, Inderjit S., Yuqiang Guan, and Brian Kulis. "Kernel k-means: spectral clustering and normalized cuts." Proceedings of the tenth ACM SIGKDD international conference on Knowledge discovery and data mining. 2004.
[2] Vankadara, Leena C., et al. "Recovery guarantees for kernel-based clustering under non-parametric mixture models." International Conference on Artificial Intelligence and Statistics. PMLR, 2021.
**One should jointly optimize clustering cost with the decision tree.**
The recent works on explainable clustering show a clear preference towards first finding the partition, and then fitting a decision tree (for two recent analysis, see [3] and [4]). Our paper demonstrates when and how worst-case guarantees can be leveraged to the kernel setting. Joint optimization of clustering cost and decision tree is a different principle. We believe it is beyond the scope of this paper.
[3] Charikar, Moses, and Lunjia Hu. "Near-optimal explainable k-means for all dimensions." Proceedings of the 2022 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA). Society for Industrial and Applied Mathematics, 2022.
[4] Esfandiari, Hossein, Vahab Mirrokni, and Shyam Narayanan. "Almost tight approximation algorithms for explainable clustering." Proceedings of the 2022 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA). Society for Industrial and Applied Mathematics, 2022.
**The paper completely discards the approach of fitting a tree on the kernel k-means clustering labels as a supervised problem.**
The (Kernel) IMM algorithm is supervised, it minimizes the number of mistakes at every iteration. Similarly, the Kernel Expand algorithm can also be run on an empty tree, choosing threshold cuts that minimize the number of points assigned to a wrong cluster. For completeness, we compared CART to Kernel IMM on three clustering datasets (checking Laplace and Gaussian kernel, as well as Taylor and kernel matrix features as in the experiments). Both are similar in terms of price of explainability.
| Datasets | Kernel | POE(Kernel IMM,KKM) | POE(CART, KKM) |
| --- | --- | --- | --- |
| Pathbased | Gaussian (kernel matrix features) | 1.06645 | 1.07004 |
| Aggregation | Laplace (kernel matrix features) | 1.00042 | 1.00065 |
| Flame | Gaussian (Taylor features) | 1.02256 | 1.02732 |
Decision tree algorithms such as CART may perform well in practice. However, we do not know of any theoretical guarantees in the clustering regime. This is precisely what has motivated the extensive works on explainable $k$-means that this paper builds on.
**The theoretical bounds on the price of explainability seem to be quite high and not really helpful in practice. For example, does the bound of mean that say for MNIST the tree will do worse by (784 * 100) times than the unconstrained kernel k-means?**
Using explainable clustering directly on images (such as MNIST) may not be meaningful because pixel values are typically not interpretable. We agree that experiments indicate that the worst-case upper bounds can be loose. It is however actually a strength of worst-case approximation guarantees that they hold for all possible data, in spite of being potentially loose for well-behaved data. Previous works on explainable $k$-means have also reported approximation ratios close to $1$ in practice (for example, see Table 1 in [5]). It is likely that algorithms tailored to specific kernels allow for tighter bounds (much in the same way as for different $\|\cdot \|_p$-norms, distinct algorithms have emerged [6]). Better guarantees could possibly be given under assumptions on the distribution, and these are interesting avenues to explore in the future. Since explainable kernel clustering has not been studied prior to our work, this paper takes a broader perspective, proposing general methods that work for a large class of kernels, regardless of the distribution.
[5] Laber, Eduardo S., and Lucas Murtinho. "On the price of explainability for some clustering problems." International Conference on Machine Learning. PMLR, 2021.
[6] Gamlath, Buddhima, et al. "Nearly-tight and oblivious algorithms for explainable clustering." Advances in Neural Information Processing Systems 34 (2021): 28929-28939.