特許 6973575 情報処理装置、情報処理方法、およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B1)

(11)【特許番号】6973575

(24)【登録日】2021年11月8日

(45)【発行日】2021年12月1日

(54)【発明の名称】情報処理装置、情報処理方法、およびプログラム

(51)【国際特許分類】

   G06Q 10/04 20120101AFI20211118BHJP

   G06N 20/00 20190101ALI20211118BHJP

【ＦＩ】

   G06Q10/04

   G06N20/00

【請求項の数】6

【全頁数】27

(21)【出願番号】特願2020-118696(P2020-118696)

(22)【出願日】2020年7月9日

【審査請求日】2020年7月9日

(73)【特許権者】

【識別番号】000004237

【氏名又は名称】日本電気株式会社

(74)【代理人】

【識別番号】100110928

【弁理士】

【氏名又は名称】速水 進治

(72)【発明者】

【氏名】中川 淳子

【審査官】 阿部 潤

(56)【参考文献】

【文献】 特開２０１２−１１７８２９（ＪＰ，Ａ）

【文献】 特開２００５−１２８８０８（ＪＰ，Ａ）

【文献】 特開２０１５−１８７７７３（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０６Ｑ １０／００ − ９９／００

Ｇ１６Ｈ １０／００ − ８０／００

Ｇ１６Ｚ ９９／００

Ｇ０６Ｎ ２０／００

(57)【特許請求の範囲】

【請求項1】

過去に発生した事象の履歴データのうち第１の期間に含まれる第１の履歴データに基づいて複数の計算手法それぞれについて計算される対象エリアの各セルにおける事象の発生し易さと、前記第１の期間よりも後の第２の期間に含まれる第２の履歴データとに基づいて、特徴量集合を前記計算手法別に生成する手法別特徴量生成手段と、
前記計算手法別に生成された特徴量集合を各次元のラベルとして有し、かつ、各行列要素に対応するセルにおける前記第２の履歴データに基づく前記事象の発生件数を行列要素として有する第１行列を生成する第１行列生成手段と、
前記第１行列と同じラベルを有し、かつ、各行列要素に対応する各次元のラベル値の合計値を行列要素として有する第２行列を生成する第２行列生成手段と、
所定の閾値に基づいて、前記第２行列において行列要素の境界を設定する境界設定手段と、
前記境界の位置と対応する位置に存在する前記第１行列の行列要素である境界行列要素のそれぞれについて、前記境界行列要素を基準として前記第１行列を複数の領域に分割し、前記複数の領域の少なくとも１つの領域の行列要素の和を用いた評価指標を生成する評価指標生成手段と、
前記境界行列要素それぞれについて生成された評価指標に基づいて、評価が最大となる前記境界行列要素に対応する各次元のラベル値の組み合わせを、前記事象の予測に用いる閾値組み合わせとして特定する閾値組み合わせ特定手段と、
前記閾値組み合わせに基づいて、所定の予測時点における前記対象エリアの各セルの事象の発生し易さの予測結果を分類し、その分類結果を含む情報を出力する予測結果出力手段と、
を備える情報処理装置。

【請求項2】

前記予測結果出力手段は、前記計算手法別に計算した前記結果を分類する処理を、前記事象の種類ごとに実行する、
請求項１に記載の情報処理装置。

【請求項3】

前記評価指標の可視化情報を表示装置に出力する評価指標可視化手段を更に備える、
請求項１または２に記載の情報処理装置。

【請求項4】

前記評価指標可視化手段は、前記評価指標の大きさに基づいて、当該評価指標の表示態様を変更する、
請求項３に記載の情報処理装置。

【請求項5】

コンピュータが、
過去に発生した事象の履歴データのうち第１の期間に含まれる第１の履歴データに基づいて複数の計算手法それぞれについて計算される対象エリアの各セルにおける事象の発生し易さと、前記第１の期間よりも後の第２の期間に含まれる第２の履歴データとに基づいて、特徴量集合を前記計算手法別に生成し、
前記計算手法別に生成された特徴量集合を各次元のラベルとして有し、かつ、各行列要素に対応するセルにおける前記第２の履歴データに基づく前記事象の発生件数を行列要素として有する第１行列を生成し、
前記第１行列と同じラベルを有し、かつ、各行列要素に対応する各次元のラベル値の合計値を行列要素として有する第２行列を生成し、
所定の閾値に基づいて、前記第２行列において行列要素の境界を設定し、
前記境界の位置と対応する位置に存在する前記第１行列の行列要素である境界行列要素のそれぞれについて、前記境界行列要素を基準として前記第１行列を複数の領域に分割し、前記複数の領域の少なくとも１つの領域の行列要素の和を用いた評価指標を生成し、
前記境界行列要素それぞれについて生成された評価指標に基づいて、評価が最大となる前記境界行列要素に対応する各次元のラベル値の組み合わせを、前記事象の予測に用いる閾値組み合わせとして特定し、
前記閾値組み合わせに基づいて、所定の予測時点における前記対象エリアの各セルの事象の発生し易さの予測結果を分類し、その分類結果を含む情報を出力する、
ことを含む情報処理方法。

【請求項6】

コンピュータに、請求項５に記載の情報処理方法を実行させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、情報処理装置、情報処理方法、およびプログラムに関する。

【背景技術】

【0002】

犯罪や疾病などの事象について、過去の発生履歴に関するデータを用いて、ある時点での当該事象の発生を予測する技術がある。例えば、このような技術の一例が、下記特許文献１乃至３、非特許文献１乃至４に開示されている。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特許第６６２７８４２号公報

【特許文献2】米国特許第９１２９２１９号明細書

【特許文献3】米国特許出願公開第２０１５／０３７９４１３号明細書

【非特許文献】

【0004】

【非特許文献1】Perry, W.L., Mclnnis, B., Price, C.C., Smith, S.C., Hollywood, J.S.: Predictive Policing: The Role of Crime Forecasting in Law Enforcement Operations, RAND Corporation (2013), available from < https://www.rand.org/content/dam/rand/pubs/research_reports/RR200/RR233/RAND_RR233.pdf> (accessed 2019-07-02).

【非特許文献2】雨宮護：潜在成長曲線を用いた地域レベルでの犯罪の時系列変化と地区環境との関連の分析，都市計画論文集, Vol.48, No.3, pp.351-356 (2013).

【非特許文献3】野貴泰，糸井川栄一：犯罪多発地点の予測に基づく防犯パトロール経路に関する提案，地域安全学会論文集, No.31, pp.195-204 (2017).

【非特許文献4】Huang, C., Zhang, J., Zheng, Y., Chawla, N.V.,: DeepCrime: Attentive Hierarchical Recurrent Networks for Crime Prediction, Proc. CIKM'18, pp.1423-1432, ACM (2018).

【発明の概要】

【発明が解決しようとする課題】

【0005】

上述の文献にあるような技術において、コンピュータによる予測がどのような判断基準でなされたか、すなわち、コンピュータによる予測の根拠（説明性）に対するニーズが高まっている。特許文献１に開示される技術において、事象を予測する場合の制約（例えば、人的リソースに基づく制約など）を考慮した予測を行うことができるが、この技術において、予測の根拠を説明する機能が備わっていない。

【0006】

本発明は、上記の課題に鑑みてなされたものである。本発明の目的の一つは、ある時点での事象の発生をコンピュータで予測するとともに、当該コンピュータによる予測の根拠を説明可能とする技術を提供することである。

【課題を解決するための手段】

【0007】

本発明の情報処理装置は、
過去に発生した事象の履歴データのうち第１の期間に含まれる第１の履歴データに基づいて複数の計算手法それぞれについて計算される対象エリアの各セルにおける事象の発生し易さと、前記第１の期間よりも後の第２の期間に含まれる第２の履歴データとに基づいて、特徴量集合を前記計算手法別に生成する手法別特徴量生成手段と、
前記計算手法別に生成された特徴量集合を各次元のラベルとして有し、かつ、各行列要素に対応するセルにおける前記第２の履歴データに基づく前記事象の発生件数を行列要素として有する第１行列を生成する第１行列生成手段と、
前記第１行列と同じラベルを有し、かつ、各行列要素に対応する各次元のラベル値の合計値を行列要素として有する第２行列を生成する第２行列生成手段と、
所定の閾値に基づいて、前記第２行列において行列要素の境界を設定する境界設定手段と、
前記境界の位置と対応する位置に存在する前記第１行列の行列要素である境界行列要素のそれぞれについて、前記境界行列要素を基準として前記第１行列を複数の領域に分割し、前記複数の領域の少なくとも１つの領域の行列要素の和を用いた評価指標を生成する評価指標生成手段と、
前記境界行列要素それぞれについて生成された評価指標に基づいて、評価が最大となる前記境界行列要素に対応する各次元のラベル値の組み合わせを、前記事象の予測に用いる閾値組み合わせとして特定する閾値組み合わせ特定手段と、
前記閾値組み合わせに基づいて、所定の予測時点における前記対象エリアの各セルの事象の発生し易さの予測結果を分類し、その分類結果を含む情報を出力する予測結果出力手段と、
を備える。

