Title: RELWALK -A LATENT VARIABLE MODEL APPROACH TO KNOWLEDGE GRAPH EMBEDDING Anonymous authors Paper under double-blind review

Abstract: Knowledge Graph Embedding (KGE) is the task of jointly learning entity and relation embeddings for a given knowledge graph. Existing methods for learning KGEs can be seen as a two-stage process where (a) entities and relations in the knowledge graph are represented using some linear algebraic structures (embeddings), and (b) a scoring function is defined that evaluates the strength of a relation that holds between two entities using the corresponding relation and entity embeddings. Unfortunately, prior proposals for the scoring functions in the first step have been heuristically motivated, and it is unclear as to how the scoring functions in KGEs relate to the generation process of the underlying knowledge graph. To address this issue, we propose a generative account of the KGE learning task. Specifically, given a knowledge graph represented by a set of relational triples (h, R, t), where the semantic relation R holds between the two entities h (head) and t (tail), we extend the random walk model (Arora et al., 2016a) of word embeddings to KGE. We derive a theoretical relationship between the joint probability p(h, R, t) and the embeddings of h, R and t. Moreover, we show that marginal loss minimisation, a popular objective used by much prior work in KGE, follows naturally from the log-likelihood ratio maximisation under the probabilities estimated from the KGEs according to our theoretical relationship. We propose a learning objective motivated by the theoretical analysis to learn KGEs from a given knowledge graph. The KGEs learnt by our proposed method obtain state-of-the-art performance on FB15K237 and WN18RR benchmark datasets, providing empirical evidence in support of the theory.

1. INTRODUCTION:

Knowledge graphs such as Freebase (Bollacker et al., 2008) organise information in the form of graphs, where entities are represented by vertices in the graph and the relation between two entities is represented by the edge that connects the corresponding two vertices. By embedding entities and relations that exist in a knowledge graph in some (possibly lower-dimensional and latent) space we can infer previously unseen relations between entities, thereby expanding a given knowledge graph (Nickel et al., 2016; Yang et al., 2015; Lin et al., 2015; Nickel et al., 2011; Trouillon et al., 2016; Wang et al., 2017; Bordes et al., 2011) .

Existing KGE methods can be seen as involving two main steps. First, given a knowledge graph represented by a set of relational triples (h, R, t), where a semantic relation R holds between a head entity h and a tail entity t, entities and relations are represented using some mathematical structures such as vectors, matrices or tensors. Second, a scoring function is proposed that evaluates the relational strength of a triple (h, R, t) and entity and relation embeddings that optimise the defined scoring function are learnt using some optimisation method. Table 1 shows some of the scoring functions proposed in prior work in KGE learning.

Despite the wide applications of entity and relation embeddings created via KGE methods, the existing scoring functions are motivated heuristically to capture some geometric requirements of the embedding space. For example, TransE (Bordes et al., 2011) assumes that the entity and relation embeddings co-exist in the same (possibly lower dimensional) vector space and translating (shifting) the head entity embedding by the relation embedding must make it closer to the tail entity embedding, whereas ComplEx (Trouillon et al., 2016) models the asymmetry in relations using

Model:

Score function f (h, R, t) Relation parameters Unstructured (Bordes et al., 2011) ht 1/2 none Structured embeddings (Bordes et al., 2011 ) (Bordes et al., 2011) h + Rt 1/2 R ∈ R d DistMult (Yang et al., 2015) h, R, t R ∈ R d RESCAL (Nickel et al., 2011) h Rt R d×d ComplEx (Trouillon et al., 2016) h, R, t R ∈ C d Table 1 : Score functions proposed in selected prior work on KGE. Entity embeddings h, t ∈ R d are vectors in all models, except in ComplEx where h, t ∈ C d . Here, x 1/2 denotes either 1 or 2 norm of the vector x. In ComplEx, x is the elementwise complex conjugate, and •, •, • denotes the component-wise multi-linear inner-product.

R 1 h -R 2 t 1,2 R 1 , R 2 ∈ R d×d TransE

the component-wise multi-linear inner-product among entity and relation embeddings. Relational triples extracted from a given knowledge graph are used as positive training instances, whereas pseudo-negative (Bordes et al., 2011) instances are automatically generated by randomly corrupting positive instances. Finally, KGE are learnt such that the prediction loss computed over the positive and negative instances is minimised.

