Title: YOU CAN TEACH AN OLD DOG NEW TRICKS! ON TRAINING KNOWLEDGE GRAPH EMBEDDINGS

Abstract: Knowledge graph embedding (KGE) models learn algebraic representations of the entities and relations in a knowledge graph. A vast number of KGE techniques for multi-relational link prediction have been proposed in the recent literature, often with state-of-the-art performance. These approaches differ along a number of dimensions, including different model architectures, different training strategies, and different approaches to hyperparameter optimization. In this paper, we take a step back and aim to summarize and quantify empirically the impact of each of these dimensions on model performance. We report on the results of an extensive experimental study with popular model architectures and training strategies across a wide range of hyperparameter settings. We found that when trained appropriately, the relative performance differences between various model architectures often shrinks and sometimes even reverses when compared to prior results. For example, RESCAL (Nickel et al., 2011) , one of the first KGE models, showed strong performance when trained with state-of-the-art techniques; it was competitive to or outperformed more recent architectures. We also found that good (and often superior to prior studies) model configurations can be found by exploring relatively few random samples from a large hyperparameter space. Our results suggest that many of the more advanced architectures and techniques proposed in the literature should be revisited to reassess their individual benefits. To foster further reproducible research, we provide all our implementations and experimental results as part of the open source LibKGE framework.

1. INTRODUCTION:

Knowledge graph embedding (KGE) models learn algebraic representations, termed embeddings, of the entities and relations in a knowledge graph. They have been successfully applied to knowledge graph completion (Nickel et al., 2015) as well as in downstream tasks and applications such as recommender systems (Wang et al., 2017) or visual relationship detection (Baier et al., 2017) .

A vast number of different KGE models for multi-relational link prediction have been proposed in the recent literature; e.g., RESCAL (Nickel et al., 2011) , TransE (Bordes et al., 2013) , DistMult, ComplEx (Trouillon et al., 2016) , ConvE (Dettmers et al., 2018) , TuckER (Balazevic et al., 2019) , RotatE (Sun et al., 2019a) , SACN (Shang et al., 2019) , and many more. Model architectures generally differ in the way the entity and relation embeddings are combined to model the presence or absence of an edge (more precisely, a subject-predicate-object triple) in the knowledge graph; they include factorization models (e.g., RESCAL, DistMult, ComplEx, TuckER), translational models (TransE, RotatE) , and more advanced models such as convolutional models (ConvE) . In many cases, the introduction of new models went along with new approaches for training these modelse.g., new training types (such as negative sampling or 1vsAll scoring), new loss functions (such as Table 1 : Selected KGE models and training strategies from the literature. Entries marked in bold were introduced (or first used) in the context of KGE in the corresponding publication.

pairwise margin ranking or binary cross entropy), new forms of regularization (such as unweighted and weighted L2), or the use of reciprocal relations (Kazemi & Poole, 2018; Lacroix et al., 2018) and ablation studies were not always performed. Table 1 shows an overview of selected models and training techniques along with the publications that introduced them.

The diversity in model training makes it difficult to compare performance results for various model architectures, especially when results are reproduced from prior studies that used a different experimental setup. Model hyperparameters are commonly tuned using grid search on a small grid involving hand-crafted parameter ranges or settings known to "work well" from prior studies. A grid suitable for one model may be suboptimal for another, however. Indeed, it has been observed that newer training strategies can considerably improve model performance (Kadlec et al., 2017; Lacroix et al., 2018; Salehi et al., 2018) .

In this paper, we take a step back and aim to summarize and quantify empirically the impact of different model architectures and different training strategies on model performance. We performed an extensive set of experiments using popular model architectures and training strategies in a common experimental setup. In contrast to most prior work, we considered many training strategies as well as a large hyperparameter space, and we performed model selection using quasi-random search (instead of grid search) followed by Bayesian optimization. We found that this approach was able to find good (and often superior to prior studies) model configurations with relatively low effort.

Our study complements and expands on the results of Kotnis & Nastase (2018) (focus on negative sampling) and Mohamed et al. (2019) (focus on loss functions) as well as similar studies in other areas, including language modeling (Melis et al., 2017) , generative adversarial networks (Lucic et al., 2018) , or sequence tagging (Reimers & Gurevych, 2017) .