【0008】

本発明の情報処理方法は、
コンピュータが、
過去に発生した事象の履歴データのうち第１の期間に含まれる第１の履歴データに基づいて複数の計算手法それぞれについて計算される対象エリアの各セルにおける事象の発生し易さと、前記第１の期間よりも後の第２の期間に含まれる第２の履歴データとに基づいて、特徴量集合を前記計算手法別に生成し、
前記計算手法別に生成された特徴量集合を各次元のラベルとして有し、かつ、各行列要素に対応するセルにおける前記第２の履歴データに基づく前記事象の発生件数を行列要素として有する第１行列を生成し、
前記第１行列と同じラベルを有し、かつ、各行列要素に対応する各次元のラベル値の合計値を行列要素として有する第２行列を生成し、
所定の閾値に基づいて、前記第２行列において行列要素の境界を設定し、
前記境界の位置と対応する位置に存在する前記第１行列の行列要素である境界行列要素のそれぞれについて、前記境界行列要素を基準として前記第１行列を複数の領域に分割し、前記複数の領域の少なくとも１つの領域の行列要素の和を用いた評価指標を生成し、
前記境界行列要素それぞれについて生成された評価指標に基づいて、評価が最大となる前記境界行列要素に対応する各次元のラベル値の組み合わせを、前記事象の予測に用いる閾値組み合わせとして特定し、
前記閾値組み合わせに基づいて、所定の予測時点における前記対象エリアの各セルの事象の発生し易さの予測結果を分類し、その分類結果を含む情報を出力する、
ことを含む。

【0009】

本発明のプログラムは、コンピュータに、上述の情報処理方法を実行させる。

【発明の効果】

【0010】

本発明によれば、ある時点での事象の発生をコンピュータで予測するとともに、当該コンピュータによる予測の根拠を説明可能とする技術が提供される。

【図面の簡単な説明】

【0011】

【図1】第１実施形態における情報処理装置の機能構成を例示する図である。

【図2】情報処理装置のハードウエア構成を例示するブロック図である。

【図3】第１実施形態における情報処理装置の処理の流れを例示するフローチャートである。

【図4】第１実施形態における情報処理装置の処理の流れを例示するフローチャートである。

【図5】事象の発生履歴データの一例を示す図である。

【図6】第１の履歴データおよび第２の履歴データの一例を示す図である。

【図7】第１行列生成部により生成される第１行列の一例を示す図である。

【図8】第２行列生成部により生成される第２行列の一例を示す図である。

【図9】第２行列において設定される境界の一例を示す図である。

【図10】評価指標生成部により生成された境界を基に決定される、第１行列における境界行列要素を例示する図である。

【図11】評価指標生成部による第１行列の分割結果の一例を示す図である。

【図12】評価指標生成部による第１行列の分割結果の他の一例を示す図である。

【図13】評価指標生成部が生成する評価情報の一例を示す情報である。

【図14】予測結果出力部が出力する情報の一例を示す図である。

【図15】予測結果出力部が出力する情報の他の一例を示す図である。

【図16】予測結果出力部が出力する情報の他の一例を示す図である。

【図17】予測結果出力部が出力する情報の他の一例を示す図である。

【図18】第２実施形態における情報処理装置の機能構成例を示す図である。

【図19】評価指標可視化部が表示する評価指標の可視化情報の一例を示す図である。

【図20】第３実施形態における情報処理装置の機能構成例を示す図である。

【図21】予測ルール出力部により出力される情報の一例を示す図である。

【発明を実施するための形態】

【0012】

以下、本発明の実施形態について、図面を用いて説明する。尚、すべての図面において、同様な構成要素には同様の符号を付し、適宜説明を省略する。また、特に説明する場合を除き、各ブロック図において、各ブロックは、ハードウエア単位の構成ではなく、機能単位の構成を表している。また、図中の矢印の向きは、情報の流れを分かり易くするためのものであり、特に説明のない限り通信の方向（一方向通信／双方向通信）を限定しない。

【0013】

［第１実施形態］
図１は、第１実施形態における情報処理装置１０の機能構成を例示する図である。図１に示されるように、本実施形態の情報処理装置１０は、手法別特徴量生成部１１０、第１行列生成部１２０、第２行列生成部１３０、境界設定部１４０、評価指標生成部１５０、および閾値組み合わせ特定部１６０、予測結果出力部１７０を備える。

【0014】

手法別特徴量生成部１１０は、過去に発生したある事象の履歴データのうち、任意の時点を基準とする第１の期間に含まれる第１の履歴データを用いて、ある対象エリアの各セルにおける事象の発生し易さを計算する。ここで、手法別特徴量生成部１１０は、当該事象の発生し易さを計算するための複数の計算手法を使って、対象エリアの各セルにおける事象の発生し易さを計算手法別に生成する。また、手法別特徴量生成部１１０は、計算手法別に計算された各セルにおける事象の発生し易さと、上述の第１の期間よりも後の第２の期間に含まれる第２履歴データとに基づいて、特徴量集合を計算手法別に生成する。手法別特徴量生成部１１０の具体的な動作例については、後述する。

【0015】

ここで、本明細書における「事象」は、事件、事故または疾病といった避けるべきリスクに関する事象のみならず、喜ばしい事象を含み得る。なお、以下においては、リスクに関する事象（主に事件）を扱うケースを例示する。

【0016】

第１行列生成部１２０は、以下のような処理を実行して、第１行列を生成する。まず、第１行列生成部１２０は、計算手法別に生成された特徴量集合に基づいて、第１行列の各次元のラベルを設定する。例えば、手法別特徴量生成部１１０により２つの計算手法（第１の計算手法、第２の計算手法）のそれぞれについて特徴量集合が生成された場合、第１行列生成部１２０は、第１の計算手法について生成された第１の特徴量集合を第１の次元（例えば、縦軸）のラベルとして設定し、第２の計算手法について生成された第２の特徴量集合を第２の次元（例えば、横軸）のラベルとして設定する。そして、第１行列生成部１２０は、各行列要素に対応するセルにおける事象の発生件数を第２の履歴データを基に計算し、第１行列の各行列要素の値として入力することによって、第１行列を生成する。第１行列生成部１２０の具体的な動作例については、後述する。

【0017】

第２行列生成部１３０は、以下のような処理を実行して、第２行列を生成する。まず、第２行列生成部１３０は、第１行列と同じラベルを、第２行列のラベルとして設定する。つまり、第１行列および第２行列について、行要素の数および列要素の数は同一となる。第２行列生成部１３０は、第２行列の各行列要素の値として、その行列要素に対応する各次元のラベル値の合計値を入力することによって、第２行列を生成する。第２行列生成部１３０の具体的な動作例については後述する。

【0018】

境界設定部１４０は、所定の閾値に基づいて、第２行列において行列要素の境界を設定する。境界設定部１４０の具体的な動作例については、後述する。

【0019】

評価指標生成部１５０は、境界設定部１４０により設定された境界に対応する位置に存在する第１行列の行列要素（以下、「境界行列要素」とも表記）のそれぞれについて、当該境界行列要素を基準として第１行列を複数の領域に分割する。そして、評価指標生成部１５０は、個々の境界行列要素を基準としたそれぞれのケースにおいて、分割した複数の領域の少なくとも１つの領域に含まれる行列要素の和を用いて評価指標を生成する。評価指標生成部１５０の具体的な動作例については、後述する。

【0020】

閾値組み合わせ特定部１６０は、評価指標生成部１５０がそれぞれのケースについて生成した評価指標に基づいて、評価が最大となる境界行列要素を特定する。そして、評価指標生成部１５０は、特定した境界行列要素に対応する各次元のラベル値の組み合わせを、事象の予測処理に用いる閾値組み合わせとして特定する。閾値組み合わせ特定部１６０の具体的な動作例については、後述する。

【0021】

予測結果出力部１７０は、所定の予測時点における対象エリアの各セルの事象の発生し易さを、複数の計算手法を用いて予測する。そして、予測結果出力部１７０は、対象エリアの各セルの事象の発生し易さに関する予測結果を、閾値組み合わせ特定部１６０により特定された閾値組み合わせに基づいて分類する。ここで、対象エリアの各セルについて、複数の計算手法それぞれの計算結果が対応付けられる。そして、予測結果出力部１７０は、複数の計算手法それぞれの計算結果を、閾値組み合わせと照らし合わせることによって、各セルの予測結果を分類することができる。そして、予測結果出力部１７０は、対象エリアの各セルの分類結果を含む情報を、例えば表示装置などに出力する。予測結果出力部１７０の具体的な動作例については、後述する。

【0022】

