# 知識グラフ埋め込みの PU 学習

垣渕太成 $* 1$ 林克彦 $\dagger 2$ 駒谷和範 $\$ 1$
1 大阪大学 産業科学研究所
2 群馬大学 社会情報学部

## 1 はじめに

知識グラフは，事実（Fact）を集積した巨大な データセットであり, Fact は, Entity $e$ とその間の Relation $r$ の三つ組 (Triple) 形式 $\left(e_{i}, r_{k}, e_{j}\right)$ で表現さ れる.ここで前者の Entity は Subject，後者は Object と呼ばれる。知識グラフの多くは Fact の欠損を含 み，これを自動予測する知識グラフ補完は実用上重要なタスクである。そして，知識グラフ埋め込みは このタスクに対する有効なアプローチの一つと考え られている。埋め込みの学習は Logistic 損失や Cross Entropy（CE）損失に基づいて，Fact とそれ以外の Triple を分類するように行うのが一般的である.

知識グラフでは Fact 以外の Triple の正負が不明で あるため，知識グラフ埋め込みの学習を行うには負例をどう定義するかが重要となる. 既存の学習法で は，Factの Subject あるいは Object を別の Entity に 入れ替えた Tripleを負例として定義する方法がよく 用いられている. しかし, この負例の定義が正確で ある保証はなく，潜在的な Factを負例として扱って しまう危険が常に伴う．また，Factの近傍のみを学習対象とするため，局所的で偏った学習を行う問題 も指摘できる．分野ではモデルの性能を Fact 近傍の ランキングで評価することが主流であるため，これ らの問題は顕在化してこなかったが，本来，知識グ ラフ補完の性能を評価するには与えられた Triple が Fact か否かを判定する分類評価が適切である。

そこで，本稿ではまず，従来の学習法によるモデ ルの分類性能を再調査する. 図 1 は Kinship データ セットに対する既存モデルの性能を PR（Precision and Recall）曲線で示しており，この結果は, Precision が早期に下がって，分類性能が極めて低いことを示 している。これは負例が正例と判定される割合が高 いことに起因するもので, この結果から, 従来の学


図 1 Kinship データセットの Fact から 3 割を除いて, CE 損失による既存の埋め込み学習（RESCAL モデルを使用） を行い，除いた Fact に対して PR 曲線を計算した。

習法における負例の定義を根本的に見直す必要性が あると結論付けられる。

よって, 本稿では, PU (Positive and Unlabeled)学習の知見を用いた知識グラフ埋め込みの新たな 学習法を提案する. PU 学習では, 正例以外全て Unlabeled であるデータセットに一部正例が含まれ ていると仮定した上で正と負への二値分類を行う ため，明示的な負例を定義する必要がない．本稿 では知識グラフ埋め込みの代表的なモデルである RESCAL を用い，Alternating Least Square（ALS）に よる最適化法を定式化した. これを既存手法と比較 したところ，分類精度の大きな向上が見られた.

## 2 関連研究

## 2.1 知識グラフ埋め込みの学習

初期には二乗誤差損失を用いて ALS で最適化す る手法 [1] が提案されたが，これには Unlabeled な データを全て負例として扱うという欠点がある.現在, 事前知識 [2] や敵対学習 [3] によって Fact 近傍で負例の可能性が高い Triple をサンプリングし, Logistic 損失や $\mathrm{CE}$ 損失に基づいて学習を行う方法 が一般的となっている。しかし，任意の知識グラフ

に対して正確な負例を定義する方法は自明でない.

## 2.2 PU 学習

PU 学習 [4] は, 正例と Unlabeled で構成される データセットに対して, Unlabeled に一部正例が含ま れている状況を仮定し，正負の二値分類を行うもの で，土地被覆分類 [5] など，負例の抽出が難しい状況で利用されている. 文献 [6] は正負のラベルの誤 りを考慮した上での学習に取り組んでいる。正例に 誤りはないが，負例のデータに正クラスが誤って混在している状況は PU 学習の仮定と等しい. 文献 [7] は，行列分解における PU 学習に取り組んでおり,本研究はそのテンソル分解への拡張とも言える.