Despite the good empirical performances of the existing KGE methods, theoretical understanding of KGE methods is comparatively under developed. For example, it is not clear how the heuristically defined KGE objectives relate to the generative process of a knowledge graph. In this paper, we attempt to fill this void by providing a theoretical analysis of KGE. Specifically, in section 2, we propose a generative process where we explain the formation of a relation R between two entities h and t using the corresponding relation and entity embeddings. Following this generative story, we derive a relationship between the probability of R holding between h and t, p(h, t | R), and the embeddings of R, h and t. Interestingly, the derived relationship is not covered by any of the previously proposed heuristically-motivated scoring functions, providing the first-ever KGE method with a provable generative explanation.

Next, in section 3, we show that the margin loss, which has been popularly used as a training objective in prior work on KGE, naturally arises as the log-likelihood ratio computed from p(h, t | R). Based on this result, we derive a training objective that we subsequently optimise for learning KGEs that satisfy our theoretical relationship. Using standard benchmark datasets proposed in prior work on KGE learning, we evaluate the learnt KGEs on a link prediction task and a triple classification task. Experimental results show that the learnt KGEs obtain state-of-the-art performance on FB15K237 and WN18RR benchmarks, thereby providing empirical evidence to support the theoretical analysis.

2. RELATIONAL WALK:

Let us consider a knowledge graph D where the knowledge is represented by relational triples (h, R, t) ∈ D. Here, R is a relational predicate of two arguments, where h (head) and t (tail) entities respectively filling the first and second arguments. We assume relations to be asymmetric in general.

In other words, if (h, R, t) ∈ D then it does not necessarily follow that (t, R, h) ∈ D. The goal of KGE is to learn embeddings (representations) for the relations and entities in the knowledge graph such that the entities that participate in similar relations are embedded closely to each other in the entity embedding space, while at the same time relations that hold between similar entities are embedded closely to each other in the relational embedding space. We call the learnt entity and relation embeddings collectively as KGEs. Following prior work on KGE (Bordes et al., 2011; Trouillon et al., 2016; Yang et al., 2015) , we assume that entities and relations are embedded in the same vector space, allowing us to perform linear algebraic operations using the embeddings in the same vector space.

Let us consider a random walk characterised by a time-dependent knowledge vector c k , where k is the current time step. The knowledge vector represents the knowledge we have about a particular group of entities and relations that express some facts about the world. For example, the knowledge that we have about people that are employed by companies can be expressed using entities of classes such as people and organisation, using relations such as CEO-of, employed-at, works-for, etc. We assume that entities h and t are represented by time-independent d-dimensional vectors, respectively h, t ∈ R d .

We assume the task of generating a relational triple (h, R, t) in a given knowledge graph to be a two-step process as described next. First, given the current knowledge vector at time k, c = c k and the relation R, we assume that the probability of an entity h satisfying the first argument of R to be given by (1).

EQUATION

Here, R 1 ∈ R d×d is a relation-specific orthogonal matrix that evaluates the appropriateness of h for the first argument of R. For example, if R is the CEO-of relation, we would require a person as the first argument and a company as the second argument of R. However, note that the role of R 1 extends beyond simply checking the types of the entities that can fill the first argument of a relation. For our example above, not all people are CEOs and R 1 evaluates the likelihood of a person to be selected as the first argument of the CEO-of relation. Z c is a normalisation coefficient such that

h∈V p(h | R, c) = 1,

where the vocabulary V is the set of all entities in the knowledge graph.1 After generating h, the state of our random walker changes to c = c k+1 , and we next generate the second argument of R with the probability given by (2).

EQUATION

Here, R 2 ∈ R d×d is a relation-specific orthogonal matrix that evaluates the appropriateness of t as the second argument of R. Z c is a normalisation coefficient such that t∈V p(t | R, c) = 1. Following our previous example of the CEO-of relation, R 2 evaluates the likelihood of an organisation to be a company with a CEO position. Importantly, R 1 and R 2 are representations of the relation R and independent of the entities. Therefore, we consider (R 1 and R 2 ) to collectively represent the embedding of R. Orthogonality of R 1 , R 2 is a requirement for the mathematical proof and also act as a regularisation constraint to prevent overfitting by restricting the relational embedding space. We first perform our mathematical analysis for relational embeddings represented by orthogonal matrices and discuss later how this requirement can be relaxed.

We assume a slow random walk where the knowledge vectors do not change significantly between consecutive time steps (c k ≈ c k+1 ). More specifically, we assume that c k -c k+1 ≤ 2 for some small 2 > 0. This is a realistic assumption for generating the two entity arguments in the same relational triple because, if the knowledge vectors were significantly different in the two generation steps, then it is likely that the corresponding relations are also different, which would not be coherent with the above-described generative process. Moreover, we assume that the knowledge vectors are distributed uniformly in the unit sphere and denote the distribution of knowledge vectors by C.

