"ユーザの移動履歴とスポットに関連付けられた属性情報からユーザの属性及び行動を推定する行動及び属性推定装置であって、スポット毎に属性情報が関連付けられた情報を格納するスポット情報記憶手段と、ユーザの位置情報収集デバイスから測位時刻、位置情報の時系列データを含む前記ユーザの移動履歴を取得する受信手段と、前記ユーザの移動履歴から前記ユーザの滞在エリアを特定する滞在エリア特定手段と、特定された前記滞在エリア内に存在するスポット数とスポットに対応するジャンルの数の比率からスコアを求めるスコア算出手段と、前記スコアと所定の条件を比較してユーザの属性及び行動を推定する推定手段と、を有し、前記スコア算出手段は、移動モデルと超幾何分布によるスポット選択モデルを用いる方法、または、訪問可能性の距離減衰を反映した確率分布推定による方法により、滞在エリアのジャンル出現希少性の値を求めるジャンル出現希少性算出手段と、前記ジャンル出現希少性の値を利用して該滞在エリア内に存在するスポット数とスポットジャンル数の比率をスコアとして計算する手段を含むことを特徴とする行動及び属性推定装置。前記ジャンル出現希少性算出手段は、一定期間の滞在エリア群から作成した滞在エリアに対応する移動範囲内のスポットの数と、該滞在エリアの移動候補エリアのジャンルの数を用いて前記ジャンル出現希少性の値を求める第1の手段、または、訪問可能性として、滞在エリアに対応する短距離移動クラスタのコアポイントと滞在エリアの距離による減少関数に従う確率分布推定により前記ジャンル出現希少性の値を求める第2の手段を有する請求項1記載の行動及び属性推定装置。前記第1の手段は、ユーザの一定期間の全滞在エリアの中心点から重心を算出し、前記移動候補エリアの中心点とし、一定量の滞在エリアが含まれる距離を半径とした円エリアを滞在エリアの移動候補エリアとして前記スポット数を求める手段、または、ユーザの一定期間の全滞在エリアの中心点に対して、長距離移動クラスタ、短距離移動クラスタ、頻繁訪問地クラスタを段階的に適用して前記移動候補エリアを求め、該移動候補エリアにより滞在エリア毎に前記スポット数を求める手段を含む請求項2記載の行動及び属性推定装置。前記スコア算出手段は、前記スコアに、滞在エリアの中心時間におけるジャンルに関連付けられた行動実施率を乗じた値をスコアとする手段を含む請求項1乃至3のいずれか1項に記載の行動及び属性推定装置。前記スコア算出手段は、前記スコアに、GPS誤差を滞在エリアを中心とした正規分布で仮定した場合のスポットの位置の上側確率を滞在エリア内における全スポット集合で除した値を乗じた値をスコアとする手段を含む請求項1乃至4のいずれか1項に記載の行動及び属性推定装置。前記スコア算出手段は、前記スコアに、ジャンル毎に、前記ユーザの移動履歴から抽出した滞在エリアと比較した場合の地理的または時間的に近いものがないという観点により求められた指標を乗じた値をスコアとする手段を含む請求項1乃至5のいずれか1項に記載の行動及び属性推定装置。ジャンルと行動が対応付けられた条件を取得し、該条件で指定されたジャンルに基づいて現在の滞在エリア内のジャンル毎、または、複数のジャンルの各スポットのスコアを合算し、第1のスポット合算値とし、該第1のスポット合算値を該条件で指定された一定期間内において合算し、第2のスポット合算値とし、該第2のスポット合算値が閾値または所定の順位の条件を満たしている場合に該第2の合算値を行動推定結果として出力するユーザ行動推定手段を更に有する請求項1乃至6のいずれか1項に記載の行動及び属性推定装置。滞在エリア毎に、前記スコア算出手段で求められたジャンル毎のスコアを合算し第1のジャンル合算値とし、該第1のジャンル合算値を一定期間内における、Stouffermethodを用いた確率結合により合算し第2のジャンル合算値とし、該第2のジャンル合算値が所定の条件を満たしている場合に該第2のジャンル合算値をユーザ属性として出力するユーザ属性推定手段を更に有する請求項1乃至7のいずれか1項に記載の行動及び属性推定装置。ユーザの移動履歴とスポットに関連付けられた属性情報からユーザの属性及び行動を推定する行動及び属性推定方法であって、受信手段が、ユーザの位置情報収集デバイスから測位時刻、位置情報の時系列データを含む前記ユーザの移動履歴を取得する受信ステップと、滞在エリア特定手段が、前記ユーザの移動履歴から前記ユーザの滞在エリアを特定する滞在エリア特定ステップと、スコア算出手段が、特定された前記滞在エリア内に存在するスポット数とスポットに対応するジャンルの数の比率をスコアとして求めるスコア算出ステップと、推定手段が、前記スコアと所定の条件を比較してユーザの属性及び行動を推定する推定ステップと、を行い、前記スコア算出ステップにおいて、前記滞在エリア特定ステップで特定された前記滞在エリアに基づいて、スポット毎に属性情報が関連付けられた情報が格納されたスポット情報記憶手段を参照して、ユーザの滞在エリア内におけるスポット及びスポットジャンルの地理分布から算出したスポット数とスポットジャンル数の比率の出現確率分布によりジャンル出現希少性の値を求め、該ジャンル出現希少性の値を利用して該滞在エリア内に存在するスポット数とスポットジャンル数の比率をスコアとして計算することを特徴とする行動及び属性推定方法。請求項1乃至8のいずれか1項に記載の行動及び属性推定装置を構成する各手段としてコンピュータを機能させるための行動及び属性推定プログラム。しかしながら、上記の従来技術では、WEB履歴を対象としており、WEB履歴には存在しない強い距離制約のような移動履歴の特徴に対応しておらず利用することができない。本発明は、上記の点に鑑みなされたもので、GPSの履歴からユーザの滞在エリア内に存在する各スポットに関連付けられた属性、ユーザの行動を精度よく推定することが可能な行動及び属性推定装置及び方法及びプログラムを提供することを目的とする。上記の課題を解決するため、本発明は、ユーザの移動履歴とスポットに関連付けられた属性情報からユーザの属性及び行動を推定する行動及び属性推定装置であって、スポット毎に属性情報が関連付けられた情報を格納するスポット情報記憶手段と、ユーザの位置情報収集デバイスから測位時刻、位置情報の時系列データを含む前記ユーザの移動履歴を取得する受信手段と、前記ユーザの移動履歴から前記ユーザの滞在エリアを特定する滞在エリア特定手段と、特定された前記滞在エリア内に存在するスポット数とスポットに対応するジャンルの数の比率からスコアを求めるスコア算出手段と、前記スコアと所定の条件を比較してユーザの属性及び行動を推定する推定手段と、を有し、前記スコア算出手段は、前記滞在エリア特定手段で特定された前記滞在エリアに基づいて前記スポット情報記憶手段を参照して、ユーザの滞在エリア内におけるスポット及びスポットジャンルの地理分布から算出したスポット数とスポットジャンル数の比率の出現確率分布によりジャンル出現希少性の値を求め、該ジャンル出現希少性の値を利用して該滞在エリア内に存在するスポット数とスポットジャンル数の比率をスコアとして計算する手段を含む。上記のように本発明によれば、ユーザの滞在エリア内に出現する希少性珍しさを考慮した、ユーザが訪問したと思われるスポットのスコアを算出することができ、このスコアを用いてジャンルつまり、ユーザの嗜好及びユーザの行動の推定が可能となる。以下図面と共に、本発明の実施の形態を説明する。まず、本発明の概要を説明する。図1は、本発明の処理のイメージを示す。ユーザがGPS機能を有する携帯端末等を操作してレストランを検索する状況を考える。ユーザからの要求が入力されると、本発明の行動属性推定装置では、GPSの移動履歴を読み込み、ユーザの嗜好属性、ユーザの訪問場所状況を推定し、ユーザの嗜好、状況に合ったレストランを検索結果として出力する。図2は、本発明の推定処理の概要フローを示す。ステップ1取得したGPS履歴から測位時刻、位置情報緯度、経度の時系列データを取得すると共に、周辺地域全てのスポットに対応するレストラン識別子、ジャンル、緯度、経度をスポットジャンル地理分布データとして読み込む。ステップ2GPS履歴から一定時間、一定の場所に滞留する滞在エリアを、既存のDBSCANアルゴリズムを用いて抽出する。データ形式は、中心点、滞在開始時刻、滞在終了時刻、滞在エリア半径からなり、このうち、滞在エリア半径は分析者が入力するパラメータである。抽出された滞在エリアの情報と予め用意されているスポットジャンルと関連付ける。なお、スポットジャンルは、スポットジャンルが格納されているデータベースまたは、インターネット等から取得するものとする。本実施の形態では、スポットデータベースを用いるものとする。ステップ3滞在エリアにおけるスポットジャンルの出現頻度をスコア化する。このとき、滞在エリア内に出現する希少性を考慮してユーザが訪問したと思われるスポットのスコアを算出する。ステップ4スコアに基づいて、ユーザの訪問スポット、行動、嗜好を推定する。図3は、本発明の一実施の形態におけるシステム構成を示す。本発明のシステムは、上記の図2の処理を行う推定処理装置100、利用者1が利用するGPS機能付き携帯端末等の位置情報収集デバイス200、サービス処理装置300から構成され、これらは通信ネットワークを介して接続されている。推定処理装置100は、位置情報収集デバイス200からGPSログを取得し、推定処理を行い、選好キーワードをサービス処理装置300に出力する。サービス処理装置300は、取得した選好キーワードに基づいて関連するコンテンツリコメンド等を位置情報収集デバイス200に送信する。以下に、上記のシステムの各装置構成を示す。図4は、本発明の一実施の形態における各装置の構成を示す。推定処理装置100は、GPS情報受信部110、GPS履歴前処理部120、スコア算出部130、訪問スポット推定部140、行動推定部150、属性嗜好推定部160、推定結果送信部170、スポットデータベース180、記憶部190から構成される。なお、記憶部190には各種テーブルが格納される。位置情報収集デバイス200は、GPSログ収集部210、GPS情報を推定処理装置100に送信するGPS情報送信部220、サービス処理装置300からのリコメンド結果を表示するコンテンツ提示部230から構成される。サービス処理装置300は、リコメンド結果を位置情報収集デバイス200に送信するコンテンツ送信処理部310、リコメンド対象コンテンツ決定処理部320、推定処理装置100から推定結果を取得する推定結果受信部330から構成される。次に、上記のシステムの全体の動作を説明する。図5は、本発明の一実施の形態におけるシステム全体の処理のフローチャートである。ステップ1110推定処理装置100のGPS情報受信部110は、位置情報収集デバイス200のGPS情報送信部220からGPS履歴情報を受信する。ステップ1120GPS履歴前処理部120は、受信したGPS履歴から位置情報収集デバイス200のユーザの滞在エリアを抽出し、新しい滞在エリアを抽出した場合はステップ1130に移行し、抽出できない場合は、当該処理を終了する。詳細な処理については後述する。ステップ1130スコア算出部130において、滞在エリア内の各スポットスコアを求める。詳細な処理については後述する。ステップ1140訪問スポット推定部140は、各スポットのスコア値からユーザの訪問スポットを推定する。ステップ1150行動推定部150は、ユーザの訪問ジャンル、行動を推定する。ステップ1160属性推定部160は、過去の一定期間のスコアを統合し、ユーザ属性を推定する。新しい滞在エリアで閾値を超える等の条件を満たした推定結果訪問スポット、行動、属性があるかを判定し、ない場合は当該処理を終了する。ステップ1170推定結果送信部170は、推定結果をサービス処理装置300に送信する。必要であれば、ユーザ属性情報をスコア算出部130に出力し、スコア式に反映する。なお、上記のステップ1140~1160処理は不要であれば省略することも可能である。以下に、上記の処理をステップ毎に詳細に説明する。<GPS情報受信部、GPS履歴前処理部の処理>まず、上記のステップ1110~1120について説明する。図6は、本発明の一実施の形態におけるGPS情報受信~滞在エリア抽出処理のフローチャートである。ステップ1111GPS受信部110は、位置情報収集デバイス200からGPS情報を受信し、GPS履歴前処理部120に出力する。図7は、受信したGPSデータ形式の例であり、経度、緯度、時間が含まれる。ステップ1121GPS履歴前処理部120は、記憶部180の図8に示すGPSデータテーブルに受信した新規GPSデータを保存し、過去時系列のGPSデータを読み出す。ステップ1122GPS履歴前処理部120は、GPSデータテーブルから読み出した時系列のGPS履歴データの緯度、経度、測定時刻から一定時間内及び一定範囲内で滞留しているデータ部分を抽出する。滞在エリアの抽出については、例えば、滞留時間閾値と、各点間距離閾値をパラメータとしてクラスタを抽出するDBSCAN法等のアルゴリズムを用いるものとする。抽出された滞在エリアのデータは、記憶部180の図9に示す滞在エリアデータテーブルに保存すると共に、スコア算出部130に送信する。ここで、抽出された滞在エリアのデータは、中心緯度、中心経度、滞在エリア半径、滞留開始時間、滞留終了時間を含むステップ1123GPS履歴前処理部120は、抽出した滞在エリア情報について、中心点から一定半径分析者が設定するパラメータ内に存在するスポットとそのスポットに関連付けられたジャンル例えば、レストランを関連付けたデータを、スポット識別子に基づいて図10に示すスポットデータベース180を参照することにより作成し、図11に示す滞在エリア内ジャンルテーブルに保存すると共に、スコア算出部130に出力する。<スコア算出部の処理>次に、上記のGPS履歴前処理部120から、滞在エリアのデータ及び、スポット及び当該スポットにジャンルを関連付けたデータを取得したスコア算出部130の処理について説明する。スコア算出部130は、滞在エリアjにおけるジャンルiを持つスポットのスコアSijxを求める。スコアSijxは、現時刻による行動確率;滞在エリアとの距離による訪問確率;周辺ジャンル分布比較による希少性;過去ユーザGPS履歴から抽出した滞在エリアとの比較による現滞在エリア希少性;の少なくとも1つ、または、いずれか、または、全てを乗じたものである。全てを乗じた場合の計算式は上記において、スポットxは滞在エリア中心点から一定範囲のスポットである。Xは、滞在エリアj内のスポット総数とする。wは滞在エリアj内の各スポットである。φ1xは累積標準正規分布の逆関数である。右辺の第1項Eは、滞在エリアの中心時間におけるジャンルに関連付けられる行動実施率であり、例えば、iがグルメジャンルの場合は、時間帯毎の食事行動を行う確率となる。但しグルメは一例であり、技術的には、カラオケや学校等の多様なジャンルがある。このとき、カラオケであれば遊び、学校であれば学習という行動が関連付けられることが考えられる。なお、この確率については、統計庁の社会生活時間調査等の外部情報を元に、時間帯毎の平均的な実施率のデータが公開されており、それら外部情報を利用することを想定している。また、上記の右辺の第2項の分子のDistはGPS誤差をy点を中心とし、μ0σGPS機器の想定誤差による値で決定した場合の正規分布によるx地点の出現確率xyの距離の出現確率である。第2項の分母のwは滞在エリアj内における全スポット集合である。locはスポットxの位置情報であり、n,mは重みパラメータ乗数である。右辺の第3項のMは滞在エリアjのジャンルiの出現数の珍しさを定量化したものである。詳細については後述する。右辺の第4項のRはユーザの過去の履歴内から抽出した滞在エリアから得られる、ユーザ嗜好に関する情報は本手法では減衰する考え方でスコアを下げる例えば、学校、会社、バイト等によるルーティン的な訪問。Rの定義例として、Rcoverageとする。ここでintervalは過去の滞在エリア内において現在滞在エリアとの中心点間距離が閾値以下の滞在エリアで、滞在エリアの滞在時刻が近いものがあるほど1に近づき、遠いほど0に近づく関数である。Coverageは過去滞在エリア内において現滞在エリアとの中心点距離が近いものがあるほど1に近づき、一定以上遠い場合は0になる関数である。上記の滞在エリアjのジャンルiの出現数の珍しさMの求め方について説明する。まず、第1のMを求める方法は、移動モデルと超機可分布によるスポット選択モデル用いるものである。当該方法は、一定期間の滞在エリア群から作成した移動モデルを利用して抽出した選択候補スポット集合から超機可分布を利用して滞在エリア毎のジャンル出現数の珍しさを以下のように定義する。上記において、Tjは滞在エリアjに対応する移動範囲内以降、移動候補エリアのスポットの数、Lijは滞在エリアjの移動候補エリア内のジャンルiの数、Nijは滞在エリアj内のジャンルiの数、Sjは滞在エリアj内のスポットの数である。このTj、Lijの決定方法は以下2つある。尚、Tjはユーザの選択候補スポット集合の総数、Lijはその中でのキーワードiを持つスポットの数であると見なせる。この選択候補スポット集合を算出するために移動モデルを定義する。移動モデルとは、ユーザの一定期間の滞在エリアの位置分布の特徴から、ユーザが各滞在エリア=訪問スポットの位置を決定する際、候補としたと考える移動候補エリアを決定するモデルである。この移動候補エリア内のスポットが、各滞在エリアについての選択候補スポット集合となる。Tj、Lij決定方法重心からの距離分布による移動モデル利用:第一の移動モデルでは、図12に示すように、各ユーザの一定期間の全滞在エリア中心点から重心を算出して移動候補エリアの中心点とする。この重心を中心として、一定量の滞在エリアが含まれる距離を半径とした円エリアを作成する。この円エリアを各滞在エリアの移動候補エリアとして選択候補スポット集合を抽出する。尚、この移動候補エリアに入らない滞在エリアがある場合は、外れ値データとして処理を行わない。Tj、Lij決定方法3層クラスタによる移動モデル利用:第二の移動モデルでは、図13に示すように、各ユーザの一定期間の全滞在エリア中心点に対して、クラスタリング例えば、DBSCAN法を3種類のパラメータで段階的に適用することで構造化して移動候補エリアを決定する。この構造化とは、全滞在エリアに対して車や電車などによる長距離移動クラスタと、車や電車を降りた後の徒歩等による短距離移動クラスタ、短距離移動クラスタ内にある自宅や職場等の頻繁訪問地クラスタとしてクラスタリングする。この構造化は、各移動種別を想定した閾値でDBSCAN法を階層的に適用して作成する。各クラスタに対する移動候補エリアは、自クラスタの下位クラスタの中で、所属する滞在エリア数最大の下位クラスタのコアポイントを中心に、自クラスタに所属する滞在エリアがすべて収まる範囲の円を作成することで決定する。コアポイントとは、各クラスタ内で滞在エリアが集中する地点を移動起点と考えて決定する位置情報である。これは、最小単位である頻繁訪問地クラスタから上位クラスタに向かって順に作成する。尚、最下位の頻繁訪問地クラスタでは、所属する滞在エリアの重心を算出することでコアポイントとする。また、頻繁訪問地クラスタ以外でも下位クラスタがない場合には、自クラスタに所属する滞在エリアの重心をコアポイントとする。各滞在エリアの移動候補エリアは、滞在エリアが所属するクラスタの移動候補エリアによって決定する。滞在エリアが長距離移動クラスタと短距離移動クラスタ階層に同時に所属する場合は、短距離移動クラスタの移動候補エリアを利用する。尚、頻繁訪問地クラスタについては、短距離移動クラスタのコアポイント算出のためだけに利用し移動候補エリア決定に利用しない。この処理で決定した移動候補エリアにより、滞在エリア毎に選択候補スポット集合を抽出する。また、Vは、滞在エリアjに対応する短距離移動クラスタのコアポイントと滞在エリアの距離による減少関数である。短距離移動クラスタに所属していない滞在エリアについてはV1とする。また、後述の第2のMを求める方法利用時は、Vは常に1とする。但し、dist=位置xと位置yの距離を返す関数、corej滞在エリアjの所属する短距離移動クラスタのコアポイント位置、centerj滞在エリアjの中心点、coreradiusj滞在エリアjの所属する短距離移動クラスタの移動候補エリア半径である。第2のMを求める方法は、訪問可能性の距離減衰を反映した確率分布推定による方法である。当該方法は、一定期間の移動範囲の中心点から各スポットに対する距離に対する減衰関数に従い、範囲内のスポットを選択する。以下に示すように、そのスポットの中心にジャンルiのスポット数比の存在率を計算し、分布を作成する。MZ上記の各変数は第1の方法と同様である。Yは一定期間の滞在エリア群の中心点から各スポットとの距離に応じて減衰する選択確率選択確率最大は中心点であり、選択確率0となる距離はKとするでスポットを多数回選択し、選択したスポットから一定範囲GPS機器誤差から想定するσの2倍範囲等のスポット群でのスポットとジャンルの比を算出することで推定するジャンル接触数確率分布に対して、引数の滞在エリアでのジャンル出現比の上側確率である。Zはif0、ifxである。当該一定期間移動で距離減衰するスポット選択モデルによる指標の例を図13に示す。上記の一定期間の滞在エリア群の中心点とKの求め方は、分析対象の全滞在エリア中心点から重心を算出して中心点とする。この重心を中心として、一定量の滞在エリアが含まれる距離に重みを掛けた値をKとする。スコア算出部130は、上記のようにして求められたスコアを訪問スポット推定部140、行動推定部150、属性推定部160に出力する。<訪問スポット推定部の処理>次に、図5のステップ1140の訪問スポット推定部140の処理を説明する。図14は、本発明の一実施の形態における訪問スポット推定部の処理のフローチャートを示す。ステップ1141訪問スポット推定部140は、履歴データが予め設定された閾値を超えているかを判定し、超えている場合はステップ1142に移行し、超えていない場合は当該処理を終了する。ステップ1142訪問スポット推定部140は、スコア算出部130から取得した各スポットxが保持するジャンルiに対するスコアSijxが予め設定された条件を超えているかを判定し、超えている場合はステップ1143に移行し、超えていない場合は当該処理を終了する。条件としては、閾値を超えている;各滞在エリアjのスポットの中で上位K番目以内;等を用いる。ステップ1143ステップ1142で条件を満たしたスポットを訪問可能性の高いスポットとして推定結果送信部170に出力する。出力結果は、例えば、図15のようになる。<行動推定部の処理>次に、図5のステップ1150の行動推定部150の処理について説明する。図16は、本発明の一実施の形態における行動推定部の処理のフローチャートである。ステップ1151行動推定部150は、記憶部180のジャンルスコア行動対応付けテーブルからスコアを算出する条件を決定する。当該ジャンルスコア行動対応付けテーブルには、ジャンルジャンル組み合わせもあり、期間等の合算条件、合算時の重み、推定結果としての閾値等の条件、条件を超えたときに推定結果とする行動が設定されている。ステップ1152ステップ1151で決定した条件に基づいて、スコア算出部130から取得した、現在の滞在エリア内のジャンル毎、もしくは、複数ジャンルの各スポットのスコアを合算する。ステップ1153ステップ1151で決定した条件に基づいて、ステップ11522で算出した滞在エリア毎のスコアを一定期間等で合算する。なお、上記のステップ1152、1153における合算については、過去期間に応じて、もしくはジャンル組み合わせでは、組み合わせジャンルに応じて重み付けしてもよい。ステップ1154合算したスコアがステップ1151で読み出した以下のような条件を満たしているかを判定し、満たしている場合はステップ1155に移行し、満たしていない場合は当該処理を終了する。条件閾値を越える、または、下回る。比較対象内で上位、または、下位になっている。ステップ1155ステップ1151で読み出した条件を満たす行動を合算値と共に推定結果送信部170に出力する。行動推定部150において、ユーザの行動の特徴を判定する方法として、ジャンルiが所定の閾値より高い場合は、ジャンルiについて訪問している可能性が高い滞在エリアが一定期間に複数回ある場合は、強く珍しいことが起こっているため、ジャンルi嗜好があると考え、また、スコアが高く寄与したスポットは訪問していると判断する。また、ジャンルiスコアが所定の閾値より低いとき、過去一定期間、ジャンルiに近接していない場合はジャンルiに関連付けられる行動をしていない、または、食事時間までにジャンルiも含む全食事系ジャンルに近づいていないのであれば、食事自体をしていないと判断する。<属性推定部の処理>次に、図5のステップ1160の属性推定部160の処理を説明する。図17は、本発明の一実施の形態における嗜好推定部の処理のフローチャートである。ステップ1161属性推定部160は、スコア算出部130から取得した滞在エリア毎、ジャンル毎のスコアを取得すると、滞在エリア毎かつジャンル毎にスコアを合算する。ステップ1162ステップ1161で合算されたスコアについて、所定の期間分のスコアを合算する。期間分の合算については、単純な合計値の算出以外に、各滞在エリアのジャンル毎スコアについて、確率として結合する方法がある。具体的にStouffermethodMetaanalyticproceduresforsocialresearchを利用した例を以下に示す。但し、Pi=分析対象とした一定期間の滞在エリアに対するキーワードiのスコア、φx:標準正規分布の累積確率密度関数、Xij=滞在エリアjにおけるキーワードiを持つ全スポット、J=一定期間の全滞在エリア、j=滞在エリアの識別子とする。ステップ1163属性推定部160は、合算したスコアが予め設定された以下のような条件を満たすかを判定する。条件閾値を超える、または、下回る。比較対象内で上位、または、下位になっている。条件を満たす場合はステップ1164に移行し、満たさない場合は当該処理を終了する。ステップ1164上記の条件を満たしたジャンルを推定結果送信部170に出力する。出力される属性リストの例を図18に示す。ステップ1165後述する訪問スポット推定補正を行うために、スコア部130に対して、フィードバックする。<推定結果送信部の処理>次に、図5のステップ1170の推定結果送信部170の処理を説明する。図19は、本発明の一実施の形態における推定結果送信部の処理のフローチャートである。ステップ1171推定結果送信部170は、訪問スポット推定部140、行動推定部150、属性推定部160から送信された訪問スポット、行動、属性の各推定結果が所定の条件を超えているかを判定し、超えている場合はステップ1172に移行し、超えていない場合は当該処理を終了する。ステップ1172条件を超えた項目をHTML、XML、CSV等の各形式に纏めたデータを生成する。ステップ1173ステップ1172で生成されたデータをサービス処理装置300に推定結果として出力する。上記のステップ1165において、属性推定部160から推定結果がスコア算出部130にフィードバックされた場合の処理について説明する。当該処理は、スポットxジャンルi保持訪問スコアの高精度化を図るものである。補正する場合の訪問スコアは、現時刻による行動確率、滞在エリアとの距離による訪問確率を乗算した値に、過去の移動履歴からのジャンルiスコアによる補正値、アンケート、Web閲覧履歴等のほかのデータからのジャンルi嗜好による補正値、Web履歴、スケジューラ、メール等の他履歴でのスポットx関連キーワード共起からの補正値、近い時間帯の滞在エリアの各ジャンルスコアからの補正値、嗜好の季節変動、時間帯例えば、昼食と夕食の違い等変動、曜日変動からの補正値、支払い、決済行動履歴による補正値のいずれかの補正値を乗算することで求められる。上記のスコア算出部130、訪問スポット推定部140、行動推定部150、属性推定部160の処理により、図20に示すような例において、ユーザの移動範囲内に、ユーザの移動範囲a;食事時刻以外の滞在エリアb;食事時刻の滞在だが情報が少ないと見做すエリアc;食事時刻の滞在エリアd;を考えるとき、エリアcのように、滞在エリアの中心点から離れすぎる、店舗が多すぎて訪問の特性が判別できない場所はスコアが小さくなる。また、エリアb、エリアc、エリアdが重複している場所は職場である可能性が高く、その中に2回も沖縄料理が入っているため、偶然とは考えにくく、ユーザの嗜好性として沖縄料理が好きであると判断できる。なお、上記の図4に示す推定処理装置の各構成要素の動作をプログラムとして構築し、推知処理装置として利用されるコンピュータにインストールして実行させる、または、ネットワークを介して流通させることが可能である。本発明は上記の実施の形態に限定されることなく、特許請求の範囲内において種々変更応用が可能である。"