(57)【要約】【課題】複数の情報源からなるデータの異常値を検出することにより、精度よく異常事例を検出する。【解決手段】入力部１０により複数の情報源からなる事例の集合を取得し、潜在ベクトル割当部３０により、事例の集合に基づいて、事例の各々について、事例に含まれる複数のデータの特徴量の組み合わせの関係性が他の事例と一貫している場合に、事例における複数の情報源に対して同一の潜在ベクトルが割り当てられ、事例に含まれる複数のデータの特徴量の組み合わせの関係性が他の事例と一貫していない場合に、事例における複数の情報源に対して異なる潜在ベクトルが割り当てられるように、事例における複数の情報源に対して潜在ベクトルを割り当て、異常スコア算出部４０により、事例の各々について、事例における潜在ベクトルの種類数に基づいて、事例の異常値スコアを算出する。【選択図】図１
2015026218.tif 000002
【特許請求の範囲】
【請求項１】
複数の情報源から得られる複数のデータからなる事例の集合を取得するデータ取得手段と、
前記データ取得手段によって取得された前記事例の集合に基づいて、前記事例の各々について、前記事例に含まれる前記複数のデータの特徴量の組み合わせの関係性が他の事例と一貫している場合に、前記事例における前記複数の情報源に対して同一の潜在ベクトルが割り当てられ、前記事例に含まれる前記複数のデータの特徴量の組み合わせの関係性が他の事例と一貫していない場合に、前記事例における前記複数の情報源に対して異なる潜在ベクトルが割り当てられるように、前記事例における前記複数の情報源に対して前記潜在ベクトルを割り当てる潜在ベクトル割当手段と、
前記事例の各々について、前記潜在ベクトル割当手段によって前記事例における前記複数の情報源に対して割り当てられた前記潜在ベクトルの種類数に基づいて、前記事例の異常値スコアを算出する異常値スコア算出手段と、
を含む異常事例検出装置。
【請求項２】
前記潜在ベクトル割当手段は、
前記事例ｎの各々について、前記事例ｎにおける前記複数の情報源ｄに対して潜在ベクトルの初期値を設定すると共に、前記複数の情報源ｄの各々について、前記潜在ベクトルを前記情報源ｄのデータの特徴量に変換するための射影行列の初期値を設定する初期値設定手段と、
前記初期値設定手段により前記事例ｎの各々について設定又は前回決定された前記複数の情報源ｄの各々に対する前記潜在ベクトルと、前記初期値設定手段により設定又は前回推定された前記複数の情報源ｄの各々の射影行列と、前記事例ｎの各々に含まれる前記複数のデータの特徴量とに基づいて、前記事例ｎの各々及び前記複数の情報源ｄの各々について、前記事例ｎにおいて前記情報源ｄに各潜在ベクトルが割り当てられる割当度、及び前記事例ｎにおいて前記情報源ｄに新規潜在ベクトルが割り当てられる割当度を算出し、前記算出した前記各潜在ベクトルの割当度及び前記新規潜在ベクトルの割当度に基づいて、各潜在ベクトル及び前記新規潜在ベクトルの中から、前記事例ｎにおいて前記情報源ｄに割り当てる潜在ベクトルを決定する割当手段と、
前記割当手段により前記事例ｎの各々について決定した前記複数の情報源ｄの各々の潜在ベクトルと、前記事例ｎの各々に含まれる前記複数のデータの特徴量とに対して、尤もらしくなるように前記複数の情報源ｄの各々の射影行列を推定する射影行列推定手段と、
予め定められた反復条件を満足するまで、前記割当手段による決定及び前記射影行列推定手段による推定を繰り返す反復手段と、を含み、
前記異常値スコア算出手段は、前記事例ｎの各々について、前記割当手段により繰り返し決定された前記事例ｎの潜在ベクトルの種類数に基づいて、前記事例ｎの異常値スコアを算出する請求項１記載の異常事例検出装置。
【請求項３】
前記異常値スコア算出手段は、前記事例ｎの各々について、前記割当手段により繰り返し決定された前記事例ｎの潜在ベクトルのうち、所定回数以上の繰り返しで決定された前記事例ｎの潜在ベクトルの種類数に基づいて、前記事例ｎの異常値スコアを算出する請求項２記載の異常事例検出装置。
【請求項４】
データ取得手段と、潜在ベクトル割当手段と、異常値スコア算出手段と、を含む異常事例検出装置における異常事例検出方法であって、
前記データ取得手段は、複数の情報源から得られる複数のデータからなる事例の集合を取得し、
前記潜在ベクトル割当手段は、前記データ取得手段によって取得された前記事例の集合に基づいて、前記事例の各々について、前記事例に含まれる前記複数のデータの特徴量の組み合わせの関係性が他の事例と一貫している場合に、前記事例における前記複数の情報源に対して同一の潜在ベクトルが割り当てられ、前記事例に含まれる前記複数のデータの特徴量の組み合わせの関係性が他の事例と一貫していない場合に、前記事例における前記複数の情報源に対して異なる潜在ベクトルが割り当てられるように、前記事例における前記複数の情報源に対して前記潜在ベクトルを割り当て、
前記異常値スコア算出手段は、前記事例の各々について、前記潜在ベクトル割当手段によって前記事例における前記複数の情報源に対して前記潜在ベクトルの種類数に基づいて、前記事例の異常値スコアを算出する
異常事例検出方法。
【請求項５】
コンピュータを、前記請求項１から請求項３の何れか１項記載の異常事例検出装置を構成する各手段として機能させるためのプログラム。
【請求項６】
コンピュータを、前記請求項１から請求項３の何れか１項記載の異常事例検出装置を構成する各手段として機能させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、異常事例検出装置、方法、プログラム、及び記録媒体に係り、特に、複数の情報源から得られる複数のデータからなる異常事例を検出する異常事例検出装置、方法、プログラム、及び記録媒体に関する。
【背景技術】
【０００２】
複数の情報源からデータが得られる場合がある。例えば、画像は色、テクスチャ、局所特徴量などの異なる情報源からのデータで表現できる。また、映画はジャンル、出演者、その映画を見たいユーザの集合で表現されたり、ある百科事典の項目が異なる言語で表現されていたりする。
【０００３】
異なる情報源のデータが一貫していない異常値を持つ事例を検出することにより、例えば、色とテクスチャの組み合わせが特殊な画像を検出できたり、特異な映画（そのジャンルからは好まれないと思われるユーザ群に好まれる映画）を発見できたり、誤情報が掲載されている多言語百科事典の項目を検出できる。
【０００４】
従来、１つの情報源からなるデータのみを用いて異常値を検出している（非特許文献１）。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Tomoharu Iwata, Kazumi Saito, "Visualization of Anomalies using Mixture Models," Journal of Intelligent Manufacturing, Vol.16, 635-643, 2005
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、非特許文献１の手法においては複数の情報源からなるデータは活用できないという問題がある。
【０００７】
本発明では、上記問題点を解決するために成されたものであり、複数の情報源から得られる複数のデータからなる事例の集合を取得し、事例の各々について、事例に含まれる複数の情報源に対して潜在ベクトルを割り当て、事例に含まれる複数の情報源に対して割り当てられた潜在ベクトルの種類数に基づいて、事例の異常値スコアを算出することにより、複数の情報源から得られる複数のデータからなる事例の集合から、精度よく異常事例を検出することができる異常事例検出装置、方法、プログラム、及び記録媒体を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するために、第１の発明の異常事例検出装置は、複数の情報源から得られる複数のデータからなる事例の集合を取得するデータ取得手段と、前記データ取得手段によって取得された前記事例の集合に基づいて、前記事例の各々について、前記事例に含まれる前記複数のデータの特徴量の組み合わせの関係性が他の事例と一貫している場合に、前記事例における前記複数の情報源に対して同一の潜在ベクトルが割り当てられ、前記事例に含まれる前記複数のデータの特徴量の組み合わせの関係性が他の事例と一貫していない場合に、前記事例における前記複数の情報源に対して異なる潜在ベクトルが割り当てられるように、前記事例における前記複数の情報源に対して前記潜在ベクトルを割り当てる潜在ベクトル割当手段と、前記事例の各々について、前記潜在ベクトル割当手段によって前記事例における前記複数の情報源に対して前記潜在ベクトルの種類数に基づいて、前記事例の異常値スコアを算出する異常値スコア算出手段と、を含んで構成されている。
【０００９】
第２の発明の異常事例検出方法は、データ取得手段と、潜在ベクトル割当手段と、異常値スコア算出手段と、を含む異常事例検出装置における異常事例検出方法であって、前記データ取得手段は、複数の情報源から得られる複数のデータからなる事例の集合を取得し、前記潜在ベクトル割当手段は、前記データ取得手段によって取得された前記事例の集合に基づいて、前記事例の各々について、前記事例に含まれる前記複数のデータの特徴量の組み合わせの関係性が他の事例と一貫している場合に、前記事例における前記複数の情報源に対して同一の潜在ベクトルが割り当てられ、前記事例に含まれる前記複数のデータの特徴量の組み合わせの関係性が他の事例と一貫していない場合に、前記事例における前記複数の情報源に対して異なる潜在ベクトルが割り当てられるように、前記事例における前記複数の情報源に対して前記潜在ベクトルを割り当て、前記異常値スコア算出手段は、前記事例の各々について、前記潜在ベクトル割当手段によって前記事例における前記複数の情報源に対して前記潜在ベクトルの種類数に基づいて、前記事例の異常値スコアを算出する。
【００１０】
このように、第１の発明及び第２の発明によれば、複数の情報源から得られる複数のデータからなる事例の集合を取得し、事例の各々について、事例に含まれる複数の情報源に対して潜在ベクトルを割り当て、事例に含まれる複数の情報源に対して割り当てられた潜在ベクトルの種類数に基づいて、事例の異常値スコアを算出することにより、複数の情報源から得られる複数のデータからなる事例の集合から、精度よく異常事例を検出することができる。
【００１１】
また、第１の発明において、前記潜在ベクトル割当手段は、前記事例ｎの各々について、前記事例ｎにおける前記複数の情報源ｄに対して潜在ベクトルの初期値を設定すると共に、前記複数の情報源ｄについて、前記潜在ベクトルを前記情報源ｄのデータの特徴量に変換するための射影行列の初期値を設定する初期値設定手段と、前記初期値設定手段により前記事例ｎの各々について設定又は前回決定された前記複数の情報源ｄの各々に対する前記潜在ベクトルと、前記初期値設定手段により設定又は前回推定された前記複数の情報源ｄの各々の射影行列と、前記事例ｎの各々に含まれる前記複数のデータの特徴量とに基づいて、前記事例ｎの各々及び前記複数の情報源ｄの各々について、前記事例ｎにおいて前記情報源ｄに各潜在ベクトルが割り当てられる割当度、及び新規潜在ベクトルが割り当てられる割当度を算出し、前記算出した前記各潜在ベクトルの割当度及び前記新規潜在ベクトルの割当度に基づいて、各潜在ベクトル及び前記新規潜在ベクトルの中から、前記事例ｎにおいて前記情報源ｄに割り当てる潜在ベクトルを決定する割当手段と、前記割当手段により前記事例ｎの各々について決定した前記複数の情報源ｄの各々の潜在ベクトルと、前記事例ｎの各々に含まれる前記複数のデータの特徴量とに対して、尤もらしくなるように前記複数の情報源ｄの各々の射影行列を推定する射影行列推定手段と、予め定められた反復条件を満足するまで、前記割当手段による決定及び前記射影行列推定手段による推定を繰り返す反復手段と、を含み、前記異常値スコア算出手段は、前記事例ｎの各々について、前記割当手段により繰り返し決定された前記事例ｎの潜在ベクトルの種類数に基づいて、前記事例ｎの異常値スコアを算出することができる。
【００１２】
また、前記異常値スコア算出手段は、前記事例ｎの各々について、前記割当手段により繰り返し決定された前記事例ｎの潜在ベクトルのうち、所定回数以上の繰り返しで決定された前記事例ｎの潜在ベクトルの種類数に基づいて、前記事例ｎの異常値スコアを算出することができる。
【００１３】
また、本発明のプログラムは、コンピュータを、上記の異常事例検出装置を構成する各手段として機能させるためのプログラムである。
【００１４】
また、本発明の記録媒体は、コンピュータを、上記の異常事例検出装置を構成する各手段として機能させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体。
【発明の効果】
【００１５】
以上説明したように、本発明の異常事例検出装置、方法、プログラム、及び記録媒体によれば、複数の情報源から得られる複数のデータからなる事例の集合を取得し、事例の各々について、事例に含まれる複数の情報源に対して潜在ベクトルを割り当て、事例に含まれる複数の情報源に対して割り当てられた潜在ベクトルの種類数に基づいて、事例の異常値スコアを算出することにより、複数の情報源から得られる複数のデータからなる事例の集合から、精度よく異常事例を検出することができる。
【図面の簡単な説明】
【００１６】
【図１】本発明の実施の形態に係る異常事例検出装置の機能的構成を示すブロック図である。
【図２】本発明の実施の形態において用いる表記の例を示す図である。
【図３】正常事例と異常事例の例を示す図である。
【図４】本発明の実施の形態に係る異常事例検出装置における異常事例検出処理ルーチンの内容を示すフローチャートである。
【図５】異なる異常値割合のデータセットでのＡＵＣの平均と標準誤差の例を示す図である。
【発明を実施するための形態】
【００１７】
以下、図面を参照して本発明の実施の形態を詳細に説明する。
【００１８】
＜異常事例検出装置の構成＞
【００１９】
本発明の実施の形態に係る異常事例検出装置について説明する。図１に示すように、本発明の実施の形態に係る異常事例検出装置１００は、ＣＰＵと、ＲＡＭと、後述する異常事例検出処理ルーチンを実行するためのプログラムや各種データを記憶したＲＯＭと、を含むコンピュータで構成することが出来る。この異常事例検出装置１００は、機能的には図１に示すように入力部１０と、演算部２０と、出力部５０とを備えている。
【００２０】
入力部１０は、キーボードなどの入力装置から、Ｄ個の情報源から得られるＤ個のデータの特徴量からなる事例をＮ個含む事例集合Ｘ（Ｘ＝｛Ｘｎ｝Ｎｎ＝１）を受け付ける。ここでＸｎ＝｛ｘｎｄ｝Ｄｄ＝１は、ｎ番目の事例に含まれるＤ個情報源のデータを表し、
【００２１】
【数１】
2015026218.tif 000003
【００２２】
はそのｄ番目の情報源から得られるデータの特徴量を表す。図２に本実施の形態で用いる表記を示す。なお、入力部１０は、ネットワーク等を介して外部から入力されたものを受け付けるようにしてもよい。
【００２３】
演算部２０は、潜在ベクトル割当集合記憶部２２、異常スコア集合記憶部２４、潜在ベクトル割当部３０、及び異常スコア算出部４０を備えている。
【００２４】
潜在ベクトル割当部３０は、初期化部３２、割当部３４、射影行列推定部３６、及び反復終了条件判定部３８を備えている。
【００２５】
潜在ベクトル割当部３０は、マルコフ連鎖モンテカルロ法（ＭＣＭＣ：Markovchain Monte Carlo methods）に従って、潜在ベクトルの現在の割り当てがどのくらい入力データを尤もらしく説明できているかを表す下記（１）式が示す尤度が収束するまで、入力部１０において受け付けたＮ個の事例の各々について、Ｄ個の情報源の各々に対する潜在ベクトルの割り当てを繰り返す。ここで、ある事例ｎの全ての情報源の特徴量の組み合わせの関係性が他の事例と一貫している場合、当該事例ｎの全情報源に１つの同一の潜在ベクトルが割り当てられていれば尤度は高くなり、ある事例ｎの全ての情報源の特徴量の組み合わせの関係性が他の事例と一貫していない場合、当該事例ｎの情報源の各々に異なる潜在ベクトルが割り当てられていれば、尤度は高くなる。なお、潜在ベクトルの次元は予め定められた次元数とする。また、潜在ベクトルの現在の割り当てがどのくらい入力データを尤もらしく説明できているかを表す尤度が収束するまでという条件が、予め定められた反復条件の一例である。また、下記（１）式は尤度の一例である。
【００２６】
図３に、全情報源の特徴量の組み合わせの関係性が他の事例と一貫している正常事例と、複数情報源の特徴量の組み合わせの関係性が他の事例と一貫していない異常事例とを示す。図３に示す例においては、６つの情報源がある。また、正常事例の場合、全ての情報源に同一の潜在ベクトルＺｎ１が割り当てられている。一方、異常事例の場合、１、２、及び５番目の情報源に１番目の潜在ベクトルＺｎ´１が、３番目の事例に２番目の潜在ベクトルＺｎ´２が、４及び６番目の事例に３番目の潜在ベクトルＺｎ´３が割り当てられている。
【００２７】
【数２】
2015026218.tif 000004
【００２８】
ここで、Ｓ＝｛｛ｓｎｄ｝Ｄｄ＝１｝Ｎｎ＝１は潜在ベクトル割当集合、ｓｎｄはｎ番目の事例のｄ番目の情報源に割り当てられた潜在ベクトルのインデックスを表し、Ｗ＝｛Ｗｄ｝Ｄｄ＝１は射影行列集合、ａ、ｂ、ｒ、及びγは予め定められたパラメータである。ｄ番目の射影行列Ｗｄは、潜在ベクトルをｄ番目の情報源の特徴量に変換するための行列である。上記（１）式の第一項は、下記（２）式で計算できる。
【００２９】
【数３】
2015026218.tif 000005
【００３０】
ここで、Ｎｎｊは、ｎ番目の事例においてｊ番目の潜在ベクトルに割り当てられた情報源の数、Ｊｎはｎ番目の事例で用いられたＮｎｊ＞０である潜在ベクトルの数を表す。また、上記（１）式の第二項は、下記（３）式で計算できる。
【００３１】
【数４】
2015026218.tif 000006
【００３２】
初期化部３２は、入力部１０において受け付けたＮ個の事例の各々の潜在ベクトル数を１に設定する。また、各事例ｎの各情報源ｄに同一の潜在ベクトルを割り当てる（ｓｎｄ＝１，ｎ＝１，…，Ｎ，ｄ＝１，…，Ｄ）。そして、情報源ｄの各々の射影行列Ｗｄからなる射影行列集合Ｗをランダムに初期化する。
【００３３】
割当部３４は、入力部１０において受け付けた各事例ｎ（ｎ＝１，…，Ｎ）の各情報源ｄ（ｄ＝１，…，Ｄ）の各々について、事例ｎにおける情報源ｄに対してどの潜在ベクトルに割り当てられるかの度合いを表す割当度を算出し、算出した割当度に基づいて、事例ｎにおける情報源ｄに対して割り当てる潜在ベクトルを決定する。
【００３４】
具体的には、ＣＲＰ（Chinese restaurant process）に従って、事例ｎ毎に、情報源ｄの各々について、予め定められた各パラメータ（ａ、ｂ、ｒ、γ）及び当該事例ｎに含まれるＤ個の情報源から得られるデータの特徴量と、当該事例ｎにおける当該情報源ｄに対する潜在ベクトルの割り当てを除いた潜在ベクトル割当集合Ｓと、初期化部３２において初期化され又は射影行列推定部３６において推定された情報源ｄの各々についての射影行列Ｗｄからなる射影行列集合Ｗとに基づいて、下記（６）式により潜在ベクトル毎の割当度を算出する。なお、下記（６）式はｊ番目の潜在ベクトルの割当度を算出するものとする。
【００３５】
【数５】
2015026218.tif 000007
【００３６】
ここで、Ｌ＝（ｎ，ｄ）はｎ番目の事例のｄ番目の情報源のインデックスを表す。＼Ｌはｎ番目の事例のｄ番目を取り除いた場合の値もしくは集合を表す。ここで、上記（６）式の第一項は下記（７）式によって計算できる。
【００３７】
【数６】
2015026218.tif 000008
【００３８】
ただし、上記（７）式はＣＲＰ（Chinese restaurant process）による情報源ｄの割り当てられる潜在ベクトルの決定式である。また、上記（６）式の第二項は下記（８）式によって計算できる。
【００３９】
【数７】
2015026218.tif 000009
【００４０】
なお、Ｉ（・）は指示関数、つまりＡが真ならばＩ（Ａ）＝１そうでなければ０を表し、ｓＬ＝ｊは、ｎ番目の事例のｄ番目の情報源がｊ番目の潜在ベクトルに割り当てられた場合の値を表し、下記（９）式及び（１０）式に示すようになる。
【数８】
2015026218.tif 000010
【００４１】
上記（６）式〜（８）式では、直観的には、当該事例ｎにおける複数のデータの特徴量の組み合わせの関係性が、他の事例と一貫しておらず、情報源が既存の潜在ベクトルでモデル化できない場合には、新規の潜在ベクトルが割り当てられるように、割当度が算出される。
【００４２】
割当部３４は、各事例ｎの各情報源ｄについて、上記（６）式から上記（１０）式に基づいて算出された潜在ベクトル毎の割当度に基づいて、当該事例ｎにおける当該情報源ｄに割り当てる潜在ベクトルを、ランダムサンプリングにより｛１，…，Ｊｎ＼Ｌ＋１｝から決定する。割当部３４による潜在ベクトルの決定は、試行毎に繰り返し行われ、試行毎に、各事例ｎの各情報源ｄに割り当てる潜在ベクトルからなる潜在ベクトル割当集合Ｓと、潜在ベクトル割当集合Ｓから求まる、各事例ｎで使用された潜在ベクトル数Ｊnとを、潜在ベクトル割当集合記憶部２２に記憶する。
【００４３】
射影行列推定部３６は、割当部３４において決定された各事例ｎ及び各情報源ｄに対する潜在ベクトル割当に基づいて、上記（１）式の尤度が高くなるように情報源ｄ毎の射影行列Ｗｄを推定しメモリ（図示省略）に記憶する。具体的には、準ニュートン法等の最適化手法を用いることにより、射影行列を推定する。
【００４４】
反復終了条件判定部３８は、割当部３４において各事例ｎの情報源ｄの各々について決定された潜在ベクトル割当からなる潜在ベクトル割当集合Ｓと、潜在ベクトル割当集合Ｓから求まる、各事例ｎで使用された潜在ベクトル数と、射影行列推定部３６において推定された情報源ｄの各々の射影行列Ｗｄからなる射影行列集合Ｗと、入力部１０において受け付けた各パラメータ（ａ，ｂ，ｒ，γ）と、当該事例ｎに含まれるＤ個の情報源から得られるデータの特徴量とに基づいて、上記（１）式に従って尤度を算出する。そして、前回求めた尤度と今回求めた尤度の差が予め定められた値の範囲内である場合には、反復終了条件を満たすと判定する。一方、前回求めた尤度と今回求めた尤度の差が予め定められた値の範囲でない場合には、割当部３４による決定処理と射影行列推定部３６による推定処理は、繰り返し反復終了条件を満足するまで、繰り返し行われる。なお、１回目の試行の場合には、繰返し終了条件は満たさないものとする。
【００４５】
潜在ベクトル割当集合記憶部２２は、潜在ベクトル割当部３０により、試行毎に決定された各事例ｎの情報源ｄの各々の潜在ベクトル割当からなる潜在ベクトル割当集合Ｓと、試行毎の潜在ベクトル割当集合Ｓから求まる試行毎の各事例ｎで使用された潜在ベクトル数とを記憶している。
【００４６】
異常スコア算出部４０は、割当読込部４２と、スコア算出部４４とを備えている。異常スコア算出部４０は、潜在ベクトル割当集合記憶部２２に記憶された、試行毎の各事例ｎで使用された潜在ベクトル数に基づいて、事例の各々について異常スコアを算出する。
【００４７】
割当読込部４２は、潜在ベクトル割当集合記憶部２２に記憶されている、試行毎の各事例ｎで使用された潜在ベクトル数を読み込む。具体的には、試行ｈ回目のｎ番目の事例で使われた潜在ベクトル数の集合｛｛Ｊ（ｈ）ｎ｝Ｈｈ＝１｝Ｎｎ＝１を読み込む。
【００４８】
スコア算出部４４は、割当読込部４２において読み込まれた試行毎の事例の各々の潜在ベクトル数に基づいて、事例の各々の異常値スコアを算出する。ここで、異常値スコアは、その事例の全情報源を表現するために１つより多くの潜在ベクトルが使われる確率であり、下記（１１）式により算出する。そして、事例の各々の異常値スコアの算出結果Ｖ＝｛ｖｎ｝Ｎｎ＝１を異常スコア集合記憶部２４に記憶すると共に、出力部５０に出力する。なお、潜在ベクトル数が潜在ベクトルの種類数の一例である。下記（１１）式のＨは潜在ベクトル割当部３０における割当部３４の決定処理の試行の回数を表す。
【００４９】
【数９】
2015026218.tif 000011
【００５０】
なお、繰り返しの初期では推定が安定していないため、所定回数未満の試行の結果を省いて、所定回数以上の試行の各々の、各事例ｎで使用された潜在ベクトル数を用いて、異常スコアを計算してもよい。
【００５１】
＜異常事例検出装置の作用＞
次に、本発明の実施の形態に係る異常事例検出装置１００の作用について説明する。まず、入力部１０によりＤ個の情報源から得られるＤ個のデータの特徴量からなる事例をＮ個含む事例集合Ｘ（Ｘ＝｛Ｘｎ｝Ｎｎ＝１）が入力される。そして、異常事例検出装置１００のＲＯＭに記憶されたプログラムを、ＣＰＵが実行することにより、図４に示す異常事例検出処理ルーチンが実行される。
【００５２】
まず、ステップＳ１００では、入力部１０により入力された事例集合Ｘ（Ｘ＝｛Ｘｎ｝Ｎｎ＝１）を受け付ける。
【００５３】
次に、ステップＳ１０２では、ステップＳ１００において取得したＮ個の事例の各々の潜在ベクトル数を１に設定する。
【００５４】
次に、ステップＳ１０４では、Ｄ個の情報源ｄの各々の射影行列Ｗｄからなる射影行列集合Ｗをランダムに初期化する。
【００５５】
次に、ステップＳ１０６では、処理対象の事例ｎにおける処理対象の情報源ｄについて、予め定められた各パラメータ（ａ、ｂ、ｒ、γ）及び当該事例ｎに含まれるＤ個の情報源から得られるデータ特徴量と、当該事例ｎにおける当該情報源ｄに対する潜在ベクトルの割り当てを除いた潜在ベクトル集合Ｓと、ステップＳ１０４において初期化され又は前回推定された情報源ｄの各々についての射影行列Ｗｄからなる射影行列集合Ｗとに基づいて、上記（６）式により潜在ベクトル毎の割当度を算出する。
【００５６】
次に、ステップＳ１１０では、処理対象の事例ｎにおける処理対象の情報源ｄについて、ステップＳ１０６において算出された潜在ベクトル毎の割当度に基づいて、処理対象の事例ｎにおける処理対象の情報源ｄに割り当てる潜在ベクトルを、ランダムサンプリングにより｛１，…，Ｊｎ＼Ｌ＋１｝から決定する。
【００５７】
次に、ステップＳ１１２では、処理対象となる事例ｎにおける全ての情報源ｄについてステップＳ１０６〜ステップＳ１１２までの処理を終了したか否かを判定する。処理対象となる事例ｎにおける全ての情報源ｄについて処理を終了している場合には、ステップＳ１１４へ移行し、処理対象となる事例ｎにおける全ての情報源ｄについて処理を終了していない場合には、処理対象となる情報源ｄを変更し、ステップＳ１０６へ移行し処理を繰り返す。
【００５８】
次に、ステップＳ１１４では、事例ｎの全てについてステップＳ１０６〜ステップＳ１１２までの処理を終了したか否かを判定する。全ての事例ｎについて処理を終了している場合には、ステップＳ１１６へ移行し、全ての事例ｎについて処理を終了していない場合には、処理対象となる事例ｎを変更し、ステップＳ１０６へ移行し処理を繰り返す。
【００５９】
次に、ステップＳ１１６では、ステップＳ１１０において決定した各事例ｎにおける各情報源ｄに割り当てた潜在ベクトルからなる潜在ベクトル割当集合Ｓと、潜在ベクトル割当集合Ｓから求まる、各事例ｎで使用された潜在ベクトル数Ｊｎとを、今回の試行の結果として潜在ベクトル割当集合記憶部２２に記憶する。
【００６０】
次に、ステップＳ１１８では、ステップＳ１１６において取得した各事例ｎ及び各情報源ｄに対する潜在ベクトル割当に基づいて、上記（１）式の尤度が高くなるように情報源ｄ毎の射影行列Ｗｄを推定し、射影行列集合Ｗを更新する。
【００６１】
次に、ステップＳ１２０では、反復終了条件を満たすか否かを判定する。反復終了条件を満たす場合には、ステップＳ１２２へ移行し、反復終了条件を満たさない場合には、ステップＳ１０６へ移行し、ステップＳ１０６〜ステップＳ１１８の処理を繰り返す。
【００６２】
次に、ステップＳ１２２では、潜在ベクトル割当集合記憶部２２に記憶されている試行毎の各事例ｎで使用された潜在ベクトル数Ｊｎを読み込む。
【００６３】
次に、ステップＳ１２４では、ステップＳ１２２において取得した試行毎の各事例ｎで使用された潜在ベクトル数Ｊｎに基づいて、上記（１１）式に従って、各事例ｎについての異常値スコアを算出する。
【００６４】
次に、ステップＳ１２６では、ステップＳ１２４において取得した各事例ｎの異常値スコアを処理結果として異常スコア集合記憶部２４に記憶すると共に、出力部５０に出力して処理を終了する。
【００６５】
＜異常事例検出装置の実施例＞
次に、本実施の形態に係る異常事例検出装置１００において１１データセットを用いて実験を行った結果を示す。本実験においては、各データセットの特徴量をランダムに分割し、複数の情報源を作成した。また、ランダムに選んだ２つの事例のある情報源の特徴量を交換することにより異常値を挿入した。
【００６６】
比較手法としてＨＯＡＤ（HOrizontal Anomaly Detection）（非特許文献２（J. Gao, W. Fan, D. Turaga, S. Parthasarathy, and J. Han. A spectral framework for detecting inconsistency across multi-source object relationships. In IEEE 11th International Conference on Data Mining (ICDM), pages 1050{1055. IEEE, 2011.）を参照）及び、Probabilistic Canonical Correlation Analysis（ＰＣＣＡ）（非特許文献３（F. R. Bach and M. I. Jordan. A probabilistic interpretation of canonical correlation analysis. Technical Report 688, Department of Statistics, University of California, Berkeley, 2005.）を参照）を用いた。異常値の割合を変化させたときの異常値検出精度を図５に示す。評価指標としてＡＵＣ（Area Under the ROC Curve）を用いた。高いＡＵＣは高い検出精度を示す。多くのデータセットで本実施の形態に係る異常事例検出装置１００（Proposed）が最も高い精度を達成しており、複数情報源から異常値を効果的に検出できることを示唆する。
【００６７】
以上、説明したように、本発明の実施の形態に係る検出装置によれば、複数の情報源から得られるデータからなる事例の異常値スコアを、各情報源に割り当てた潜在ベクトルの数に基づいて算出することにより、精度よく異常事例を検出することができる。
【００６８】
なお、本発明は、上記実施の形態に限定されるものではなく、この発明の要旨を逸脱しない範囲内で様々な変形や応用が可能である。
【００６９】
本実施の形態においては、上記（１）式を尤度として求める場合について説明したがこの限りでない。例えば、上記（１）式の代わりに事後確率など同等の性質をもつ値を用いてもよい。
【００７０】
また、本実施の形態においては、マルコフ連鎖モンテカルロ法に基づいて潜在ベクトルを割り当てる場合について説明したがこの限りでない。例えば、最尤推定法や変分ベイズ法などを用いてもよい。
【００７１】
また、本実施の形態においては、潜在ベクトルの次元は予め定められた次元数としている場合について説明したがこの限りでない。例えば、潜在ベクトルの次元を交差検定等を用いて推定してもよい。
【００７２】
また、本実施の形態においては、反復終了条件として尤度の変化の大きさを用いる場合について説明したがこの限りではない。例えば、反復終了条件として予め定めた繰り返し回数、処理を繰り返したか否かを用いてもよい。
【符号の説明】
【００７３】
１０   入力部
２０   演算部
２２   潜在ベクトル割当集合記憶部
２４   異常スコア集合記憶部
３０   潜在ベクトル割当部
３２   初期化部
３４   割当部
３６   射影行列推定部
３８   反復終了条件判定部
４０   異常スコア算出部
４２   割当読込部
４４   スコア算出部
５０   出力部
１００ 異常事例検出装置
【図１】
2015026218.tif 000012
【図４】
2015026218.tif 000013
【図２】
2015026218.tif 000014
【図３】
2015026218.tif 000015
【図５】
2015026218.tif 000016