We thank all the reviewers for their careful reading of the paper and helpful comments. We will take these comments into account when improving the paper.
Many reviewers mention things not included in the paper. Unfortunately we could not provide a detailed descriptions of geometries used, or the connectome datasets, due to space limits. These issues were described well in our references. Our dataset is the same as used in Allard & Serrano (2020), limited to the subset that is publicly available. Instead, we have focused more on extensive statistical analysis and testing, which we believe is lacking in most papers in the field.
We agree that visualizations of our embeddings would be very helpful. In particular, such a visual analysis could also help us to determine the reasons why particular geometries are better or worse for each connectome. The framework we are using (supplemental material) does have visualization capabilities. Unfortunately, we focus on 3D embeddings, which are difficult to render in paper, and again, this uses up a space, which is limited. We will consider this addition for the future version of the paper.
Answers to questions:
- HRG is the Hyperbolic Random Graph model, explained in Section 2. We will explain this acronym, thanks for catching this! In this model, the edges are indeed assumed to be independent (once the points in the hyperbolic disk are assigned).
- We are not sure how ablation study would be helpful for determining M. Our M=20000 is large enough to improve on the state-of-the-art methods in most cases. We also include an analysis comparing M=20000 to M=100000.
- We have indeed given an incorrect formula for the Minkowski inner product. Thanks for catching this! We do not understand the comment about "it is unclear what hyperbolic model the paper talks about", in all places where the model is relevant, we state the model used (Minkowski hyperboloid model to explain the basics, Poincare disk for visualization, horocyclic coordinate system to explain Solv); in the core paper, it does not matter which model is used underneath (we only care about distances between points; in some applications, the model may matter due to numerical precision issues, but our disks are small enough that these issues do not matter). (The Minkowski hyperboloid model is sometimes called the Lorentz model in ML, but this does not seem to be a popular name in general, e.g. Wikipedia does not include that name.)
- d(a1,a2) is the geometric distance between point a1 and a2. We need to approximate like that because we have found out that the method used in (Kopczyński & Celińska-Kopczyńska, 2022) does not yield usable results once the two points are too far away.
Other questions of Reviewer Qkfi: (e) in our experiments, we use Dijkstra's search just once, to construct a MxM array, which is then used for all iterations and all connectomes. This computation takes noticeable time, but it is a negligible part of the whole setup. In the last line of Page 8 we mention that the embedding time is 220 seconds per run on Mouse3 in H3 (other times can be found in the supplemental material). The dominating part is computing the new loglikelihood after the change, which is done in time O(n) where n is the number of the nodes; this is done in every iteration, so computational complexity is O(n*N_s). In our experiments, the number of iterations N_s is also linearly dependent on n. (f) we use the whole dataset to compute embeddings.
We mention "universal cover" / $\text{SL}(2,\mathcal{R}) because this name is used in some sources (e.g. Thurston 1982), so we include it so readers seeing that name in other sources know what we are talking about ("also known as"). This name is not intuitive, and explaining it would not fit the scope of the paper. Twisted product of H2xR works just like Nil (twisted product of E2xR) but the base is the hyperbolic plane.
Reviewer 7zXB: Our assumption is that the geometry of the network is different from the geometry of the Euclidean space that the connectome physically is in. This is the topic of Allard & Serrano (2020) where comparisons to physical embeddings are included, our purpose it to determine what the best network geometry it is. The paper cited by the reviewer seems to only study the physical geometry, and also we see no mention of Voronoi tessellations in it.