### W2: Synthetic misspecifications
We agree with the reviewer that our experiments showcasing the (lack of) robustness to assumption violations should be more extensive.
We appreciate the suggestion to explore the Mixture of Gaussians case (or, similarly, the Mixture of vMFs). In this setup, each cluster is represented by a set of fixed, proximal vectors, and samples are drawn from one of these vectors according to a Gaussian distribution (or a vMF, provided the proximal vectors lie on the hypersphere). The issue with this setup is that it is also identifiable up to a linear transformation. **The Mixture of vMFs case is a direct application of the results already proven in our paper**, with the fixed proximal vectors serving as the new cluster vectors. While the Gaussian case (when samples may fall of the hypersphere) is not directly covered by our theorem, it follows a similar reasoning and has an analogous proof to ours.
Instead, **we explore the amount of violation caused by distributions of varying scales and shapes**. We do this by utilising a _generalised Normal distribution_ (where the shape parameter $\beta=1$ corresponds to a Laplace and $\beta=2$ to a Normal distribution) and, furthermore, we vary the concentration parameter alpha (that corresponds to what we previously also called concentration parameter in the case of the vMF conditional).
In Appx. D and Fig. 4, the left plot reports the recovery of latent $z$, the middle the recovery of cluster vectors $v_c$, and the right the relationship between the concentration parameter and the average norm of representations (in a 1D case). These scores are the average across 5 seeds.
In Appx. D, Fig. 5, we also show that **the fat tails of the Laplace distribution make it significantly deviate from the vMF case**, leading to a much lower identifiability scores. To show this this, we use a _truncated Laplace_ conditional, where the truncation parameter restricts the distribution’s support to the symmetrical interval [-Truncation;Truncation]. As Fig. 5 Right demonstrates, a higher truncation parameter (i.e., a distribution with **fatter tails**) **corresponds to lower $R^2$ scores**, this relationship holds over multiple concentration parameters and seeds.
From these figures, we conclude:
- [Fig. 4 Left and Middle] we see that the **distribution's shape and concentration plays a critical role**. Moving towards a Laplace distribution leads to more performance degradation, which is due to its fatter tails (shown in Fig. 5 Right). This difference is even more pronounced for high concentration values.
- [Fig. 4 Right] We show that there is a **non-trivial relationship between the norm of the representation and the concentration as a function of the shape parameter $\beta$**. For a given concentration value $\alpha$, decreasing beta (i.e., moving towards a Laplace distribution) leads to a bigger average norm.
### W3: Weak ImageNet-X results
We re-evaluated our results on ImageNet-X and now provide statistical tests to have a better overview on the performance of trained models in recovering (proxy) latent features. As our **revised Fig. 3** suggests, **almost all features are recovered better than chance performance** (except “brighter” - second subplot in the second row) both for trained ViTs and ResNet50. More importantly, as **there are no true latents available for ImageNet-X, the ones we evaluate against are only binary proxies defined by human evaluators**. We hypothesize that the proxy nature and especially the binary latents are the cause of degraded performance
### Q1: The impact of repeated elements
The questions is interesting and relevant. We distinguish the theoretical and more practical cases.
In our theoretical setup, the data generating process (DGP) is defined in a way that the **occurrence of repeated elements is bound to happen**. Namely, the vMF conditional distribution has full support over the hypersphere and, hence, any sample could have resulted from any instance. **Nonetheless, the prediction of the theorem holds**: the inferred parameter vectors $w_i$’s pick up the values of the cluster vectors that correspond to the instance label $i$. The latents are reconstructed up to linear transformations and the score is unchanged.
The phenomenon may seem more perplexing **in practical, finite dataset cases. Here, seemingly the same sample is mapped to two different indices**. Yet, this _does not change the learned representation or the corresponding parameter vector $w$_. The two parameter vectors $w_i$ and $w_j$ will, in fact, match. **What changes is the value of the instance discrimination objective (DIET)**, as the model is now uncertain in assigning the index. The logit values of index $i$ and $j$ will match. However, again, this does not affect the learned representation.