We found that when trained appropriately, the performance of a particular model architecture can by far exceed the performance observed in older studies. For example, RESCAL (Nickel et al., 2011) , which constitutes one of the first KGE models but is rarely considered in newer work, showed very strong performance in our study: it was competitive to or outperformed more recent architectures such as ConvE (Dettmers et al., 2018) and TuckER (Balazevic et al., 2019) . More generally, we found that the relative performance differences between various model architectures often shrunk and sometimes even reversed when compared to prior results. This suggests that (at least currently) training strategies have a significant impact on model performance and may account for a substantial fraction of the progress made in recent years. We also found that suitable training strategies and hyperparameter settings vary significantly across models and datasets, indicating that a small search grid may bias results on model performance. Fortunately, as indicated above, large hyperparameter spaces can be (and should be) used with little additional training effort. To facilitate such efforts, we provide implementations of relevant training strategies, models, and evaluation methods as part of the open source LibKGE framework,1 which emphasizes reproducibility and extensibility.

Our study focuses solely on pure KGE models, which do not exploit auxiliary information such as textual data or logical rules (Wang et al., 2017) . Since many of the studies on these non-pure models did not (and, to be fair, could not) use current training strategies and consequently underestimated the performance of pure KGE models, their results and conclusions need to be revisited.

2. KNOWLEDGE GRAPH EMBEDDINGS: MODELS, TRAINING, EVALUATION:

The literature on KGE models is expanding rapidly. We review selected architectures, training methods, and evaluation protocols; see Table 1 . The table examplarily indicates that new model architectures are sometimes introduced along with new training strategies (marked bold). Reasonably recent survey articles about KGE models include Nickel et al. (2015) and Wang et al. (2017) .

Multi-relational link prediction. KGE models are typically trained and evaluated in the context of multi-relational link prediction for knowledge graphs (KG). Generally, given a set E of entities and a set R of relations, a knowledge graph K ⊆ E × R × E is a set of subject-predicate-object (spo) triples. The goal of multi-relational link prediction is to "complete the KG", i.e., to predict true but unobserved triples based on the information in K. Common approaches include rule-based methods (Galarraga et al., 2013; Meilicke et al., 2019) , KGE methods (Nickel et al., 2011; Bordes et al., 2013; Trouillon et al., 2016; Dettmers et al., 2018) , and hybrid methods (Guo et al., 2018) .

Knowledge graph embeddings (KGE).:

A KGE model associates with each entity i ∈ E and relation k ∈ R an embedding e i ∈ R de and r k ∈ R dr in a low-dimensional vector space, respectively; here d e , d r ∈ N + are hyperparameters for the embedding size. Each particular model uses a scoring function s : E × R × E → R to associate a score s(i, k, j) with each potential triple (i, k, j) ∈ E × R × E. The higher the score of a triple, the more likely it is considered to be true by the model. The scoring function takes form s(i, k, j) = f (e i , r k , e j ), i.e., depends on i, k, and j only through their respective embeddings. Here f , which represents the model architecture, may be either a fixed function or learned (e.g., f may be a parameterized function).

Evaluation. The most common evaluation task for KGE methods is entity ranking, which is a form of question answering. The available data is partitioned into a set of training, validation, and test triples. Given a test triple (i, k, j) (unseen during training), the entity ranking task is to determine the correct answer-i.e., the missing entity j and i, resp.-to questions (i, k, ?) and (?, k, j). To do so, potential answer triples are first ranked by their score in descending order. All triples but (i, k, j) that occur in the training, validation, or test data are subsequently filtered out so that other triples known to be true do not affect the ranking. Finally, metrics such as the mean reciprocal rank (MRR) of the true answer or the average HITS@k are computed; see the appendix. KGE models. KGE model architectures differ in their scoring function. We can roughly classify models as decomposable or monolithic: the former only allow element-wise interactions between (relation-specific) subject and object embeddings, whereas the latter allow arbitrary interactions. More specifically, decomposable models use scoring functions of form s(i, k, j) = f ( z g([h 1 (e i , e r ) • h 2 (e r , e j )] z )), where • is any element-wise function (e.g., multiplication), h 1 and h 2 are functions that obtain relation-specific subject and object embeddings, resp., and g and f are scalar functions (e.g., identity or sigmoid). The most popular models in this category are perhaps RESCAL (Nickel et al., 2011 ), TransE (Bordes et al., 2013) , DistMult (Yang et al., 2015) , ComplEx (Trouillon et al., 2016) , and ConvE (Dettmers et al., 2018) . RESCAL's scoring function is bilinear in the entity embeddings: it uses s(i, k, j) = e T i R k e j , where R k ∈ R de×de is a matrix formed from the entries of r k ∈ R dr (where d r = d 2 e ). DistMult and ComplEx can be seen as constrained variants of RESCAL with smaller relation embeddings (d r = d e ). TransE is a translation-based model and uses negative distance -e i + r k -e j p between e i + r k and e j as score, commonly using the L1 norm (p = 1) or the L2 norm (p = 2). Finally, ConvE uses a 2D convolutional layer and a large fully connected layer to obtain relation-specific entity embeddings (i.e., in h 1 above). Other recent examples for decomposable models include TuckER (Balazevic et al., 2019) , RotatE (Sun et al., 2019a) , and SACN (Shang et al., 2019) . Decomposable models are generally fast to use: once the relation-specific embeddings are (pre-)computed, score computations are cheap. Monolithic models-such as ConvKB or KBGAT-do not decompose into relation-specific embeddings: they take form s(i, k, j) = f (e i , r k , e j ). Such models are more flexible, but they are also considerably more costly to train and use. It is currently unclear whether monolithic models can achieve comparable or better performance than decomposable models Sun et al. (2019b) .