To learn KGEs, we must estimate the probability that h and t satisfy the relation R, p(h, t | R), which can be obtained by taking the expectation of

EQUATION

Here, partition functions are given by 4) follows from our two-step generative process where the generation of h and t in each step is independent given the relation and the corresponding knowledge vectors.

Z c = h∈V c∈C exp h R 1 c and Z c = t∈V c ∈C exp t R 2 c . (

Computing the expectation in ( 5) is generally difficult because of the two partition functions Z c and Z c . However, Lemma 1 shows that the partition functions are narrowly distributed around a constant value for all c (or c ) values with high probability.

Lemma 1 (Concentration Lemma). If the entity embedding vectors satisfy the Bayesian prior v = sv, where v is from the spherical Gaussian distribution, and s is a scalar random variable, which is always bounded by a constant κ, then the entire ensemble of entity embeddings satisfies that

EQUATION

for z = O(1/ √ n), and δ = exp(-Ω(log 2 n)),

where n ≥ d is the number of words and Z c is the partition function for c given by c∈V exp h R 1 c . proof: To prove the concentration lemma, we show that the mean E h [Z c ] of Z c is concentrated around a constant for all knowledge vectors c and its variance is bounded. Recall that

EQUATION

)

If P is an orthogonal matrix and x is a vector, then P x 2 2 = (P x) (P x) = x PP x = x 2 2 , because P P = I. Therefore, from ( 7) and the orthogonality of the relational embeddings, we see that R 1 c is a simple rotation of c and does not alter the length of c. We represent h = s h ĥ, where s h = h and ĥ is a unit vector (i.e. ĥ 2 = 1) distributed on the spherical Gaussian with zero mean and unit covariance matrix I d ∈ R d×d . Let s be a random variable that has the same distribution as s h . Moreover, let us assume that s is upper bounded by a constant κ such that s ≤ κ. From the assumption of the knowledge vector c, it is on the unit sphere as well, which is then rotated by R 1 .

We can write the partition function using the inner-product between two vectors h and R 1 c, Z c = h∈V exp h (R 1 c) . Arora et al. (2016a) showed that (Lemma 2.1 in their paper) the expectation of a partition function of this form can be approximated as follows:

EQUATION

EQUATION

where n = |V| is the number of entities in the vocabulary. ( 8) follows from the expectation of a sum and the independence of h and R 1 from c. The inequality of ( 9) is obtained by applying the Taylor expansion of the exponential series and the final equality is due to the symmetry of the spherical Gaussian. From the law of total expectation, we can write

EQUATION

)

where, x = h R 1 c. Note that conditioned on s h , h is a Gaussian random variable with variance σ 2 = s 2 h . Therefore, conditioned on s h , x is a random variable with variance σ 2 = σ 2 h . Using this distribution, we can evaluate E x|s h exp h R 1 c as follows:

EQUATION

Therefore, it follows that

EQUATION

)

where s is the variance of the 2 norms of the entity embeddings. Because the set of entities is given and fixed, both n and σ are constants, proving that E[Z c ] does not depend on c.

Next, we calculate the variance V c [Z c ] as follows:

EQUATION

)

Because 2h R 1 t is a Gaussian random variable with variance 4σ 2 = 4s 2 h from a similar calculation as in (11) we obtain,

E x|s h exp 2h R 1 t | s h = exp(2σ 2 ).

(16) By substituting ( 16) in ( 15) we have that

V c [Z c ] ≤ nE s h exp 2σ 2 = nE s h exp(2s 2 ) ≤ Λn (17)

for Λ = exp(8κ 2 ) a constant bounding s ≤ κ as stated.

From above, we have bounded both the mean and variance of the partition function by constants that are independent of the knowledge vector. Note that neither exp h R 1 c nor exp t R 2 c are sub-Gaussian nor sub-exponential. Therefore, standard concentration bounds derived for sub-Gaussian or sub-exponential random variables cannot be used in our analysis. However, the argument given in Appendix A.1 in Arora et al. (2016b) for a partition function with bounded mean and variance can be directly applied to Z c in our case, which completes the proof of the concentration lemma.

From the symmetry between h and t, Lemma 1 also applies for the partition function t∈V t R 2 c . Under the conditions required to satisfy Lemma 1, the following main theorem of this paper holds: Theorem 1. Suppose that the entity embeddings satisfy (1). Then, we have