## 3 知識グラフ埋め込みの PU 学習

## 3.1 記法

ベクトル，行列，テンソルをそれぞれ小文字の太文字 $\boldsymbol{a}$, 大文字の太文字 $\boldsymbol{A}$, カリグラフィック体 $\mathcal{A}$ で表す. 行列 $\boldsymbol{A}$ の $\mathrm{i}$ 番目の行ベクトルを $\boldsymbol{a}_{i}^{T}$, 行列 の各要素を $a_{i j}$, テンソルの各要素を $x_{i j k}$ で表す. $X_{[i::]}, X_{[: j:]}, X_{[:: k]}$ はそれぞれテンソル $1,2,3$ 次元の $i, j, k$ 番目のスライスを表す. また $\operatorname{vec}(\boldsymbol{A})$ は行列 $\boldsymbol{A}$ のベクトル表示 $\left(a_{11}, \ldots, a_{1 n}, \ldots, a_{m 1}, \ldots, a_{m n}\right)^{T}$ であ る. 行列のフロベニウスノルムを $\|\cdot\|_{F}, \otimes$ はクロ ネッカー積, $x_{n}$ はテンソルの $n$ モード積を表す.

## 3.2 従来法

知識グラフ埋め込みでは Fact の分布を近似する 潜在的な特徴表現の獲得を目指す、以下では, Fact の集合を $\mathcal{F}$, Entity の集合を $\mathcal{E}$ ，Relation の集合を $\mathcal{R}$ で表す.このとき，知識グラフは二値のテンソル $X \in\{0,1\}^{|\mathcal{E}| \times|\mathcal{R}| \times|\mathcal{E}|}$ の形で表現でき, 各要素 $x_{i j k}$ は $\left(e_{i}, r_{k}, e_{j}\right) \in \mathcal{F}$ なら 1 , そうでなければ 0 をる.

RESCAL モデルでは Entity $e_{i}$ を $d$ 次元の縦ベクト ル $\boldsymbol{e}_{i}$, Relation $r_{k}$ を $d \times d$ 次元の行列 $\boldsymbol{R}_{k} \in \mathbb{R}^{d \times d}$ で 表現し，隣接行列 $X_{[:: k]}$ を以下のように近似する.

$
\chi_{[\because k]} \approx \boldsymbol{E} \boldsymbol{R}_{k} \boldsymbol{E}^{T}(k=1,2, \ldots,|\mathcal{R}|) .
$

Triple に対するスコア関数は以下となる.

$
\theta_{i j k}=\boldsymbol{e}_{i}^{T} \boldsymbol{R}_{k} \boldsymbol{e}_{j}
$

他にも様々なモデルが提案されているが，RESCAL は表現力が完全であり [8], 現在でも標準的なモデ ルとして使用されている [9].一般的な知識グラフ埋め込みでは，スコア関数に 対する損失を $\ell\left(\theta_{i j k}, x_{i j k}\right)$ として目的関数は以下の ように定義される.

$L=\sum_{\left(e_{i}, r_{k}, e_{j}\right) \in \mathcal{F}} \ell\left(\theta_{i j k}, 1\right)+\sum_{\left(e_{i}, r_{k}, e_{j}\right) \in \mathcal{D}^{-}} \ell\left(\theta_{i j k}, 0\right)$

where

$\mathcal{D}^{-}=\left.\{\left(e_{i^{\prime}}, r_{k}, e_{j}\right) \mid e_{i^{\prime}} \in \mathcal{E} \wedge e_{i^{\prime}} \neq e_{i} \wedge\left(e_{i}, r_{k}, e_{j}\right) \in \mathcal{F}\right.\}$

