(57)【要約】【課題】クラスタブロックの重要度をより適切に特定することが可能な共クラスタリング装置を提供する。【解決手段】行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリング処理を行うことで、関係データを複数のクラスタブロックに分割する共クラスタリング装置１００であって、複数のクラスタブロックのそれぞれにおいて生成される関係の傾向を示す統計量の、関係データ全体における分布傾向を生成する分布傾向生成部１３０と、複数のクラスタブロックのそれぞれの統計量と、分布傾向生成部１３０が生成した分布傾向とに基づいて、分布傾向に応じて算出結果が変化する算出方法により、複数のクラスタブロックのそれぞれに対して重要度を算出する算出部１４０と、複数のクラスタブロックの少なくとも１つを示す情報と、当該少なくとも１つに対して算出部１４０が算出した重要度を示す情報とを出力する出力部１５０とを備える。【選択図】図１
2014081899.tif 000002
【特許請求の範囲】
【請求項１】
行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリング処理を行うことによって、前記関係データを複数のクラスタブロックに分割する共クラスタリング装置であって、
前記複数のクラスタブロックのそれぞれにおいて生成される関係の傾向を示す統計量の、前記関係データ全体における分布傾向を生成する分布傾向生成部と、
前記複数のクラスタブロックのそれぞれの統計量と、前記分布傾向生成部が生成した前記分布傾向とに基づいて、前記分布傾向に応じて算出結果が変化する算出方法により、前記複数のクラスタブロックのそれぞれに対して重要度を算出する算出部と、
前記複数のクラスタブロックの少なくとも１つを示す情報と、当該少なくとも１つに対して前記算出部が算出した前記重要度を示す情報とを出力する出力部とを備える
共クラスタリング装置。
【請求項２】
前記分布傾向生成部は、
前記関係データ全体の統計量を、前記分布傾向として生成する
請求項１に記載の共クラスタリング装置。
【請求項３】
前記算出部は、
前記複数のクラスタブロックのそれぞれについて、当該クラスタブロックにおける前記分布傾向により示される値と、当該クラスタブロックの統計量との距離が大きいほど、大きくなるように前記重要度を算出する
請求項２に記載の共クラスタリング装置。
【請求項４】
前記算出部は、
前記複数のクラスタブロックのそれぞれについて、前記分布傾向と、当該クラスタブロックの統計量と、当該クラスタブロックの大きさとを用いて前記重要度を算出する
請求項２に記載の共クラスタリング装置。
【請求項５】
前記分布傾向生成部は、
前記複数のクラスタブロックのそれぞれの統計量を要素とする統計量データに対してクラスタリング処理を行うことによって、前記統計量データを複数のクラスタに分割し、分割することで得られた前記複数のクラスタの情報を、前記分布傾向として生成する
請求項１に記載の共クラスタリング装置。
【請求項６】
前記算出部は、
前記複数のクラスタのそれぞれについて、当該クラスタ内の要素の数が少ないほど、当該クラスタに要素として含まれるクラスタブロックに対して大きくなるように前記重要度を算出する
請求項５に記載の共クラスタリング装置。
【請求項７】
前記算出部は、
前記複数のクラスタのそれぞれについて、当該クラスタ内の要素の数と、前記複数のクラスタの要素に対応する１以上の前記クラスタブロックの大きさとに基づいて、当該クラスタに要素として含まれる前記クラスタブロックのそれぞれに対して前記重要度を算出する
請求項５に記載の共クラスタリング装置。
【請求項８】
行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリング処理を行うことによって、前記関係データを複数のクラスタブロックに分割する共クラスタリング装置における共クラスタリング方法であって、
前記複数のクラスタブロックのそれぞれにおいて生成される関係の傾向を示す統計量の、前記関係データ全体における分布傾向を生成する分布傾向生成ステップと、
前記複数のクラスタブロックのそれぞれの統計量と、前記分布傾向生成ステップで生成した前記分布傾向とに基づいて、前記分布傾向に応じて算出結果が変化する算出方法により、前記複数のクラスタブロックのそれぞれに対して重要度を算出する算出ステップと、
前記複数のクラスタブロックの少なくとも１つを示す情報と、当該少なくとも１つに対して前記算出ステップで算出した前記重要度を示す情報とを出力する出力ステップとを含む
共クラスタリング方法。
【請求項９】
請求項８に記載の共クラスタリング方法をコンピュータに実行させるためのプログラム。
【請求項１０】
行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリング処理を行うことによって、前記関係データを複数のクラスタブロックに分割する集積回路であって、
前記複数のクラスタブロックのそれぞれにおいて生成される関係の傾向を示す統計量の、前記関係データ全体における分布傾向を生成する分布傾向生成部と、
前記複数のクラスタブロックのそれぞれの統計量と、前記分布傾向生成部が生成した前記分布傾向とに基づいて、前記分布傾向に応じて算出結果が変化する算出方法により、前記複数のクラスタブロックのそれぞれに対して重要度を算出する算出部と、
前記複数のクラスタブロックの少なくとも１つを示す情報と、当該少なくとも１つに対して前記算出部が算出した前記重要度を示す情報とを出力する出力部とを備える
集積回路。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリングを行う共クラスタリング装置、共クラスタリング方法、プログラム及び集積回路に関する。
【背景技術】
【０００２】
関係データを分析するための有効な手段の一つがクラスタリングである。関係データにおいて、オブジェクトの集合（以降、ドメインと呼ぶ）が複数存在する場合、個々のドメインに対するクラスタリングを同時に行うことができる。個々のドメインに対するクラスタリングを同時に行うことを、特に、共クラスタリング（co-clustering）と呼び、様々な研究がなされている。
【０００３】
このような従来の技術としては、例えば非特許文献１に記載されたものが知られている。非特許文献１で提案されているInfinite Relational Model（以降、ＩＲＭと呼ぶ）は、関係データの生成過程（generative process）を表すノンパラメトリックベイズモデル（non-parametric bayes model)であり、行列又は３次元以上のテンソル（tensor)の形式で表現される関係データに対して、関係の類似性に基づき共クラスタリングを行うことができる。
【０００４】
また、従来の技術として、例えば特許文献１に記載されたものが知られている。特許文献１は、関係データに対して関係の類似性に基づいた共クラスタリングを行い、入力された関係データを複数のクラスタブロックに分割する。分割する際には、複数のクラスタブロックについて、それぞれ統計量（相関強度）を算出して、算出された統計量を重要度とみなして複数のクラスタブロックを重要度が高い順が分かる態様になるよう並び替えて表示する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許第４６９０１９９号公報
【非特許文献】
【０００６】
【非特許文献１】C. Kemp, J. Tenenbaum, T. Griffiths, T. Yamada and U. Naonori: “Learning systems of concepts with an infinite relational model”, in Proceedings of the 21st national conference on Artificial intelligence - Volume 1, ser. AAAI’06. AAAI Press, 2006, pp. 381-388.
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、従来技術では、クラスタブロックの適切な重要度を特定することができないという問題がある。
【０００８】
そこで、本発明は、クラスタブロックのより適切な重要度を特定することが可能な共クラスタリング装置等を提供する。
【課題を解決するための手段】
【０００９】
本発明の一態様に係る共クラスタリング装置は、行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリング処理を行うことによって、前記関係データを複数のクラスタブロックに分割する共クラスタリング装置であって、前記複数のクラスタブロックのそれぞれにおいて生成される関係の傾向を示す統計量の、前記関係データ全体における分布傾向を生成する分布傾向生成部と、前記複数のクラスタブロックのそれぞれの統計量と、前記分布傾向生成部が生成した前記分布傾向とに基づいて、前記分布傾向に応じて算出結果が変化する算出方法により、前記複数のクラスタブロックのそれぞれに対して重要度を算出する算出部と、前記複数のクラスタブロックの少なくとも１つを示す情報と、当該少なくとも１つに対して前記算出部が算出した前記重要度を示す情報とを出力する出力部とを備える。
【００１０】
なお、これらの包括的または具体的な態様は、システム、方法、集積回路、コンピュータプログラムまたはコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラムおよび記録媒体の任意な組み合わせで実現されてもよい。
【発明の効果】
【００１１】
本発明の共クラスタリング装置は、クラスタブロックのより適切な重要度を特定することができる。
【図面の簡単な説明】
【００１２】
【図１】図１は、実施の形態１に係る共クラスタリング装置の構成の例を示すブロック図である。
【図２】図２は、実施の形態１に係る関係データの例を示す図である。
【図３】図３は、実施の形態１に係る関係データの他の例を示す図である。
【図４】図４は、実施の形態１に係る共クラスタリングの説明図である。
【図５】図５は、実施の形態１に係る共クラスタリング装置の動作の例を示すフローチャートである。
【図６】図６は、実施の形態１に係る共クラスタリング装置の処理の例を示す図である。
【図７】図７は、実施の形態１に係る共クラスタリング装置の構成の他の例を示すブロック図である。
【図８】図８は、実施の形態２に係る共クラスタリング装置の構成の例を示すブロック図である。
【発明を実施するための形態】
【００１３】
（本発明の基礎となった知見）
本発明者は、「背景技術」の欄において記載した、関係データの分析方法に関し、以下の問題が生じることを見出した。
【００１４】
昨今では、生活やビジネスなどの様々な状況において、インターネットの利用が必要不可欠なものになった。それに伴い、個人が「どんな商品を買ったか」、あるいは、「誰と知り合いか」など、個人の社会活動において必ず生まれる他者（或いは物）との関係が、電子化情報として蓄積されるようになり、これら関係を表す情報（以降、関係データ（relational data）と呼ぶ）からニーズや嗜好などの潜在的な傾向を分析することがますます重要になってきている。
【００１５】
関係データを分析するための有効な手段の一つがクラスタリングである。関係データのクラスタリングでは、人や物など、関係を生成する主体（以降、オブジェクトと呼ぶ）は、自分が所属するクラスタに依存して、特徴的な傾向にもとづいて他のオブジェクトとの関係を生成するという仮説に基づき、類似したオブジェクトをグループ化する。
【００１６】
関係データにおいて、例えば、購買履歴における人の集合と商品の集合とのように、オブジェクトの集合（以降、ドメインと呼ぶ）が複数存在することが一般的であり、個々のドメインに対するクラスタリングを同時に行うことができる。個々のドメインに対するクラスタリングを同時に行うことを、特に、共クラスタリング(co-clustering)と呼び様々な研究がなされている。
【００１７】
このような従来の技術としては、例えば非特許文献１に記載されたものが知られている。非特許文献１で提案されているInfinite Relational Model（以降、ＩＲＭと呼ぶ）は、関係データの生成過程（generative process）を表すノンパラメトリックベイズモデル（non-parametric bayes model）であり、行列又は３次元以上のテンソル（tensor）の形式で表現される関係データに対して、関係の類似性に基づき共クラスタリングを行うことができる。関係データを共クラスタリングすると、各ドメインは複数のクラスタにクラスタ分割され、異なるドメインのクラスタの組み合わせごとに、ブロック状の領域（以降、クラスタブロックと呼ぶ）に分割される。個々のクラスタブロックは、関係の生成しやすさ（し難さ）が類似した単位と解釈することができる。例えば、人が商品を購入する場合における、人ごとの商品の購買履歴を表す関係データに対して、共クラスタリングを行うことで得られた個々のクラスタブロックを確認することで、特定の人のクラスタと特定のアイテムのクラスタとの間にある、購入しやすい、又は、購入しにくい、といった傾向を知ることができる。しかしながら、これらの手法においては、全てのクラスタブロックを確認しないと、どのクラスタブロックが重要であるかを知ることができない。そのため、クラスタブロックの数が非常に多い場合にどのクラスタブロックが注目すべき重要なクラスタブロックであるかを判断することが困難である。
【００１８】
上記の問題を解決する技術として、例えば、特許文献１に開示されたものが知られている。特許文献１に開示される技術では、関係データに対して関係の類似性に基づいた共クラスタリングを行い、入力された関係データを複数のクラスタブロックに分割する。分割する際には、複数のクラスタブロックについて、それぞれの統計量として相関強度を算出し、算出した相関強度を重要度とみなし、クラスタブロックの重要度の順序が分かるような態様で、複数のクラスタブロックを並び替えて表示する。
【００１９】
しかしながら、入力された関係データの性質によっては、算出された相関強度を、重要度と見なすことが適切ではない場合がある。
【００２０】
例えば、関係データ全体を一つのクラスタブロックとみなして算出された相関強度が高い値を取る場合、相関強度が低いクラスタブロックが重要度の高いクラスタブロックとなりうる。なぜならば、そのような場合においては、関係データ全体と異なる性質を有するクラスタブロック、つまり、相関強度が低いクラスタブロックが注目すべき重要なクラスタブロックであるからである。
【００２１】
一方、関係データ全体を一つのクラスタブロックとみなして算出された相関強度が低い値を取る場合、相関強度が高いクラスタブロックが重要度の高いクラスタブロックとなりうる。なぜならば、そのような場合においては、関係データ全体と異なる性質を有するクラスタブロック、つまり、相関強度が高いクラスタブロックが注目すべき重要なクラスタブロックであるためである。
【００２２】
上記の２つの場合、従来技術では、クラスタブロックの重要度を特定することが困難である。なぜならば、そのような状況では、関係データ全体の相関強度の値によって、各クラスタブロックの重要度が変化するため、ただ単に各クラスタブロックについて相関強度を計算するだけではクラスタブロックの重要度を特定することができないからである。
【００２３】
つまり、従来技術のように各クラスタブロックについて統計量を計算するだけでは、関係データ全体の分布の傾向によって各クラスタブロックの重要度が変化する状況において、クラスタブロックの重要度を特定することが困難であるという問題がある。
【００２４】
そこで、本発明は、クラスタブロックの重要度を特定することが可能な共クラスタリング装置を提供する。
【００２５】
つまり、本発明は、行列や３次元以上のテンソルの形式で表現される関係データに対して、関係データ全体の分布の傾向を考慮して、クラスタブロックの重要度を特定することが可能な共クラスタリング装置を提供することを目的とする。
【００２６】
上記目的を達成するために、本発明の一態様に係る共クラスタリング装置は、行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリング処理を行うことによって、前記関係データを複数のクラスタブロックに分割する共クラスタリング装置であって、前記複数のクラスタブロックのそれぞれにおいて生成される関係の傾向を示す統計量の、前記関係データ全体における分布傾向を生成する分布傾向生成部と、前記複数のクラスタブロックのそれぞれの統計量と、前記分布傾向生成部が生成した前記分布傾向とに基づいて、前記分布傾向に応じて算出結果が変化する算出方法により、前記複数のクラスタブロックのそれぞれに対して重要度を算出する算出部と、前記複数のクラスタブロックの少なくとも１つを示す情報と、当該少なくとも１つに対して前記算出部が算出した前記重要度を示す情報とを出力する出力部とを備える。
【００２７】
これにより、共クラスタリング装置は、行列又は３次元以上のテンソルの形式で表現される関係データに対して共クラスタリング処理を行う場合に、各クラスタブロックの統計量の分布傾向を考慮したクラスタブロックの重要度を出力する。ここで出力されるクラスタブロックの重要度は、関係データ全体の統計量と、当該クラスタブロックの統計量とを考慮した結果である。よって、たとえ当該クラスタブロック内の要素が同一であったとしても、関係データ全体の統計量が異なる場合には、異なる重要度が出力される。つまり、分布傾向を用いることで、入力された関係データ全体の傾向を考慮して、各クラスタブロックの重要度を算出することができるので、関係データの性質に応じたクラスタブロックの重要度を特定することができる。
【００２８】
例えば、前記分布傾向生成部は、前記関係データ全体の統計量を、前記分布傾向として生成する。
【００２９】
これにより、共クラスタリングを行うことで得られた各クラスタブロックの統計量を、それぞれ、共クラスタリングを行う前の関係データ全体の統計量と比較することができるため、各クラスタブロックが、入力された関係データにおいて、どれぐらい希少であるかを評価して、評価を重要度に反映させることができる。
【００３０】
例えば、前記算出部は、前記複数のクラスタブロックのそれぞれについて、当該クラスタブロックにおける前記分布傾向により示される値と、当該クラスタブロックの統計量との距離が大きいほど、大きくなるように前記重要度を算出する。
【００３１】
これにより、共クラスタリングを行うことで得られた各クラスタブロックの統計量を、それぞれ、共クラスタリングを行う前の関係データ全体の統計量と比較して、差異が大きいクラスタブロックほど、相対的に重要度が高いと判定することができる。
【００３２】
例えば、前記算出部は、前記複数のクラスタブロックのそれぞれについて、前記分布傾向と、当該クラスタブロックの統計量と、当該クラスタブロックの大きさとを用いて前記重要度を算出する。
【００３３】
これにより、関係データ全体を一つのクラスタブロックと見なした場合の統計量と、各クラスタブロックの統計量との比較に加えて、クラスタブロックの大きさを考慮した重要度を算出することができる。
【００３４】
例えば、前記分布傾向生成部は、前記複数のクラスタブロックのそれぞれの統計量を要素とする統計量データに対してクラスタリング処理を行うことによって、前記統計量データを複数のクラスタに分割し、分割することで得られた前記複数のクラスタの情報を、前記分布傾向として生成する。
【００３５】
これにより、各クラスタブロックの統計量の分布傾向が複雑な関係データにおいても、各クラスタブロックの統計量の分布傾向を考慮して、重要度の高いクラスタブロックを特定することができる。
【００３６】
例えば、前記算出部は、前記複数のクラスタのそれぞれについて、当該クラスタ内の要素の数が少ないほど、当該クラスタに要素として含まれるクラスタブロックに対して大きくなるように前記重要度を算出する。
【００３７】
これにより、各クラスタブロックの統計量の分布傾向が複雑な関係データにおいて、共クラスタリングを行うことで得られた各クラスタブロックが、それぞれ、入力された関係データにおいて、どれぐらい希少であるかを評価して、評価を重要度に反映させることができる。
【００３８】
例えば、前記算出部は、前記複数のクラスタのそれぞれについて、当該クラスタ内の要素の数と、前記複数のクラスタの要素に対応する１以上の前記クラスタブロックの大きさとに基づいて、当該クラスタに要素として含まれる前記クラスタブロックのそれぞれに対して前記重要度を算出する。
【００３９】
これにより、所属するクラスタブロックの数と、各クラスタブロックの統計量と、に加えて、クラスタブロックの大きさを考慮した重要度を算出することができる。
【００４０】
なお、本発明は、共クラスタリング装置として実現できるだけではなく、当該共クラスタリング装置を構成する処理手段をステップとする方法として実現することもできる。また、これらステップをコンピュータに実行させるプログラムとして実現してもよい。さらに、当該プログラムを記録したコンピュータ読み取り可能なＣＤ−ＲＯＭ（Ｃｏｍｐａｃｔ  Ｄｉｓｃ−Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）などの記録媒体、並びに、当該プログラムを示す情報、データ又は信号として実現してもよい。そして、それらプログラム、情報、データ及び信号は、インターネットなどの通信ネットワークを介して配信してもよい。
【００４１】
また、上記の各分類装置を構成する構成要素の一部又は全部は、１個のシステムＬＳＩ（Ｌａｒｇｅ  Ｓｃａｌｅ  Ｉｎｔｅｇｒａｔｉｏｎ：大規模集積回路）から構成されていてもよい。システムＬＳＩは、複数の構成部を１個のチップ上に集積して製造された超多機能ＬＳＩであり、具体的には、マイクロプロセッサ、ＲＯＭ及びＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）などを含んで構成されるコンピュータシステムである。
【００４２】
なお、これらの包括的または具体的な態様は、システム、方法、集積回路、コンピュータプログラムまたはコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラムまたは記録媒体の任意な組み合わせで実現されてもよい。
【００４３】
以下、実施の形態について、図面を参照しながら具体的に説明する。
【００４４】
なお、以下で説明する実施の形態は、いずれも包括的または具体的な例を示すものである。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、本発明を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、最上位概念を示す独立請求項に記載されていない構成要素については、任意の構成要素として説明される。
【００４５】
（実施の形態１）
まず、本発明の実施の形態１に係る共クラスタリング装置の概要について説明する。本発明の実施の形態１に係る共クラスタリング装置は、行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリング処理を行うことによって、前記関係データを複数のクラスタブロックに分割する共クラスタリング装置であって、前記複数のクラスタブロックのそれぞれにおいて生成される関係の傾向を示す統計量の、前記関係データ全体における分布傾向を生成する分布傾向生成部と、前記複数のクラスタブロックのそれぞれの統計量と、前記分布傾向生成部が生成した前記分布傾向とに基づいて、前記分布傾向に応じて算出結果が変化する算出方法により、前記複数のクラスタブロックのそれぞれに対して重要度を算出する算出部と、前記複数のクラスタブロックの少なくとも１つを示す情報と、当該少なくとも１つに対して前記算出部が算出した前記重要度を示す情報とを出力する出力部とを備える。
【００４６】
これにより、共クラスタリング装置は、行列又は３次元以上のテンソルの形式で表現される関係データに対して共クラスタリング処理を行う場合に、各クラスタブロックの統計量の分布傾向を考慮したクラスタブロックの重要度を出力する。ここで出力されるクラスタブロックの重要度は、関係データ全体の統計量と、当該クラスタブロックの統計量とを考慮した結果である。よって、たとえ当該クラスタブロック内の要素が同一であったとしても、関係データ全体の統計量が異なる場合には、異なる重要度が出力される。つまり、分布傾向を用いることで、入力された関係データ全体の傾向を考慮して、各クラスタブロックの重要度を算出することができるので、関係データの性質に応じたクラスタブロックの重要度を特定することができる。
【００４７】
以下では、まず、本実施の形態に係る共クラスタリング装置の構成について説明する。図１は、本実施の形態に係る共クラスタリング装置１００の構成の例を示すブロック図である。図１に示すように、本実施の形態に係る共クラスタリング装置１００は、データ入力部１１０と、共クラスタリング部１２０と、分布傾向生成部１３０と、算出部１４０と、出力部１５０とを備える。
【００４８】
データ入力部１１０は、行列又は３次元以上のテンソルの形式で表現される（表現可能な）関係データを共クラスタリング装置１００に入力する。なお、データ入力部１１０を介して入力される関係データは、ＨＤＤ（Hard disk drive）などの磁気ディスク装置又はメモリカードなどから読み出される構成であってもよいし、ユーザインターフェイスを介して入力されてもよいし、インターネット上のデータのうち、ユーザが探索して収集したデータが入力されてもよい。
【００４９】
ここで、関係データの定義について説明する。
【００５０】
関係データは、１つ以上のドメイン情報と、オブジェクト間関係情報とを含んでいる。なお、ドメイン情報は、ドメインを構成する複数のオブジェクトを識別する情報を含んでいる。例えば、インターネットショッピングサービスにおける購買履歴を表す関係データの例を考える。この場合、関係データには、「Ｔ１:ユーザ集合」と「Ｔ２:アイテム集合」との２つのドメインが含まれる。ユーザ集合は、当該インターネットショッピングサービスを利用可能なユーザの全集合を表す。また、アイテム集合は、当該インターネットショッピングサービスでユーザが購入することができるアイテムの全集合を表す。このとき、ユーザ集合に関するドメイン情報は、ユーザ集合に含まれる各ユーザを特定する情報のことを意味する。また、アイテム集合に関するドメイン情報は、アイテム集合に含まれる各アイテムを特定する情報のことを意味する。また、オブジェクト間関係情報は、オブジェクト同士の関係を表す情報である。例えば、購買履歴を表す関係データの場合、オブジェクト間関係情報は、「ユーザ集合」に含まれるある任意のユーザと、「アイテム集合」に含まれるある任意のアイテムの組に対して、「購入」又は「未購入」の二値の関係のうちのいずれであるかを特定可能にする情報である。この購買履歴の例における関係データの形式を、
【数１】
2014081899.tif 000003
と表す。この意味するところは、関係データＲは２つのドメイン情報Ｔ１とＴ２とを含み、オブジェクト間関係情報は、Ｔ１に含まれるオブジェクトと、Ｔ２に含まれるオブジェクトとの間に、二値の関係｛０，１｝が定義されていることを意味する。先に述べた購買履歴の例では、Ｔ１はユーザの集合、Ｔ２はアイテムの集合、二値の値｛０，１｝はそれぞれ購入と未購入とを表現している。今、Ｔ１がＮ１人のユーザで構成され、Ｔ２がＮ２個のアイテムで構成されている場合、関係データは、Ｎ１行Ｎ２列の行列の形で図示できる。
【００５１】
図２は、本実施の形態に係る関係データの例を示す図である。図２に示される関係データは、Ｎ１行Ｎ２列の行列の形で示される関係データの例である。図２の（ａ）は、ユーザと、当該ユーザが購入したアイテムとの対応を示す表である。また、図２の（ｂ）は、当該関係データを、Ｔ１（ユーザ集合、縦軸）及びＴ２（アイテム集合、横軸）に対して、白及び黒で示す図である。
【００５２】
今、ｉをＴ１に含まれるオブジェクトのインデックス、ｊをＴ２に含まれるオブジェクトのインデックスとすると、ｉ行ｊ列目の要素Ｒ（ｉ，ｊ）は、ｉ番目のユーザ
【数２】
2014081899.tif 000004
がｊ番目のアイテム
【数３】
2014081899.tif 000005
を購入したか否かを表している。図２の（ｂ）では、未購入であること（０）を白で表現し、購入したこと（１）を黒で表現している。
【００５３】
なお、関係データには様々なバリエーションがある。図３は、本実施の形態に係る関係データの他の例を示す図である。
【００５４】
図３の（ａ）は、複数の質問項目を有するアンケートに対して、ユーザが５段階で回答する場合のアンケート結果を示す図である。これは、ユーザ集合、及び、質問項目集合に対して、複数の段階の関係（アンケート回答）を有する関係データの例である。
【００５５】
図３の（ｂ）は、３つのドメインに対して、関係が多値になる場合の関係データの例である。
【００５６】
例えば、ソーシャルネットワーク（ＳＮＳ（Social network service））上の友人関係は、
【数４】
2014081899.tif 000006
で表される関係データである。また、関係が二値ではなく、多値の場合、
【数５】
2014081899.tif 000007
連続値の場合、
【数６】
2014081899.tif 000008
なども考えられる。さらに、３つ以上のドメインの関係を表す関係データ、
【数７】
2014081899.tif 000009
などを考える事もできる。この場合は、関係データは行列ではなく、行列を一般化した概念であるテンソルと見なせる。
【００５７】
このような、様々な関係データのバリエーションは全て、本発明の実施の形態１に係る共クラスタリング装置における関係データの範疇である。以降の説明では、簡単のため、２つのドメインの間の二値関係を表す関係データ
【数８】
2014081899.tif 000010
を具体例として説明するが、本発明がこれに限定されない。
【００５８】
以上が、関係データの定義についての説明である。
【００５９】
共クラスタリング部１２０は、関係データＲを入力として、共クラスタリングを行い、共クラスタリング結果としてクラスタブロック（又は、クラスタブロックを示す情報）を出力する。共クラスタリングとは、クラスタリングの一種であり、関係データに含まれる個々のドメインに対するクラスタリングを同時に行うことである。クラスタリング結果は、各ドメインに含まれるオブジェクトが所属するクラスタを特定する情報を少なくとも含む。具体的には、２つのドメインで構成される関係データ
【数９】
2014081899.tif 000011
の場合、共クラスタリング装置１００は、関係データＲに対して、関係の類似性に基づきＴ１のクラスタ割り当て
【数１０】
2014081899.tif 000012
と、Ｔ２のクラスタ割り当て
【数１１】
2014081899.tif 000013
とを求め、ｚ１及びｚ２をクラスタリング結果として出力する。ただし、
【数１２】
2014081899.tif 000014
は、Ｔ１に対するクラスタのカテゴリの集合であり、
【数１３】
2014081899.tif 000015
は、Ｔ２に対するクラスタカテゴリの集合である。
【００６０】
実際に共クラスタリングを実行するアルゴリズムは様々なものが挙げられるが、ここでは、非特許文献１として引用したＩＲＭを用いて、共クラスタリングを実行する手順を詳細に説明する。なお、ここで説明する共クラスタリングにより、図４の（ａ）に示される関係データが、図４の（ｂ）に示されるような共クラスタリング結果データとなる。
【００６１】
Ｋｅｍｐらによって提案されたＩＲＭは、関係データの生成過程を表現した確率モデルであり、関係データ
【数１４】
2014081899.tif 000016
の場合の生成モデル（generative model）は、（式１−１）〜（式１−４）で与えられる。
【００６２】
【数１５】
2014081899.tif 000017
【００６３】
ここで、ＣＲＰ（・）は中華料理店過程（Chinese Restaurant Process）を意味し、Ｂｅｔａ（・，・）はベータ分布を意味し、Ｂｅｒｎｏｕｌｌｉ（・，・）はベルヌーイ分布を意味する。また、γは中華料理店過程のパラメータ、βはベータ分布のパラメータを表す。
【００６４】
（式１−１）〜（式１−４）の生成モデルの簡単な説明をする。まず初めに、各ドメインごとにクラスタ割り当てが生成される（式１−１及び１−２）。次に、各クラスタブロック（ｋ，ｌ）に対して、クラスタブロック内に関係が生成される確率η（ｋ，ｌ）が、ベータ分布に従って生成される（式１−３）。最後に、関係データを構成する各関係Ｒ（ｉ，ｊ）は、オブジェクトｉが所属するクラスタ
【数１６】
2014081899.tif 000018
と、オブジェクトｊが所属するクラスタ
【数１７】
2014081899.tif 000019
との組で特定されるη（ｋ，ｌ）をパラメータとするベルヌーイ分布に従って生成される。
【００６５】
（式１−１）〜（式１−４）の生成モデルにおいて、関係データＲが生成される確率は、（式２）で計算される。
【００６６】
【数１８】
2014081899.tif 000020
【００６７】
ここで、ベータ分布は、ベルヌーイ分布の自然共役事前分布（natural conjugate prior distribution）なので、（式２）は、（式３）のように、η積分消去（integrated out）した形に書ける。
【００６８】
【数１９】
2014081899.tif 000021
【００６９】
今、クラスタ割り当てｚ１及びｚ２が得られている場合、（式３）を計算することにより、関係データＲが生成される確率を求めることができる。つまり、以下の最適化問題を解くことにより、共クラスタリング部１２０の出力であるクラスタ割り当てｚ１及びｚ２が得られる。
【００７０】
【数２０】
2014081899.tif 000022
【００７１】
（式４）を実際に解くために、様々な方法が提案されている。ここでは、一例として、ギブスサンプリング（Gibbs sampling）を用いた推定方法を説明する。ギブスサンプリングは、マルコフ連鎖モンテカルロ法（Markov Chain Monte Carlo method）と呼ばれる方法群の一つであり、適当な初期値から確率分布空間の探索を開始して、確率密度が高い箇所を推定することができる。つまり、（式４）の場合、ギブスサンプリングを用いることで、ｚ１及びｚ２を変数とする確率分布空間
【数２１】
2014081899.tif 000023
を探索して、尤度（likelihood）が最大となるときのｚ１及びｚ２の推定値を得ることができる。ここでは、理論的な説明は省略して結論のみを説明する。（式４）の問題を解くためのギブスサンプリングの手順は以下で与えられる。
（手順１）ｚ１及びｚ２の初期値を適当に決める。
（手順２）ｉ＝１，２，・・・，Ｎ１に対して以下の処理を行う。
（手順２−１）
【数２２】
2014081899.tif 000024
に従う確率で、
【数２３】
2014081899.tif 000025
の値を更新する。
（手順３）ｊ＝１，２，・・・，Ｎ２に対して以下の処理を行う。
（手順３−１）
【数２４】
2014081899.tif 000026
に従う確率で、
【数２５】
2014081899.tif 000027
の値を更新する。
（手順４）
【数２６】
2014081899.tif 000028
の値を計算し、収束していなければ（手順２）の処理を実行する。収束した場合は、終了する。ただし、
【数２７】
2014081899.tif 000029
である。ここで、
【数２８】
2014081899.tif 000030
は、現時点でドメインＴ１において、ｉ番目のオブジェクトを無視した条件下で、クラスタｋ＊に割り当てられているオブジェクトの数である。また、Ｌは、現時点でのドメインＴ２に関するクラスタの数である。
【００７２】
【数２９】
2014081899.tif 000031
は、関係データＲのｉ行目を無視して数えた、クラスタブロック（ｋ＊，ｌ）内のリンク（Ｒ（ｉ，ｊ）＝１）の数である。
【００７３】
【数３０】
2014081899.tif 000032
は、同様に数えた非リンク（Ｒ（ｉ，ｊ）＝０）の数である。
【００７４】
【数３１】
2014081899.tif 000033
は、関係データＲのｉ行目のクラスタ割り当て
【数３２】
2014081899.tif 000034
が、ｋ＊であると仮定して数えたリンクの数である。
【００７５】
【数３３】
2014081899.tif 000035
は同様にして数えた非リンクの数である。なお、
【数３４】
2014081899.tif 000036
は、同様に導出ができるために説明を省略する。 
【００７６】
以上の手続きにより、図４に示されるような関係データの共クラスタリングが実行される。なお、以上で説明した共クラスタリング手順は一例に過ぎず、本発明はこれに限定されるものではない。入力される関係データＲによっては、３つ以上のドメインを扱うような生成モデルであってもよいし、各ドメインに含まれるオブジェクトが所属するクラスタを特定する情報を少なくとも含む全く異なる共クラスタリング方式であってもよい。また、生成モデルの推定にギブスサンプリングを用いることも一例に過ぎず、本発明はこれに限定されるものではない。変分ベイズ法（Variational Bayes Inference）などの任意の生成モデルの推定法を用いてもよい。
【００７７】
分布傾向生成部１３０は、関係データＲを入力して、関係データＲを共クラスタリングして生成される複数クラスタブロックにおいて、個々のクラスタブロックを特徴付ける統計量の分布傾向情報を生成する。ここで、クラスタブロックを特徴付ける統計量とは、あるクラスタブロックに含まれる個々の関係が取る値の傾向を示す情報であり、例えば、クラスタブロック内の関係が取る値の平均若しくは分散などの数値、又は、クラスタブロック内の関係に対して任意の確率分布を当てはめて得られるパラメータを表す数値の集合を用いることができる。分布傾向情報は少なくとも、関係データＲを共クラスタリングして生成される複数クラスタブロックに対応する統計量がどのようにばらつくかを示す情報を含む。例えば、分布傾向情報の一例として、関係データＲ全体を一つのクラスタブロックと見なしたときの個々の関係の平均値が考えられる。例えば、２つのドメインに関する二値の関係データ
【数３５】
2014081899.tif 000037
の例を考えると、関係データＲ全体を一つのクラスタブロックと見なしたときの個々の関係の平均値は、
【数３６】
2014081899.tif 000038
で計算できる。この値は、二値の関係データにおいて、オブジェクト間の関係が１を取る割合を意味する。そのため、
【数３７】
2014081899.tif 000039
が０．０に近い場合、関係データＲはほとんどの関係の値が０である疎（sparse）なデータであるため、関係データＲを共クラスタリングして生成される複数クラスタブロックの統計量も０．０に近い値の付近に集まる可能性が高いことが分かる。一方で、
【数３８】
2014081899.tif 000040
が１．０に近い場合、関係データＲはほとんどの関係の値が１である密（dense)なデータであるため、関係データＲを共クラスタリングして生成される複数クラスタブロックの統計量も１．０に近い値の付近に集まる可能性が高いことが分かる。なお、ここでは例として、関係の平均値を分布傾向情報としたが、これはあくまでも一例であり本発明はこれに限定されるものではない。分布傾向情報は分散であってもよいし、その他の統計的な量、又は、その統計的な量の集合であってもよい。
【００７８】
算出部１４０は、関係データＲと、共クラスタリング結果ｚ１及びｚ２と、分布傾向情報とを入力として、個々のクラスタブロックの重要度情報を出力する。重要度情報は、個々のクラスタブロックがどれぐらい注目すべきであるかを表す数値であり、少なくとも分布傾向情報に応じて変化するという特徴を有する。例えば、分布傾向情報が関係データＲ全体を一つのクラスタブロックと見なしたときの個々の関係の平均値
【数３９】
2014081899.tif 000041
である場合、関係データＲと、共クラスタリング結果ｚ１及びｚ２とから個々のクラスタブロックの統計量
【数４０】
2014081899.tif 000042
を求めて、
【数４１】
2014081899.tif 000043
と
【数４２】
2014081899.tif 000044
とを引数とする関数
【数４３】
2014081899.tif 000045
によりクラスタブロック（ｋ，ｌ）の重要度を算出することが考えられる。個々のクラスタブロックの統計量
【数４４】
2014081899.tif 000046
は、例えば、クラスタブロック内の関係の平均値として、
【数４５】
2014081899.tif 000047
で計算し、関数Ｄ（・，・）はユークリッド距離を返す距離関数として、クラスタブロック（ｋ，ｌ）の重要度Ｉ（ｋ，ｌ）を
【数４６】
2014081899.tif 000048
で計算してもよい。
【００７９】
なお、本例では、分布傾向情報が、全体を一つのクラスタブロックと見なしたときの個々の関係の平均値
【数４７】
2014081899.tif 000049
である場合の例に合わせて、個々のクラスタブロックの統計量
【数４８】
2014081899.tif 000050
を、クラスタブロック内の関係の平均値としているが、これはあくまで一例であり、本発明はこれに限定されるものではない。例えば、分散であってもよいし、その他の任意の統計的指標であってもよい。また、本発明では、重要度Ｉ（ｋ，ｌ）を、
【数４９】
2014081899.tif 000051
と
【数５０】
2014081899.tif 000052
との間のユークリッド距離（Euclidean Distance）として定義しているが、これはあくまで一例であり、本発明はこれに限定されるものではない。重要度Ｉ（ｋ，ｌ）は、少なくとも、分布傾向情報と、個々のクラスタブロックの統計量とに依存するように算出される値であればよい。
【００８０】
出力部１５０は、関係データＲと、共クラスタリング結果ｚ１及びｚ２と、重要度情報とを入力として、クラスタブロックの重要度を示す情報を出力する。クラスタブロックの重要度を示す情報は、共クラスタリング部１２０により生成された複数のクラスタブロックのうちの少なくとも１つを示す情報と、当該クラスタブロックを示す情報のことである。例えば、各クラスタブロックの重要度、各クラスタブロックに含まれるオブジェクトを識別する情報の集合を出力する。なお、重要クラスタブロック情報を出力する先は、ＨＤＤ又はメモリカードなどのストレージであってもよいし、ネットワークを介して配信されてもよいし、モニタ等のディスプレイデバイスへ表示されてもよい。
【００８１】
次に、本実施の形態に係る共クラスタリング装置１００の動作の一例について説明する。図５は、本実施の形態に係る共クラスタリング装置１００の動作の例を示すフローチャートである。
【００８２】
まず、データ入力部１１０により、関係データが入力される（Ｓ１１０）。
【００８３】
次に、共クラスタリング部１２０は、入力された関係データに対して、共クラスタリングを行い、共クラスタリング結果を出力する（Ｓ１２０）。
【００８４】
次に、分布傾向生成部１３０は、関係データを入力して、分布傾向情報を生成する（Ｓ１３０）。
【００８５】
次に、算出部１４０は、関係データと、共クラスタリング結果と、分布傾向情報とを入力として、個々のクラスタブロックの重要度を出力する（Ｓ１４０）。
【００８６】
最後に、出力部１５０は、クラスタブロックの重要度を示す情報を出力する（Ｓ１５０）。
【００８７】
次に、共クラスタリング装置１００によるクラスタリング処理の例を示す。図６は、実施の形態１に係る共クラスタリング装置１００の処理の例を示す図である。
【００８８】
図６において、関係データ全体の統計量が比較的大きい場合（入力データが密（dense）な場合）の共クラスタリング装置１００の処理を、（ａ）から（ｅ）に示す。また、関係データ全体の統計量が比較的小さい場合（入力データが疎（sparce）な場合）の共クラスタリング装置１００の処理を、（ｋ）から（ｏ）に示す。
【００８９】
図６の（ａ）は、データ入力部１１０が入力する、全体の統計量が比較的大きい関係データを示す図である。二値の関係（｛０，１｝）を、黒及び白で示す。
【００９０】
図６の（ｂ）は、関係データを共クラスタリングした結果である。
【００９１】
図６の（ｃ）は、共クラスタリングした結果のデータにおける、各クラスタブロックの統計量を示す図である。ここでは、統計量は、クラスタブロック内の全要素の数に対する、二値の関係が１である要素の数の割合（充填率）として算出しており、各クラスタブロックに示されている。
【００９２】
図６の（ｄ）は、図６の（ｃ）に示される各クラスタブロックの統計量の、関係データ全体における分布傾向を示す図である。ここでは、関係データ全体における分布傾向は、各クラスタブロックの統計量の関係データ全体における平均値として算出している。
【００９３】
図６の（ｅ）は、各クラスタブロックの重要度を示す図である。ここでは、各クラスタブロックの重要度は、各クラスタブロックの統計量（図６の（ｃ））と、関係データ全体の統計量（図６の（ｃ））との差の絶対値として算出している。図６の（ｅ）において、重要度が最大であるクラスタブロックは、クラスタブロック６０１であることがわかる。つまり、関係データ全体の統計量が比較的大きい場合には、統計量が比較的小さいクラスタブロックに対して大きい重要度が算出される。
【００９４】
図６の（ｋ）は、データ入力部１１０が入力する、全体の統計量が比較的小さい関係データを示す図である。図６の（ｌ）から（ｏ）までは、それぞれ、上記の図６の（ｂ）から（ｅ）に対応する。
【００９５】
図６の（ｏ）において、重要度が最大であるクラスタブロックは、クラスタブロック６０２であることがわかる。つまり、関係データ全体の統計量が比較的小さい場合には、統計量が比較的大きいクラスタブロックに対して大きい重要度が算出される。
【００９６】
以上のように、共クラスタリング装置１００は、関係データ全体の統計量に基づいて、各クラスタブロックの重要度を算出し出力する。なお、この算出方法は、関係データ全体の統計量に応じて算出結果が変化すると表現することもできる。関係データ全体の統計量に応じて算出結果が変化するので、たとえ、当該クラスタブロック内の要素が同一であったとしても、関係データ全体の統計量が異なる場合には、異なる重要度が出力される。
【００９７】
以上のように、本実施の形態に係る共クラスタリング装置は、行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリング処理を行うことによって、前記関係データを複数のクラスタブロックに分割する共クラスタリング装置であって、前記複数のクラスタブロックのそれぞれにおいて生成される関係の傾向を示す統計量の、前記関係データ全体における分布傾向を生成する分布傾向生成部と、前記複数のクラスタブロックのそれぞれの統計量と、前記分布傾向生成部が生成した前記分布傾向とに基づいて、前記分布傾向に応じて算出結果が変化する算出方法により、前記複数のクラスタブロックのそれぞれに対して重要度を算出する算出部と、前記複数のクラスタブロックの少なくとも１つを示す情報と、当該少なくとも１つに対して前記算出部が算出した前記重要度を示す情報とを出力する出力部とを備える。
【００９８】
これにより、共クラスタリング装置は、行列又は３次元以上のテンソルの形式で表現される関係データに対して共クラスタリング処理を行う場合に、各クラスタブロックの統計量の分布傾向を考慮したクラスタブロックの重要度を出力する。ここで出力されるクラスタブロックの重要度は、関係データ全体の統計量と、当該クラスタブロックの統計量とを考慮した結果である。よって、たとえ当該クラスタブロック内の要素が同一であったとしても、関係データ全体の統計量が異なる場合には、異なる重要度が出力される。つまり、分布傾向を用いることで、入力された関係データ全体の傾向を考慮して、各クラスタブロックの重要度を算出することができるので、関係データの性質に応じたクラスタブロックの重要度を特定することができる。
【００９９】
本実施の形態に係る共クラスタリング装置は、様々な用途が考えられる。例えば、本実施の形態に係る共クラスタリング装置は、関係データを分析するためのソフトウェアとして、実装されることが考えられる。具体的には、ソーシャルネットワークサービス上の人間関係の分析、インターネットショッピングにおける商品購買履歴やコンテンツ配信サービスにおけるコンテンツ視聴履歴からの嗜好若しくは傾向の分析、又は、バイオ技術分野における関係の分析などの用途が考えられる。また、本実施の形態に係る共クラスタリング装置をシステムの一部に組み込むことで、レコメンド（recommendation）などのサービスを実現することも考えられる。
【０１００】
なお、本実施の形態に係る共クラスタリング装置１００において、算出部１４０は、重要度を算出するために個々のクラスタブロックの大きさを表す情報を加味してもよい。例えば、共クラスタリング結果ｚ１及びｚ２から各クラスタブロックの面積が分かるので、
【数５１】
2014081899.tif 000053
が同じ、あるいは近接しているクラスタブロックが複数存在した場合に、クラスタブロックの面積が大きい場合に、重要度Ｉ（ｋ，ｌ）が大きくなるように算出する。これにより、より多くのオブジェクトが所属するクラスタブロックが相対的に高い重要度を示すようになる。
【０１０１】
また、本実施の形態では、共クラスタリング部１２０と、分布傾向生成部１３０と、算出部１４０とはそれぞれ、入力と、処理と、出力とが、定義された独立な手続き（アルゴリズム）であったが、必ずしもそれぞれの機能ブロック（構成要素）が独立したアルゴリズムでなくてもよい。例えば、関係データの生成モデルとして例示したＩＲＭを拡張し、分布傾向生成部１３０と、算出部１４０とに相当する構成を生成モデルのレベルで含めてもよい。このような構成の共クラスタリング装置について、本実施の形態の変形例において詳細に説明する。
【０１０２】
（実施の形態１の変形例）
本発明の実施の形態１の変形例について説明する。
【０１０３】
図７は、本実施の形態に係る共クラスタリング装置１００Ａの構成を示すブロック図である。図７に示すように、共クラスタリング装置１００Ａは、データ入力部１１０と、共クラスタリング部１２０Ａと、出力部１５０とを備える。共クラスタリング部１２０Ａは、その内部機能として、分布傾向生成部１３０と、算出部１４０とを有する。
【０１０４】
データ入力部１１０と、出力部１５０とは、共クラスタリング装置１００と同様であるので、説明を省略する。
【０１０５】
共クラスタリング部１２０Ａは、関係データＲを入力として共クラスタリングを行い、共クラスタリング結果を出力する。また、共クラスタリングを行うと同時に、又は、平行して、分布傾向生成部１３０がクラスタブロックの統計量の分布傾向を生成し、算出部１４０が、クラスタブロックの重要度を算出する。つまり、共クラスタリング処理と重要度算出処理とを、同時に、又は、平行して、実行することができる。
【０１０６】
具体的には、関係データ
【数５２】
2014081899.tif 000054
の場合、例えば、以下のような生成モデルを考えることができる。
【０１０７】
【数５３】
2014081899.tif 000055
【０１０８】
（式９−１）〜（式９−６）の生成モデルの簡単な説明をする。まず初めに、ＩＲＭと同様に、各ドメインに毎にクラスタ割り当てが生成し（式９−１及び９−２）、各クラスタブロック（ｋ，ｌ）に対して、クラスタブロックの重要度Ｉ（ｋ，ｌ）が、ベータ分布に従って生成される（式９−３）。次に、関係データ全体を一つのクラスタブロックと見なして、関係データ全体に渡る関係生成確率η０を、ベータ分布に従って生成する（式９−４）。次に、クラスタブロック固有の関係生成確率η（ｋ，ｌ）が、クラスタブロックの重要度Ｉ（ｋ，ｌ）と、関係データ全体に渡る関係生成確率η０とから計算される（式９−５）。ここでσは混合率（mixture rate）を表す値であり、０より大きく１以下の所定の値をとる。最後に、関係生成確率η（ｋ，ｌ）をパラメータとするベルヌーイ分布に従って、関係データを構成する各関係Ｒ（ｉ，ｊ）が生成される（式９−６）。
【０１０９】
上記の生成モデルにおいて、関係データＲが生成される確率は、
【数５４】
2014081899.tif 000056
で計算される。つまり、ＩＲＭのときに説明したように、ギブスサンプリングや変分ベイズ法などの任意のパラメータ推定方法を使うことで、未知のパラメータである、ｚ１、ｚ２、Ｉ、η０、σを推定することが可能である。ここで、（式９−１）と（式９−２）とは、クラスタ割り当てｚ１及びｚ２を未知のパラメータとして生成モデルに組み込む役割を果たしている。また、（式９−６）は、クラスタ割り当てｚ１及びｚ２の結果に依存して、関係データＲが生成されることを表している。つまり、（式９−１）、（式９−２）及び（式９−６）の式は、本実施の形態における共クラスタリング部１２０に相当する。また、η０は関係データ全体に渡る関係生成確率であったことに注意すると、η０は分布傾向情報の一例と見なせる。つまり、（式９−４）は、本実施の形態に係る共クラスタリング装置における分布傾向生成部１３０に相当することが分かる。また、（式９−５）は、クラスタブロックの重要度Ｉ（ｋ，ｌ）と、関係データ全体に渡る関係生成確率η０とを用いて、クラスタブロック固有の関係生成確率η（ｋ，ｌ）を算出する式になっていることに注意すると、（式９−５）は、
【数５５】
2014081899.tif 000057
と等価である。この式は、クラスタブロック固有の関係生成確率η（ｋ，ｌ）と、関係データ全体に渡る関係生成確率η０とを用いて、クラスタブロックの重要度Ｉ（ｋ，ｌ）を算出していると見なすことができるため、（式９−５）は、本実施の形態に係る共クラスタリング装置における算出部１４０に相当する。
【０１１０】
以上のように、（式９−１）〜（式９−６）は、生成モデルのレベルで、本実施の形態に係る共クラスタリング装置１００を構成する構成要素が含まれていることになる。
【０１１１】
（実施の形態２）
次に、本実施の形態２に係る共クラスタリング装置２００の概要について説明する。本実施の形態に係る共クラスタリング装置では、前記分布傾向生成部は、前記複数のクラスタブロックのそれぞれの統計量を要素とする統計量データに対してクラスタリング処理を行うことによって、前記統計量データを複数のクラスタに分割し、分割することで得られた前記複数のクラスタの情報を、前記分布傾向として生成する。
【０１１２】
これにより、各クラスタブロックの統計量の分布傾向が複雑な関係データにおいても、各クラスタブロックの統計量の分布傾向を考慮して、重要度の高いクラスタブロックを特定することができる。
【０１１３】
図８は、本実施の形態に係る共クラスタリング装置２００の構成の例を示すブロック図である。図８に示すように、本実施の形態に係る共クラスタリング装置２００は、共クラスタリング装置１００（図１）における分布傾向生成部１３０と算出部１４０との代わりに、分布傾向生成部２３０と、算出部２４０とを備える点が異なる。以下では、本実施の形態に係る共クラスタリング装置１００と異なる点を説明し、同じ点は説明を省略する。
【０１１４】
分布傾向生成部２３０は、関係データＲと、共クラスタリング結果ｚ１及びｚ２とを入力して、複数のクラスタブロックにおいて、個々のクラスタブロックを特徴付ける統計量が類似したクラスタブロックをクラスタリングによりグループ化して、グループ化した結果を分布傾向情報として生成する。個々でグループ化した結果は、個々のクラスタブロックがどのグループに所属するかを示す情報である。つまり、実施の形態１における関係データを共クラスタリングして得られるクラスタブロックのそれぞれの統計量を要素とする統計量データをクラスタリングすることで、関係データの全体傾向を得る。
【０１１５】
例えば、実施の形態１に係る共クラスタリング装置１００の説明のときと同様に、２つのドメインに関する二値の関係データ
【数５６】
2014081899.tif 000058
の例を考える。今、個々のクラスタブロックを特徴付ける統計量が、関係の値の平均である（式７）であり、共クラスタリング結果ｚ１及びｚ２にそれぞれ、Ｋ個及びＬ個のクラスタが含まれる場合、分布傾向生成部２３０は、Ｋ×Ｌ個のクラスタブロックを、個々のクラスタブロックの統計量
【数５７】
2014081899.tif 000059
の類似性に基づいて任意の数Ｍ（＜Ｋ×Ｌ）個のグループにクラスタリングする。このクラスタリングには、ｋ−ｍｅａｎｓなどの著名なクラスタリングアルゴリズムを用いてもよいし、所定の閾値を設けて、クラスタブロックの統計量
【数５８】
2014081899.tif 000060
が上記所定の閾値の範囲内であるときに同じグループにする、といった単純な方法でもよい。以上により、分布傾向生成部２３０は、Ｋ×Ｌ個のクラスタブロックのそれぞれがどのグループに所属するかを示す情報を、グループ化した結果として出力する。
【０１１６】
算出部２４０は、グループ化した結果を入力として、グループ化した結果に応じて変化するように、個々のクラスタブロックの重要度を算出する。例えば、グループ化した結果において、同一のグループに所属するクラスタブロックの数が小さいクラスタブロックほど、相対的に高い値になるように重要度を算出することが考えられる。具体的には、今、Ｋ×Ｌ個のクラスタブロックの、Ｍ（＜Ｋ×Ｌ）個のグループへのクラスタ割り当てが、
【数５９】
2014081899.tif 000061
であったとき、個々のクラスタブロック（ｋ，ｌ）に対応する同一のグループに所属するクラスタブロックの数Δ（ｋ，ｌ）は、（式１２）で計算できる。
【０１１７】
【数６０】
2014081899.tif 000062
【０１１８】
ただし、（式１２）において、δ（・）は、括弧の中の式の評価結果が真（ｔｒｕｅ）の時に１、偽（ｆａｌｓｅ)の時に０を返す関数である。そして、個々のクラスタブロックの重要度Ｉ（ｋ，ｌ）を、（式１３）で算出する。
【０１１９】
【数６１】
2014081899.tif 000063
【０１２０】
（式１３）により重要度を算出することで、個々のクラスタブロックの重要度Ｉ（ｋ，ｌ）は、統計量
【数６２】
2014081899.tif 000064
が自分と類似しているクラスタブロックの数が少ないほど、大きな値を示す。つまり、希少なクラスタブロックであるほど、相対的に重要度が高くなる。このように個々のクラスタブロックの重要度を計算することで、実施の形態１に係る共クラスタリング装置１００の場合のように、分布傾向情報を唯一の統計量
【数６３】
2014081899.tif 000065
で表すことができない複雑な関係データが与えられた場合にでも、希少性が高い重要なクラスタブロックを特定することができる。
【０１２１】
以上のように、本実施の形態に係る共クラスタリング装置では、前記分布傾向生成部は、前記複数のクラスタブロックのそれぞれの統計量を要素とする統計量データに対してクラスタリング処理を行うことによって、前記統計量データを複数のクラスタに分割し、分割することで得られた前記複数のクラスタの情報を、前記分布傾向として生成する。
【０１２２】
これにより、各クラスタブロックの統計量の分布傾向が複雑な関係データにおいても、各クラスタブロックの統計量の分布傾向を考慮して、重要度の高いクラスタブロックを特定することができる。
【０１２３】
本実施の形態に係る共クラスタリング装置は、様々な用途が考えられる。最も基本的な例としては、関係データを分析するためのソフトウェアとして、本実施の形態に係る共クラスタリング装置を実装することが考えられる。具体的には、ソーシャルネットワークサービス上の人間関係の分析、インターネットショッピングにおける商品購買履歴やコンテンツ配信サービスにおけるコンテンツ視聴履歴からの嗜好や傾向の分析、バイオ技術分野における関係の分析、等の用途が考えられる。また、本実施の形態に係る共クラスタリング装置をシステムの一部に組み込むことで、レコメンド（recommendation）などのサービスを実現することも考えられる。
【０１２４】
なお、本実施の形態に係る共クラスタリング装置２００によれば、算出部２４０は、重要度を算出するために個々のクラスタブロックの大きさを表す情報を加味してもよい。例えば、共クラスタリング結果ｚ１及びｚ２から各クラスタブロックの面積が分かるので、
【数６４】
2014081899.tif 000066
が同じ、あるいは近接しているクラスタブロックが複数存在した場合に、同一のグループに所属するクラスタブロックの面積の和が大きいほど、重要度Ｉ（ｋ，ｌ）が大きくなるように（式１３）を補正した式で重要度Ｉ（ｋ，ｌ）を算出する。これにより、面積の大きなクラスタブロックが所属するグループが相対的に高い重要度を示すようになる。
【０１２５】
また、本実施の形態では、共クラスタリング部１２０と、分布傾向生成部２３０と、算出部２４０とはそれぞれ、入力と、処理と、出力とが、定義された独立な手続き（アルゴリズム）であったが、必ずしもそれぞれの機能ブロック（構成要素）が独立したアルゴリズムでなくてもよい。 
【０１２６】
例えば、関係データの生成モデルとして例示したＩＲＭを拡張し、分布傾向生成部２３０又は算出部２４０に相当する構成の一部又は全部を生成モデルのレベルで含めてもよい。
【０１２７】
具体的には、関係データ
【数６５】
2014081899.tif 000067
の場合、例えば、以下のような分布傾向生成部を含めた生成モデルを考えることができる。
【０１２８】
【数６６】
2014081899.tif 000068
【０１２９】
（式１４−１）〜（式１４−６）の生成モデルの簡単な説明をする。まず初めに、ＩＲＭと同様に、各ドメインごとにクラスタ割り当てを生成する（式１４−１及び１４−２）。次に、各クラスタブロック（ｋ，ｌ）をグループ化した結果ｚＣＢを生成する（式１４−３）。次に、それぞれのクラスタブロックのグループｕに対して、固有の関係生成確率θｕが生成される（式１４−４）。次に、各クラスタブロック固有の関係生成確率η（ｋ，ｌ）が、自分の所属するグループ
【数６７】
2014081899.tif 000069
に依存して、θの中から選ばれる（式１４−５）。最後に、関係生成確率η（ｋ，ｌ）をパラメータとするベルヌーイ分布に従って、関係データを構成する各関係Ｒ（ｉ，ｊ）が生成される（式１４−６）。
【０１３０】
上記の生成モデルにおいて、関係データＲが生成される確率は、
【数６８】
2014081899.tif 000070
で計算される。つまり、ＩＲＭのときに説明したように、ギブスサンプリング又は変分ベイズ法などの任意のパラメータ推定方法を使うことで、未知のパラメータである、ｚ１、ｚ２、ｚＣＢ、ηを推定することが可能である。ここで、（式１４−１）と（式１４−２）とは、クラスタ割り当てｚ１及びｚ２を未知のパラメータとして生成モデルに組み込む役割を果たしている。また、（式１４−６）は、クラスタ割り当てｚ１及びｚ２の結果に依存して、関係データＲが生成されることを表している。つまり、（式９−１）、（式９−２）及び（式９−６）の各式は、本実施の形態における共クラスタリング部１２０に相当する。また、ｚＣＢは、クラスタ割り当てｚ１及びｚ２により特定されるＫ×Ｌ個のクラスタブロックを、Ｍ（＜Ｋ×Ｌ）個のグループへクラスタリングすることを表しているため、分布傾向情報の一例と見なせる。つまり、（式１４−４）は、本実施の形態に係る共クラスタリング装置における分布傾向生成部２３０に相当することが分かる。以上に説明したモデルの未知のパラメータを推定することで、共クラスタリング部１２０の出力である共クラスタリング結果ｚ１及びｚ２と、分布傾向生成部２３０の出力である分布傾向情報ｚＣＢを同時に取得することができる。ｚ１、ｚ２、ｚＣＢが得られた場合、（式１２）及び（式１３）を用いることで重要度も計算することができる。
【０１３１】
以上のように、（式１４−１）〜（式１４−６）には、生成モデルのレベルで、本実施の形態に係る共クラスタリング装置２００を構成する構成要素が含まれていることになる。
【０１３２】
（その他の変形例）
なお、上記の実施の形態の共クラスタリング装置は、典型的には半導体集積回路であるＬＳＩとして実現される。これらは個別に１チップ化されてもよいし、一部またはすべてを含むように１チップ化されてもよい。ここではＬＳＩとしたが、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。
【０１３３】
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサで実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。
【０１３４】
さらには、半導体技術の進歩または派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。バイオ技術の適応などが可能性として有り得る。
【０１３５】
さらに加えて、本実施の形態の共クラスタリング装置を集積化した半導体チップと、画像を描画するためのディスプレイとを組み合せて、様々な用途に応じた描画機器を構成することができる。携帯電話やテレビ、デジタルビデオレコーダ、デジタルビデオカメラ、カーナビゲーション等における情報描画手段として、本発明を利用することが可能である。ディスプレイとしては、ブラウン管（ＣＲＴ）の他、液晶やＰＤＰ（プラズマディスプレイパネル）、有機ＥＬなどのフラットディスプレイ、プロジェクターを代表とする投射型ディスプレイなどと組み合わせることが可能である。
【０１３６】
なお、上記各実施の形態において、各構成要素は、専用のハードウェア（電子回路）で構成されるか、各構成要素に適したソフトウェアプログラムを実行することによって実現されてもよい。各構成要素は、ＣＰＵまたはプロセッサなどのプログラム実行部が、ハードディスクまたは半導体メモリなどの記録媒体に記録されたソフトウェアプログラムを読み出して実行することによって実現されてもよい。ここで、上記各実施の形態の共クラスタリング装置などを実現するソフトウェアは、次のようなプログラムである。
【０１３７】
すなわち、このプログラムは、コンピュータに、行列又は３次元以上のテンソルの形式で表現可能な関係データに対して共クラスタリング処理を行うことによって、前記関係データを複数のクラスタブロックに分割する共クラスタリング装置における共クラスタリング方法であって、前記複数のクラスタブロックのそれぞれにおいて生成される関係の傾向を示す統計量の、前記関係データ全体における分布傾向を生成する分布傾向生成ステップと、前記複数のクラスタブロックのそれぞれの統計量と、前記分布傾向生成ステップで生成した前記分布傾向とに基づいて、前記分布傾向に応じて算出結果が変化する算出方法により、前記複数のクラスタブロックのそれぞれに対して重要度を算出する算出ステップと、前記複数のクラスタブロックの少なくとも１つを示す情報と、当該少なくとも１つに対して前記算出ステップで算出した前記重要度を示す情報とを出力する出力ステップとを含む共クラスタリング方法を実行させる。
【０１３８】
以上、一つまたは複数の態様に係る共クラスタリング装置について、実施の形態に基づいて説明したが、本発明は、この実施の形態に限定されるものではない。本発明の趣旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したものや、異なる実施の形態における構成要素を組み合わせて構築される形態も、一つまたは複数の態様の範囲内に含まれてもよい。
【産業上の利用可能性】
【０１３９】
本発明の共クラスタリング装置は、様々な用途に利用可能である。例えば、携帯電話、携帯音楽プレーヤー、デジタルカメラ若しくはデジタルビデオカメラ等の電池駆動の携帯表示端末、又は、テレビ、デジタルビデオレコーダ若しくはカーナビゲーション等の高解像度の情報表示機器におけるメニュー表示、又は、Ｗｅｂブラウザ、エディタ、ＥＰＧ若しくは地図表示等における情報表示手段として利用価値が高い。
【符号の説明】
【０１４０】
１００、１００Ａ、２００    共クラスタリング装置
１１０    データ入力部
１２０、１２０Ａ    共クラスタリング部
１３０、１３０Ａ、２３０    分布傾向生成部
１４０、１４０Ａ、２４０    算出部
１５０    出力部
【図１】
2014081899.tif 000071
【図２】
2014081899.tif 000072
【図３】
2014081899.tif 000073
【図４】
2014081899.tif 000074
【図５】
2014081899.tif 000075
【図６】
2014081899.tif 000076
【図７】
2014081899.tif 000077
【図８】
2014081899.tif 000078