EQUATION

)

for = O(1/ √ n) + O(1/d),

where

EQUATION

The complete proof of Theorem 1 is given in Appendix A. Below we briefly sketch the main steps.

Proof sketch: Let F be the event that both c and c are within (1 ± z )Z. Then, from Lemma 1 and the union bound, event F happens with probability at least 1 -2 exp(-Ω(log 2 n)). The R.H.S. of (5) can be split into two parts T 1 and T 2 according to whether F happens or not.

EQUATION

)

T 1 can be approximated as given by ( 21).

T 1 = 1 ± O( z ) Z 2 E c,c exp h R 1 c exp t R 2 c (21)

On the other hand, T 2 can be shown to be a constant, independent of d, given by ( 22).

EQUATION

The vocabulary size n of real-world knowledge graphs is typically over 10 5 , for which T 2 becomes negligibly small. Therefore, it suffices to consider only T 1 . Because of the slowness of the random walk we have c ≈ c

Using the law of total expectation we can write T 1 as follows:

EQUATION

where A(c) := E c |c exp t R 2 c . Doing some further evaluations we show that The relationship given by ( 18) indicates that head and tail entity embeddings are first transformed respectively by R 1 and R 2 , and the squared 2 norm of the sum of the transformed vectors is proportional to the probability p(h, t | R).

EQUATION

3. LEARNING KNOWLEDGE GRAPH EMBEDDINGS:

In this section, we derive a training objective from Theorem 1 that we can then optimise to learn KGE. The goal is to empirically validate the theoretical result by evaluating the learnt KGEs. Knowledge graphs represent information about relations between two entities in the form of relational triples. The joint probability p(h, R, t) given by Theorem 1 is useful for determining whether a relation R exists between two given entities h and t. For example, if we know that with a high probability that R holds between h and t, then we can append (h, R, t) to the knowledge graph. The task of expanding knowledge graphs by predicting missing links between entities or relations is known as the link prediction problem (Trouillon et al., 2016) . In particular, if we can automatically append such previously unknown knowledge to the knowledge graph, we can expand the knowledge graph and address the knowledge acquisition bottleneck.

To derive a criteria for determining whether a link must be predicted among entities and relations, let us consider a relational triple (h, R, t) ∈ D that exists in a given knowledge graph D. We call such relational triples as positive triples because from the assumption it is known that R holds between h and t. On the other hand, consider a negative relational triple (h , R, t ) ∈ D formed by, for example, randomly perturbing a positive triple. A popular technique for generating such (pseudo) negative triples is to replace h or t with a randomly selected different instance of the same entity type. As an alternative for random perturbation, Cai and Wang (2018) proposed a method for generating negative instances using adversarial learning. Here, we are not concerned about the actual method used for generating the negative triples but assume a set of negative triples, D, generated using some method, to be given.

Given a positive triple (h, R, t) ∈ D and a negative triple (h , R, t ) ∈ D, we would like to learn KGEs such that a higher probability is assigned to (h, R, t) than that assigned to (h , R, t ). We can formalise this requirement using the likelihood ratio given by (25).

EQUATION

Here, η > 1 is a threshold that determines how higher we would like to set the probabilities for the positive triples compares to that of the negative triples.

By taking the logarithm of both sides in (25) we obtain

EQUATION

If a positive triple (h, R, t) is correctly assigned a higher probability than a negative triple p(h , R, t ), then the left hand side of (26) will be negative, indicating that there is no loss incurred during this classification task. Therefore, we can re-write (26) to obtain the marginal loss Bordes et al. (2013; 2011) , L(D, D), a popular choice as a learning objective in prior work in KGE, as shown in (27) .

L(D, D) = (h,R,t)∈D (h ,R,t )∈ D max (0, log η + log p(h , R, t ) -log p(h, R, t)) = max 0, 2d log η + R 1 h + R 2 t 2 2 -R 1 h + R 2 t 2 2 (27)

We can assume 2d log η to be the margin for the constraint violation.

Theorem 1 requires R 1 and R 2 to be orthogonal. To reflect this requirement, we add two 2 regularisation terms R 1 R 1 -I 2 2 and R 2 R 2 -I 2 2 respectively with regularisation coefficients λ 1 and λ 2 to the objective function given by ( 27). In our experiments, we compute the gradients (27) w.r.t. each of the parameters h, t, R 1 and R 2 and use stochastic gradient descent (SGD) for optimisation. This approach can be easily extended to learn from multiple negative triples as shown in Appendix B.