＜ハードウエア構成例＞
情報処理装置１０の各機能構成部は、各機能構成部を実現するハードウエア（例：ハードワイヤードされた電子回路など）で実現されてもよいし、ハードウエアとソフトウエアとの組み合わせ（例：電子回路とそれを制御するプログラムの組み合わせなど）で実現されてもよい。以下、情報処理装置１０が、ハードウエアとソフトウエアとの組み合わせで各機能構成部を実現する場合について、図２を用いてさらに説明する。図２は、情報処理装置１０のハードウエア構成を例示するブロック図である。

【0023】

情報処理装置１０は、バス１０１０、プロセッサ１０２０、メモリ１０３０、ストレージデバイス１０４０、入出力インタフェース１０５０、及びネットワークインタフェース１０６０を有する。

【0024】

バス１０１０は、プロセッサ１０２０、メモリ１０３０、ストレージデバイス１０４０、入出力インタフェース１０５０、及びネットワークインタフェース１０６０が、相互にデータを送受信するためのデータ伝送路である。ただし、プロセッサ１０２０などを互いに接続する方法は、バス接続に限定されない。

【0025】

プロセッサ１０２０は、ＣＰＵ（Central Processing Unit） やＧＰＵ（Graphics Processing Unit）などで実現されるプロセッサである。

【0026】

メモリ１０３０は、ＲＡＭ（Random Access Memory）などで実現される主記憶装置である。

【0027】

ストレージデバイス１０４０は、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）、メモリカード、又はＲＯＭ（Read Only Memory）などで実現される補助記憶装置である。ストレージデバイス１０４０は情報処理装置１０の各機能（手法別特徴量生成部１１０、１２０、第２行列生成部１３０、境界設定部１４０、評価指標生成部１５０、閾値組み合わせ特定部１６０、予測結果出力部１７０など）を実現するプログラムモジュールを記憶している。プロセッサ１０２０がこれら各プログラムモジュールをメモリ１０３０上に読み込んで実行することで、各プログラムモジュールに対応する各機能が実現される。

【0028】

入出力インタフェース１０５０は、情報処理装置１０と各種入出力デバイスとを接続するためのインタフェースである。入出力インタフェース１０５０には、キーボードやマウスといった入力装置（図示せず）、ディスプレイやスピーカーといった出力装置（図示せず）などが接続され得る。

【0029】

ネットワークインタフェース１０６０は、情報処理装置１０をネットワークに接続するためのインタフェースである。このネットワークは、例えばＬＡＮ（Local Area Network）やＷＡＮ（Wide Area Network）である。ネットワークインタフェース１０６０がネットワークに接続する方法は、無線接続であってもよいし、有線接続であってもよい。

【0030】

なお、情報処理装置１０のハードウエア構成は図２に例示される構成に限定されない。

【0031】

＜処理の流れ＞
図３および図４は、第１実施形態における情報処理装置１０の処理の流れを例示するフローチャートである。図３では、事象の発生予測を行うにあたり前もって実行される準備フェーズの処理が示される。また、図４は準備フェーズの処理結果に基づく予測フェーズの処理が示される。

【0032】

＜＜準備フェーズの処理＞＞
まず、手法別特徴量生成部１１０は、分析（予測）の対象エリアに関する情報、および、分析（予測）に利用する事象の履歴データを取得する（Ｓ１０２）。対象エリアに関する情報は、例えば、オペレータの入力操作（住所の入力・指定操作や地図上でのエリア指定操作など）によって生成される。また、履歴データは、例えば、図５に示すような形式で、所定の記憶領域（例えば、ストレージデバイス１０４０など）に記憶されている。図５は、事象の発生履歴データの一例を示す図である。図５に示すように、発生履歴データは、少なくとも、事象の発生時刻を示す情報（発生年月日時刻）と当該事象の発生位置を示す情報（発生位置）とを含む。発生履歴データは、事象の種類や当該事象に関連する人物の属性といった情報を含んでいてもよい。そして、手法別特徴量生成部１１０は、分析（予測）すべき事象の種類（例えば、犯罪の種類など）ごとの履歴データを取得することができる。例えば、分析（予測）すべき事象の種類をオペレータが入力装置を介して入力し、手法別特徴量生成部１１０は、その入力に従って読み込むべき履歴データを選択するように構成されていてもよい。

【0033】

そして、手法別特徴量生成部１１０は、分析の対象となるエリア（対象エリア）を複数の小領域（以下、「セル」と表記）に分割する（Ｓ１０４）。手法別特徴量生成部１１０は、複数のセルの大きさや形状を任意に決定することができる。例えば、手法別特徴量生成部１１０は、予め記憶されているパラメータに従って各セルの大きさや形状を決定してもよいし、オペレータが入力装置を用いて入力内容に従って各セルの大きさや形状を決定してもよい。また、複数のセルそれぞれの大きさや形状は、均一であっても異なっていてもよい。一例として、同一サイズの矩形のセルによって対象エリアが分割されてもよい。他の一例として、住所の区分（都道府県、区市町村、町丁目など）といった不均一なサイズおよび形状のセルによって対象エリアが分割されてもよい。他の一例として、住所と市街化地域の対応データがあれば、当該対応データを参照して、市街化地域とそれ以外の地域に対象エリアが分割されてもよい。また、例えば対象エリアが複数の都道府県を含む場合に、市街化地域が占める面積割合の指定入力を受け付け、市街化地域の面積割合が指定された割合以上の都道府県と、それ以外の都道府県とに対象エリアが分割されてもよい。なお、対象エリアが既に複数のセルに分割されている場合（対象エリアに関するセルの情報が既に与えられている場合）、手法別特徴量生成部１１０は、Ｓ１０４の処理を実行しなくてもよい。

【0034】

対象エリアが複数のセルに分割された後、手法別特徴量生成部１１０は、事象の発生し易さを計算する複数の計算手法それぞれについて、ある時点を基準とする第１の期間に含まれる事象の履歴データ（第１の履歴データ）を基に各セルの事象の発生し易さを計算する（Ｓ１０６）。なお、手法別特徴量生成部１１０は、第１の期間に含まれる履歴データの全てを第１の履歴データとして用いて各セルの事象の発生し易さを計算してもよいし、第１の期間に含まれる履歴データの一部を第１の履歴データとして用いて各セルの事象の発生し易さを計算してもよい。ここで、手法別特徴量生成部１１０は、既知の計算手法の中から任意の複数の計算手法を選択して利用することができる。既知の計算手法の例としては、特に限定されないが、例えば、上述の特許文献１に開示されるカーネル密度推定による分布計算手法および時空間Ｋ関数による分布計算手法や、上述の非特許文献１に開示されるグリッドマッピングによる分布計算手法、上述の非特許文献２乃至４に開示されるリスク値に関係する環境要因を用いた計算手法などがある。その他にも、ベテラン担当者の経験に基づく各セルのリスクの発生し易さや対象エリアの各セルのパトロールのし易さといった、人的知見による手法も既知の計算手法に含まれる。

【0035】

また、手法別特徴量生成部１１０は、事象の履歴データと他のデータとを組み合わせて、各セルの事象の発生し易さを計算する手法を選択してもよい。例えば、手法別特徴量生成部１１０は、リスク値に関係する様々な環境要因を用いて、各セルの事象の発生し易さを計算する手法を選択するように構成され得る。第１の具体例として、手法別特徴量生成部１１０は、第１の履歴データ（計算に使用する履歴データ）が示す事象の発生時刻および発生位置を用いて特定できる天候データに基づいて、各セルの事象の発生し易さを計算する計算手法を選択するように構成され得る。第２の具体例として、手法別特徴量生成部１１０は、第１の履歴データ（計算に使用する履歴データ）が示す事象の発生位置を用いて特定できるエリアの人口密度に基づいて、各セルの事象の発生し易さを計算する計算手法を選択するように構成され得る。なお、この場合において、各計算手法を定義するデータは、例えば、その計算手法を選択する基準となる情報（環境要因に関する情報）と対応付けた状態で、手法別特徴量生成部１１０がアクセス可能な所定の記憶領域に予め記憶される。例えば、第１の具体例のケースでは、各計算手法のデータは、選択時の基準となる天候に関する情報（例：晴れ、くもり、雨、雪などを示す情報）と対応付けて記憶される。また例えば、第２の具体例のケースでは、各計算手法のデータは、選択時の基準となる人口密度に関する情報（例：人口密度の大きさを段階的な数値（ランク）で示す情報や、人口密度の数値を示す情報）と対応付けて記憶される。

【0036】

ここで、第１の履歴データと第２の履歴データについて説明する。図６は、第１の履歴データおよび第２の履歴データの一例を示す図である。図６において、バツ印は、後述の｛Ｉ^ｋ｝（ｋ＝１，２，・・・，Ｋ）に含まれる履歴データを例示している。また、図６において、三角の印は、後述の｛Ｅ^ｋ｝（ｋ＝１，２，・・・，Ｋ）に含まれる履歴データを例示している。

【0037】