Training type. There are three commonly used approaches to train KGE models, which differ mainly in the way negative examples are generated. First, training with negative sampling (NegSamp) (Bordes et al., 2013) obtains for each positive triple t = (i, k, j) from the training data a set of (pseudo-)negative triples obtained by randomly perturbing the subject, relation, or object position in t (and optionally verifying that the so-obtained triples do not exist in the KG). An alternative approach (Lacroix et al., 2018) , which we term 1vsAll, is to omit sampling and take all triples that can be obtained by perturbing the subject and object positions as negative examples for t (even if these tuples exist in the KG). 1vsAll is generally more expensive than NegSamp, but it is feasible (and even surprisingly fast in practice) if the number of entities is not excessively large. Finally, Dettmers et al. (2018) proposed a training type that we term KvsAll2 : this approach (i) constructs batches from non-empty rows (i, k, * ) or ( * , k, j) instead of from individual triples, and (ii) labels all such triples as either positive (occurs in training data) or negative (otherwise).

Loss functions. Several loss functions for training KGEs have been introduced so far. RESCAL originally used squared error between the score of each triple and its label (positive or negative). TransE used pairwise margin ranking with hinge loss (MR), where each pair consists of a positive triple and one of its negative triples (only applicable to NegSamp and 1vsAll) and the margin η is a hyperparameter. Trouillon et al. (2016) proposed to use binary cross entropy (BCE) loss: it applies a sigmoid to the score of each (positive or negative) triple and uses the cross entropy between the resulting probability and that triple's label as loss. BCE is suitable for multi-class and multi-label classification. Finally, Kadlec et al. (2017) used cross entropy (CE) between the model distribution (softmax distribution over scores) and the data distribution (labels of corresponding triples, normalized to sum to 1). CE is more suitable for multi-class classification (as in NegSamp and 1vsAll), but it has also been used in the multi-label setting (KvsAll). Mohamed et al. (2019) found that the choice of loss function can have a significant impact on model performance, and that the best choice is data and model dependent. Our experimental study provides additional evidence for this finding.

Reciprocal relations. Kazemi & Poole (2018) and Lacroix et al. (2018) introduced the technique of reciprocal relations into KGE training. Observe that during evaluation and also most training methods discussed above, the model is solely asked to score subjects (for questions of form (?, k, j)) or objects (for questions of form (i, k, ?)). The idea of reciprocal relations is to use separate scoring functions s sub and s obj for each of these tasks, resp. Both scoring functions share entity embeddings, but they do not share relation embeddings: each relation thus has two embeddings. 3 The use of reciprocal relations may decrease the computational cost (as in the case of ConvE), and it may also lead to better model performance Lacroix et al. ( 2018) (e.g., for relations in which one direction is easier to predict). On the downside, the use of reciprocal relations means that a model does not provide a single triple score s(i, k, j) anymore; generally, s sub (i, k, j) = s obj (i, k, j). Kazemi & Poole (2018) proposed to take the average of the two triple scores and explored the resulting models.