4. RELATED WORK:

At a high-level of abstraction, KGE methods can be seen as differing in their design choices for the following two main problems: (a) how to represent entities and relations, and (b) how to model the interaction between two entities and a relation that holds between them. Next, we briefly discuss prior proposals to those two problems (refer (Wang et al., 2017; Nickel et al., 2015; Nguyen, 2017) for an extended survey on KGE).

A popular choice for representing entities is to use vectors, whereas relations have been represented by vectors, matrices or tensors. For example, TransE (Bordes et al., 2011) , TransH (Wang et al., 2014) , TransD (Ji et al., 2015) , TransG (Xiao et al., 2016) , TransR (Lin et al., 2015) , lppTransD (Yoon et al., 2016) , DistMult (Yang et al., 2015) , HolE (Nickel et al., 2016) and ComplEx (Trouillon et al., 2016) represent relations by vectors, whereas Structured Embeddings (Bordes et al., 2011) , TranSparse (Ji et al., 2016 ), STransE (Nguyen et al., 2016) , RESCAL (Nickel et al., 2011) use matrices and Neural Tensor Network (NTN) (Socher et al., 2013) uses 3D tensors. ComplEx (Trouillon et al., 2016) introduced complex vectors for KGEs to capture the asymmetry in semantic relations. (Ding et al., 2018) obtained state-of-the-art performance for KGE by imposing non-negativity and entailment constraints to ComplEx.

Given entity and relation embeddings, a scoring function is defined that evaluates the strength of a relation R between two entities h and t in a triple (h, R, t). The scoring functions that encode various intuitions have been proposed such as the 1 or 2 norms of the vector formed by a translation of the head entity embedding by the relation embedding over the target embedding, or by first performing a projection from the entity embedding space to the relation embedding space (Yoon et al., 2016) As an alternative to using vector norms as scoring functions, DistMult and ComplEx use the component-wise multi-linear dot product.

Once a scoring function is defined, KGEs are learnt that assign better scores to relational triples in existing knowledge graphs (positive triples) over triples where the relation does not hold (negative triples) by minimising a loss function such as the logistic loss (RESCAL, DistMult, ComplEx) or marginal loss (TransE, TransH, TransD, TransD) . Because knowledge graphs record only positive triples, a popular method to generate pseudo negative triples is to perturb a positive instance by replacing its head or tail entity by an entity selected uniformly at random from the vocabulary of the entities. However, uniformly sampled negative triples are likely to be obvious examples that do not provide much information to the learning process and can be detected by simply checking for the type of the entities in a triple. Cai and Wang (2018) proposed an adversarial learning approach where a generator assigns a probability to each relation triple and negative instances are sampled according to this probability distribution to train a discriminator that discriminates between positive and negative instances. (Xiao et al., 2016) proposed TransG, a generative model based on the Chinese restaurant process, to model multiple relations that exist between a pair of entities. However, their relation embeddings are designed to satisfy vector translation similar to TransE.

As an alternative to directly learning embeddings from a graph, several methods (Grover and Leskovec, 2016; Perozzi et al., 2014; Ristoski et al., 2018) have considered the vertices visited during truncated random walks over the graph as pseudo sentences, and have applied popular word embedding learning algorithms such as skip-gram with negative sampling or continuous bag-of-words model (Mikolov et al., 2013) to learn vertex embeddings. However, pseudo sentences generated this way are syntactically very different from sentences in natural languages.

On the other hand, our work extends the random walk analysis by Arora et al. (2016a) that derives a useful connection between the joint co-occurrence probability of two words and the 2 norm of the sum of the corresponding word embeddings. Specifically, they proposed a latent variable model where the words in a corpus are generated by a probabilistic model parametrised by a time-dependent discourse vector that performs a random walk. However, unlike in our work, they do not consider the relations between two co-occurring words in a corpus. Bollegala et al. (2018) extended the model proposed by Arora et al. (2016a) to capture co-occurrences involving more than two words. They defined the co-occurrence of k unique words in a given context as a k-way co-occurrence, where Arora et al. (2016a) 's result could be seen as a special case corersponding to k = 2. Moreover, Bollegala et al. (2018) showed that it is possible to learn word embeddings that capture some types of semantic relations such as antonymy and collocation using 3-way co-occurrences more accurately than using 2-way co-occurrences. However, their model does not explicitly consider the relations between words/entities and uses only a corpus for learning the word embeddings. Cai and Wang (2018) . All other results for the baselines are taken from their original papers.

FB15K237 WN18RR:

Method MRR MR H@1 H@3 H@10 MRR MR H@1 H@3 H@10 

5. EMPIRICAL VALIDATION:

To empirically evaluate the theoretical result stated in Theorem 1, we learn KGEs (denoted by RelWalk) by minimising the marginal loss objective derived in section 3. We use the FB15k237, FB13 (subsets of Freebase) and WN11, WN18RR (subsets of WordNet) datasets, which are standard benchmarks for KGE. We use the standard training, validation and test splits as detailed in Table 4 . We generate negative triples by replacing a head or a tail entity in a positive triple by a randomly selected different entity and learn KGEs. We train the model until convergence or at most 1000 epochs over the training data where each epoch is divided into 100 mini-batches. The best model is selected by early stopping based on the performance of the learnt embeddings on the validation set (evaluated after each 20 epochs). The training details and hyperparameter settings are detailed in Appendix C. RelWalk is implemented in the open-source toolkit OpenKE (Han et al., 2018) .2 

We conduct two evaluation tasks: link prediction (predict the missing head or tail entity in a given triple (h, R, ?) or (?, R, t)) (Bordes et al., 2011) and triple classification (predict whether a relation R holds between h and t in a given triple (h, R, t)) (Socher et al., 2013) . We evaluate the performance in the link prediction task using mean reciprocal rank (MRR), mean rank (MR (the average of the rank assigned to the original head or tail entity in a corrupted triple) and hits at ranks 1, 3 and 10 (H@1,3,10), whereas in the triple classification task we use accuracy (percentage of the correctly classified test triples). We only report scores under the filtered setting Bordes et al. (2013) , which removes all triples appeared in training, validating and testing sets from candidate triples before obtaining the rank of the ground truth triple. In link prediction, we consider all entities that appear in the corresponding argument in the entire knowledge graph as candidates.

In Tables 2 and 3 we compare the KGEs learnt by RelWalk against prior work using the published results. For link prediction, RelWalk reports SoTA on both WN18RR and FB15K237 in all evaluation measures, except against ConvE in WN18RR measured by MRR. WN18RR excludes triples from WN18 that are simply inverted between train and test partitions (Toutanova and Chen, 2015; Dettmers et al., 2017) . RelWalk's consistently good performance on both versions of this dataset shows that it is considering the global structure in the knowledge graph when learning KGEs. For triple classification, RelWalk reports the best performance on FB13, whereas TransG reports the best performance on WN11. Considering that both TransG and RelWalk are generative models, it would be interesting to further investigate generative approaches for KGE in the future. Overall, the experimental results support our theoretical claim and emphasise the importance of theoretically motivating the scoring function design process.

APPENDIX A PROOF OF THEOREM 1:

Let us consider the probabilistic event that (1

-z )Z ≤ Z c ≤ (1 + z )Z to be F c and (1 -z )Z ≤ Z c ≤ (1 + z )Z to be F c . From Lemma 1 we have Pr c [F c ] ≥ 1 -δ.

Then from the union bound we have,

EQUATION

Moreover, let F be the probabilistic event that both F c and F c being True. Then from

Pr[F ] = 1 -Pr[ Fc ∨ Fc ] we have, Pr[F ] ≥ 1 -2δ.

We can decompose the expectation in the R.H.S. in (5) into two terms T 1 and T 2 depending on whether respectively F is True or False as follows:

EQUATION

Here, 1 F and 1 F are indicator functions given by:

EQUATION

and

EQUATION

)

Let us first show that T 2 is negligibly small.

For two real integrable functions ψ 1 (x) and ψ 2 (x) in [a, b], the Cauchy-Schwarz's inequality states that b a

EQUATION

Applying (32) to T 2 in (29) we have:

EQUATION

Note that Z c ≥ 1 because Z c is the sum of positive numbers and if h R 1 c ≥ 0 for at least one of the h ∈ V, then the total sum will be greater than 1. Therefore, by dropping Z c term from the denominator we can further increase the first term in (33) as given by ( 34).

EQUATION

Let us split the expectation on the R.H.S. of (34) into two cases depending on whether h R 1 c > 0 or otherwise, indicated respectively by 1 (h R1c>0) and 1 (h R1c≤0) .

E c exp h R 1 c 2 E c |c [1 F ] = E c exp h R 1 c 2 1 (h R1c>0) E c |c [1 F ] + E c exp h R 1 c 2 1 (h R1c≤0) E c |c [1 F ] (35)

The second term of ( 35) is upper bounded by