図６に示されるように、手法別特徴量生成部１１０は、少なくとも１つの評価時点｛ｔ^ｋ｝（ｋ＝１，２，・・・，Ｋ）を設定する。またここで、各評価時点ｔ^ｋを基準として決まる過去評価期間ｈ_ｔは、それぞれ、上述の「第１の期間」に該当する。また、各評価時点ｔ^ｋを基準として決まる各評価期間Δｔは、それぞれ、上述の「第２の期間」に該当する。そして、手法別特徴量生成部１１０は、各評価時点ｔ^ｋを基準として決まる各期間ｈ_ｔに含まれる履歴データ｛Ｉ^ｋ｝（ｋ＝１，２，・・・，Ｋ）を用いて、各セルにおける事象の発生し易さを計算手法別に計算する。

【0038】

そして、手法別特徴量生成部１１０は、複数の計算手法それぞれについて計算された各セルにおける事象の発生し易さと、各評価時点ｔ^ｋを基準として決まる各評価期間Δｔに含まれる履歴データ｛Ｅ^ｋ｝（ｋ＝１，２，・・・，Ｋ）とを用いて、計算手法別の特徴量集合を生成する（Ｓ１０８）。具体的には、手法別特徴量生成部１１０は、以下のように動作する。

【0039】

手法別特徴量生成部１１０は、各セルにおける事象の発生し易さに基づいて、計算手法別に各セルを順位付けする。そして、手法別特徴量生成部１１０は、事象の発生し易さに基づいて決まる各セルの順位を、対象エリアの全セル数で割ることによって得られる評価値を算出する。そして、手法別特徴量生成部１１０は、第２の期間に含まれる履歴データ｛Ｅ^ｋ｝（ｋ＝１，２，・・・，Ｋ）に基づいて、対象エリアにおいて事象が発生したセルを識別する。ここで、各セルには、計算手法別に計算された評価値が対応付けられている。手法別特徴量生成部１１０は、計算手法別に、識別した各セルの評価値を比較する処理を行い、当該評価値を大きさ順に並べる。これにより、計算手法別の特徴量集合が生成される。

【0040】

なお、手法別特徴量生成部１１０が算出する評価値は、複数の計算手法に対して共通に定義でき、かつ、事象の発生し易さ、つまり発生確率に準じた意味を持たせることができる何らかの尺度であれば、特に制限されない。例えば、複数の計算手法の計算結果のすべてが確率（事象の発生確率）として解釈可能であれば、手法別特徴量生成部１１０は、それら複数の計算手法の計算結果に基づく事象の発生確率を評価値として用いてもよい。また、複数の計算手法の計算結果として、特定のスポット（事象に関連する施設など）からの距離が利用される場合、手法別特徴量生成部１１０は、それら複数の計算手法から得られる距離を評価値として用いてもよい。

【0041】

複数の計算手法それぞれの特徴量集合が生成された後、第１行列生成部１２０は、第１行列を生成する（Ｓ１１０）。第１行列は、例えば、次のように生成される。

【0042】

まず、第１行列生成部１２０は、計算手法別の特徴量集合を、行列の各次元のラベルに設定する。例えば、使用された計算手法の数が２つである場合、第１行列生成部１２０は、第１の計算手法について生成された特徴量集合を第１の次元（例えば、縦軸）のラベルＲ１_ｉ（ｉ＝１，２，・・・，Ｉ）_、第２の計算手法について生成された特徴量集合を第２の次元（例えば、横軸）のラベルＲ２_ｊ（ｊ＝１，２，・・・，Ｊ）として設定する。

【0043】

そして、第１行列生成部１２０は、第２の期間に含まれる履歴データを用いて、第１行列の各行列要素に値を入力していく。なお、第１行列の各行列要素の初期値はいずれも０である。第１行列生成部１２０は、第２の期間に含まれる履歴データ｛Ｅ^ｋ｝（ｋ＝１，２，・・・，Ｋ）に基づいて、第２の期間において事象が発生したセルを識別する。そして、第１行列生成部１２０は、識別したセルに対応付けられた計算手法別の評価値を用いて、値を入力すべき行列要素を特定する。例えば、２つの計算手法（第１の計算手法、第２の計算手法）を使って、２次元の第１行列が生成されるケースを考える。第１行列生成部１２０は、識別したセルそれぞれについて、第１の計算手法に基づく評価値および第２の計算手法に基づく評価値の組み合わせを取得する。そして、第１行列生成部１２０は、第１行列の各次元に設定されたラベルの値を参照し、当該評価値の組み合わせと同じ組み合わせとなる行列要素を特定する。そして、第１行列生成部１２０は、特定した行列要素に１を加算する。第１行列生成部１２０は、履歴データ｛Ｅ^ｋ｝（ｋ＝１，２，・・・，Ｋ）の全てについて、上述の行列要素を特定して値を加算する処理を繰り返すことによって、第１行列を生成することができる。

【0044】

図７は、第１行列生成部１２０により生成される第１行列の一例を示す図である。図７には、２次元の第１行列の一例が示されている。なお図７では、図の見易さを考慮して、初期値（＝０）のままの行列要素は空欄で示されている。図７の例において、第１の計算手法の特徴量集合および第２の計算手法の特徴量集合が、それぞれ、縦軸方向のラベルＲ１および横方向のラベルＲ２として設定されている。このような第１行列において、ある行列要素（Ｒ１_ｉ，Ｒ２_ｊ）の値が「２」であったとする。これは、第１の計算手法での評価値が「Ｒ１_ｉ」であり、かつ、第２の計算手法での評価値が「Ｒ２_ｊ」であるセルについては、累計２件の事象が発生したことを意味している。なお、各計算手法において、複数の異なるセルに同じ評価値が付与される場合もある。そのため、第１行列の各行列要素に格納される値は、ある特定の１つのセルで発生した事象の件数を示すとは限らない。また、第１行列の次元数は、２次元に限られない。第１行列は、計算手法の数に応じた次元を有する。例えば、３種類の計算手法を用いた場合、第１行列は３次元行列となる。

【0045】

次に、第２行列生成部１３０は、第２行列を生成する（Ｓ１１２）。第２行列は、例えば、次のように生成される。

【0046】

まず、第２行列生成部１３０は、第１行列と同じラベルを、第２行列のラベルとして設定する。つまり、第１行列および第２行列について、行要素の数および列要素の数は同じとなる。第２行列生成部１３０は、各行列要素に、その行列要素に対応する各次元のラベルの値の合計値を入力する。図８は、第２行列生成部１３０により生成される第２行列の一例を示す図である。図８で示されるように、第２行列の各行列要素の値は、対応する各次元のラベル値によって決定される。具体的には、第２行列のある行列要素（Ｒ１_ｉ，Ｒ２_ｊ）の値は、当該行列要素の行ラベルの値Ｒ１_ｉと列ラベルの値Ｒ２_ｊとの合計値となる。

【0047】

次に、境界設定部１４０は、メモリ１０３０などの記憶領域に予め用意された或いはオペレータの入力により与えられた所定の閾値を用いて、第２行列において境界を設定する（Ｓ１１４）。本例において、この閾値は、対象エリアの全セルに対する、"事象が発生する可能性の高いセル"として予測するセル数の占める割合を示す。例えば、所定の閾値として「０．０５」という値が与えられた場合、"事象が発生する可能性の高いセル"として予測するセルの数が全セル数の５％に制約されることを意味する。この制約は、例えば、事象に関連する業務（例：パトロールなど）に対応可能なリソース（人的リソース、機材リソース）の量に基づいて決定される。なお、所定の閾値は、時間的リソースおよび費用などに基づいて決定されてもよい。この場合、境界設定部１４０は、例えば、時間的リソースまたは費用に関する数値をパラメータとして所定の関数に代入することによって、所定の閾値を導出することができる。この場合の所定の関数は、例えば、メモリ１０３０などの記憶領域に予め用意される。

【0048】

そして、境界設定部１４０は、第２行列の各行列要素の値と所定の閾値とを比較し、例えば図９に示されるように、各行列要素の値が所定の閾値を超える位置に境界を設定する。

【0049】

図９は、第２行列において設定される境界の一例を示す図である。図９において、網掛けで示す行列要素の位置が境界の位置を示している。境界設定部１４０は、所定の閾値β（本図では「０．０５」）と、行列要素の行ラベルの値Ｒ１_ｉおよび列ラベルの値Ｒ２_ｊとの関係において、「Ｒ１_ｉ＋Ｒ２_ｊ≦β」が成り立つ行列要素と、「Ｒ１_ｉ＋Ｒ２_ｊ＞β」が成り立つ行列要素とが隣接する位置を特定し、前者が成り立つ行列要素のうち最大の行列要素の位置を、境界の位置として設定する。これら境界に位置する行列要素に対応するラベル値の組み合わせの集合を｛ＴＨ^ｌ｝（ｌ＝１，２，・・・，Ｌ）と表記する。ｌは、例えば、各行および各列についての探索順に付与される。また、｛ＴＨ^ｌ｝の１つの要素ＴＨ^ｌの値を、行ラベル値および列ラベル値をそれぞれｔｈ１^ｌおよびｔｈ２^ｌとして、ＴＨ^ｌ（ｔｈ１^ｌ，ｔｈ２^ｌ）と表す。するとすべてのｌについて，「ｔｈ１^ｌ＋ｔｈ２^ｌ≦β」が成り立つ。