$\cup\left.\{\left(e_{i}, r_{k}, e_{j^{\prime}}\right) \mid e_{j^{\prime}} \in \mathcal{E} \wedge e_{j^{\prime}} \neq e_{j} \wedge\left(e_{i}, r_{k}, e_{j}\right) \in \mathcal{F}\right.\}$.

$\mathcal{D}^{-}$は学習で負に分類される Triple の集合であ る.これは正例以外の定義可能な Triple 全て を含むわけではなく, $\operatorname{Fact}\left(e_{i}, r_{k}, e_{j}\right)$ に対して, $\left(e_{i}, r_{k}, ?\right),\left(?, r_{k}, e_{j}\right)$ となるものだけが含まれる.

## 3.3 知識グラフ埋め込みの PU 定式化

$\mathrm{PU}$ 学習 $[6,7]$ の知見を用いて, 知識グラフ埋め込 みの新たな学習法を定式化する．知識グラフに本来存在する Fact数から $\rho$ の割合でランダムに欠損が起 こると仮定する。 すなわち $\rho$ は，真の Fact の集合を $\hat{\mathcal{F}}$, 現在の集積されている Fact の集合を $\mathcal{F}$ として $\rho=1-|\mathcal{F}| /|\hat{\mathcal{F}}|$ と定義される. $\hat{\mathcal{F}}$ に基づいて損失関数 $\ell$ で分類した場合の損失を近似するような損失関数 $\tilde{\ell}$ 定義し, 目的関数は以下のように修正できる.