E c,c [1 F ] ≤ exp -Ω(log 2 n) (36)

The first term of (35) can be bounded as follows:

EQUATION

where α > 1. Therefore, it is sufficient to bound

E c exp(αh R 1 c) 2 E c |c [1 F ] when h = Ω( √ d).

Let us denote by z the random variable 2h R 1 c. Moreover, let r 

(z) = E c |z [1 F ],

EQUATION

where

t satisfies that E c [1 [t,+∞] (z)] = Pr[z ≥ t] = E c [r(z)] ≤ exp(-Ω(log 2 n)).

Here, 1 [t,+∞] (z) is a function that takes the value 1 when z ≥ t and zero elsewhere. Then, we claim Pr c [z ≥ t] ≤ exp(-Ω(log 2 n)) implies that t ≥ Ω(log .9 n).

If c was distributed as N (0, 1 d I), this would be a simple tail bound. However, as c is distributed uniformly on the sphere, this requires special care, and the claim follows by applying the tail bound for the spherical distribution given by Lemma A.1 in (Arora et al., 2016a) instead. Finally, applying Corollary A.3 in (Arora et al., 2016a) , we have:

EQUATION

From a similar argument as above we can obtain the same bound for c as well. Therefore, T 2 in (29) can be upper bounded as follows:

EQUATION

Because n = |V|, the size of the entity vocabulary, is large (ca. n > 10 5 ) in most knowledge graphs, we can ignore the T 2 term in (29). Combining this with (29) we obtain an upper bound for p(h, t | R) given by (41).

EQUATION

)

where |D| is the number of relational tuples (h, R, t) in the KB D and δ 0 = |D| exp(-Ω(log 1.8 n)) ≤ exp(-Ω(log 1.8 n)) by the fact that Z ≤ exp(2κ)n = O(n), where κ is the upper bound on h R 1 c and t R 2 c , which is regarded as a constant.

On the other hand, we can lower bound p(h, t | R) as given by (42).

p(h, t | R) ≥ (1 -z ) 2 1 Z 2 E c,c exp h R 1 c exp t R 2 c 1 F ≥ (1 -z ) 2 1 Z 2 E c,c exp h R 1 c exp t R 2 c -|D| exp(-Ω(log 1.8 n)) ≥ (1 -z ) 2 1 Z 2 E c,c exp h R 1 c exp t R 2 c -δ 0 (42)

Taking the logarithm of both sides, from ( 41) and ( 42), the multiplicative error translates to an additive error given by (43).

EQUATION

where A(c) := E c |c exp t R 2 c .

We assumed that c and c are on the unit sphere and R 1 and R 2 to be orthogonal matrices. Therefore, R 1 c and R 2 c are also on the unit sphere. Moreover, if we let the upper bound of the 2 norm of the entity embeddings to be κ √ d, then we have h ≤ κ √ d and t ≤ κ √ d. Therefore, we have

EQUATION

Then we can lower bound A(c) as follows:

EQUATION

)

For some 2 > 0. The last inequality holds because To obtain a lower bound on A(c) from the first-order Taylor approximation of exp(x) ≥ 1 + x we observe that

EQUATION

)

Therefore, from our model assumptions we have

EQUATION

Hence,

EQUATION

Therefore, from ( 46) and (49) we have

EQUATION

)

Plugging A(c) back in (43) we obtain 

log p(h, t | R) = log E c exp h R 1 c A(c) ± δ 0 -2 log Z + 2 log(1 ± z ) = log E c exp h R 1 c (1 ± 2 ) exp t R 2 c ± δ 0 -2 log Z + 2 log(1 ± z ) (51) = log E c exp h R 1 c exp t R 2 c ± δ 0 -2 log Z + 2 log(1 ± z ) + log(1 ± 2 ) = log E c exp h R 1 c + t R 2 c ± δ 0 -2 log Z + 2 log(1 ± z ) + log(1 ± 2 ) = log E c exp R 1 h + R 2 t c ± δ 0 -2 log Z + 2 log(1 ± z ) + log(1 ± 2 )

C TRAINING DETAILS:

The statistics of the benchmark datasets are show in Table 4 .

We selected the initial learning rate (α) for SGD in {0.01, 0.001}, the regularisation coefficients (λ 1 , λ 2 ) for the orthogonality constraints of relation matrices in {0, 1, 10, 100}. The number of randomly generated negative triples n neg for each positive example is varied in {1, 10, 20, 50, 100} and d ∈ {50, 100}. Optimal hyperparameter settings were: λ 1 = λ 2 = 10, n neg = 100 for all the datasets, α = 0.001 for FB15K, FB15K237 and FB13, α = 0.01 for WN18, WN18RR and WN11. For FB15K237 and WN18RR d = 100 was the best, whereas for all other datasets d = 50 performed best.

