R E: K G E R - R C S: Otat Nowledge Raph Mbedding by ELA Tional Otation in Omplex Pace
R E: K G E R - R C S: Otat Nowledge Raph Mbedding by ELA Tional Otation in Omplex Pace
[email protected]
[email protected]
A BSTRACT
1 I NTRODUCTION
Knowledge graphs are collections of factual triplets, where each triplet (h, r, t) represents a rela-
tion r between a head entity h and a tail entity t. Examples of real-world knowledge graphs include
Freebase (Bollacker et al., 2008), Yago (Suchanek et al., 2007), and WordNet (Miller, 1995). Knowl-
edge graphs are potentially useful to a variety of applications such as question-answering (Hao et al.,
2017), information retrieval (Xiong et al., 2017), recommender systems (Zhang et al., 2016), and
natural language processing (Yang & Mitchell, 2017). Research on knowledge graphs is attracting
growing interests in both academia and industry communities.
Since knowledge graphs are usually incomplete, a fundamental problem for knowledge graph is pre-
dicting the missing links. Recently, extensive studies have been done on learning low-dimensional
representations of entities and relations for missing link prediction (a.k.a., knowledge graph embed-
ding) (Bordes et al., 2013; Trouillon et al., 2016; Dettmers et al., 2017). These methods have been
shown to be scalable and effective. The general intuition of these methods is to model and infer the
connectivity patterns in knowledge graphs according to the observed knowledge facts. For example,
some relations are symmetric (e.g., marriage) while others are antisymmetric (e.g., filiation); some
relations are the inverse of other relations (e.g., hypernym and hyponym); and some relations may be
composed by others (e.g., my mother’s husband is my father). It is critical to find ways to model and
infer these patterns, i.e., symmetry/antisymmetry, inversion, and composition, from the observed
facts in order to predict missing links.
Indeed, many existing approaches have been trying to either implicitly or explicitly model one or
a few of the above relation patterns (Bordes et al., 2013; Wang et al., 2014; Lin et al., 2015b;
∗
This work was done when the first author was visiting Mila and Université de Montréal.
1
Published as a conference paper at ICLR 2019
Table 1: The score functions fr (h, t) of several knowledge graph embedding models, where h·i
denotes the generalized dot product, ◦ denotes the Hadamard product, ⊗ denotes circular correlation,
σ denotes activation function and ∗ denotes 2D convolution. · denotes conjugate for complex
vectors, and 2D reshaping for real vectors in ConvE model. TransX represents a wide range of
TransE’s variants, such as TransH (Wang et al., 2014), TransR (Lin et al., 2015b), and STransE
(Nguyen et al., 2016), where gr,i (·) denotes a matrix multiplication with respect to relation r.
Yang et al., 2014; Trouillon et al., 2016). For example, the TransE model (Bordes et al., 2011),
which represents relations as translations, aims to model the inversion and composition patterns; the
DisMult model (Yang et al., 2014), which models the three-way interactions between head entities,
relations, and tail entities, aims to model the symmetry pattern. However, none of existing models is
capable of modeling and inferring all the above patterns. Therefore, we are looking for an approach
that is able to model and infer all the three types of relation patterns.
In this paper, we propose such an approach called RotatE for knowledge graph embedding. Our
motivation is from Euler’s identity eiθ = cos θ + i sin θ, which indicates that a unitary complex
number can be regarded as a rotation in the complex plane. Specifically, the RotatE model maps
the entities and relations to the complex vector space and defines each relation as a rotation from
the source entity to the target entity. Given a triplet (h, r, t), we expect that t = h ◦ r, where
h, r, t ∈ Ck are the embeddings, the modulus |ri | = 1 and ◦ denotes the Hadamard (element-wise)
product. Specifically, for each dimension in the complex space, we expect that:
ti = hi ri , where hi , ri , ti ∈ C and |ri | = 1. (1)
It turns out that such a simple operation can effectively model all the three relation patterns: sym-
metric/antisymmetric, inversion, and composition. For example, a relation r is symmetric if and
only if each element of its embedding r, i.e. ri , satisfies ri = e0/iπ = ±1; two relations r1 and r2
are inverse if and only if their embeddings are conjugates: r2 = r̄1 ; a relation r3 = eiθ3 is a combi-
nation of other two relations r1 = eiθ1 and r2 = eiθ2 if and only if r3 = r1 ◦ r2 (i.e. θ3 = θ1 + θ2 ).
Moreover, the RotatE model is scalable to large knowledge graphs as it remains linear in both time
and memory.
To effectively optimizing the RotatE, we further propose a novel self-adversarial negative sampling
technique, which generates negative samples according to the current entity and relation embed-
dings. The proposed technique is very general and can be applied to many existing knowledge graph
embedding models. We evaluate the RotatE on four large knowledge graph benchmark datasets in-
cluding FB15k (Bordes et al., 2013), WN18 (Bordes et al., 2013), FB15k-237 (Toutanova & Chen,
2015) and WN18RR (Dettmers et al., 2017). Experimental results show that the RotatE model sig-
nificantly outperforms existing state-of-the-art approaches. In addition, RotatE also outperforms
state-of-the-art models on Countries (Bouchard et al., 2015), a benchmark explicitly designed for
composition pattern inference and modeling. To the best of our knowledge, RotatE is the first model
that achieves state-of-the-art performance on all the benchmarks.1
2 R ELATED W ORK
1
The codes of our paper are available online: https://github.com/DeepGraphLearning/
KnowledgeGraphEmbedding. pP
2
The p-norm of a complex vector v is defined as kvkp = p |vi |p . We use L1-norm for all distance-
based models in this paper and drop the subscript of k·k1 for brevity.
2
Published as a conference paper at ICLR 2019
Predicting missing links with knowledge graph embedding (KGE) methods has been extensively
investigated in recent years. The general methodology is to define a score function for the triplets.
Formally, let E denote the set of entities and R denote the set of relations, then a knowledge graph
is a collection of factual triplets (h, r, t), where h, t ∈ E and r ∈ R. Since entity embeddings are
usually represented as vectors, the score function usually takes the form fr (h, t), where h and t are
head and tail entity embeddings. The score function fr (h, t) measures the salience of a candidate
triplet (h, r, t). The goal of the optimization is usually to score true triplet (h, r, t) higher than the
corrupted false triplets (h0 , r, t) or (h, r, t0 ). Table 1 summarizes different score functions fr (h, t)
in previous state-of-the-art methods as well as the model proposed in this paper. These models gen-
erally capture only a portion of the relation patterns. For example, TransE represents each relation
as a bijection between source entities and target entities, and thus implicitly models inversion and
composition of relations, but it cannot model symmetric relations; ComplEx extends DistMult by
introducing complex embeddings so as to better model asymmetric relations, but it cannot infer the
composition pattern. The proposed RotatE model leverages the advantages of both.
A relevant and concurrent work to our work is the TorusE (Ebisu & Ichise, 2018) model, which
defines knowledge graph embedding as translations on a compact Lie group. The TorusE model
can be regarded as a special case of RotatE, where the modulus of embeddings are set fixed; our
RotatE is defined on the entire complex space, which has much more representation capacity. Our
experiments show that this is very critical for modeling and inferring the composition patterns.
Moreover, TorusE focuses on the problem of regularization in TransE while this paper focuses on
modeling and inferring multiple types of relation patterns.
There are also a large body of relational approaches for modeling the relational patterns on knowl-
edge graphs (Lao et al., 2011; Neelakantan et al., 2015; Das et al., 2016; Rocktäschel & Riedel,
2017; Yang et al., 2017). However, these approaches mainly focus on explicitly modeling the rela-
tional paths while our proposed RotatE model implicitly learns the relation patterns, which is not
only much more scalable but also provides meaningful embeddings for both entities and relations.
Another related problem is how to effectively draw negative samples for training knowledge graph
embeddings. This problem has been explicitly studied by Cai & Wang (2017), which proposed a
generative adversarial learning framework to draw negative samples. However, such a framework
requires simultaneously training the embedding model and a discrete negative sample generator,
which are difficult to optimize and also computationally expensive. We propose a self-adversarial
sampling scheme which only relies on the current model. It does require any additional optimization
component, which make it much more efficient.
The key of link prediction in knowledge graph is to infer the connection patterns, e.g., relation
patterns, with observed facts. According to the existing literature (Trouillon et al., 2016; Toutanova
& Chen, 2015; Guu et al., 2015; Lin et al., 2015a), three types of relation patterns are very important
3
Published as a conference paper at ICLR 2019
and widely spread in knowledge graphs: symmetry, inversion and composition. We give their formal
definition here:
Definition 1. A relation r is symmetric (antisymmetric) if ∀x, y
r(x, y) ⇒ r(y, x) ( r(x, y) ⇒ ¬r(y, x) )
A clause with such form is a symmetry (antisymmetry) pattern.
Definition 2. Relation r1 is inverse to relation r2 if ∀x, y
r2 (x, y) ⇒ r1 (y, x)
A clause with such form is a inversion pattern.
Definition 3. Relation r1 is composed of relation r2 and relation r3 if ∀x, y, z
r2 (x, y) ∧ r3 (y, z) ⇒ r1 (x, z)
A clause with such form is a composition pattern.
According to the definition of the above three types of relation patterns, we provide an analysis of
existing models on their abilities in inferring and modeling these patterns. Specifically, we provide
an analysis on TransE, TransX, DistMult, and ComplEx.3 We did not include the analysis on HolE
and ConvE since HolE is equivalent to ComplEx (Hayashi & Shimbo, 2017), and ConvE is a black
box that involves two-layer neural networks and convolution operations, which are hard to analyze.
The results are summarized into Table 2. We can see that no existing approaches are capable of
modeling all the three relation patterns.
In this part, we introduce our proposed model that is able to model and infer all the three types of
relation patterns. Inspired by Euler’s identity, we map the head and tail entities h, t to the complex
embeddings, i.e., h, t ∈ Ck ; then we define the functional mapping induced by each relation r as
an element-wise rotation from the head entity h to the tail entity t. In other words, given a triple
(h, r, t), we expect that:
t = h ◦ r, where |ri | = 1, (2)
and ◦ is the Hadmard (or element-wise) product. Specifically, for each element in the embeddings,
we have ti = hi ri . Here, we constrain the modulus of each element of r ∈ Ck , i.e., ri ∈ C,
to be |ri | = 1. By doing this, ri is of the form eiθr,i , which corresponds to a counterclockwise
rotation by θr,i radians about the origin of the complex plane, and only affects the phases of the
entity embeddings in the complex vector space. We refer to the proposed model as RotatE due to its
rotational nature. According to the above definition, for each triple (h, r, t), we define the distance
function of RotatE as:
dr (h, t) = kh ◦ r − tk (3)
By defining each relation as a rotation in the complex vector spaces, RotatE can model and infer all
the three types of relation patterns introduced above. Formally, we have following results4 :
3
See discussion at Appendix A
4
We relegate all proofs to the appendix.
4
Published as a conference paper at ICLR 2019
Table 3: Number of entities, relations, and observed triples in each split for four benchmarks.
Lemma 1. RotatE can infer the symmetry/antisymmetry pattern. (See proof in Appendix B)
Lemma 2. RotatE can infer the inversion pattern. (See proof in Appendix C)
Lemma 3. RotatE can infer the composition pattern. (See proof in Appendix D)
These results are also summarized into Table 2. We can see that the RotatE model is the only model
that can model and infer all the three types of relation patterns.
Connection to TransE. From Table 2, we can see that TransE is able to infer and model all the
other relation patterns except the symmetry pattern. The reason is that in TransE, any symmetric
relation will be represented by a 0 translation vector. As a result, this will push the entities with
symmetric relations to be close to each other in the embedding space. RotatE solves this problem
and is a able to model and infer the symmetry pattern. An arbitrary vector r that satisfies ri = ±1
can be used for representing a symmetric relation in RotatE, and thus the entities having symmetric
relations can be distinguished. Different symmetric relations can be also represented with different
embedding vectors. Figure 1 provides illustrations of TransE and RotatE with only 1-dimensional
embedding and shows how RotatE models a symmetric relation.
3.3 O PTIMIZATION
Negative sampling has been proved quite effective for both learning knowledge graph embedding
(Trouillon et al., 2016) and word embedding (Mikolov et al., 2013). Here we use a loss function
similar to the negative sampling loss (Mikolov et al., 2013) for effectively optimizing distance-based
models:
n
X 1
L = − log σ(γ − dr (h, t)) − log σ(dr (h0i , t0i ) − γ), (4)
i=1
k
where γ is a fixed margin, σ is the sigmoid function, and (h0i , r, t0i ) is the i-th negative triplet.
We also propose a new approach for drawing negative samples. The negative sampling loss samples
the negative triplets in a uniform way. Such a uniform negative sampling suffers the problem of
inefficiency since many samples are obviously false as training goes on, which does not provide
any meaningful information. Therefore, we propose an approach called self-adversarial negative
sampling, which samples negative triples according to the current embedding model. Specifically,
we sample negative triples from the following distribution:
exp αfr (h0j , t0j )
p(h0j , r, t0j |{(hi , ri , ti )}) = P 0 0 (5)
i exp αfr (hi , ti )
where α is the temperature of sampling. Moreover, since the sampling procedure may be costly,
we treat the above probability as the weight of the negative sample. Therefore, the final negative
sampling loss with self-adversarial training takes the following form:
X n
L = − log σ(γ − dr (h, t)) − p(h0i , r, t0i ) log σ(dr (h0i , t0i ) − γ) (6)
i=1
5
Published as a conference paper at ICLR 2019
4 E XPERIMENTS
We evaluate our proposed model on four widely used knowledge graphs. The statistics of these
knowledge graphs are summarized into Table 3.
• FB15k (Bordes et al., 2013) is a subset of Freebase (Bollacker et al., 2008), a large-scale
knowledge graph containing general knowledge facts. Toutanova & Chen (2015) showed
that almost 81% of the test triplets (x, r, y) can be inferred via a directly linked triplet
(x, r0 , y) or (y, r0 , x). Therefore, the key of link prediction on FB15k is to model and infer
the symmetry/antisymmetry and inversion patterns.
• WN18 (Bordes et al., 2013) is a subset of WordNet (Miller, 1995), a database featuring
lexical relations between words. This dataset also has many inverse relations. So the main
relation patterns in WN18 are also symmetry/antisymmetry and inversion.
• FB15k-237 (Toutanova & Chen, 2015) is a subset of FB15k, where inverse relations are
deleted. Therefore, the key of link prediction on FB15k-237 boils down to model and infer
the symmetry/antisymmetry and composition patterns.
• WN18RR (Dettmers et al., 2017) is a subset of WN18. The inverse relations are deleted,
and the main relation patterns are symmetry/antisymmetry and composition.
Hyperparameter Settings. We use Adam (Kingma & Ba, 2014) as the optimizer and fine-tune the
hyperparameters on the validation dataset. The ranges of the hyperparameters for the grid search are
set as follows: embedding dimension5 k ∈ {125, 250, 500, 1000}, batch size b ∈ {512, 1024, 2048},
self-adversarial sampling temperature α ∈ {0.5, 1.0}, and fixed margin γ ∈ {3, 6, 9, 12, 18, 24, 30}.
Both the real and imaginary parts of the entity embeddings are uniformly initialized, and the phases
of the relation embeddings are uniformly initialized between 0 and 2π. No regularization is used
since we find that the fixed margin γ could prevent our model from over-fitting.
Evaluation Settings. We evaluate the performance of link prediction in the filtered setting: we
rank test triples against all other candidate triples not appearing in the training, validation, or test
set, where candidates are generated by corrupting subjects or objects: (h0 , r, t) or (h, r, t0 ). Mean
Rank (MR), Mean Reciprocal Rank (MRR) and Hits at N (H@N) are standard evaluation measures
for these datasets and are evaluated in our experiments.
Baseline. Apart from RotatE, we propose a variant of RotatE as baseline, where the modulus of
the entity embeddings are also constrained: |hi | = |ti | = C, and the distance function is thus
2C sin θh +θ2r −θt (See Equation 17 at Appendix F for a detailed derivation). In this way, we can
investigate how RotatE works without modulus information and with only phase information. We
refer to the baseline as pRotatE. It is obvious to see that pRotatE can also model and infer all the
three relation patterns.
We compare RotatE to several state-of-the-art models, including TransE (Bordes et al., 2013), Dist-
Mult (Yang et al., 2014), ComplEx (Trouillon et al., 2016), HolE (Nickel et al., 2016), and ConvE
(Dettmers et al., 2017), as well as our baseline model pRotatE, to empirically show the importance
of modeling and inferring the relation patterns for the task of predicting missing links.
Table 4 summarizes our results on FB15k and WN18. We can see that RotatE outperforms all the
state-of-the-art models. The performance of pRotatE and RotatE are similar on these two datasets.
Table 5 summarizes our results on FB15k-237 and WN18RR, where the improvement is much
more significant. The difference between RotatE and pRotatE is much larger on FB15k-237 and
5
Following Trouillon et al. (2016), we treat complex number as the same as real number with regard to the
embedding dimension. If the same number of dimension is used for both the real and imaginary parts of the
complex number as the real number, the number of parameters for the complex embedding would be twice the
number of parameters for the embeddings in the real space.
6
Published as a conference paper at ICLR 2019
FB15k WN18
MR MRR H@1 H@3 H@10 MR MRR H@1 H@3 H@10
TransE [♥] - .463 .297 .578 .749 - .495 .113 .888 .943
DistMult [♦] 42 .798 - - .893 655 .797 - - .946
HolE - .524 .402 .613 .739 - .938 .930 .945 .949
ComplEx - .692 .599 .759 .840 - .941 .936 .945 .947
ConvE 51 .657 .558 .723 .831 374 .943 .935 .946 .956
pRotatE 43 .799 .750 .829 .884 254 .947 .942 .950 .957
RotatE 40 .797 .746 .830 .884 309 .949 .944 .952 .959
Table 4: Results of several models evaluated on the FB15K and WN18 datasets. Results of [♥] are
taken from (Nickel et al., 2016) and results of [♦] are taken from (Kadlec et al., 2017). Other results
are taken from the corresponding original papers.
FB15k-237 WN18RR
MR MRR H@1 H@3 H@10 MR MRR H@1 H@3 H@10
TransE [♥] 357 .294 - - .465 3384 .226 - - .501
DistMult 254 .241 .155 .263 .419 5110 .43 .39 .44 .49
ComplEx 339 .247 .158 .275 .428 5261 .44 .41 .46 .51
ConvE 244 .325 .237 .356 .501 4187 .43 .40 .44 .52
pRotatE 178 .328 .230 .365 .524 2923 .462 .417 .479 .552
RotatE 177 .338 .241 .375 .533 3340 .476 .428 .492 .571
Table 5: Results of several models evaluated on the FB15k-237 and WN18RR datasets. Results of
[♥] are taken from (Nguyen et al., 2017). Other results are taken from (Dettmers et al., 2017).
WN18RR, where there are a lot of composition patterns. This indicates that modulus is very impor-
tant for modeling and inferring the composition pattern.
Moreover, the performance of these models on different datasets is consistent with our analysis on
the three relation patterns (Table 2):
7
Published as a conference paper at ICLR 2019
Countries (AUC-PR)
DistMult ComplEx ConvE RotatE
S1 1.00 ± 0.00 0.97 ± 0.02 1.00 ± 0.00 1.00 ± 0.00
S2 0.72 ± 0.12 0.57 ± 0.10 0.99 ± 0.01 1.00 ± 0.00
S3 0.52 ± 0.07 0.43 ± 0.07 0.86 ± 0.05 0.95 ± 0.00
Table 6: Results on the Countries datasets. Other results are taken from (Dettmers et al., 2017).
150
20
200
100 15
10
100
50
5
0 0 0
0 π 1 π 2 π 0 π 1 π 2 π 0 π 1 π 2 π
30 15
15
300
20 10
10 200
10 5
5 100
0 0 0 0
0 π 1 π 2 π 0 π 1 π 2 π 0 π 1 π 2 π 0 π 1 π 2 π
We also evaluate our model on the Countries dataset (Bouchard et al., 2015; Nickel et al., 2016),
which is carefully designed to explicitly test the capabilities of the link prediction models for com-
position pattern modeling and inferring. It contains 2 relations and 272 entities (244 countries, 5
regions and 23 subregions). Unlike link prediction on general knowledge graphs, the queries in
Countries are of the form locatedIn(c, ?), and the answer is one of the five regions. The Coun-
tries dataset has 3 tasks, each requiring inferring a composition pattern with increasing length and
difficulty. For example, task S2 requires inferring a relatively simpler composition pattern:
neighborOf(c1 , c2 ) ∧ locatedIn(c2 , r) ⇒ locatedIn(c1 , r),
while task S3 requires inferring the most complex composition pattern:
neighborOf(c1 , c2 ) ∧ locatedIn(c2 , s) ∧ locatedIn(s, r) ⇒ locatedIn(c1 , r).
In Table 6, we report the results with respect to the AUC-PR metric, which is commonly used in the
literature. We can see that RotatE outperforms all the previous models. The performance of RotatE
is significantly better than other methods on S3, which is the most difficult task.
In this section, we verify whether the relation patterns are implicitly represented by RotatE relation
embeddings. We ignore the specific positions in the relation embedding θr and plot the histogram
of the phase of each element in the relation embedding, i.e., {θr,i }.
Symmetry pattern requires the symmetric relations to have property r ◦ r = 1, and the solution is
ri = ±1. We investigate the relation embeddings from a 500-dimensional RotatE trained on WN18.
Figure 2a gives the histogram of the embedding phases of a symmetric relation similar to. We can
find that the embedding phases are either π (ri = −1) or 0, 2π (ri = 1). It indicates that the RotatE
8
Published as a conference paper at ICLR 2019
Table 7: TransE with different negative sampling techniques. The results in first 2 rows are taken
from (Cai & Wang, 2017), where KBGAN uses a ComplEx negative sample generator.
Table 8: Results of TransE and ComplEx with self-adversarial sampling and negative sampling loss
on FB15k, FB15k-237 and Countries datasets.
model does infer and model the symmetry pattern. Figure 2b is the histogram of relation hypernym,
which shows that the embedding of a general relation does not have such a ±1 pattern.
Inversion pattern requires the embeddings of a pair of inverse relations to be conjugate. We use the
same RotatE model trained on WN18 for an analysis. Figure 2c illustrates the element-wise addition
of the embedding phases from relation r1 = hypernym and its inversed relation r2 = hyponym. All
the additive embedding phases are 0 or 2π, which represents that r1 = r−1 2 . This case shows that
the inversion pattern is also inferred and modeled in the RotatE model.
Composition pattern requires the embedding phases of the composed relation to be the addition of
the other two relations. Since there is no significant composition pattern in WN18, we study the
inference of the composition patterns on FB15k-237, where a 1000-dimensional RotatE is trained.
Figure 2d - 2g illustrate such a r1 = r2 ◦ r3 case, where θ2,i + θ3,i = θ1,i or θ2,i + θ3,i = θ1,i + 2π.
More results of implicitly inferring basic patterns are presented in the appendix.
In this part, we compare different negative sampling techniques including uniform sampling, our
proposed self-adversarial technique, and the KBGAN model (Cai & Wang, 2017), which aims to
optimize a generative adversarial network to generate the negative samples. We re-implement a 50-
dimension TransE model with the margin-based ranking criterion that was used in (Cai & Wang,
2017), and evaluate its performance on FB15k-237, WN18RR and WN18 with self-adversarial neg-
ative sampling. Table 7 summarizes our results. We can see that self-adversarial sampling is the
most effective negative sampling technique.
One may argue that the contribution of RotatE comes from the self-adversarial negative sampling
technique. In this part, we conduct further experiments on TransE and ComplEx in the same setting
as RotatE to make a fair comparison among the three models. Table 8 shows the results of TransE
and ComplEx trained with the self-adversarial negative sampling technique on FB15k and FB15k-
237 datasets, where a large number of relations are available. In addition, we evaluate these three
models on the Countries dataset, which explicitly requires inferring the composition pattern. We
also provide a detailed ablation study on TransE and RotatE in the appendix.
From Table 8, we can see that similar results are observed as Table 4 and 5. The RotatE model
achieves the best performance on both FB15k and FB15k-237, as it is able to model all the three
relation patterns. The TransE model does not work well on the FB15k datasets, which requires
modeling the symmetric pattern; the ComplEx model does not work well on FB15k-237, which
requires modeling the composition pattern. The results on the Countries dataset are a little bit
different, where the TransE model slightly outperforms RoateE on the S3 task. The reason is that
9
Published as a conference paper at ICLR 2019
Relation Category 1-to-1 1-to-N N-to-1 N-to-N 1-to-1 1-to-N N-to-1 N-to-N
Tasks Prediction Head (Hits@10) Prediction Tail (Hits@10)
TransE .437 .657 .182 .472 .437 .197 .667 .500
TransH (bern) .668 .876 .287 .645 .655 .398 .833 .672
KG2E KL (bern) .923 .946 .660 .696 .926 .679 .944 .734
TransE .894 .972 .567 .880 .879 .671 .964 .910
ComplEx .939 .969 .692 .893 .938 .823 .952 .910
RotatE .922 .967 .602 .893 .923 .713 .961 .922
Tasks Prediction Head (MRR) Prediction Tail (MRR)
TransE .701 .912 .424 .737 .701 .561 .894 .761
ComplEx .832 .914 .543 .787 .826 .661 .869 .800
RotatE .878 .934 .465 .803 .872 .611 .909 .832
Table 9: Experimental results on FB15k by relation category. The first three rows are taken from
(He et al., 2015). The rest of the results are from RotatE trained with the self-adversarial negative
sampling technique.
the Countries datasets do not have the symmetric relation between different regions, and all the
three tasks in the Countries datasets only require inferring the region for a given city. Therefore, the
TransE model does not suffer from its inability of modeling symmetric relations. For ComplEx, we
can see that it does not perform well on Countries since it cannot infer the composition pattern.
We also did some further investigation on the performance of RotatE on different relation categories:
one-to-many, many-to-one, and many-to-many relations6 . The results of RotatE on different relation
categories on the data set FB15k are summarized into Table 9. We also compare an additional
approach KG2E KL (He et al., 2015), which is a probabilistic framework for knowledge graph
embedding methods and aims to model the uncertainties of the entities and relations in knowledge
graphs with the TransE model. We also summarize the statistics of different relation categories into
Table 10 in the appendix.
We can see that besides the one-to-one relation, the RotatE model also performs quite well on the
non-injective relations, especially on many-to-many relations. We also notice that the probabilistic
framework KG2E KL(bern) (He et al., 2015) is quite powerful, which consistently outperforms
its corresponding knowledge graph embedding model, showing the importance of modeling the
uncertainties in knowledge graphs. We leave the work of modeling the uncertainties in knowledge
graphs with RotatE as our future work.
5 C ONCLUSION
We have proposed a new knowledge graph embedding method called RotatE, which represents enti-
ties as complex vectors and relations as rotations in complex vector space. In addition, we propose a
novel self-adversarial negative sampling technique for efficiently and effectively training the RotatE
model. Our experimental results show that the RotatE model outperforms all existing state-of-the-
art models on four large-scale benchmarks. Moreover, RotatE also achieves state-of-the-art results
on a benchmark that is explicitly designed for composition pattern inference and modeling. A deep
investigation into RotatE relation embeddings shows that the three relation patterns are implicitly
represented in the relation embeddings. In the future, we plan to evaluate the RotatE model on more
datasets and leverage a probabilistic framework to model the uncertainties of entities and relations.
6
Following Wang et al. (2014), for each relation r, we compute the average number of tails per head (tphr)
and the average number of head per tail (hptr). If tphr < 1.5 and hptr < 1.5, r is treated as one-to-one; if
tphr ≥ 1.5 and hptr ≥ 1.5, r is treated as a many-to-many; if tphr < 1.5 and hptr ≥ 1.5, r is treated as
one-to-many.
10
Published as a conference paper at ICLR 2019
R EFERENCES
Kurt Bollacker, Colin Evans, Praveen Paritosh, Tim Sturge, and Jamie Taylor. Freebase: a collab-
oratively created graph database for structuring human knowledge. In Proceedings of the 2008
ACM SIGMOD international conference on Management of data, pp. 1247–1250. AcM, 2008.
Antoine Bordes, Jason Weston, Ronan Collobert, Yoshua Bengio, et al. Learning structured embed-
dings of knowledge bases. In AAAI, volume 6, pp. 6, 2011.
Antoine Bordes, Nicolas Usunier, Alberto Garcia-Duran, Jason Weston, and Oksana Yakhnenko.
Translating embeddings for modeling multi-relational data. In Advances in neural information
processing systems, pp. 2787–2795, 2013.
Guillaume Bouchard, Sameer Singh, and Theo Trouillon. On approximate reasoning capabilities
of low-rank vector spaces. AAAI Spring Syposium on Knowledge Representation and Reasoning
(KRR): Integrating Symbolic and Neural Approaches, 2015.
Liwei Cai and William Yang Wang. Kbgan: Adversarial learning for knowledge graph embeddings.
arXiv preprint arXiv:1711.04071, 2017.
Rajarshi Das, Arvind Neelakantan, David Belanger, and Andrew McCallum. Chains of rea-
soning over entities, relations, and text using recurrent neural networks. arXiv preprint
arXiv:1607.01426, 2016.
Tim Dettmers, Pasquale Minervini, Pontus Stenetorp, and Sebastian Riedel. Convolutional 2d
knowledge graph embeddings. arXiv preprint arXiv:1707.01476, 2017.
Takuma Ebisu and Ryutaro Ichise. Toruse: Knowledge graph embedding on a lie group. In Pro-
ceedings of the Thirty-Second AAAI Conference on Artificial Intelligence, pp. 1819–1826. AAAI
Press, 2018.
Kelvin Guu, John Miller, and Percy Liang. Traversing knowledge graphs in vector space. arXiv
preprint arXiv:1506.01094, 2015.
Yanchao Hao, Yuanzhe Zhang, Kang Liu, Shizhu He, Zhanyi Liu, Hua Wu, and Jun Zhao. An end-
to-end model for question answering over knowledge base with cross-attention combining global
knowledge. In Proceedings of the 55th Annual Meeting of the Association for Computational
Linguistics (Volume 1: Long Papers), volume 1, pp. 221–231, 2017.
Katsuhiko Hayashi and Masashi Shimbo. On the equivalence of holographic and complex embed-
dings for link prediction. arXiv preprint arXiv:1702.05563, 2017.
Shizhu He, Kang Liu, Guoliang Ji, and Jun Zhao. Learning to represent knowledge graphs with
gaussian embedding. In Proceedings of the 24th ACM International on Conference on Information
and Knowledge Management, pp. 623–632. ACM, 2015.
Rudolf Kadlec, Ondrej Bajgar, and Jan Kleindienst. Knowledge base completion: Baselines strike
back. arXiv preprint arXiv:1705.10744, 2017.
Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. arXiv preprint
arXiv:1412.6980, 2014.
Ni Lao, Tom Mitchell, and William W Cohen. Random walk inference and learning in a large scale
knowledge base. In Proceedings of the Conference on Empirical Methods in Natural Language
Processing, pp. 529–539. Association for Computational Linguistics, 2011.
Yankai Lin, Zhiyuan Liu, Huanbo Luan, Maosong Sun, Siwei Rao, and Song Liu. Modeling relation
paths for representation learning of knowledge bases. arXiv preprint arXiv:1506.00379, 2015a.
Yankai Lin, Zhiyuan Liu, Maosong Sun, Yang Liu, and Xuan Zhu. Learning entity and relation
embeddings for knowledge graph completion. In AAAI, volume 15, pp. 2181–2187, 2015b.
Farzaneh Mahdisoltani, Joanna Biega, and Fabian M Suchanek. Yago3: A knowledge base from
multilingual wikipedias. In CIDR, 2013.
11
Published as a conference paper at ICLR 2019
Tomas Mikolov, Ilya Sutskever, Kai Chen, Greg S Corrado, and Jeff Dean. Distributed represen-
tations of words and phrases and their compositionality. In Advances in neural information pro-
cessing systems, pp. 3111–3119, 2013.
George A Miller. Wordnet: a lexical database for english. Communications of the ACM, 38(11):
39–41, 1995.
Arvind Neelakantan, Benjamin Roth, and Andrew McCallum. Compositional vector space models
for knowledge base completion. arXiv preprint arXiv:1504.06662, 2015.
Dai Quoc Nguyen, Tu Dinh Nguyen, Dat Quoc Nguyen, and Dinh Phung. A novel embedding
model for knowledge base completion based on convolutional neural network. arXiv preprint
arXiv:1712.02121, 2017.
Dat Quoc Nguyen, Kairit Sirts, Lizhen Qu, and Mark Johnson. Stranse: a novel embedding model
of entities and relationships in knowledge bases. arXiv preprint arXiv:1606.08140, 2016.
Maximilian Nickel, Lorenzo Rosasco, Tomaso A Poggio, et al. Holographic embeddings of knowl-
edge graphs. In AAAI, volume 2, pp. 3–2, 2016.
Tim Rocktäschel and Sebastian Riedel. End-to-end differentiable proving. In Advances in Neural
Information Processing Systems, pp. 3788–3800, 2017.
Fabian M Suchanek, Gjergji Kasneci, and Gerhard Weikum. Yago: a core of semantic knowledge.
In Proceedings of the 16th international conference on World Wide Web, pp. 697–706. ACM,
2007.
Kristina Toutanova and Danqi Chen. Observed versus latent features for knowledge base and text
inference. In Proceedings of the 3rd Workshop on Continuous Vector Space Models and their
Compositionality, pp. 57–66, 2015.
Théo Trouillon, Johannes Welbl, Sebastian Riedel, Éric Gaussier, and Guillaume Bouchard. Com-
plex embeddings for simple link prediction. In International Conference on Machine Learning,
pp. 2071–2080, 2016.
Zhen Wang, Jianwen Zhang, Jianlin Feng, and Zheng Chen. Knowledge graph embedding by trans-
lating on hyperplanes. In AAAI, volume 14, pp. 1112–1119, 2014.
Chenyan Xiong, Russell Power, and Jamie Callan. Explicit semantic ranking for academic search
via knowledge graph embedding. In Proceedings of the 26th international conference on world
wide web, pp. 1271–1279. International World Wide Web Conferences Steering Committee, 2017.
Bishan Yang and Tom Mitchell. Leveraging knowledge bases in lstms for improving machine read-
ing. In Proceedings of the 55th Annual Meeting of the Association for Computational Linguistics
(Volume 1: Long Papers), volume 1, pp. 1436–1446, 2017.
Bishan Yang, Wen-tau Yih, Xiaodong He, Jianfeng Gao, and Li Deng. Embedding entities and
relations for learning and inference in knowledge bases. arXiv preprint arXiv:1412.6575, 2014.
Fan Yang, Zhilin Yang, and William W Cohen. Differentiable learning of logical rules for knowledge
base completion. CoRR, abs/1702.08367, 2017.
Fuzheng Zhang, Nicholas Jing Yuan, Defu Lian, Xing Xie, and Wei-Ying Ma. Collaborative knowl-
edge base embedding for recommender systems. In Proceedings of the 22nd ACM SIGKDD
international conference on knowledge discovery and data mining, pp. 353–362. ACM, 2016.
12
Published as a conference paper at ICLR 2019
A PPENDIX
B P ROOF OF L EMMA 1
Proof. if r(x, y) and r(y, x) hold, we have
y =r◦x∧x=r◦y ⇒r◦r=1
Otherwise, if r(x, y) and ¬r(y, x) hold, we have
y = r ◦ x ∧ x 6= r ◦ y ⇒ r ◦ r 6= 1
C P ROOF OF L EMMA 2
Proof. if r1 (x, y) and r2 (y, x) hold, we have
y = r1 ◦ x ∧ x = r2 ◦ y ⇒ r1 = r−1
2
D P ROOF OF L EMMA 3
Proof. if r1 (x, z), r2 (x, y) and r3 (y, z) hold, we have
z = r1 ◦ x ∧ y = r2 ◦ x ∧ z = r3 ◦ y ⇒ r1 = r2 ◦ r3
E P ROPERTIES OF ROTAT E
A useful property for RotatE is that the inverse of a relation can be easily acquired by complex
conjugate. In this way, the RotatE model treats head and tail entities in a uniform way, which is
potentially useful for efficient 1-N scoring (Dettmers et al., 2017):
kh ◦ r − tk = k(h ◦ r − t) ◦ rk = kt ◦ r − hk (7)
Moreover, considering the embeddings in the polar form, i.e., hi = mh,i eiθh,i , ri = eiθr,i , ti =
mt,i eiθt,i , we can rewrite the RotatE distance function as:
k r
X θh,i + θr,i − θt,i
kh ◦ r − tk = (mh,i − mt,i )2 + 4mh,i mt,i sin2 (8)
i=1
2
13
Published as a conference paper at ICLR 2019
YAGO3-10
MR MRR H@1 H@3 H@10
DistMult 5926 .34 .24 .38 .54
ComplEx 6351 .36 .26 .40 .55
ConvE 1671 .44 .35 .49 .62
RotatE 1767 .495 .402 .550 .670
Table 11: Results of several models evaluated on the YAGO3-10 datasets. Other results are taken
from (Dettmers et al., 2017).
F P ROOF OF T HEOREM 4
YAGO3-10 is a subset of YAGO3 (Mahdisoltani et al., 2013), which consists of entities that have a
minimum of 10 relations each. It has 123,182 entities and 37 relations. Most of the triples deal with
descriptive attributes of people, such as citizenship, gender, profession and marital status.
Table 11 shows that the RotatE model also outperforms state-of-the-art models on YAGO3-10.
H H YPERPARAMETERS
We list the best hyperparameter setting of RotatE w.r.t the validation dataset on several benchmarks
in Table 12.
14
Published as a conference paper at ICLR 2019
Table 13: Results of ablation study on FB15k-237, where “adv” represents “self-adversarial”.
I A BLATION S TUDY
Table 13 shows our ablation study of self-adversarial sampling and negative sampling loss on
FB15k-237. We also re-implement a 1000-dimension TransE and do ablation study on it. From
the table, We can find that self-adversarial sampling boosts the performance for both models, while
negative sampling loss is only effective on RotatE; in addition, our re-implementation of TransE
also outperforms all the state-of-the-art models on FB15k-237.
15
Published as a conference paper at ICLR 2019
Table 14: The average and variance of the MRR results of RotatE on FB15k, WN18, FB15k-237
and WN18RR.
20
10
15
10
5
0 0
0 π 1 π 2 π 0 π 1 π 2 π
(a) has part (b) instance hypernym
80
100
60
40
50
20
0 0
0 π 1 π 2 π 0 π 1 π 2 π
(c) also see (d) verb group
150
100
100
50
50
0 0
0 π 1 π 2 π 0 π 1 π 2 π
(e) derivationally related form (f) similar to
Figure 3: Histograms of embedding phases from two general relations and four symmetric relations
on WN18. ( k = 500 )
16
Published as a conference paper at ICLR 2019
300
300
200 200
100 100
0 0
0 π 1 π 2 π 0 π 1 π 2 π
(a) /celebrities/celebrity- (b) /award/award winner-
/celebrity friends./celebrities- /awards won./award/award honor-
/friendship/friend /award winner
300 300
200 200
100 100
0 0
0 π 1 π 2 π 0 π 1 π 2 π
(c) /location/statistical region- (d) /base/popstra/celebrity-
/places exported to./location- /breakup./base/popstra-
/imports and exports/exported to /breakup/participant
300
400
300
200
200
100
100
0 0
0 π 1 π 2 π 0 π 1 π 2 π
(e) /base/popstra/celebrity- (f) /government/legislative session-
/dated./base/popstra- /members./government-
/dated/participant /government position held/legislative sessions
Figure 4: Histograms of embedding phases from six symmetric relations on FB15k-237. (k = 1000)
17
Published as a conference paper at ICLR 2019
200
200
100
100
0
0 π 1 π 2 π 0
0 π 1 π 2 π
(a) member of domain topic
(b) has part ◦ part of
◦ synset domain topic of
150
150
100
100
50
50
0 0
0 π 1 π 2 π 0 π 1 π 2 π
(c) synset domain usage of (d) synset domain region of
◦ member of domain usage ◦ member of domain region
200
100
100 50
0 0
0 π 1 π 2 π 0 π 1 π 2 π
(e) member meronym ◦ member holonym (f) instance hypernym ◦ instance hyponym
Figure 5: Histograms of element-wise additions of inversed relation embedding phases on WN18.
(k = 500)
18