$
\begin{aligned}
& L=\sum_{\left(e_{i}, r_{k}, e_{j}\right) \in \mathcal{F}} \tilde{\ell}\left(\theta_{i j k}, 1\right)+\sum_{\left(e_{i}, r_{k}, e_{j}\right) \notin \mathcal{F}} \tilde{\ell}\left(\theta_{i j k}, 0\right) \\
& \text { where }\left.\{\begin{array}{c}
\tilde{\ell}\left(\theta_{i j k}, 1\right)=\frac{\ell\left(\theta_{i j k}, 1\right)-\rho \ell\left(\theta_{i j k}, 0\right)}{1-\rho} \\
\tilde{\ell}\left(\theta_{i j k}, 0\right)=\ell\left(\theta_{i j k}, 0\right)
\end{array}\right.
\end{aligned}
$

知識グラフ $\hat{\mathcal{F}}$ に対応するテンソルが $\hat{X} ， そ の$ 要素を $\hat{x}_{i j k}$ とすると, 損失について以下の等式が成り立つ.

$
\mathbb{E}\left[\tilde{\ell}\left(\theta_{i j k}, x_{i j k}\right)\right]=\mathbb{E}\left[\ell\left(\theta_{i j k}, \hat{x}_{i j k}\right)\right]
$

知識グラフにおいて $\rho$ の値は不明だが， ハイパーパ ラメータとして与えた $\rho$ が実際に Fact の欠損してい る割合と等しい時，欠損を含まない知識グラフを学習した場合の損失に近似できる。

## 3.4 ALS による最適化

本研究では, 損失 $\ell$ に二乗誤差 $\left(x_{i j k}-\theta_{i j k}\right)^{2}$ を採用し，ALSによる最適化を行う。このとき，

$
\tilde{\ell}\left(\theta_{i j k}, 1\right)=\left(\theta_{i j k}-\frac{1}{1-\rho}\right)^{2}-\frac{\rho}{(1-\rho)^{2}}
$

と書ける。第二項目は定数項であるから，

$\operatorname{argmin} L=\operatorname{argmin}\left[\sum_{x_{i j k}=1}\left(\theta_{i j k}-\frac{1}{1-\rho}\right)^{2}+\sum_{x_{i j k}=0} \theta_{i j k}^{2}\right]$

図 2 AUC-PR と $\rho$ の関係 (Kinship).

が成立する. よって，要素が 0 と $1 /(1-\rho)$ のテンソ ル分解の問題に帰着する. 以下では，この二値のテ ンソル $X^{\prime}$ を $X^{\prime}:=1 /(1-\rho) X$ と定義する.

ALS は行列・テンソル分解において使用されるア ルゴリズムであり，分解する行列の片方を固定した 上でもう一方の行列をコスト関数の偏微分が 0 にな るよう更新し，これを交互に繰り返す.

本研究では RESCAL モデルを使用した. 正則化項を加えたコスト関数と ALS 更新式を以下に示す.

$
L=\frac{1}{2} \sum_{k}\left.\|X_{[:: k]}^{\prime}-\boldsymbol{E} \boldsymbol{R}_{k} \boldsymbol{E}^{T}\right.\|_{F}^{2}+\lambda\left(\|\boldsymbol{E}\|_{F}^{2}+\sum_{k}\left.\|\boldsymbol{R}_{k}\right.\|_{F}^{2}\right)
$

更新式

$
\begin{aligned}
& \boldsymbol{E} \leftarrow\left[\sum_{k} \mathcal{X}_{[:: k]}^{\prime} \boldsymbol{E}_{k}+\mathcal{X}_{[\because: k]}^{\prime T} \boldsymbol{E}_{k}^{T}\right]\left[\sum_{k} \boldsymbol{B}_{k}+\boldsymbol{C}_{k}+\lambda \boldsymbol{I}\right]^{-1} \\
& \quad \text { where } \boldsymbol{B}_{k}=\boldsymbol{R}_{k} \boldsymbol{E}^{T} \boldsymbol{E} \boldsymbol{R}_{k}^{T}, \boldsymbol{C}_{k}=\boldsymbol{R}_{k}^{T} \boldsymbol{E}^{T} \boldsymbol{E} \boldsymbol{R}_{k} \\
& \boldsymbol{R}_{k} \leftarrow\left(\boldsymbol{Z}^{T} \boldsymbol{Z}+\lambda \boldsymbol{I}\right)^{-1} \boldsymbol{Z}_{\operatorname{vec}}\left(\mathcal{X}_{[\because k]}^{\prime}\right)
\end{aligned}
$

where $\boldsymbol{Z}=\boldsymbol{E} \otimes \boldsymbol{E}$.

ALS を用いる利点として，大規模な知識グラフへの 対応が挙げられる. 文献 [10] では，知識グラフのス パース性を利用し 300 万の Entity，40の Relation，お よび 7000 万の Fact で構成される知識グラフの埋め 込みに成功している。

## 4 実験

## 4.1 分類精度の評価

3 節で提案した PU 学習に基づく定式化を実験的に検証する. 使用したデータセットは Kinship と

図 3 AUC-PR と $\rho$ の関係 (UMLS).

UMLS で，共に Fact が全て集積されておりそれ以外 の Triple は全て負例である性質を持つ [11].

データセットから Fact をランダムに一定割合 $\{0.1,0.3,0.5,0.7\}$ 除き残りを学習に与え，それ以外 の全 Triple に対し分類精度をテストする．検証の目的は以下の三つである.

1. PU 学習により分類精度は向上するか.
2. $\rho$ と除いた正例の割合の関係. 両者が等しい時 に分類精度が最も高くなるか.

3. 従来法との比較.

目的 1,2 の検証のため，ハイパーパラメータ $\rho$ ごと に AUC-PR の值を調べた（図 2,3）。ここで，次元数 100 , 正則化項 $\lambda=0.0001$ で固定し 10 分割交差検証 を行った. Factを除いた割合によらず， $\rho>0.0$ の適当な值で $\rho=0.0$ （PUなし）と比較して分類精度が 向上している。一方, 理論的には除いた正例の割合 に対し $\rho$ が等しい場合（図赤点）で AUC が最も高く なることが予想されるが，そのような傾向は得られ なかった。一般的な分類問題と異なり知識グラフ埋 め込みでは潜在的な特徴表現を通じて各 Triple のス コアが相互に影響を及ぼすため，理論通りの学習が 行われてない可能性があり，知識グラフ埋め込みに より適した PU 学習の定式化は今後の課題である.

続いて従来法と提案手法を比較する。従来法にはRESCAL モデル，目的関数に式 (1)，損失関数に Sigmoid Cross Entoropyを用いて SGD （Stochastic Gradient Descent）で最適化を行った。 ハイパーパラメータは次元数 $d \in\{30,50,100\}$,学習率は $\eta \in\{0.005,0.003,0.001\}$ から Grid Search を行った. Epochは 50 で固定したが，モデル
表 1 AUC-PR の手法ごとの比較 (Kinship).

表 2 AUC-PRの手法ごとの比較（UMLS）.

は，学習データ以外の一部 Factを対象として Epoch ごとにランキングによる評価を行い，最 も精度が高いものを使用した. 実装はhttps: //github.com/ibalazevic/TuckERにRESCAL を加えた1). 実験結果を表 1,2 に示す. どの設定で も提案手法が従来法より優れた分類精度を示した。

## 4.2 学習範囲による影響の分析

表 1,2 では従来法と比較して PUなしの ALS $(\rho=0.0)$ も高い分類精度を示している. 全 Tripleを 学習する ALS に対し, 従来法は Fact と負例の一部の み学習している. この違いが分類に与える影響を調査するため, 従来法が学習において負に分類する範囲を Negative Area，学習を行わない範囲を Unlearned Area として，それぞれの分類精度を AUC-PR で評価 した (図 4) ${ }^{2)}$. 二つの領域における正例と負例の比率を合わせるため, Unlearned Area で負例をランダ ムにサンプリングしテストした．従来法では，限ら れた範囲の学習によって Triple 全体における Fact の 補完が可能であるとされてきたが，Unlearned Area ではほとんど分類ができていない. ALS での分類精度はこれと比較して高い。これは ALS が Unlearned Area の学習を行っているので直感的だが, それで も値が高くないのは, これらの Triple は近傍に正例 を持たないため学習が困難だからだと考えられる。 また Negative Area でも ALS の方が分類精度が高く, Unlearned Area の学習が特徴表現を通じ Fact の補完 に重要な情報を与えていると考えられる。