【0050】

そして、評価指標生成部１５０は、境界設定部１４０により設定された境界の位置に対応する位置に存在する第１行列の行列要素（境界行列要素）を基に、複数の計算手法それぞれの評価値の組み合わせに関する評価指標を生成する（Ｓ１１６）。以下では、計算手法が２つである場合の評価指標生成部１５０の動作例について説明する。

【0051】

評価指標生成部１５０は、次の処理を実行する。
（１）ラベル値組み合わせ集合｛ＴＨ^ｌ｝の中から１つの要素ＴＨ^ｌを取り出す。
（２）取り出した要素ＴＨ^ｌのラベル値（ｔｈ１^ｌ，ｔｈ２^ｌ）を用いて、第１行列を複数の領域に分割する。
（３）複数の領域の少なくとも１つについて、行列要素の和を算出する。
（４）（３）で算出した行列要素の和を基に、取り出した要素ＴＨ^ｌの評価指標を生成する。
（５）｛ＴＨ^ｌ｝の全ての要素を取り出すまで、（１）〜（４）の処理を繰り返す。

【0052】

上述の動作を、図１０乃至図１２を用いて説明する。図１０は、評価指標生成部１５０により生成された境界を基に決定される、第１行列における境界行列要素を例示する図である。図１０に示されるように、評価指標生成部１５０は、境界設定部１４０により第２行列において設定された境界に基づいて、当該境界の位置と対応する位置に存在する第１行列の行列要素を境界行列要素として特定する。評価指標生成部１５０は、境界行列要素の１つを取り出し、当該取り出した境界行列要素を基準として、第１行列を複数の領域に分割する。ここで、第１行列の分割結果は、取り出された境界行列要素に応じて変わる。例えば、評価指標生成部１５０が複数の境界行列要素の中から境界行列要素ＴＨ^７を取り出した場合、第１行列は図１１に示されるように分割される。図１１は、評価指標生成部１５０による第１行列の分割結果の一例を示す図である。また、評価指標生成部１５０が複数の境界行列要素の中から境界行列要素ＴＨ^１８を取り出した場合、第１行列は図１２に示されるように分割される。図１２は、評価指標生成部１５０による第１行列の分割結果の他の一例を示す図である。

【0053】

図１１および図１２に示されるように、２次元の第１行列は、取り出した境界行列要素ＴＨ^ｌ（ｔｈ１^ｌ，ｔｈ２^ｌ）により、４つの領域Ａ〜Ｄに分割される。ここで、領域Ａに含まれる行列要素の和は、第１の計算手法での評価値がｔｈ１^ｌ以下であり、かつ、第２の計算手法での評価値がｔｈ２^ｌ以下であるという条件を満たすセルで発生した事象の件数を意味する。言い換えると、領域Ａには、先に述べたリソースに関する制約下での予測において、第１の計算手法および第２の計算手法の双方で"事象が発生する可能性の高いセル"として予測されたセルで発生した事象が含まれる。同様に、領域Ｂ、Ｃ、Ｄには、それぞれ、第１の計算手法だけで"事象が発生する可能性の高いセル"として予測されたセルで発生した事象、第２の計算手法だけで"事象が発生する可能性の高いセル"として予測されたセルで発生した事象、第１の計算手法および第２の計算手法のいずれにおいても"事象が発生する可能性の高いセル"として予測されなかったセルで発生した事象が含まれる。なお，予測されたセルが含まれる領域Ａ〜Ｄを，予測されたセルの属性とよぶ。

【0054】

そして、評価指標生成部１５０は、図１１や図１２に示されるように分割された複数の領域（領域Ａ〜Ｄ）の少なくともいずれか１つの領域に含まれる行列要素の和を算出する。複数の領域のうちどの領域について行列要素の和を算出するかは、あらかじめ決められた評価指標の算出の仕方によって決まる。なお、評価指標の生成方法は様々である。例えば、評価指標生成部１５０は、領域Ｄに含まれる行列要素の和に関する以下の式（１）に基づいて、評価指標を生成することができる。

【0055】

【数1】

【0056】

また例えば、評価指標生成部１５０は、上述の履歴データ｛Ｅ^ｋ｝の総数に対して上述の領域Ａに含まれる行列要素の和が占める割合に基づく評価指標を生成してもよい。また例えば、評価指標生成部１５０は、上述の履歴データ｛Ｅ^ｋ｝の総数に対して上述の領域Ｂ、あるいは、領域Ｃに含まれる行列要素の和が占める割合に基づく評価指標を生成してもよい。また、評価指標生成部１５０は、上述の履歴データ｛Ｅ^ｋ｝の総数に対して上述の領域Ａ〜Ｄの少なくとも２つ以上の領域に含まれる行列要素の和が占める割合に基づく評価指標を生成してもよい。また例えば、評価指標生成部１５０は、事象が発生する可能性の高いセル"として予測されたセル数に対する、評価期間Δｔにおいて実際に事象が発生したセル数の割合を評価指標として生成してもよい。また例えば、評価指標生成部１５０は、評価期間Δｔにおいて事象が発生したセルの総数に対する、"事象が発生する可能性の高いセル"として予測されたセルの中で当該評価期間Δｔにおいて実際に事象が発生したセル数の割合を評価指標として生成してもよい。

【0057】

上述の式（１）に基づいて、評価指標生成部１５０が生成する情報の一例を図１３に示す。図１３は、評価指標生成部１５０が生成する評価情報の一例を示す情報である。図１３では、境界行列要素｛ＴＨ^ｌ｝（ｌ＝１，２，・・・，３２）における行ラベル値ｔｈ１および列ラベル値ｔｈ２の組み合わせと、領域Ｄに含まれる行列要素の和と、上記式（１）に基づく評価指標とを対応付けた情報が例示されている。

【0058】

閾値組み合わせ特定部１６０は、｛ＴＨ^ｌ｝の全ての要素それぞれの評価指標を比較し、評価が最大となる要素ＴＨ^ｌの値の組み合わせ（ｔｈ１^ｌ，ｔｈ２^ｌ）を、閾値組み合わせとして特定する（Ｓ１１８）。一例として、評価指標生成部１５０が、上述の式（１）を用いて評価指標として生成した場合を考える。この場合、評価指標の値が大きいほど、複数の計算手法のいずれでも予測できなかった事象の数が少なくなる。図１３の例では、ＴＨ^３１（０．０００７５、０．０４７３６）に対応する評価指標値「０．３４８」が最大となる。この「０．３４８」という数値は、すべての評価時点｛ｔ^ｋ｝（ｋ＝１，２，・・・，Ｋ）において「第１の計算手法による評価値Ｒ１≦０．０００７５」および「第２の計算手法による評価値Ｒ２≦０．０４７３６」の少なくともいずれか一方の条件が成り立つセルを"事象が発生する可能性の高いセル"として予測していた場合に、対応する評価期間Δｔ内に発生した全事象｛Ｅ^ｋ｝（ｋ＝１，２，・・・，Ｋ）のうちの３４．８％について予測できていたことを意味する。つまり、この場合、閾値組み合わせ特定部１６０は、ＴＨ^３１（０．０００７５、０．０４７３６）の組み合わせを、評価を最大とする組み合わせとして特定する。他の一例として、評価指標生成部１５０が、「履歴データ｛Ｅ^ｋ｝の総数に対して領域Ａに含まれる行列要素の和が占める割合」を評価指標として生成した場合を考える。この場合、評価指標の値が大きいほど、複数の計算手法の双方で予測された事象の数が多くなる。つまり、この場合、閾値組み合わせ特定部１６０は、領域Ａの行列要素の和が最大となる（ｔｈ１^ｌ，ｔｈ２^ｌ）の組み合わせを、評価を最大とする組み合わせとして特定する。

【0059】

そして、閾値組み合わせ特定部１６０は、特定した要素ＴＨ^ｌの値の組み合わせ（ｔｈ１^ｌ，ｔｈ２^ｌ）を、閾値組み合わせとして、メモリ１０３０やストレージデバイス１０４０等の記憶領域に記憶する（Ｓ１２０）。このとき、閾値組み合わせ特定部１６０は、複数の計算手法の組み合わせを示す情報および制約条件（所定の閾値）を示す情報と共に、特定した閾値組み合わせをメモリ等の記憶領域に記憶する。ここで記憶領域に記憶された閾値組み合わせは、後述の予測フェーズの処理で利用される。

【0060】