Regularization. The most popular form of regularization in the literature is L2 regularization on the embedding vectors, either unweighted or weighted by the frequency of the corresponding entity or relation (Yang et al., 2015) . Lacroix et al. (2018) proposed to use L3 regularization. TransE normalized the embeddings to unit norm after each update. ConvE used dropout (Srivastava et al., 2014) in its hidden layers (and only in those). In our study, we additionally consider L1 regularization and the use of dropout on the entity and/or relation embeddings.

Hyperparameters. Many more hyperparameters have been used in prior work. This includes, for example, different methods to initialize model parameters, different optimizers, different optimizer parameters such as the learning rate or batch size, the number of negative examples for NegSamp, the regularization weights for entities and relations, and so on. To deal with such a large search space, most prior studies favor grid search over a small grid where most of these hyperparameters remain fixed. As discussed before, this approach may lead to bias in the results, however.

3. EXPERIMENTAL STUDY:

In this section, we report on the design and results of our experimental study. We focus on the most salient points here; more details can be found in the appendix.

3.1. EXPERIMENTAL SETUP:

Datasets. We used the FB15K-237 (Toutanova & Chen, 2015) (extracted from Freebase) and WNRR (Dettmers et al., 2018) (extracted from WordNet) datasets in our study. We chose these datasets because (i) they are frequently used in prior studies, (ii) they are "hard" datasets that have been designed specifically to evaluate multi-relational link prediction techniques, (iii) they are diverse in that relative model performance often differs, and (iv) they are of reasonable size for a large study. Dataset statistics are given in Table 4 in the appendix.

Models. We selected RESCAL (Nickel et al., 2011) , TransE (Bordes et al., 2013) , DistMult (Yang et al., 2015) , ComplEx (Trouillon et al., 2016) and ConvE (Dettmers et al., 2018) for our study. These models are perhaps the most well-known KGE models and include both early and more recent models. We did not consider monolithic models due to their excessive training cost.

Evaluation. We report filtered MRR (%) and filtered HITS@10 (%); see Sec. 2 and the appendix for details. We use test data to compare final model performance (Table 2 ) and validation data for our more detailed analysis.

Hyperparameters. We used a large hyperparameter space to ensure sure that suitable hyperparameters for each model are not excluded a priori. We included all major training types (NegSamp, 1vsAll, KvsAll), use of reciprocal relations, loss functions (MR, BCE, CE), regularization techniques (none/L1/L2/L3, dropout), optimizers (Adam, Adagrad), and initialization methods (4 in total) used in the KGE community as hyperparameters. We considered three embeddings sizes (128, 256, 512) and used separate weights for dropout/regularization for entity and relation embeddings. Table 5 in the appendix provides additional details. To the best of our knowledge, no prior study used such a large hyperparameter search space.

Training. All models were trained for a maximum of 400 epochs. We validated models using filtered MRR (on validation data) every five epochs and performed early stopping with a patience of 50 epochs. To keep search tractable, we stopped training on models that did not reach ≥ 5% filtered MRR after 50 epochs; in our experience, such configurations did not produce good models.