図 4 目的関数式 (1) が負例として学習する Triple (Negative Area）と学習をしない Triple (Unlearned Area) を分け, AUC-PR で評価. 左は(1) に基づく従来法, RESCAL-ALS $(\rho=0.0)$ による学習.

図 5 Factを 0.8 除いた Kinship における正例負例のスコ アの分布. PU なし $\rho=0.0$ (左) と PU あり $\rho=0.3$ (右).

## 4.3 PU 学習が与える影響の分析

PU 学習なしの ALS に対して PU ありの ALS が高 い分類精度を示しており，特に正例を除いた割合が 高い時その差が大きい. PU の影響を調査するため,正例を 0.8 の割合除いて，学習に与えてない正例と 負例のスコアの分布を $\rho=0.0$ と $\rho=0.3$ で比較した (図 5 $)^{3)} \cdot \rho=0.0$ では Fact 以外を全て負に分類する ためスコアが非常に低い値に密集しており正例と負例のスコアに差を持たせることができておらず分類精度が低い（AUC-PR=0.075），一方 $\rho=0.3$ ではスコ アの密集が防がれ正例と負例のスコアに差が生まれ ており，分類精度が向上している（AUC-PR=0.378）。

## 4.4 実験まとめ

従来法は学習範囲が局所的という問題があり，分類精度が非常に低いのに対し，ALS はこの問題を回避している。また PU 学習を用いることで ALS によ る埋め込みの分類精度が向上した. 一方 PU 学習の 理論に適合しない結果もあり，より知識グラフに適 した PU 定式化は今後の課題である.

3）引き続き次元数 100 , 正則化項 $\lambda=0.0001$ を使用.