＜＜予測フェーズの処理＞＞
予測結果出力部１７０は、複数の計算手法それぞれについて、事象の履歴データを基にある予測時点での事象の発生し易さを計算し、複数の計算手法それぞれにおける計算結果を基に"事象が発生する可能性の高いセル"を予測する（Ｓ１２２）。例えば、予測結果出力部１７０は、予測時点（予測年月日時刻ｔ^ｐ）を取得し、その予測年月日時刻ｔ^ｐについて、複数の計算手法それぞれについて、事象の履歴データを基に、事象の発生し易さを計算する。そして、予測結果出力部１７０は、閾値組み合わせ特定部１６０により特定された閾値組み合わせを使って、"事象が発生する可能性が高いセル"を対象エリアのセルの中から選択することにより予測する。また、予測結果出力部１７０は、対象エリアにおいて"事象が発生する可能性の高いセル"と予測されたセルを、閾値組み合わせ特定部１６０により特定された閾値組み合わせに基づいて分類する（Ｓ１２４）。例えば、予測結果出力部１７０は、"事象が発生する可能性の高い"と予測されたあるセルについて、第１の計算手法によるそのセルの評価値Ｒ１および第２の計算手法によるそのセルの評価値Ｒ２の組み合わせと、メモリ等に記憶された閾値組み合わせ（ｔｈ１^ｌ，ｔｈ２^ｌ）とを比較する。そして、予測結果出力部１７０は、評価値Ｒ１と閾値ｔｈ１^ｌとの比較結果、および、評価値Ｒ２と閾値ｔｈ２^ｌとの比較結果に基づいて、予測されたセルの属性（当該セルが領域Ａ〜Ｄのどの領域に含まれるか）を特定する。そして、予測結果出力部１７０は、Ｓ１２２の処理による予測結果とＳ１２４の処理による分類結果とを含む情報（例：図１４）を出力する（Ｓ１２６）。

【0061】

図１４は、予測結果出力部１７０が出力する情報の一例を示す図である。図１４に示される例では、予測結果出力部１７０は、"事象が発生する可能性の高いセル"として予測されたセルと、当該セルが閾値組み合わせ特定部１６０により特定された組み合わせによって分類される領域のいずれに含まれるかを示す属性情報とを含む一覧を出力している。図１４のような情報によれば、"事象が発生する可能性の高いセル"として予測されたセルが、どの計算手法を根拠として選択されたかが容易に判断できる。例えば、セル（１２，４）およびセル（２３，９）は、第１の計算手法および第２の計算手法の双方において"事象が発生する可能性の高いセル"として予測されたことがわかる。なお、セル（１２，４）は、手法別特徴量生成部１１０において対象エリアを複数セルに分割（Ｓ１０４）した際に付与したセル識別ＩＤが（１２，４）であるセルを示す。また、セル（９，３）、セル（１４，８）、セル（２７，４３）、セル（３５，７）およびセル（４２，３５）は、第１の計算手法だけで"事象が発生する可能性の高いセル"として予測されたことがわかる。また、セル（１１，４０）、セル（１４，８）およびセル（２７，３５）は、第２の計算手法だけで"事象が発生する可能性の高いセル"として予測されたことがわかる。

【0062】

なお、予測結果出力部１７０が出力する情報は、図１４の例に制限されない。例えば、予測結果出力部１７０は、図１５や図１６に示されるような情報を出力するように構成されていてもよい。図１５及び図１６は、予測結果出力部１７０が出力する情報の他の一例を示す図である。図１５の例では、予測結果出力部１７０は、「閾値組み合わせに基づく分類」の代わりに、予測に用いた閾値組み合わせと予測されたセルの計算手法別の評価値との大小関係を示す情報を出力している。また、図１６の例では、予測結果出力部１７０は、予測に用いた閾値組み合わせに基づく条件の成立有無を示す情報を出力している。これらの図に示すような情報によっても、"事象が発生する可能性の高いセル"として予測されたセルが、どの計算手法を根拠として選択されたかが容易に判断できる。

【0063】

また、予測結果出力部１７０は、上述したような各セルの分類結果と、ＧＩＳ（Geographic Information System）などの地図情報を提供するシステムから得られるデータとを組み合わせた情報を表示装置に出力するように構成されていてもよい。この場合、予測結果出力部１７０は、例えば図１７に例示されるような情報を出力する。図１７は、予測結果出力部１７０が出力する情報の他の一例を示す図である。図１７の例において、水平線のハッチングパターンを有するセルは「第１の計算手法で計算された評価値Ｒ１_ｉ≦閾値組み合わせにおけるｔｈ１^ｌ」という条件を満たすセルである。また、垂直線のハッチングパターンを有するセルは、「第２の計算手法で計算された評価値Ｒ２_ｊ≦閾値組み合わせにおけるｔｈ２^ｌ」という条件を満たすセルである。また、斜線のハッチングパターンを有するセルは、「第１の計算手法で計算された評価値Ｒ１_ｉ≦閾値組み合わせにおけるｔｈ１^ｌ」かつ「第２の計算手法で計算された評価値Ｒ２_ｊ≦閾値組み合わせにおけるｔｈ２^ｌ」という条件を満たすセルである。図１７に例示されるような情報によっても、各セルがどの計算手法によって"事象が発生する可能性の高いセル"と予測されたかが容易に判断できる。なお、予測結果出力部１７０は、計算手法それぞれに割り当てられた固有色（例えば、第１の計算手法は赤、第２の計算手法は青、双方が重なる場合は紫など）によってセルを塗り分けるように構成されていてもよい。またこの場合、予測結果出力部１７０は、セルの色の濃さを、そのセルにおける評価値の大きさに応じて変更してもよい。このような情報によれば、「どの計算手法で予測されたか」ということと共に、その予測の根拠の度合いを容易に判断することができる。

【0064】

また、予測結果出力部１７０は、図１７に例示されるような情報を表示する場合に、セルの属性による表示の絞り込み機能を提供するように構成されていてもよい。ここで、セルの属性とは、例えば、第１の計算手法だけで予測されたセル、第２の計算手法だけで予測されたセル、第１の計算手法および第２の計算手法の双方で予測されたセル、などである。この場合、予測結果出力部１７０は、例えばオペレータによる属性の選択入力結果に応じて、選択された属性を有するセルのみにハッチング表示を行うように構成される。

【0065】

なお、予測結果出力部１７０は、事象の種類（例えば、「ひったくり」や「自転車盗」などの犯罪の種類）毎の発生し易さについて複数の計算手法を使った予測を行い、その結果として予測されたセルにつき、閾値組み合わせに基づく分類処理を行ってもよい。この場合、予測結果出力部１７０は、例えば、図１４乃至図１６に例示されるような一覧情報を、事象の種類ごとに生成して表示装置に表示する。また、予測結果出力部１７０は、例えば、図１７に例示されるような地図情報を事象の種類ごとに生成し、それらを統合した状態で表示装置に表示してもよい。また、この場合に、予測結果出力部１７０は、事象の種類を指定する入力を受け付け、当該入力によって指定された種類に関する情報を表示装置に表示させるように構成されていてもよい。

【0066】

以上、本実施形態によれば、コンピュータを用いて事象の発生予測を行った場合に、その予測がどのような判断基準でなされたかを簡単に説明できるようになる。

【0067】

＜変形例＞
上記では、情報処理装置１０が１つの対象エリアについて閾値組み合わせを特定する構成例について説明した。ここで、情報処理装置１０は、複数の対象エリアについて共通の閾値組み合わせを特定するように構成されていてもよい。ここでは、複数の対象エリアについて共通の閾値組み合わせを特定する場合の情報処理装置１０の構成および動作の例について説明する。

【0068】

本変形例のケースでは、まず、複数の対象エリアを示す情報と、各々の対象エリアにおける事象の履歴データとが情報処理装置１０に対して与えられる。このとき、情報処理装置１０は、任意の基準に基づいて、複数の対象エリアを１以上の組に分類する。一例として、情報処理装置１０は、「単位面積あたりの事象の発生件数（発生件数密度）」という基準を用いて、複数の対象エリアを１以上の組に分類することができる。この場合、情報処理装置１０は、次のように動作する。まず、情報処理装置１０は、各対象エリアについて、履歴データの件数（事象の発生件数）を面積で割った値（発生件数密度）を算出する。そして、情報処理装置１０は、各対象エリアの分類先（組）を、各対象エリアについて算出された発生件数密度に基づいて特定する。これにより、互いに発生密度件数（特徴）の類似する複数の対象エリアが、同一の組に分類される。この処理は、手法別特徴量生成部１１０によって実行されてもよいし、図示しない新たな処理部によって実行されてもよい。

【0069】

