We sincerely thank the reviewer for the thoughtful and constructive feedback, which has helped us improve the clarity, rigor, and quality of our manuscript. Below, we provide point-by-point responses to the comments and outline the corresponding changes made to the paper.
## Weaknesses 1
> The paper is not organized clearly, which is not friendly for understanding. For example, there is a lack of preliminary details for how to model MB and other module in 3.1 GENERATIVE MODEL.
### Response
We thank the reviewer for highlighting the issues regarding organization and clarity.
First of all, we would like to clarify the modeling described in Section 3.1 GENERATIVE MODEL, where we present a **probabilistic framework**. The “modules” described in this section correspond to probabilistic distributions assumed for the variables. For instance, Equation 4, $w_{hr} \sim \mathcal{N}(w_{hr} | \mathbf{0}, \sigma^2 \mathbf{I})$ means that $w_{hr}$ follows a Gaussian distribution. Similarly, the "$\mathcal{MB}$ module" refers to a multivariate Bernoulli distribution, as defined in the Preliminary section (Equation 1)
We have structured the modeling and implementation details as follows:
- **In Section 3.1**, we describe the **probabilistic distributions for latent variables**, such as $p(z_{hr}) and p(w_{hr})$, as well as for observed triples $A_{hr, t}$ in Equation 7.
- **In Section 3.2**, we explain how to **learn the posterior distributions** of the latent variables $\mathcal{H}$ (Equation 11 to 15) using a VAE encoder (Equations 11–15)
- **In Section 3.3**, we describe the decoding process, which **uses the latent variables to complete missing triples**.
Thus, Section 3.1 introduces the core probabilistic modeling, while Sections 3.2 and 3.3 detail the implementation steps for encoding and decoding. We hope this clarifies the structure and flow of the manuscript.
## Weaknesses 2
> Figure 2 lacks of explanation, \textit{e.g.,} how the modules work together and match the equations in the main paper. The paper lacks the necessary reproduction file for the results.
### Response
We agree that Figure 2 would benefit from more detailed explanations. We have **updated the caption of Figure 2** to include a comprehensive description of how the modules work together.
We have prepared a detailed supplementary material file, including all the necessary reproduction files, such as datasets, code, and instructions for running the experiments.
## Weaknesses 3
> The paper lacks the analysis of time complexity as well as space complexity, which is necessary to study the efficiency of the model.
### Response
We appreciate the reviewer’s emphasis on the importance of time andspace complexity in model evaluation. In the revised manuscript, we have made the following enhancements to address this concern:
1. We provide a detailed **analysis of the time and space complexity** involved in computing the ELBO (Lines 281, Page 6 and Lines 796, Page 15).
2. We have also added a **new table (Table 7, Line 930, Page 18) comparing the parameter count and training time** of our method against the baseline model, SimKGC. The results indicate that introducing SBM and VAE **only increases the parameter count by 1.4M (0.8%)**. While the ELBO training requires **additional epochs to converge**, resulting in a slight increase in total training time, we find this increase to be **within an managable range**.
## Weaknesses 4
> The authors do not compare the model with other SOTA KGE methods, e.g.,[1][2][3]. The performance of, MRR in FB15K-237 is 0.36 while that of the proposed paper is 0.355. In this way, the performance of the proposed paper is not significant and the authors may better give a reasonable explanation.
### Response
Thank you for pointing out the need for comparison with SOTA methods.
1. **Baseline Inclusion**. In the revised manuscript, we **have included the SOTA results from KRACL [3]** for its best performance on WN18RR, as shown in Table 1. However, we would like to clarify that the **original submission already included several SOTA KGE methods, such as HittER and N-Former**, which achieve strong performance on FB15k-237, with MRR scores of 37.3 and 37.2, respectively. The inclusion of the additional baseline **does not alter the overall comparative assessment presented in our paper**.
2. **Performance on FB15k-237**. Regarding the relatively poorer performance on FB15k-237, we have dedicated substantial effort to **analyzing this issue in the original submission**, covering the discussion in detail **from Section 5.2 (MAIN RESULTS) to Section 6 (ANALYSIS)**. In summary, the dense connections and highly correlated relations in FB15k-237 lead to **fewer distinguishable clustering patterns**, making it challenging for our model to extract cluster-based features. Despite this, our method demonstrates **significant improvements on WN18RR and Wikidata5M and achieves a 2-3 point gain over text-based methods on FB15k-237**.