:

We can consider different vocabularies for the entities that can fill the first argument and second argument of relations in a knowledge graph. However, for simplicity, we use a common vocabulary here.

To facilitate the double blind policy, the source code for RelWalk will be released upon paper acceptance

CONCLUSIONWe proposed RelWalk, a generative model of KGE and derived a theoretical relationship between the probability of a triple and entity, relation embeddings. We then proposed a learning objective based on the theoretical relationship we derived. Experimental results on a link prediction and a triple classification tasks show that RelWalk obtains strong performances in multiple benchmark datasets.

Figures and Tables:

24)Plugging (50) back in (23) provides the claim of the theorem.

Triple classification.

Table content: <table><tr><td/><td colspan="2">Accuracy</td></tr><tr><td>Method</td><td colspan="2">WN11 FB13</td></tr><tr><td>SE</td><td>53.0</td><td>75.2</td></tr><tr><td>TransE</td><td>75.9</td><td>81.5</td></tr><tr><td>TransR</td><td>85.9</td><td>82.5</td></tr><tr><td>TransG</td><td>87.4</td><td>87.3</td></tr><tr><td>NTN</td><td>70.4</td><td>87.1</td></tr><tr><td colspan="2">RelWalk 75.48</td><td>87.5</td></tr></table>

Link prediction. Results marked with [ ] are taken from Dettmers et al. (2017), [•] from Nguyen et al. (2017), [ ] from and

Table content: <table/>

which is a function of z between [0, 1]. We wish to upper bound E c [exp(z)r(z)]. The worst-case r(z) can be quantified using a continuous version of Abel's inequality (proved as Lemma A.4 in Arora et al. (2016b)), we can upper bound E c [exp(z)r(z)] as follows:

Table content: <table/>

Statistics of the datasets

Table content: <table><tr><td>Dataset</td><td colspan="3">Relations Entities Train</td><td>Test</td><td>Validation</td></tr><tr><td>FB15K</td><td>1,345</td><td>14,951</td><td colspan="3">483,142 59,071 50,000</td></tr><tr><td colspan="2">FB15K237 237</td><td>14,541</td><td colspan="3">272,115 17,535 20,466</td></tr><tr><td>WN18</td><td>18</td><td>40,943</td><td colspan="2">141,442 5,000</td><td>5,000</td></tr><tr><td>WN18RR</td><td>11</td><td>40,943</td><td>86,835</td><td>3,134</td><td>3,034</td></tr><tr><td>WN11</td><td>11</td><td>38,588</td><td colspan="3">112,581 10,544 2,609</td></tr><tr><td>FB13</td><td>13</td><td>75,043</td><td colspan="3">316,232 23,733 5,908</td></tr></table>

Peer review (iclr2019):
 This paper proposes to perform the link prediction in knowledge bases by introducing a new scoring function and theoretically motivating their method. The authors validate their proposed approach through several experiments. This paper reads well and the results appear sound. I personally find the theoretical argument behind the proposed scoring function very interesting. Unfortunately, the contribution seems rather small to be accepted for ICLR. This is a straight application and combination of existing pieces with not much originality and without being backed up by very strong experimental results. My concerns are as follows: - Having only results on two flawed datasets (considering the inverse relations in them) makes it hard to evaluate the quality of the method. I suggest conducting the experiments on the FB15K-237 and WN18RR from [1] instead. - You only evaluate on MR and Hits@10, but it is standard to include metrics like MRR and Hits@1 and 3 also, since no metric is perfect for this task. - Since the goal of the work is not providing a state of the art method, and focus on the theoretical understanding of their scoring function, it is of high importance to assess the characteristics of their embeddings and scoring function through designing other experiments. As a result, I suggest to study the geometric behavior of their embeddings and compare it to the other methods. Further, investigating the semantic purity of the embeddings by calculating the entropy of the type distribution of the entities, similar as [2], can shed more light on the significance of their method. On overall, although the proposed method seems a direct application of Arora et al.,2016a, I find their extension novel and quite interesting, But the paper needs more experimental results to validate the idea. [1] Dettmers, Tim, et al. "Convolutional 2d knowledge graph embeddings.", AAAI-18. [2] Ding, Boyang, et al. "Improving Knowledge Graph Embedding Using Simple Constraints.", ACL-18.