## 参考文献

[1] Maximilian Nickel, Volker Tresp, and Hans-Peter Kriegel. A three-way model for collective learning on multirelational data. In Proceedings of the 28th International Conference on International Conference on Machine Learning, ICML'11, p. 809-816, Madison, WI, USA, 2011. Omnipress.

[2] Antoine Bordes, Nicolas Usunier, Alberto Garcia-Duran, Jason Weston, and Oksana Yakhnenko. Translating embeddings for modeling multi-relational data. In C. J. C. Burges, L. Bottou, M. Welling, Z. Ghahramani, and K. Q. Weinberger, editors, Advances in Neural Information Processing Systems, Vol. 26, pp. 2787-2795. Curran Associates, Inc., 2013.

[3] Zhiqing Sun, Zhi-Hong Deng, Jian-Yun Nie, and Jian Tang. Rotate: Knowledge graph embedding by relational rotation in complex space. In 7th International Conference on Learning Representations, ICLR 2019, New Orleans, LA, USA, May 6-9, 2019. OpenReview.net, 2019.

[4] Charles Elkan and Keith Noto. Learning classifiers from only positive and unlabeled data. In Proceedings of the 14th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, KDD ’08, p. 213-220, New York, 2008. Association for Computing Machinery.

[5] Qinghua Guo Wenkai Li and Charles Elkan. A positive and unlabeled learning algorithm for one-class classification of remote-sensing data. IEEE Transactions on Geoscience and Remote Sensing, Vol. 49, No. 2, pp. 717-725, 2011.

[6] Nagarajan Natarajan, Inderjit S Dhillon, Pradeep K Ravikumar, and Ambuj Tewari. Learning with noisy labels. In C. J. C. Burges, L. Bottou, M. Welling, Z. Ghahramani, and K. Q. Weinberger, editors, Advances in Neural Information Processing Systems 26, pp. 1196-1204. Curran Associates, Inc., 2013.

[7] Cho-Jui Hsieh, Nagarajan Natarajan, and Inderjit $S$ Dhillon. Pu learning for matrix completion. In Proceedings of the 32nd International Conference on International Conference on Machine Learning - Volume 37, ICML'15, p. 2445-2453. JMLR.org, 2015.

[8] Yanjie Wang, Rainer Gemulla, and Hui Li. On multirelational link prediction with bilinear models. In Sheila A. McIlraith and Kilian Q. Weinberger, editors, Proceedings of the Thirty-Second AAAI Conference on Artificial Intelligence, (AAAI-18), the 30th innovative Applications of Artificial Intelligence (IAAI-18), and the 8th AAAI Symposium on Educational Advances in Artificial Intelligence (EAAI-18), New Orleans, Louisiana, USA, February 2-7, 2018, pp. 4227-4234. AAAI Press, 2018.

[9] Daniel Ruffinelli, Samuel Broscheit, and Rainer Gemulla. You CAN teach an old dog new tricks! on training knowledge graph embeddings. In 8th International Conference on Learning Representations, ICLR 2020, Addis Ababa, Ethiopia, April 26-30, 2020. OpenReview.net, 2020.

[10] Maximilian Nickel, Volker Tresp, and Hans-Peter Kriegel. Factorizing yago: Scalable machine learning for linked data. In Proceedings of the 21st International Conference on World Wide Web, WWW '12, p. 271-280, New York,
NY, USA, 2012. Association for Computing Machinery.

[11] Théo Trouillon, Christopher R. Dance, Johannes Welbl, Sebastian Riedel, Éric Gaussier, and Guillaume Bouchard. Knowledge graph completion via complex tensor factorization. CoRR, Vol. abs/1702.06879, , 2017.
 垣渕太成 $* 1$ 林克彦 $\dagger 2$ 駒谷和範 $\$ 1$ 1 大阪大学 産業科学研究所 2 群馬大学 社会情報学部  