複数の対象エリアの分類が完了した後、手法別特徴量生成部１１０は、複数の対象エリアそれぞれについて、上述したように、複数の計算手法を用いて各セルの事象の発生し易さを計算する。また、手法別特徴量生成部１１０は、対象エリアの組毎に、計算手法別の特徴量集合を生成する。具体的な例として、３つの対象エリアが所定の基準に基づいて同一の組に分類された場合を考える。この場合、手法別特徴量生成部１１０は、これら３つの対象エリアを１つのエリアＡ１と見做し、上述したように、計算手法別の特徴量集合を生成する。具体的には、手法別特徴量生成部１１０は、エリアＡ１中の各セルについて、第１の計算手法を用いて事象の発生し易さを計算し、その値に基づいた各セルの評価値を計算して、エリアＡ１に関する一つの特徴量集合を生成する。手法別特徴量生成部１１０は、第１の計算手法以外に用いた他の計算手法についても、同様にして、特徴量集合を生成する。なお、上述の分類処理によって対象エリアの組が複数得られた場合、手法別特徴量生成部１１０は、ここで説明したような処理を対象エリアの組毎に実行する。このようにして、手法別特徴量生成部１１０は、対象エリアの組毎に、計算手法別の特徴量集合を生成する。

【0070】

次に、第１行列生成部１２０は、対象エリアの組毎に、第１行列を生成する処理を実行する。例えば上述の分類処理によって対象エリアの組が複数得られた場合、第１行列生成部１２０は、対象エリアの複数の組の１つを選択する。そして、第１行列生成部１２０は、その選択された組について生成された計算手法別の特徴量集合を、その組の第１行列の各次元のラベルに設定する。そして、第１行列生成部１２０は、選択された組に分類されている各対象エリアの履歴データのうち、先に説明した「第２の期間」に含まれる履歴データを用いて、第１行列の行列要素に値を入力していく。これにより、選択された組の第１行列が生成される。第１行列生成部１２０は、複数の組それぞれについて同様の処理を繰り返し、複数の組それぞれの第１行列を生成する。

【0071】

次に、第２行列生成部１３０は、対象エリアの組毎に、第２行列を生成する処理を実行する。例えば上述の分類処理によって対象エリアの組が複数得られた場合、第２行列生成部１３０は、まず、対象エリアの複数の組の１つを選択する。そして、第２行列生成部１３０は、その選択された組について生成された第１行列と同じラベルを、第２行列のラベルとして設定する。そして、第２行列生成部１３０は、選択された組の第２行列の各行列要素に、その行列要素に対応する各次元のラベルの値の合計値を入力する。これにより、選択された組の第２行列が生成される。第２行列生成部１３０は、複数の組それぞれについて同様の処理を繰り返し、複数の組それぞれの第２行列を生成する。

【0072】

次に、境界設定部１４０は、所定の閾値を用いて、複数の組それぞれの第２行列において境界を設定する。例えば上述の分類処理によって対象エリアの組が複数得られた場合、境界設定部１４０は、まず、対象エリアの複数の組の１つを選択する。そして、境界設定部１４０は、その選択された組について生成された第２行列の各行列要素の値と所定の閾値とを比較し、その組の第２行列において境界を設定する。境界設定部１４０は、複数の組それぞれについて同様の処理を繰り返し、複数の組それぞれの第２行列において境界を設定する。

【0073】

次に、評価指標生成部１５０は、対象エリアの組毎に、境界設定部１４０により設定された境界を用いて、上述の「評価指標」を算出する。例えば上述の分類処理によって対象エリアの組が複数得られた場合、評価指標生成部１５０は、まず、対象エリアの複数の組の１つを選択する。そして、評価指標生成部１５０は、その選択された組の第２行列において設定された境界の位置に対応する位置に存在する、その組の第１行列の行列要素（境界行列要素）を特定する。そして、評価指標生成部１５０は、境界行列要素毎に第１行列を複数の領域に分割する。そして、評価指標生成部１５０は、各分割領域に含まれる行列要素を基に、複数の計算手法それぞれの評価値の組み合わせに関する評価指標を生成する。評価指標生成部１５０は、複数の組それぞれについて同様の処理を繰り返し、複数の組それぞれについて評価指標を生成する。

【0074】

次に、閾値組み合わせ特定部１６０は、対象エリアの組毎に、評価指標生成部１５０により生成された評価指標を基に、上述の「閾値組み合わせ」を特定する。例えば上述の分類処理によって対象エリアの組が複数得られた場合、閾値組み合わせ特定部１６０は、まず、対象エリアの複数の組の１つを選択する。そして、閾値組み合わせ特定部１６０は、選択された組について、上述の境界行列要素毎に第１行列を分割した結果から得られた各評価指標のうち、所定の基準による評価が最大となる評価指標を特定する。そして、閾値組み合わせ特定部１６０は、特定した評価指標に対応する境界行列要素の各ラベル値の組み合わせを、閾値組み合わせとして特定する。閾値組み合わせ特定部１６０は、複数の組それぞれについて同様の処理を繰り返し、複数の組それぞれについて閾値組み合わせを生成する。複数の組それぞれについて生成された閾値組み合わせは、当該組を識別する組識別情報（例：組毎に一意な識別子など）と対応付けて、所定の記憶領域（例：メモリ１０３０やストレージデバイス１０４０など）に記憶される。

【0075】

本変形例では、対象エリアの組毎に閾値組み合わせが生成される。これにより、ある組に含まれる全ての対象エリアの予測において共通で利用可能な閾値組み合わせを特定することができる。

【0076】

なお、予測フェーズにおいて、予測結果出力部１７０は、対象エリアが属する組を示す情報を入力情報として更に取得することによって、その組に対応する閾値組み合わせを読み出して予測処理を実行することができる。

【0077】

このように、本変形例では、所定の基準（既存知識など）を基に類似性が想定されるデータをまとめて、学習が行われる。これにより、ある対象エリアで学習用のデータとして用いる事象の履歴データの数が少ないケースであっても、当該対象エリアと類似する特徴を有する他のエリアの履歴データによってデータの不足を補うことができる。その結果、学習によって得られるモデルの予測精度を高める効果が期待できる。

【0078】

［第２実施形態］
本実施形態は、以下の点を除き、第１実施形態と同様の構成を有する。

【0079】

＜機能構成例＞
図１８は、第２実施形態における情報処理装置１０の機能構成例を示す図である。図１８に示されるように、本実施形態の情報処理装置１０は、評価指標可視化部１５２を更に有する。評価指標可視化部１５２は、評価指標生成部１５０により生成された評価指標の可視化情報を表示装置に出力する。本実施形態において、評価指標可視化部１５２は、例えば、評価指標生成部１５０により生成された評価指標の大きさに基づいて、可視化情報の表示態様を変更する。

【0080】

評価指標可視化部１５２は、例えば、図１９に示されるような評価指標の可視化情報を表示装置に表示する。図１９は、評価指標可視化部１５２が表示する評価指標の可視化情報の一例を示す図である。図１９の例において、評価指標可視化部１５２は、ラベル値組み合わせ集合の要素ＴＨ^ｌ（ｔｈ１^ｌ，ｔｈ２^ｌ）（ｌ＝１，２，・・・，Ｌ）のそれぞれのｔｈ１^ｌとｔｈ２^ｌとの組み合わせに対応する位置に、ｔｈ１^ｌとｔｈ２^ｌとの組み合わせについて評価指標生成部１５０により生成された評価指標値に関連する数値を半径とする、円形のマーカーを表示した散布図を、可視化情報として生成している。図１９に例示される可視化情報によれば、例えば、以下が把握される。
（１）第１の計算手法の評価値ｔｈ１^ｌが小さく、かつ、第２の計算手法の評価値ｔｈ２^ｌが大きいほど、評価指標は大きくなる。
（２）図１９の右上の範囲での組み合わせであれば、評価指標はそれほど大きく変化しない。

【0081】

このように、評価指標を可視化することによって、事象の発生を予測する際に利用する閾値組み合わせに対する評価指標の変化傾向を定性的に把握できるようになる。このような情報は例えば、以下のように利用できる。一例として、評価指標が大きく変化しない閾値組み合わせが複数存在する場合、閾値組み合わせ特定部１６０は、それらの組み合わせに係る閾値に基づく統計値（例えば、平均値など）を計算し、その計算結果を予測に利用する閾値組み合わせとして特定してもよい。そのほかにも、閾値組み合わせ特定部１６０は、事象の発生を予測する計算手法がｍ個存在する場合に、図１９に例示されるような可視化情報を使って、評価指標に対する寄与の小さい計算手法に関する閾値を閾値組み合わせに含めないように設定する、あるいは、寄与傾向が類似している複数の計算手法については、その中のいずれか１つを選択して閾値組み合わせを特定するようにしてもよい。また、予測結果出力部１７０は、事象の発生を予測する計算手法がｍ個存在する場合に、図１９に例示されるような可視化情報を使って、評価指標に対する寄与の小さい計算手法を予測に使用しないように設定する、あるいは、寄与傾向が類似している複数の計算手法については、その中のいずれか１つを選択するようにしてもよい。このようにすることで、行列の次元を削減して、計算処理量の削減効果が期待できる。

【0082】

