(57)【特許請求の範囲】
【請求項１】
ユーザの移動履歴とスポットに関連付けられた属性情報からユーザの属性及び行動を推定する行動及び属性推定装置であって、
スポット毎に属性(ジャンル)情報が関連付けられた情報を格納するスポット情報記憶手段と、
ユーザの位置情報収集デバイスから測位時刻、位置情報の時系列データを含む前記ユーザの移動履歴を取得する受信手段と、
前記ユーザの移動履歴から前記ユーザの滞在エリアを特定する滞在エリア特定手段と、
特定された前記滞在エリア内に存在するスポット数とスポットに対応するジャンルの数(スポットジャンル数)の比率からスコアを求めるスコア算出手段と、
前記スコアと所定の条件を比較してユーザの属性及び行動を推定する推定手段と、
を有し、
前記スコア算出手段は、
移動モデルと超幾何分布によるスポット選択モデルを用いる方法、または、訪問可能性の距離減衰を反映した確率分布推定による方法により、滞在エリアのジャンル出現希少性の値を求めるジャンル出現希少性算出手段と、
前記ジャンル出現希少性の値を利用して該滞在エリア内に存在するスポット数とスポットジャンル数の比率をスコアとして計算する手段を含む
ことを特徴とする行動及び属性推定装置。
【請求項２】
前記ジャンル出現希少性算出手段は、
一定期間の滞在エリア群から作成した滞在エリアに対応する移動範囲内(移動候補エリア)のスポットの数と、該滞在エリアの移動候補エリアのジャンルの数を用いて前記ジャンル出現希少性の値を求める第１の手段、
または、
訪問可能性として、滞在エリアに対応する短距離移動クラスタのコアポイントと滞在エリアの距離による減少関数に従う確率分布推定により前記ジャンル出現希少性の値を求める第２の手段を有する
請求項１記載の行動及び属性推定装置。
【請求項３】
前記第１の手段は、
ユーザの一定期間の全滞在エリアの中心点から重心を算出し、前記移動候補エリアの中心点とし、一定量の滞在エリアが含まれる距離を半径とした円エリアを滞在エリアの移動候補エリアとして前記スポット数を求める手段、
または、
ユーザの一定期間の全滞在エリアの中心点に対して、長距離移動クラスタ、短距離移動クラスタ、頻繁訪問地クラスタを段階的に適用して前記移動候補エリアを求め、該移動候補エリアにより滞在エリア毎に前記スポット数を求める手段
を含む請求項２記載の行動及び属性推定装置。
【請求項４】
前記スコア算出手段は、
前記スコアに、滞在エリアの中心時間におけるジャンルに関連付けられた行動実施率を乗じた値をスコアとする手段を含む
請求項１乃至３のいずれか１項に記載の行動及び属性推定装置。
【請求項５】
前記スコア算出手段は、
前記スコアに、ＧＰＳ誤差を滞在エリアを中心とした正規分布で仮定した場合のスポットの位置の上側確率を滞在エリア内における全スポット集合で除した値を乗じた値をスコアとする手段を含む
請求項1乃至４のいずれか１項に記載の行動及び属性推定装置。
【請求項６】
前記スコア算出手段は、
前記スコアに、ジャンル毎に、前記ユーザの移動履歴から抽出した滞在エリアと比較した場合の地理的または時間的に近いものがないという観点により求められた指標を乗じた値をスコアとする手段を含む
請求項１乃至５のいずれか１項に記載の行動及び属性推定装置。
【請求項７】
ジャンルと行動が対応付けられた条件を取得し、該条件で指定されたジャンルに基づいて現在の滞在エリア内のジャンル毎、または、複数のジャンルの各スポットのスコアを合算し、第１のスポット合算値とし、該第１のスポット合算値を該条件で指定された一定期間内において合算し、第２のスポット合算値とし、該第２のスポット合算値が閾値または所定の順位の条件を満たしている場合に該第２の合算値を行動推定結果として出力するユーザ行動推定手段を更に有する
請求項１乃至６のいずれか１項に記載の行動及び属性推定装置。
【請求項８】
滞在エリア毎に、前記スコア算出手段で求められたジャンル毎のスコアを合算し第１のジャンル合算値とし、該第１のジャンル合算値を一定期間内における、Stouffer-methodを用いた確率結合により合算し第２のジャンル合算値とし、該第２のジャンル合算値が所定の条件を満たしている場合に該第２のジャンル合算値をユーザ属性として出力するユーザ属性推定手段を更に有する
請求項１乃至７のいずれか１項に記載の行動及び属性推定装置。
【請求項９】
ユーザの移動履歴とスポットに関連付けられた属性情報からユーザの属性及び行動を推定する行動及び属性推定方法であって、
受信手段が、ユーザの位置情報収集デバイスから測位時刻、位置情報の時系列データを含む前記ユーザの移動履歴を取得する受信ステップと、
滞在エリア特定手段が、前記ユーザの移動履歴から前記ユーザの滞在エリアを特定する滞在エリア特定ステップと、
スコア算出手段が、特定された前記滞在エリア内に存在するスポット数とスポットに対応するジャンルの数(スポットジャンル数)の比率をスコアとして求めるスコア算出ステップと、
推定手段が、前記スコアと所定の条件を比較してユーザの属性及び行動を推定する推定ステップと、
を行い、
前記スコア算出ステップにおいて、
前記滞在エリア特定ステップで特定された前記滞在エリアに基づいて、スポット毎に属性(ジャンル)情報が関連付けられた情報が格納されたスポット情報記憶手段を参照して、ユーザの滞在エリア内におけるスポット及びスポットジャンルの地理分布から算出したスポット数とスポットジャンル数の比率の出現確率分布によりジャンル出現希少性の値を求め、該ジャンル出現希少性の値を利用して該滞在エリア内に存在するスポット数とスポットジャンル数の比率をスコアとして計算する
ことを特徴とする行動及び属性推定方法。
【請求項１０】
請求項１乃至８のいずれか１項に記載の行動及び属性推定装置を構成する各手段としてコンピュータを機能させるための行動及び属性推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、行動及び属性推定装置及び方法及びプログラムに係り、特に、ユーザの移動履歴（ＧＰＳ(Global Positioning system)履歴）と、店舗等のスポットに関連付いた属性情報と組み合わせた情報から、ユーザの属性、状況を推定する行動及び属性推定装置及び方法及びプログラムに関する。
【背景技術】
【０００２】
ユーザの行動や状況に合わせて適切なサービス・コンテンツをリコメンドする技術が望まれている。このため、書籍通販サイトにて、サイト内での書籍情報閲覧履歴からユーザの興味情報を推定して書籍をリコメンドする等、多様な履歴情報からユーザの趣味嗜好等の属性や状況を推定する技術が提案されている。ここで、このような取り組みにおいてWEB閲覧履歴については入手性が高いこともあり検討があるが、移動履歴については、昨今、携帯電話への搭載等にてＧＰＳ移動履歴の入手性は高まっているが比較的取り組みが少ない。ここで、移動履歴を利用する場合については昨今スマートフォン等の普及により日中全時間帯の履歴をとることが安易となっているため、ユーザの状況に適応するコンテキストアウェアなサービスが実現できる可能性もあり潜在的な有用性は高い。
【０００３】
移動履歴での取り組みが少ない理由は、ユーザ属性、状況等を推定するためには、店舗/スポットの訪問を把握することが有効であるが、ＧＰＳの誤差により、どの店舗/スポットに訪問したか把握することは技術的難易度が高いためであると考えられる。尚、このＧＰＳ誤差は、見通しの良い屋外では1m以下を期待できるが、市街地での店舗への訪問（屋内）を識別したい場合には、数百mの誤差を想定する必要がある。
【０００４】
この課題を解決する方法として、ＧＰＳ履歴から店舗訪問の可能性が高い履歴として、一定時間の滞留を抽出しエリア化（以降、「滞在エリア」と記す）し、エリア内に存在するスポットに設定されているスポットの属性（キーワード）の出現頻度を元に訪問可能性を考えて統計的にユーザ属性を推定する方法が考えられる。尚、この方法で利用するスポットの情報は、ＷＥＢ等に広く存在するデータから、各スポットの位置情報（緯度・経度等）とスポット属性（ラーメン屋、焼き肉屋等）を収集する、もしくは市販の店舗情報ＤＢを購入することを想定する。
【０００５】
このようなユーザ履歴における「キーワード出現頻度」から、ユーザ嗜好や属性を推定する点で共通する関連研究が、ＷＥＢ閲覧履歴に対して行われている（例えば、非特許文献１参照）。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】松尾豊、福田隼人、石塚満  ユーザ個人の閲覧履歴からのキーワード抽出によるブラウジング支援  人工知能学会論文誌,Vol 18 No4 pp 203-211,2003
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、上記の従来技術では、ＷＥＢ履歴を対象としており、ＷＥＢ履歴には存在しない強い距離制約のような移動履歴の特徴に対応しておらず利用することができない。
【０００８】
本発明は、上記の点に鑑みなされたもので、GPSの履歴からユーザの滞在エリア内に存在する各スポットに関連付けられた属性、ユーザの行動を精度よく推定することが可能な行動及び属性推定装置及び方法及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記の課題を解決するため、本発明は、  ユーザの移動履歴とスポットに関連付けられた属性情報からユーザの属性及び行動を推定する行動及び属性推定装置であって、
スポット毎に属性(ジャンル)情報が関連付けられた情報を格納するスポット情報記憶手段と、
ユーザの位置情報収集デバイスから測位時刻、位置情報の時系列データを含む前記ユーザの移動履歴を取得する受信手段と、
前記ユーザの移動履歴から前記ユーザの滞在エリアを特定する滞在エリア特定手段と、
特定された前記滞在エリア内に存在するスポット数とスポットに対応するジャンルの数(スポットジャンル数)の比率からスコアを求めるスコア算出手段と、
前記スコアと所定の条件を比較してユーザの属性及び行動を推定する推定手段と、
を有し、前記スコア算出手段は、
前記滞在エリア特定手段で特定された前記滞在エリアに基づいて前記スポット情報記憶手段を参照して、ユーザの滞在エリア内におけるスポット及びスポットジャンルの地理分布から算出したスポット数とスポットジャンル数の比率の出現確率分布によりジャンル出現希少性の値を求め、該ジャンル出現希少性の値を利用して該滞在エリア内に存在するスポット数とスポットジャンル数の比率をスコアとして計算する手段を含む。
【発明の効果】
【００１０】
上記のように本発明によれば、ユーザの滞在エリア内に出現する希少性（珍しさ）を考慮した、ユーザが訪問したと思われるスポットのスコアを算出することができ、このスコアを用いてジャンル（つまり、ユーザの嗜好）及びユーザの行動の推定が可能となる。
【図面の簡単な説明】
【００１１】
【図１】本発明の処理のイメージである。
【図２】本発明の推定処理の概要フローである。
【図３】本発明の一実施の形態におけるシステム構成図である。
【図４】本発明の一実施の形態における各装置の構成図である。
【図５】本発明の一実施の形態におけるシステム全体の処理のフローチャートである。
【図６】本発明の一実施の形態におけるＧＰＳ情報受信〜滞在エリア抽出処理のフローチャートである。
【図７】本発明の一実施の形態におけるＧＰＳログソースデータの例である。
【図８】本発明の一実施の形態におけるＧＰＳデータテーブルの例である。
【図９】本発明の一実施の形態における滞在エリアデータテーブルの例である。
【図１０】本発明の一実施の形態におけるスポットデータベースの例である。
【図１１】本発明の一実施の形態における滞在エリア内ジャンルテーブルの例である。
【図１２】本発明の一実施の形態における重心からの距離分布による移動モデルを示す図である。
【図１３】本発明の一実施の形態における３層クラスタによる移動モデルを示す図である。
【図１４】本発明の一実施の形態におけるユーザ訪問スポット推定結果の例である。
【図１５】本発明の一実施の形態における訪問スポット推定部の処理のフローチャートである。
【図１６】本発明の一実施の形態における行動推定部の処理のフローチャートである。
【図１７】本発明の一実施の形態における属性（嗜好）推定部の処理のフローチャートである。
【図１８】本発明の一実施の形態における出力される属性リストの例である。
【図１９】本発明の一実施の形態における推定結果送信部の処理のフローチャートである。
【図２０】本発明の一実施の形態におけるユーザ促成・状況推定の例である。
【発明を実施するための形態】
【００１２】
以下図面と共に、本発明の実施の形態を説明する。
【００１３】
まず、本発明の概要を説明する。
【００１４】
図１は、本発明の処理のイメージを示す。
【００１５】
ユーザがＧＰＳ機能を有する携帯端末等を操作してレストランを検索する状況を考える。
【００１６】
ユーザからの要求が入力されると、本発明の行動・属性推定装置では、ＧＰＳの移動履歴を読み込み、ユーザの嗜好（属性）、ユーザの訪問場所（状況)を推定し、ユーザの嗜好、状況に合ったレストランを検索結果として出力する。
【００１７】
図２は、本発明の推定処理の概要フローを示す。
【００１８】
ステップ１）  取得したＧＰＳ履歴から測位時刻、位置情報（緯度、経度）の時系列データ+を取得すると共に、周辺地域全てのスポット(例えば、レストラン店舗)に対応するレストラン識別子、ジャンル、緯度、経度をスポットジャンル地理分布データとして読み込む。
【００１９】
ステップ２）  ＧＰＳ履歴から一定時間、一定の場所に滞留する滞在エリアを、既存のDBSCANアルゴリズムを用いて抽出する。データ形式は、中心点、滞在開始時刻、滞在終了時刻、滞在エリア半径からなり、このうち、滞在エリア半径は分析者が入力するパラメータである。抽出された滞在エリアの情報と予め用意されているスポットジャンルと関連付ける。
【００２０】
なお、スポットジャンルは、スポットジャンルが格納されているデータベースまたは、インターネット等から取得するものとする。本実施の形態では、スポットデータベースを用いるものとする。
【００２１】
ステップ３）  滞在エリアにおけるスポットジャンルの出現頻度をスコア化する。このとき、滞在エリア内に出現する希少性(珍しさ)を考慮してユーザが訪問したと思われるスポットのスコアを算出する。
【００２２】
ステップ４）  スコアに基づいて、ユーザの訪問スポット、行動、嗜好を推定する。
【００２３】
図３は、本発明の一実施の形態におけるシステム構成を示す。
【００２４】
本発明のシステムは、上記の図２の処理を行う推定処理装置１００、利用者1が利用するＧＰＳ機能付き携帯端末等の位置情報収集デバイス２００、サービス処理装置３００から構成され、これらは通信ネットワークを介して接続されている。
【００２５】
推定処理装置１００は、位置情報収集デバイス２００からＧＰＳログを取得し、推定処理を行い、選好キーワードをサービス処理装置３００に出力する。サービス処理装置３００は、取得した選好キーワードに基づいて関連するコンテンツリコメンド等を位置情報収集デバイス２００に送信する。
【００２６】
以下に、上記のシステムの各装置構成を示す。
【００２７】
図４は、本発明の一実施の形態における各装置の構成を示す。
【００２８】
推定処理装置１００は、ＧＰＳ情報受信部１１０、ＧＰＳ履歴前処理部１２０、スコア算出部１３０、訪問スポット推定部１４０、行動推定部１５０、属性（嗜好）推定部１６０、推定結果送信部１７０、スポットデータベース１８０、記憶部１９０から構成される。
【００２９】
なお、記憶部１９０には各種テーブルが格納される。
【００３０】
位置情報収集デバイス２００は、ＧＰＳログ収集部２１０、ＧＰＳ情報を推定処理装置１００に送信するＧＰＳ情報送信部２２０、サービス処理装置３００からのリコメンド結果を表示するコンテンツ提示部２３０から構成される。
【００３１】
サービス処理装置３００は、リコメンド結果を位置情報収集デバイス２００に送信するコンテンツ送信処理部３１０、リコメンド対象コンテンツ決定処理部３２０、推定処理装置１００から推定結果を取得する推定結果受信部３３０から構成される。
【００３２】
次に、上記のシステムの全体の動作を説明する。
【００３３】
図５は、本発明の一実施の形態におけるシステム全体の処理のフローチャートである。
【００３４】
ステップ１１１０）  推定処理装置１００のＧＰＳ情報受信部１１０は、位置情報収集デバイス２００のＧＰＳ情報送信部２２０からＧＰＳ履歴情報を受信する。
【００３５】
ステップ１１２０）  ＧＰＳ履歴前処理部１２０は、受信したＧＰＳ履歴から位置情報収集デバイス２００のユーザの滞在エリアを抽出し、新しい滞在エリアを抽出した場合はステップ１１３０に移行し、抽出できない場合は、当該処理を終了する。詳細な処理については後述する。
【００３６】
ステップ１１３０）  スコア算出部１３０において、滞在エリア内の各スポットスコアを求める。詳細な処理については後述する。
【００３７】
ステップ１１４０）  訪問スポット推定部１４０は、各スポットのスコア値からユーザの訪問スポットを推定する。
【００３８】
ステップ１１５０）  行動推定部１５０は、ユーザの訪問ジャンル、行動を推定する。
【００３９】
ステップ１１６０）  属性(嗜好)推定部１６０は、過去の一定期間のスコアを統合し、ユーザ属性を推定する。新しい滞在エリアで閾値を超える等の条件を満たした推定結果（訪問スポット、行動、属性）があるかを判定し、ない場合は当該処理を終了する。
【００４０】
ステップ１１７０）  推定結果送信部１７０は、推定結果をサービス処理装置３００に送信する。必要であれば、ユーザ属性情報をスコア算出部１３０に出力し、スコア式に反映する。
【００４１】
なお、上記のステップ１１４０〜１１６０処理は不要であれば省略することも可能である。
【００４２】
以下に、上記の処理をステップ毎に詳細に説明する。
【００４３】
＜GPS情報受信部、GPS履歴前処理部の処理＞
まず、上記のステップ１１１０〜１１２０について説明する。
【００４４】
図６は、本発明の一実施の形態におけるＧＰＳ情報受信〜滞在エリア抽出処理のフローチャートである。
【００４５】
ステップ１１１１）  ＧＰＳ受信部１１０は、位置情報収集デバイス２００からＧＰＳ情報を受信し、ＧＰＳ履歴前処理部１２０に出力する。図７は、受信したＧＰＳデータ形式の例であり、経度、緯度、時間が含まれる。
【００４６】
ステップ１１２１）  ＧＰＳ履歴前処理部１２０は、記憶部１８０の図８に示すＧＰＳデータテーブルに受信した新規ＧＰＳデータを保存し、過去時系列のＧＰＳデータを読み出す。
【００４７】
ステップ１１２２）  ＧＰＳ履歴前処理部１２０は、ＧＰＳデータテーブルから読み出した時系列のＧＰＳ履歴データの緯度、経度、測定時刻から一定時間内及び一定範囲内で滞留しているデータ部分を抽出する。滞在エリアの抽出については、例えば、滞留時間閾値と、各点間距離閾値をパラメータとしてクラスタを抽出するDBSCAN法等のアルゴリズムを用いるものとする。抽出された滞在エリアのデータは、記憶部１８０の図９に示す滞在エリアデータテーブルに保存すると共に、スコア算出部１３０に送信する。ここで、抽出された滞在エリアのデータは、中心緯度、中心経度、滞在エリア半径、滞留開始時間、滞留終了時間を含む
ステップ１１２３）  ＧＰＳ履歴前処理部１２０は、抽出した滞在エリア情報について、中心点から一定半径（分析者が設定するパラメータ）内に存在するスポットとそのスポットに関連付けられたジャンル（例えば、「レストラン」）を関連付けたデータを、スポット識別子に基づいて図１０に示すスポットデータベース１８０を参照することにより作成し、図１１に示す滞在エリア内ジャンルテーブルに保存すると共に、スコア算出部１３０に出力する。
【００４８】
＜スコア算出部の処理＞
次に、上記のＧＰＳ履歴前処理部１２０から、滞在エリアのデータ及び、スポット及び当該スポットにジャンルを関連付けたデータを取得したスコア算出部１３０の処理について説明する。
【００４９】
スコア算出部１３０は、滞在エリアｊにおけるジャンルｉを持つスポットのスコア（Ｓijx）を求める。
【００５０】
スコアＳijxは、
・現時刻による行動確率；
・滞在エリアとの距離による訪問確率；
・周辺ジャンル分布比較による希少性；
・過去ユーザＧＰＳ履歴から抽出した滞在エリアとの比較による現滞在エリア希少性；
の少なくとも１つ、または、いずれか、または、全てを乗じたものである。全てを乗じた場合の計算式は
【００５１】
【数１】
0005410462.tif 000002
上記において、スポットｘは滞在エリア中心点から一定範囲(滞在エリア内)のスポットである。Xは、滞在エリアj内のスポット総数とする。wは滞在エリアj内の各スポットである。φ-1（x）は累積標準正規分布の逆関数である。右辺の第１項E(t)は、滞在エリアの中心時間におけるジャンルに関連付けられる行動実施率であり、例えば、ｉが「グルメジャンル」の場合は、時間帯毎の食事行動を行う確率となる。但し「グルメ」は一例であり、技術的には、カラオケや学校等の多様なジャンルがある。このとき、「カラオケ」であれば遊び(交友)、「学校」であれば学習という行動が関連付けられることが考えられる。なお、この確率については、統計庁の社会生活時間調査等の外部情報を元に、時間帯毎の平均的な実施率のデータが公開されており、それら（外部情報)を利用することを想定している。
【００５２】
また、上記の右辺の第２項の分子のDist(x,y)はＧＰＳ誤差をｙ点を中心とし、μ=0,σ=（ＧＰＳ機器の想定誤差による値）で決定した場合の正規分布によるｘ地点の出現確率（x,yの距離の出現確率）である。第２項の分母のｗは滞在エリアｊ内における全スポット集合である。loc(x)はスポットｘの位置情報であり、ｎ，ｍは重みパラメータ（乗数）である。右辺の第３項のＭ(i,j)は滞在エリアｊのジャンルｉの出現数の珍しさを定量化したものである。詳細については後述する。右辺の第４項のR(i,j)はユーザの過去の履歴内から抽出した滞在エリアから得られる、ユーザ嗜好(意図、属性、行動等)に関する情報は本手法では減衰する考え方でスコアを下げる（例えば、学校、会社、バイト等によるルーティン的な訪問）。Ｒ(j)の定義例として、
R(j)=(1-(interval(j)*coverage(j))
とする。ここでinterval(j)は過去の滞在エリア内において現在滞在エリアとの中心点間距離が閾値以下の滞在エリアで、滞在エリアの滞在時刻が近いものがあるほど１に近づき、遠いほど(例えば、２週間以上)０に近づく関数である。Coverage(j)は過去滞在エリア内において現滞在エリアとの中心点距離が近いものがあるほど１に近づき、一定以上遠い場合は０になる関数である。
【００５３】
上記の滞在エリアｊのジャンルｉの出現数の珍しさＭ(i,j)の求め方について説明する。
【００５４】
まず、第１のＭ(i,j)を求める方法は、移動モデルと超機可分布によるスポット選択モデル用いるものである。当該方法は、一定期間の滞在エリア群から作成した移動モデルを利用して抽出した選択候補スポット集合から超機可分布を利用して滞在エリア毎のジャンル出現数の珍しさを以下のように定義する。
【００５５】
【数２】
0005410462.tif 000003
上記において、Tjは滞在エリアjに対応する移動範囲内（以降、移動候補エリア）のスポットの数、Lijは滞在エリアjの移動候補エリア内のジャンルｉの数、Ｎijは滞在エリアｊ内のジャンルｉの数、Sjは滞在エリアｊ内のスポットの数である。
【００５６】
このTj、Lijの決定方法は以下2つある。尚、Tjはユーザの選択候補スポット集合の総数、Lijはその中でのキーワードiを持つスポットの数であると見なせる。この選択候補スポット集合を算出するために移動モデルを定義する。移動モデルとは、ユーザの一定期間の滞在エリアの位置分布の特徴から、ユーザが各滞在エリア（＝訪問スポット）の位置を決定する際、候補としたと考える移動候補エリアを決定するモデルである。この移動候補エリア内のスポットが、各滞在エリアについての選択候補スポット集合となる。
【００５７】
Tj、Lij決定方法(1)  重心からの距離分布による移動モデル利用：
第一の移動モデルでは、図１２に示すように、各ユーザの一定期間の全滞在エリア中心点から重心を算出して移動候補エリアの中心点とする。この重心を中心として、一定量(例95%)の滞在エリアが含まれる距離を半径とした円エリアを作成する。この円エリアを各滞在エリアの移動候補エリアとして選択候補スポット集合を抽出する。尚、この移動候補エリアに入らない滞在エリアがある場合は、外れ値データとして処理を行わない。
【００５８】
Tj、Lij決定方法(2)  ３層クラスタによる移動モデル利用：
第二の移動モデルでは、図１３に示すように、各ユーザの一定期間の全滞在エリア中心点に対して、クラスタリング（例えば、DBSCAN法）を3種類のパラメータで段階的に適用することで構造化して移動候補エリアを決定する。この構造化とは、全滞在エリアに対して車や電車などによる「長距離移動クラスタ」と、車や電車を降りた後の徒歩等による「短距離移動クラスタ」、短距離移動クラスタ内にある自宅や職場等の「頻繁訪問地クラスタ」としてクラスタリングする。この構造化は、各移動種別を想定した閾値でDBSCAN法を階層的に適用して作成する。
【００５９】
各クラスタに対する移動候補エリアは、自クラスタの下位クラスタの中で、所属する滞在エリア数最大の下位クラスタの「コアポイント」を中心に、自クラスタに所属する滞在エリアがすべて収まる範囲の円を作成することで決定する。コアポイントとは、各クラスタ内で滞在エリアが集中する地点を移動起点と考えて決定する位置情報である。これは、最小単位である頻繁訪問地クラスタから上位クラスタに向かって順に作成する。尚、最下位の頻繁訪問地クラスタでは、所属する滞在エリアの重心を算出することでコアポイントとする。また、頻繁訪問地クラスタ以外でも下位クラスタがない場合には、自クラスタに所属する滞在エリアの重心をコアポイントとする。
【００６０】
各滞在エリアの移動候補エリアは、滞在エリアが所属するクラスタの移動候補エリアによって決定する。滞在エリアが長距離移動クラスタと短距離移動クラスタ階層に同時に所属する場合は、短距離移動クラスタの移動候補エリアを利用する。尚、頻繁訪問地クラスタについては、短距離移動クラスタのコアポイント算出のためだけに利用し移動候補エリア決定に利用しない。この処理で決定した移動候補エリアにより、滞在エリア毎に選択候補スポット集合を抽出する。
【００６１】
また、Ｖ(j)は、滞在エリアjに対応する短距離移動クラスタのコアポイントと滞在エリアの距離による減少関数である。短距離移動クラスタに所属していない滞在エリアについてはV(j)=1とする。また、後述の第２のＭ(i,j)を求める方法利用時は、V(j)は常に1とする。
【００６２】
【数３】
0005410462.tif 000004
但し、dist(x,y)＝位置ｘと位置ｙの距離を返す関数、corej=滞在エリアｊの所属する短距離移動クラスタのコアポイント位置、centerj=滞在エリアｊの中心点、core_radiusj=滞在エリアｊの所属する短距離移動クラスタの移動候補エリア半径である。
【００６３】
第２のＭ(i,j)を求める方法は、訪問可能性の距離減衰を反映した確率分布推定による方法である。当該方法は、一定期間の移動範囲の中心点から各スポットに対する距離に対する減衰関数に従い、範囲内のスポットを選択する。以下に示すように、そのスポットの中心にジャンルiのスポット数比の存在率を計算し、分布を作成する。
【００６４】
M(i,j)=Z(Y(Nij,Sj))
上記の各変数は第１の方法と同様である。Y(x,y)は一定期間の滞在エリア群の中心点から各スポットとの距離に応じて減衰する選択確率（選択確率最大は中心点であり、選択確率0となる距離はKとする）でスポットを多数回選択し、選択したスポットから一定範囲（ＧＰＳ機器誤差から想定するσの2倍範囲等）のスポット群でのスポットとジャンルの比を算出することで推定するジャンル接触数確率分布に対して、引数の滞在エリアでのジャンル出現比(x/y)の上側確率である。Z(x)は
if(x≧0.5) 0、if(x<0.5) x
である。当該一定期間移動で距離減衰するスポット選択モデルによる指標の例を図１３に示す。上記の一定期間の滞在エリア群の中心点とKの求め方は、分析対象の全滞在エリア中心点から重心を算出して中心点とする。この重心を中心として、一定量(例95%)の滞在エリアが含まれる距離に重みを掛けた値をKとする。
【００６５】
スコア算出部１３０は、上記のようにして求められたスコアを訪問スポット推定部１４０、行動推定部１５０、属性推定部１６０に出力する。
【００６６】
＜訪問スポット推定部の処理＞
次に、図５のステップ１１４０の訪問スポット推定部１４０の処理を説明する。
【００６７】
図１４は、本発明の一実施の形態における訪問スポット推定部の処理のフローチャートを示す。
【００６８】
ステップ１１４１）  訪問スポット推定部１４０は、履歴データ(GPSポイント、抽出滞在エリア数、データ収集期間)が予め設定された閾値を超えているかを判定し、超えている場合はステップ１１４２に移行し、超えていない場合は当該処理を終了する。
【００６９】
ステップ１１４２）  訪問スポット推定部１４０は、スコア算出部１３０から取得した各スポットxが保持するジャンルiに対するスコアSijxが予め設定された条件を超えているかを判定し、超えている場合はステップ１１４３に移行し、超えていない場合は当該処理を終了する。条件としては、
・閾値を超えている；
・各滞在エリアjのスポットの中で上位K番目以内；
等を用いる。
【００７０】
ステップ１１４３）  ステップ１１４２で条件を満たしたスポットを訪問可能性の高いスポットとして推定結果送信部１７０に出力する。
【００７１】
出力結果は、例えば、図１５のようになる。
【００７２】
＜行動推定部の処理＞
次に、図５のステップ１１５０の行動推定部１５０の処理について説明する。
【００７３】
図１６は、本発明の一実施の形態における行動推定部の処理のフローチャートである。
【００７４】
ステップ１１５１）  行動推定部１５０は、記憶部１８０のジャンルスコア行動対応付けテーブル(図示せず)からスコアを算出する条件を決定する。当該ジャンルスコア行動対応付けテーブルには、ジャンル（ジャンル組み合わせもあり）、期間等の合算条件、合算時の重み、推定結果としての閾値等の条件、条件を超えたときに推定結果とする行動が設定されている。
【００７５】
ステップ１１５２）  ステップ１１５１で決定した条件に基づいて、スコア算出部１３０から取得した、現在の滞在エリア内のジャンル毎、もしくは、複数ジャンルの各スポットのスコアを合算する。
【００７６】
ステップ１１５３）  ステップ１１５１で決定した条件に基づいて、ステップ１１５２２で算出した滞在エリア毎のスコアを一定期間等で合算する。
【００７７】
なお、上記のステップ１１５２、１１５３における合算については、過去期間に応じて、もしくはジャンル組み合わせでは、組み合わせジャンルに応じて重み付けしてもよい。
【００７８】
ステップ１１５４）  合算したスコアがステップ１１５１で読み出した以下のような条件を満たしているかを判定し、満たしている場合はステップ１１５５に移行し、満たしていない場合は当該処理を終了する。
【００７９】
条件）
・閾値を越える、または、下回る。
【００８０】
・比較対象内で上位、または、下位になっている。
【００８１】
ステップ１１５５）  ステップ１１５１で読み出した条件を満たす行動を合算値と共に推定結果送信部１７０に出力する。
【００８２】
行動推定部１５０において、ユーザの行動の特徴を判定する方法として、ジャンルiが所定の閾値より高い場合は、ジャンルiについて訪問している可能性が高い滞在エリアが一定期間に複数回ある場合は、強く珍しいことが起こっているため、ジャンルi嗜好があると考え、また、スコアが高く寄与したスポットは訪問していると判断する。
【００８３】
また、ジャンルiスコアが所定の閾値より低いとき、過去一定期間、ジャンルiに近接していない場合はジャンルiに関連付けられる行動(例えば、飲食)をしていない、または、食事時間(例えば、昼食や夕食時間帯)までにジャンルiも含む全食事系ジャンルに近づいていないのであれば、食事自体をしていないと判断する。
【００８４】
＜属性(嗜好)推定部の処理＞
次に、図５のステップ１１６０の属性(嗜好)推定部１６０の処理を説明する。
【００８５】
図１７は、本発明の一実施の形態における嗜好推定部の処理のフローチャートである。
【００８６】
ステップ１１６１）  属性推定部１６０は、スコア算出部１３０から取得した滞在エリア(スポット)毎、ジャンル毎のスコアを取得すると、滞在エリア毎かつジャンル毎にスコアを合算する。
【００８７】
ステップ１１６２）  ステップ１１６１で合算されたスコアについて、所定の期間分のスコアを合算する。期間分の合算については、単純な合計値の算出以外に、各滞在エリアのジャンル毎スコアについて、確率として結合する方法がある。具体的にStouffer-method(Rosental,R. (1984). Meta-analytic procedures for social research.)を利用した例を以下に示す。
【００８８】
【数４】
0005410462.tif 000005
但し、Pi＝分析対象とした一定期間の滞在エリアに対するキーワードｉのスコア、φ（ｘ）：標準正規分布の累積確率密度関数、Ｘij＝滞在エリアｊにおけるキーワードｉを持つ全スポット、Ｊ＝一定期間の全滞在エリア、ｊ＝滞在エリアの識別子とする。
【００８９】
ステップ１１６３）  属性推定部１６０は、合算したスコアが予め設定された以下のような条件を満たすかを判定する。
【００９０】
条件）
・閾値を超える、または、下回る。
【００９１】
・比較対象内で上位、または、下位になっている。
【００９２】
条件を満たす場合はステップ１１６４に移行し、満たさない場合は当該処理を終了する。
【００９３】
ステップ１１６４）  上記の条件を満たしたジャンルを推定結果送信部１７０に出力する。出力される属性リストの例を図１８に示す。
【００９４】
ステップ１１６５）  後述する訪問スポット推定補正を行うために、スコア部１３０に対して、フィードバックする。
【００９５】
＜推定結果送信部の処理＞
次に、図５のステップ１１７０の推定結果送信部１７０の処理を説明する。
【００９６】
図１９は、本発明の一実施の形態における推定結果送信部の処理のフローチャートである。
【００９７】
ステップ１１７１）  推定結果送信部１７０は、訪問スポット推定部１４０、行動推定部１５０、属性推定部１６０から送信された訪問スポット、行動、属性の各推定結果が所定の条件を超えているかを判定し、超えている場合はステップ１１７２に移行し、超えていない場合は当該処理を終了する。
【００９８】
ステップ１１７２）  条件を超えた項目をHTML、XML、CSV等の各形式に纏めたデータを生成する。
【００９９】
ステップ１１７３）  ステップ１１７２で生成されたデータをサービス処理装置３００に推定結果として出力する。
【０１００】
上記のステップ１１６５において、属性推定部１６０から推定結果がスコア算出部１３０にフィードバックされた場合の処理について説明する。
【０１０１】
当該処理は、スポットｘ（ジャンルi保持）訪問スコアの高精度化を図るものである。
【０１０２】
補正する場合の訪問スコアは、現時刻による行動確率、滞在エリアとの距離による訪問確率を乗算した値に、過去の移動履歴からのジャンルiスコアによる補正値、アンケート、Web閲覧履歴等のほかのデータからのジャンルi嗜好による補正値、Web履歴、スケジューラ、メール等の他履歴でのスポットｘ関連キーワード共起からの補正値、近い時間帯の滞在エリアの各ジャンルスコアからの補正値、嗜好の季節変動、時間帯（例えば、昼食と夕食の違い等）変動、曜日変動からの補正値、支払い、決済行動履歴による補正値のいずれかの補正値を乗算することで求められる。
【０１０３】
上記のスコア算出部１３０、訪問スポット推定部１４０、行動推定部１５０、属性(嗜好)推定部１６０の処理により、図２０に示すような例において、ユーザの移動範囲内に、
・ユーザの移動範囲ａ；
・食事時刻以外の滞在エリアｂ；
・食事時刻の滞在だが情報が少ないと見做すエリアｃ；
・食事時刻の滞在エリアｄ；
を考えるとき、エリアｃのように、滞在エリアの中心点から離れすぎる、店舗が多すぎて訪問の特性が判別できない場所はスコアが小さくなる。また、エリアｂ、エリアｃ、エリアｄが重複している場所は職場である可能性が高く、その中に２回も沖縄料理が入っているため、偶然とは考えにくく、ユーザの嗜好性として沖縄料理が好きであると判断できる。
【０１０４】
なお、上記の図４に示す推定処理装置の各構成要素の動作をプログラムとして構築し、推知処理装置として利用されるコンピュータにインストールして実行させる、または、ネットワークを介して流通させることが可能である。
【０１０５】
本発明は上記の実施の形態に限定されることなく、特許請求の範囲内において種々変更・応用が可能である。
【符号の説明】
【０１０６】
１００  推定処理装置
１１０  ＧＰＳ情報受信部
１２０  ＧＰＳ履歴前処理部
１３０  スコア算出部
１４０  訪問スポット推定部
１５０  行動推定部
１６０  属性(嗜好)推定部
１７０  推定結果送信部
１８０  スポットデータベース
１９０  記憶部
２００  位置情報収集デバイス
２１０  ＧＰＳログ収集部
２２０  ＧＰＳ情報送信部
２３０  コンテンツ提示部
３００  サービス処理装置
３１０  コンテンツ送信処理部
３２０  リコメンド対象コンテンツ決定処理部
３３０  推定結果受信部
【図１】
0005410462.tif 000006
【図２】
0005410462.tif 000007
【図３】
0005410462.tif 000008
【図４】
0005410462.tif 000009
【図５】
0005410462.tif 000010
【図６】
0005410462.tif 000011
【図７】
0005410462.tif 000012
【図８】
0005410462.tif 000013
【図９】
0005410462.tif 000014
【図１０】
0005410462.tif 000015
【図１１】
0005410462.tif 000016
【図１２】
0005410462.tif 000017
【図１３】
0005410462.tif 000018
【図１４】
0005410462.tif 000019
【図１５】
0005410462.tif 000020
【図１６】
0005410462.tif 000021
【図１７】
0005410462.tif 000022
【図１８】
0005410462.tif 000023
【図１９】
0005410462.tif 000024
【図２０】
0005410462.tif 000025