Model selection. Model selection was performed using filtered MRR on validation data. We used the Ax framework (https://ax.dev/) to conduct quasi-random hyperparameter search via a Sobol sequence. Quasi-random search methods aim to distribute hyperparameter settings evenly and try avoid "clumping" effects (Bergstra & Bengio, 2012) . More specifically, for each dataset and model, we generated 30 different configurations per valid combination of training type and loss function (2 for TransE, which only supports NegSamp+MR and NegSamp+CE; 7 for all other models). After the quasi-random hyperparameter search, we added a short Bayesian optimization phase (best configuration so far + 20 new trials, using expected improvement; also provided by Ax) to tune the numerical hyperparameters further. Finally, we trained five models with the so-obtained hyperparameter configuration and selected the best-performing model according to validation MRR as the final model. The standard deviation of the validation MRR was relatively low; see Table 9 .

Reproducibility. We implemented all models, training strategies, evaluation, and hyperparameter search in the LibKGE framework. The framework emphasizes reproducibility and extensibility and is highly configurable. We provide all configurations used in this study as well as the detailed data of our experiments for further analysis. 4 We hope that these resources facilitate the evaluation of KGE models in a comparable environment.

3.2. COMPARISON OF MODEL PERFORMANCE:

Table 2 shows the filtered MRR and filtered Hits@10 on test data of various models both from prior studies and the best models (according filtered validation MRR) found in our study.

First reported performance vs. our observed performance. We compared the first reported performance on FB15K-237 and WNRR ("First" block of Table 2 ) with the performance obtained in our study ("Ours" block). We found that the performance of a single model can vary wildly across studies. For example, ComplEx was first run on FB15K-237 by Dettmers et al. (2018) , where it achieved a filtered MRR of 24.7%. This is a relatively low number by today's standards. In our FB15K-237 WNRR MRR Hits@10 MRR Hits@10 First RESCAL (Wang et al., 2019) 27.0 42.7 42.0 44.7 TransE (Nguyen et al., 2018) 29.4 46.5 22.6 50.1 DistMult (Dettmers et al., 2018) 24.1 41.9 43.0 49.0 ComplEx (Dettmers et al., 2018) 24.7 42.8 44.0 51.0 ConvE (Dettmers et al., 2018) 32 Relative model performance (our study). Next, we compared the performance across the models models used in our study ("Ours" block). We found that the relative performance differences between model architectures often shrunk and sometimes reversed when compared to prior results ("First" block). For example, ConvE showed the best overall performance in prior studies, but is consistently outperformed by RESCAL, DistMult, and ComplEx in our study. As another example, RESCAL (Nickel et al., 2011) , which constitutes one of the first KGE models but is rarely considered in newer work, showed strong performance and outperformed all models except ComplEx. This suggests that (at least currently) training strategies have a significant impact on model performance and may account for a large fraction of the progress made in recent years.

Relative model performance (overall). Table 2 additionally shows the best performance results obtained in prior studies of selected recent models ("Recent" block) and very large models with very good performance ("Large" block). We found that TuckER (Balazevic et al., 2019) , RotatE (Sun et al., 2019a) , and SACN (Shang et al., 2019) all achieve state-of-the-art performance, but the performance difference to the best prior models ("Ours" block) in terms of MRR is small or even vanishes. Even for HITS@10, which we did not consider for model selection, the advantage of more recent models is often small, if present. The models in the last block ("Large") show the best performance numbers we have seen in the literature. For DistMult and ComplEx, these numbers have been obtained using very large embedding sizes (up to 4000). Such large models were not part of our search space.

Limitations. We note that all models used in our study are likely to yield better performance when hyperparameter tuning is further improved. For example, we found a ComplEx configuration (of size 100) that achieved 49.0% MRR and 56.5% Hits@10 on WNRR, and a RESCAL configuration that achieved 36.3% MRR and 54.6% Hits@10 on FB15K-237 in preliminary experiments. Sun et al. (2019a) report on a TransE configuration that achieved 33.3% MRR (52.2% Hits@10) on FB15K-237. Since the focus of this study is to compare models in a common experimental setup and without manual intervention, we do not further report on these results.

3.3. IMPACT OF HYPERPARAMETERS:

Anatomy of search space. Figure 1 shows the distribution of filtered MRR for each model and dataset (on validation data, Sobol configurations only). Each distribution consists of roughly 200 different quasi-random hyperparameter configurations (except TransE, for which 60 configurations have been used). Perhaps unsurprisingly, we found that all models showed a wide dispersion on both datasets and only very few configurations achieved state-of-the-art results. There are, however, notable differences across datasets and models. On FB15K-237, for example, the median ConvE configuration is best; on WNRR, DistMult and Complex have a much higher median MRR. Generally, the impact of the hyperparameter choice is more pronounced on WNRR (higher variance) than on FB15K-237.

Best configurations (quasi-random search). The hyperparameters of the best performing models during our quasi-random search are reported in Table 3 (selected hyperparameters) and Tables 6 and 7 (in the appendix, all hyperparameters). First, we found that the optimum choice of hyperparameters is often model and dataset dependent: with the exception of the loss function (discussed below), almost no hyperparameter setting was consistent across datasets and models. For example, the NegSample training type was frequently selected on FB15K-237 but not on WNRR, where KvsAll was more prominent. Our findings provide further evidence that grid search on a small search space is not suitable to compare model performance because the result may be heavily influenced by the specific grid points being used. Moreover, the budget that we used for quasi-random search was comparable to a grid search over a small (or at least not very large) grid: we tested merely 30 different model configurations for every combination of training type and loss function.

It is therefore both feasible and beneficial to work with a large search space.

Best configurations (Bayesian optimization). We already obtained good configurations solely from quasi-random search. After Bayesian optimization (which we used to tune the numerical hyperparameters), the performance on test data almost always improved only marginally. The best hyperparameters after the Bayesian optimization phase are reported in Table 8 in the appendix.

Impact of specific hyperparameters. It is difficult to asses the impact of each hyperparameter individually. As a proxy for the importance of each hyperparameter setting, we report in Table 3 and Tables 6 and 7 (appendix) the best performance in terms of filtered MRR on validation data for a different choice of value for each hyperparameter (difference in parenthesis). For example, the best configuration during quasi-random search for RESCAL on FB15K-237 did not use reciprocal relations. The best configuration that did use reciprocal relations had an MRR that was 0.5 points lower (35.6 instead of 36.1). This does not mean that the gap is explained by the use of the reciprocal relations alone, as other hyperparameters may also be different, but it does show the best performance we obtained when enforcing reciprocal relations. We can see that most hyperparameters appear to have moderate impact on model performance: many hyperparameter settings did not produce a significant performance improvement that others did not also produce. A clear exception here is the use of the loss function, which is consistent across models and datasets (always CE), partly with large performance improvements over alternative loss functions. Impact of loss function. In order to further explore the impact of the loss function, we show the distribution of the filtered MRR on validation data for every (valid) combination of training type and loss function in Figure 10 . We found that the use of CE (the three rightmost bars) generally lead to better configurations than using other losses. This is surprising (and somewhat unsatisfying) in that we use the CE loss, which is a multi-class loss, for a multi-label task. Similar observations have been made for computer vision tasks (Joulin et al., 2016; Mahajan et al., 2018) though. Note that the combination of KvsAll with CE (on the very right in figure) has not been explored previously.

3.4. MODEL SELECTION:

We briefly summarize the behavior of various models during model selection. Figure 3 shows the best validation MRR (over all quasi-random hyperparameter configurations) achieved by each model when trained for the specified number of epochs. The corresponding hyperparameter configuration may change as models were trained longer; in the figure, we marked such changes with a dot. We found that a suitable hyperparameter configuration as well as a good model can be found using significantly less than the 400 epochs used in our study. Nevertheless, longer training did help in many settings and some models may benefit when trained for even more than 400 epochs.

4. CONCLUSION:

We reported the results of an extensive experimental study with popular KGE model architectures and training strategies across a wide range of hyperparameter settings. We found that when trained appropriately, the relative performance differences between various model architectures often shrunk and sometimes even reversed when compared to prior results. This suggests that (at least currently) training strategies have a significant impact on model performance and may account for a substantial fraction of the progress made in recent years. To enable an insightful model comparison, we recommend that new studies compare models in a common experimental setup using a sufficiently large hyperparameter search space. Our study shows that such an approach is possible with reasonable computational cost, and we provide our implementation as part of the LibKGE framework to facilitate such model comparison. Moreover, we feel that many of the more advanced architectures and techniques proposed in the literature need to be revisited to reassess their individual benefits. Validation MRR WNRR 1 2 8 2 5 6 5 1 2 1 2 8 2 5 6 5 1 2 1 2 8 2 5 6 5 1 2 1 2 8 2 5 6 5 1 2 Validation MRR WNRR 1 2 8 2 5 6 5 1 2 1 0 2 4 1 2 8 2 5 6 5 1 2 1 0 2 4 1 2 8 2 5 6 5 1 2 1 0 2 4 1 2 8 2 5 6 5 1 2 1 0 2 4 

:

https://github.com/uma-pi1/kge

Note that the KvsAll strategy is called 1-N scoring inDettmers et al. (2018).

Alternatively, we can view such an approach as only predicting objects, but doing so for both each original relation (k) and a new reciprocal relation formed by switching subject and object (k -1 ).

https://github.com/uma-pi1/kge-iclr20

Figures and Tables:

Figure 1: Distribution of filtered MRR (%) on validation data over the quasi-random hyperparameter configurations explored in our study.

Figure 2: Distribution of filtered MRR (%) on validation data over the quasi-random hyperparameter configurations for different training type and loss functions.

Figure 3: Best filtered MRR (%) on validation data achieved during quasi-random search as a function of the number of training epochs. Changes in the corresponding "best-so-far" hyperparameter configuration are marked with a dot.

Figure 4: Distribution of filtered MRR (on validation data, quasi-random search only) for different embedding sizes (top row: FB15K-237, bottom row: WNRR)

Figure 5: Distribution of filtered MRR (on validation data, quasi-random search only) for different batch sizes (top row: FB15K-237, bottom row: WNRR)

Figure 9: Distribution of filtered MRR (on validation data, quasi-random search only) for different penalty terms in the loss function (top row: FB15K-237, bottom row: WNRR)

Figure 10: Distribution of filtered MRR (on validation data, quasi-random search only) with and without dropout (top row: FB15K-237, bottom row: WNRR)

Model performance in prior studies and our study (as percentages, on test data).

Table content: <table><tr><td>.5</td><td>50.1</td><td>43.0</td><td>52.0</td></tr></table>

Hyperparameters of best performing models after quasi-random hyperparameter search and Bayesian optimization w.r.t. filtered MRR (on validation data). For each hyperparameter, we also give the reduction in filtered MRR for the best configuration that does not use this value of the hyperparameter (in parenthesis).

Table content: <table><tr><td/><td/><td/><td>RESCAL</td><td>TransE</td><td>DistMult</td><td>ComplEx</td><td>ConvE</td></tr><tr><td>FB15K-237</td><td>Validation MRR</td><td>10 20 30 40</td><td/><td/><td/><td/></tr><tr><td/><td/><td>0</td><td/><td/><td/><td/></tr><tr><td>WNRR</td><td>Validation MRR</td><td>10 20 30 40</td><td/><td/><td/><td/></tr><tr><td/><td/><td>0</td><td/><td/><td/><td/></tr><tr><td/><td colspan="2">NegSamp+MR</td><td>NegSamp+BCE</td><td>1vsAll+BCE</td><td>KvsAll+BCE</td><td>NegSamp+CE</td><td>1vsAll+CE</td><td>KvsAll+CE</td></tr></table>

Peer review (iclr2020):
The paper presents an experimental study about some KGE methods. It argues that papers often propose changes in several different dimensions, such as model, loss, training, regularizer, etc., at once without providing a sufficient investigation about the individual components* contributions. The experimental study considers two datasets (FB15k-237 and WNRR) and five different models (RESCAL, TransE, DistMult, ComplEx, ConvE). The models were selected using a quasi-random hyperparameter search, followed by a short Bayesian optimization phase to fine-tune the parameters. The performance of the best models found during this hyperparameter search are compared to first published results for the same model, as well as to a small selection of recent papers. To analyse the influence of single hyperparameters, the best found configuration is compared to the best configuration which does not use this specific value for the given hyperparameter. Overall, the paper adresses an important problem, as papers about new KGE methods often lack a clear separation of the individual changes* contribution. The experimental results show that older, simpler can compete with recently proposed models when trained properly. The intra-model comparison lacks statistical rigorousity, yet hints a few directions to further explore. The experiments are based on a quasi-random hyperparameter search. While it is necessary for efficient exploration of larger search spaces [1], and should be the standard methodology for hyperparameter search of a new method, the interpretability of the comparison of two runs suffers. For the comparison between the trained models and previously published results, the sample size might be sufficient to draw the conclusions. However, the intra-model comparison, e.g. in Figure 2, are now comparing subsets of the runs which only comprise approx. (200/6) runs. Furthermore, the influence of random initialization is not accounted for. Another place where this can be witnessed is Table 3. Here, for some ablations, e.g. TransE + Reciprocal, no reduction is given. If I understood it correctly, this is due to not having a configuration which uses TransE and reciprocal relations. Also for the other ablations, it is unclear how statistically significant the reduction is. Further Comments: 1. Please add the best published results for a specific model-dataset combination to table 2. 2. Do the plots in Figure 1 include the runs which were stopped after 50 epochs due to insufficient MRR? 3. Could you elaborate on the combination of KvsAll and CE? 4. The combination of subject and object triple scores has for instance been used in SimplE [2]. [1] Bergstra, James, and Yoshua Bengio. *Random search for hyper-parameter optimization.* Journal of Machine Learning Research 13.Feb (2012): 281-305. [2] Kazemi, Seyed Mehran, and David Poole. *Simple embedding for link prediction in knowledge graphs.* Advances in Neural Information Processing Systems. 2018.