なお、評価指標可視化部１５２が生成する可視化情報は、組み合わせ毎の評価指標の大小関係を把握し易い情報であればよく、図１９の例に制限されない。例えば、評価指標可視化部１５２は、すべての組み合わせの中で最大の評価指標値に対する割合やその自然対数をとった値に基づくグラフを可視化情報として表示してもよい。

【0083】

［第３実施形態］
本実施形態は、以下の点を除き、第１実施形態または第２実施形態と同様の構成を有する。以下では、第１実施形態の構成をベースとする構成について説明する。

【0084】

＜機能構成例＞
図２０は、第３実施形態における情報処理装置１０の機能構成例を示す図である。図２０に示されるように、本実施形態の情報処理装置１０は、予測ルール出力部１８０をさらに備える。予測ルール出力部１８０は、閾値組み合わせ特定部１６０によって特定された閾値組み合わせに基づいて、セルの属性を分類する処理（例：図４のＳ１２４の処理）における分類ルールに関する情報を表示装置に出力する。

【0085】

予測ルール出力部１８０は、例えば、図２１に例示されるような情報を表示装置に出力する。図２１は、予測ルール出力部１８０により出力される情報の一例を示す図である。図２１の例において、予測ルール出力部１８０は、閾値組み合わせ特定部１６０により特定された閾値組み合わせに基づいて、セルの属性を分類するルール（決定木）を示す情報を生成し、表示装置に出力している。なお、図２１においてＲ１（ｐ，ｇ）とは、予測時点ｔ^ｐにおけるセルｇの事象の発生し易さを第１の計算手法で計算した値の評価値を意味する。また、図２１においてＲ２（ｐ，ｇ）とは、予測時点ｔ^ｐにおけるセルｇの事象の発生し易さを第２の計算手法で計算した値の評価値を意味する。予測結果出力部１７０による予測では、この決定木に従って、各セルがＡ〜Ｄのいずれかに分類される。なお、一般的に、計算手法の数がｍ個の場合、予測ルール出力部１８０により表示される決定木は、ｍ層の決定木となる。この場合、最上位層に各セルのｍ次元の評価値ベクトルが入力され、最終的に、各セルが２^ｍの分類のいずれかに振り分けられる。

【0086】

本実施形態の構成によれば、予測結果出力部１７０での予測ルールが可視化される。本発明の情報処理装置１０を使って事象の発生を予測する人物は、このように可視化された予測ルールを見ることによって、予測結果がどのように分類されるかについて簡単に理解できるようになる。

【0087】

以上、図面を参照して本発明の実施の形態について述べたが、本発明はこれらに限定されて解釈されるべきものではなく、本発明の要旨を逸脱しない限りにおいて、当業者の知識に基づいて、種々の変更、改良等を行うことができる。実施形態に開示されている複数の構成要素は、適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素からいくつかの構成要素を削除してもよいし、異なる実施形態の構成要素を適宜組み合わせてもよい。

【0088】

また、上述の説明で用いた複数のフローチャートでは、複数の工程（処理）が順番に記載されているが、各実施形態で実行される工程の実行順序は、その記載の順番に制限されない。各実施形態では、図示される工程の順番を内容的に支障のない範囲で変更することができる。また、上述の各実施形態は、内容が相反しない範囲で組み合わせることができる。

【0089】

上記の実施形態の一部または全部は、以下の付記のようにも記載されうるが、以下に限られない。
１．
過去に発生した事象の履歴データのうち第１の期間に含まれる第１の履歴データに基づいて複数の計算手法それぞれについて計算される対象エリアの各セルにおける事象の発生し易さと、前記第１の期間よりも後の第２の期間に含まれる第２の履歴データとに基づいて、特徴量集合を前記計算手法別に生成する手法別特徴量生成手段と、
前記計算手法別に生成された特徴量集合を各次元のラベルとして有し、かつ、各行列要素に対応するセルにおける前記第２の履歴データに基づく前記事象の発生件数を行列要素として有する第１行列を生成する第１行列生成手段と、
前記第１行列と同じラベルを有し、かつ、各行列要素に対応する各次元のラベル値の合計値を行列要素として有する第２行列を生成する第２行列生成手段と、
所定の閾値に基づいて、前記第２行列において行列要素の境界を設定する境界設定手段と、
前記境界の位置と対応する位置に存在する前記第１行列の行列要素である境界行列要素のそれぞれについて、前記境界行列要素を基準として前記第１行列を複数の領域に分割し、前記複数の領域の少なくとも１つの領域の行列要素の和を用いた評価指標を生成する評価指標生成手段と、
前記境界行列要素それぞれについて生成された評価指標に基づいて、評価が最大となる前記境界行列要素に対応する各次元のラベル値の組み合わせを、前記事象の予測に用いる閾値組み合わせとして特定する閾値組み合わせ特定手段と、
前記閾値組み合わせに基づいて、所定の予測時点における前記対象エリアの各セルの事象の発生し易さの予測結果を分類し、その分類結果を含む情報を出力する予測結果出力手段と、
を備える情報処理装置。
２．
前記予測結果出力手段は、前記計算手法別に計算した前記結果を分類する処理を、前記事象の種類ごとに実行する、
１．に記載の情報処理装置。
３．
前記評価指標の可視化情報を表示装置に出力する評価指標可視化手段を更に備える、
１．または２．に記載の情報処理装置。
４．
前記評価指標可視化手段は、前記評価指標の大きさに基づいて、当該評価指標の表示態様を変更する、
３．に記載の情報処理装置。
５．
コンピュータが、
過去に発生した事象の履歴データのうち第１の期間に含まれる第１の履歴データに基づいて複数の計算手法それぞれについて計算される対象エリアの各セルにおける事象の発生し易さと、前記第１の期間よりも後の第２の期間に含まれる第２の履歴データとに基づいて、特徴量集合を前記計算手法別に生成し、
前記計算手法別に生成された特徴量集合を各次元のラベルとして有し、かつ、各行列要素に対応するセルにおける前記第２の履歴データに基づく前記事象の発生件数を行列要素として有する第１行列を生成し、
前記第１行列と同じラベルを有し、かつ、各行列要素に対応する各次元のラベル値の合計値を行列要素として有する第２行列を生成し、
所定の閾値に基づいて、前記第２行列において行列要素の境界を設定し、
前記境界の位置と対応する位置に存在する前記第１行列の行列要素である境界行列要素のそれぞれについて、前記境界行列要素を基準として前記第１行列を複数の領域に分割し、前記複数の領域の少なくとも１つの領域の行列要素の和を用いた評価指標を生成し、
前記境界行列要素それぞれについて生成された評価指標に基づいて、評価が最大となる前記境界行列要素に対応する各次元のラベル値の組み合わせを、前記事象の予測に用いる閾値組み合わせとして特定し、
前記閾値組み合わせに基づいて、所定の予測時点における前記対象エリアの各セルの事象の発生し易さの予測結果を分類し、その分類結果を含む情報を出力する、
ことを含む情報処理方法。
６．
前記コンピュータが、前記計算手法別に計算した前記結果を分類する処理を、前記事象の種類ごとに実行する、
ことを含む５．に記載の情報処理方法。
７．
前記コンピュータが、前記評価指標の可視化情報を表示装置に出力する、
ことを含む５．または６．に記載の情報処理方法。
８．
前記コンピュータが、前記評価指標の大きさに基づいて、当該評価指標の表示態様を変更する、
ことを含む７．に記載の情報処理方法。
９．
コンピュータに、５．から８．のいずれか１つに記載の情報処理方法を実行させるためのプログラム。

【符号の説明】

【0090】

１０ 情報処理装置
１１０ 手法別特徴量生成部
１２０ 第１行列生成部
１３０ 第２行列生成部
１４０ 境界設定部
１５０ 評価指標生成部
１５２ 評価指標可視化部
１６０ 閾値組み合わせ特定部
１７０ 予測結果出力部
１８０ 予測ルール出力部
１０１０ バス
１０２０ プロセッサ
１０３０ メモリ
１０３０ 例：メモリ
１０４０ ストレージデバイス
１０５０ 入出力インタフェース
１０６０ ネットワークインタフェース

【要約】

【課題】コンピュータによる事象の予測の根拠を説明可能とする。
【解決手段】情報処理装置１０は、事象の履歴データを基に、事象の発生し易さを計算する複数の計算手法別に特徴量集合を生成し、複数の計算手法別の特徴量集合を各次元のラベルとして有し、かつ、各行列要素に対応するセルにおける事象の発生件数を行列要素として有する第１行列を生成し、第１行列と同じラベルを有し、かつ、各行列要素に対応するラベル値の合計値を行列要素として有する第２行列を生成し、所定の閾値を基に第２行列において設定された境界に対応する位置に存在する境界行列要素を基に第１行列を複数の領域に分割した結果それぞれについて評価指標を生成し、事象の予測に用いる閾値組み合わせを各評価指標を基に特定し、ある予測時点での対象エリアの各セルの事象の予測結果を閾値組み合わせを基に分類し、その結果を含む情報を出力する。
【選択図】図１

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】

【図20】

【図21】