(57)【要約】【課題】Ｗｅｂ関係のデータの予測値を正確に求める。【解決手段】時系列分析装置は、Ｗｅｂ関係のデータに特有の増減傾向を表し、パラメータが未定の関数をトレンド関数として、特定の事柄について集計されたＷｅｂ関係のデータの時系列値に基づいて該トレンド関数のパラメータをトレンドパラメータとして算出し、Ｗｅｂ関係のデータに特有の周期性を表し、パラメータが未定の関数を周期関数として、前記トレンドパラメータを適用した前記トレンド関数を除去した前記時系列値に基づいて、該周期関数のパラメータを周期パラメータとして算出し、算出された前記トレンドパラメータを適用した前記トレンド関数と、算出された前記周期パラメータを適用した前記周期関数とを加算することにより、前記特定の事柄に関するＷｅｂ関係のデータの推移を予測するための予測式を算出する。【選択図】図１
2011081491.tif 000002
【特許請求の範囲】
【請求項１】
Ｗｅｂ関係のデータに特有の増減傾向を表し、パラメータが未定の関数をトレンド関数として予め定義しておき、特定の事柄について集計されたＷｅｂ関係のデータの時系列値に基づいて該トレンド関数のパラメータをトレンドパラメータとして算出するトレンドパラメータ算出手段と、
Ｗｅｂ関係のデータに特有の周期性を表し、パラメータが未定の関数を周期関数として予め定義しておき、前記トレンドパラメータを適用した前記トレンド関数を除去した前記時系列値に基づいて、該周期関数のパラメータを周期パラメータとして算出する周期パラメータ算出手段と、
前記トレンドパラメータ算出手段により算出された前記トレンドパラメータを適用した前記トレンド関数と、前記周期パラメータ算出手段により算出された前記周期パラメータを適用した前記周期関数とを加算することにより、前記特定の事柄に関するＷｅｂ関係のデータの推移を予測するための予測式を算出する予測式算出手段と、
を有する時系列分析装置。
【請求項２】
前記周期パラメータ算出手段は、周期が異なる複数の周期関数を定義しておき、それぞれの周期関数に対応する複数の周期パラメータを算出する、請求項１に記載の時系列分析装置。
【請求項３】
前記時系列値のうち、突発的に増減した値である外れ値を検出し、該外れ値を補正する補正手段を更に有し、
前記トレンドパラメータ算出手段は、前記補正手段により外れ値が検出された場合、該補正手段により該外れ値が補正された前記時系列値から前記トレンドパラメータを算出する、請求項１又は２に記載の時系列分析装置。
【請求項４】
前記補正手段は、前記時系列値が集計された期間の各時期を最初の時期とする複数のバースト期間において、それぞれの該バースト期間内の最初の時期をmintime、該バースト期間内の最後の時期をmaxtime、正の実数をα、時系列値を集計した時期をｔ、時期ｔにおける前記時系列値をｙt、自然数をｍ、該バースト期間におけるバースト値をＢＵＲＳＴとして、
【数１】
2011081491.tif 000019
【数２】
2011081491.tif 000020
【数３】
2011081491.tif 000021
の式により各バースト期間のバースト値を算出し、該バースト値の絶対値が高い順に所定個のバースト期間を抽出し、抽出した該バースト期間内の時系列値を外れ値として補正する、請求項３に記載の時系列分析装置。
【請求項５】
前記補正手段は、前記バースト期間の長さを、最小値から最大値まで変化させて、各バースト値を算出する、請求項４に記載の時系列分析装置。
【請求項６】
前記補正手段は、バースト値の絶対値が高いバースト期間を優先して順に抽出し、既に抽出したバースト期間と時期が重複するバースト期間は抽出しない、請求項５に記載の時系列分析装置。
【請求項７】
時系列値の値を予測すべき予測時期が入力されたとき、前記予測式算出手段により算出された前記予測式に、該予測時期を代入することにより、該予測時期における時系列値を予測値として算出する予測値算出手段を更に有する、請求項１乃至６のいずれか１項に記載の時系列分析装置。
【請求項８】
イベントが開催された時期から所定の期間が経過するまでに集計された時系列値と、前記予測値算出手段により該期間内で算出された予測値とに基づいて該イベントの影響度を算出する影響度算出手段を更に有する、請求項７に記載の時系列分析装置。
【請求項９】
前記予測値算出手段は、第１の集計期間と、第２の集計期間と、予測時期とが入力されたとき、該第１の集計期間内に集計された時系列値に基づいて第１の予測式を算出し、該第２の集計期間内に集計された時系列値に基づいて第２の予測式を算出し、該第１の予測式に該予測時期を代入して得た予測値に対する、該第２の予測式に該予測時期を代入して得た予測値の度合を成長度として算出する、請求項７に記載の時系列分析装置。
【請求項１０】
前記予測式算出手段により算出された前記予測式と、該時系列値が集計された条件とを対応付けて記憶する記憶手段を更に有し、
前記予測値算出手段は、前記記憶手段に記憶された前記条件と類似する条件で集計された時系列値と予測時期とが入力されたとき、前記記憶手段から該類似する条件に対応する前記予測式を読み出し、該予測式に該予測時期を代入することにより予測値を算出する、請求項７に記載の時系列分析装置。
【請求項１１】
前記トレンド関数は、対数関数、一次関数、累乗関数、及び指数関数のうち、いずれかの関数である、請求項１乃至１０のいずれか１項に記載の時系列分析装置。
【請求項１２】
前記トレンドパラメータ算出手段は、最小ニ乗法を使用して前記時系列値と、該時系列値が集計された時期をトレンド関数に代入して得た値との差異が最小となる前記トレンドパラメータを算出する、請求項１乃至１１のいずれか１項に記載の時系列分析装置。
【請求項１３】
前記周期パラメータ算出手段は、前記減算後の前記時系列値をフーリエ解析することにより、前記周期パラメータを算出する、請求項１乃至１２のいずれか１項に記載の時系列分析装置。
【請求項１４】
前記周期パラメータ算出手段は、前記減算後の前記時系列値をＡＲモデル解析することにより、前記周期パラメータを算出する、請求項１乃至１３のいずれか１項に記載の時系列分析装置。
【請求項１５】
Ｗｅｂ関係のデータに特有の増減傾向を表し、パラメータが未定の関数をトレンド関数として予め保持したトレンドパラメータ算出手段が、特定の事柄について集計されたＷｅｂ関係のデータの時系列値に基づいて該トレンド関数のパラメータをトレンドパラメータとして算出し、
Ｗｅｂ関係のデータに特有の周期性を表し、パラメータが未定の関数を周期関数として予め保持した周期パラメータ算出手段が、前記トレンドパラメータを適用した前記トレンド関数を除去した前記時系列値に基づいて、該周期関数のパラメータを周期パラメータとして算出し、
予測式算出手段が、前記トレンドパラメータ算出手段により算出された前記トレンドパラメータを適用した前記トレンド関数と、前記周期パラメータ算出手段により算出された前記周期パラメータを適用した前記周期関数とを加算することにより、前記特定の事柄に関するＷｅｂ関係のデータの推移を予測するための予測式を算出する、時系列分析方法。
【請求項１６】
コンピュータに、
Ｗｅｂ関係のデータに特有の増減傾向を表し、パラメータが未定の関数をトレンド関数として予め定義しておき、特定の事柄について集計されたＷｅｂ関係のデータの時系列値に基づいて該トレンド関数のパラメータをトレンドパラメータとして算出するトレンドパラメータ算出手順、
Ｗｅｂ関係のデータに特有の周期性を表し、パラメータが未定の関数を周期関数として予め定義しておき、前記トレンドパラメータを適用した前記トレンド関数を除去した前記時系列値に基づいて、該周期関数のパラメータを周期パラメータとして算出する周期パラメータ算出手順、及び
前記トレンドパラメータ算出手順で算出された前記トレンドパラメータを適用した前記トレンド関数と、前記周期パラメータ算出手順で算出された前記周期パラメータを適用した前記周期関数とを加算することにより、前記特定の事柄に関するＷｅｂ関係のデータの推移を予測するための予測式を算出する予測式算出手順、
を実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、時系列データの推移を推測する技術に関する。
【背景技術】
【０００２】
非特許文献１に示すように、ブログや掲示板上のクチコミ情報をリアルタイムに収集、分析するＡＳＰ型サービスが各社から提供されている。但し、これらのサービスでは、将来のブログの書き込み数がどのように推移していくかは予測されない。
【０００３】
ここで、経済や交通など、他の分野であれば、株価や交通情報の推移を予測する技術が用いられている。
【０００４】
例えば、特許文献１に記載された学習装置は、ブログの書き込み数を予測するものでないが、交通情報に関する時系列データを取得し、ＡＲ（Auto Regression）モデルのクラスを使用して、周期的なパターンを付加した、将来の予測値を算出することができる。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】[online]、[平成２１年９月２５日検索]、インターネット<URL：http://kandoreport.jp/>
【特許文献】
【０００６】
【特許文献１】特開２００６−１１５７２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかし、特許文献１に記載された学習装置では、ブログの書き込み数など、Ｗｅｂ関係のデータの推移を正確に予測できない場合があった。
【０００８】
Ｗｅｂ関係のデータの変動には、周期性のほか、特有の傾向も存在する。これに対し、特許文献１の学習装置は、周期的な変動を考慮して、予測値を算出するが、データが上昇傾向、下降傾向にあるといった特有の傾向については考慮していない。
【０００９】
このため、特許文献１に記載された技術をＷｅｂ関係のデータに適用すると、予測値と実際の値との間のずれが大きくなってしまうという問題があった。
【００１０】
本発明は、Ｗｅｂ関係のデータの予測値を正確に求めることを目的とする。
【課題を解決するための手段】
【００１１】
上記目的を達成するために、本発明の時系列分析装置は、Ｗｅｂ関係のデータに特有の増減傾向を表し、パラメータが未定の関数をトレンド関数として予め定義しておき、特定の事柄について集計されたＷｅｂ関係のデータの時系列値に基づいて該トレンド関数のパラメータをトレンドパラメータとして算出するトレンドパラメータ算出手段と、Ｗｅｂ関係のデータに特有の周期性を表し、パラメータが未定の関数を周期関数として予め定義しておき、前記トレンドパラメータを適用した前記トレンド関数を除去した前記時系列値に基づいて、該周期関数のパラメータを周期パラメータとして算出する周期パラメータ算出手段と、前記トレンドパラメータ算出手段により算出された前記トレンドパラメータを適用した前記トレンド関数と、前記周期パラメータ算出手段により算出された前記周期パラメータを適用した前記周期関数とを加算することにより、前記特定の事柄に関するＷｅｂ関係のデータの推移を予測するための予測式を算出する予測式算出手段と、を有する。
【００１２】
本発明の時系列分析方法は、Ｗｅｂ関係のデータに特有の増減傾向を表し、パラメータが未定の関数をトレンド関数として予め保持したトレンドパラメータ算出手段が、特定の事柄について集計されたＷｅｂ関係のデータの時系列値に基づいて該トレンド関数のパラメータをトレンドパラメータとして算出し、Ｗｅｂ関係のデータに特有の周期性を表し、パラメータが未定の関数を周期関数として予め保持した周期パラメータ算出手段が、前記トレンドパラメータを適用した前記トレンド関数を除去した前記時系列値に基づいて、該周期関数のパラメータを周期パラメータとして算出し、予測式算出手段が、前記トレンドパラメータ算出手段により算出された前記トレンドパラメータを適用した前記トレンド関数と、前記周期パラメータ算出手段により算出された前記周期パラメータを適用した前記周期関数とを加算することにより、前記特定の事柄に関するＷｅｂ関係のデータの推移を予測するための予測式を算出する、時系列分析方法である。
【００１３】
本発明のプログラムは、コンピュータに、Ｗｅｂ関係のデータに特有の増減傾向を表し、パラメータが未定の関数をトレンド関数として予め定義しておき、特定の事柄について集計されたＷｅｂ関係のデータの時系列値に基づいて該トレンド関数のパラメータをトレンドパラメータとして算出するトレンドパラメータ算出手順、Ｗｅｂ関係のデータに特有の周期性を表し、パラメータが未定の関数を周期関数として予め定義しておき、前記トレンドパラメータを適用した前記トレンド関数を除去した前記時系列値に基づいて、該周期関数のパラメータを周期パラメータとして算出する周期パラメータ算出手順、及び前記トレンドパラメータ算出手順で算出された前記トレンドパラメータを適用した前記トレンド関数と、前記周期パラメータ算出手順で算出された前記周期パラメータを適用した前記周期関数とを加算することにより、前記特定の事柄に関するＷｅｂ関係のデータの推移を予測するための予測式を算出する予測式算出手順、を実行させるためのプログラムである。
【発明の効果】
【００１４】
本発明によれば、時系列分析装置は、Ｗｅｂ関係のデータに特有の増減傾向や周期性を表す関数のパラメータを算出するので、周期性のほか、増減傾向も存在するブログなどに関する時系列値について、予測値を正確に算出できる。
【図面の簡単な説明】
【００１５】
【図１】本発明の第１の実施形態の時系列データ分析システムの一構成例を示す全体図である。
【図２】本発明の第１の実施形態のブログ抽出・集計部の一構成例を示すブロック図である。
【図３】本発明の第１の実施形態の時系列分析要求データに含まれる情報を示す図である。
【図４】本発明の第１の実施形態の学習データの示す内容をまとめた表である。
【図５】本発明の第１の実施形態のバースト期間を可変して算出されたバースト値を示す表である。
【図６】（ａ）本発明の第１の実施形態の整列したバースト値を示す図である。（ｂ）本発明の第１の実施形態の整列したバースト値を示す図である。
【図７】本発明の第１の実施形態の予測式データの示す内容をまとめた表である。
【図８】本発明の第１の実施形態の分析・可視化サーバの時系列分析処理を示すフローチャートである。
【図９】本発明の第１の実施形態のブログ抽出・集計処理を示すフローチャートである。
【図１０】本発明の第１の実施形態の周期成分抽出処理を示すフローチャートである。
【図１１】外れ値の補正前における、本発明の第１の実施形態の時系列値をプロットした折れ線グラフである。
【図１２】外れ値の補正後における、本発明の第１の実施形態の時系列値をプロットした折れ線グラフである。
【図１３】トレンド関数算出後における、本発明の第１の実施形態の時系列値をプロットした折れ線グラフである。
【図１４】トレンド成分除去後における、本発明の第１の実施形態の時系列値をプロットした折れ線グラフである。
【図１５】本発明の第１の実施形態の予測値をプロットした折れ線グラフである。
【図１６】本発明の第１の実施形態の予測値をプロットした折れ線グラフである。
【図１７】本発明の第２の実施形態の影響度分析要求データに含まれる情報を示す図である。
【図１８】本発明の第２の実施形態の予測式データの示す内容をまとめた表である。
【図１９】本発明の第２の実施形態の分析・可視化サーバの影響度分析処理を示すフローチャートである。
【図２０】本発明の第２の実施形態の予測値をプロットした折れ線グラフである。
【図２１】本発明の第３の実施形態の成長度スコア算出要求データに含まれる情報を示す図である。
【図２２】本発明の第３の実施形態の分析・可視化サーバの成長度算出処理を示すフローチャートである。
【図２３】本発明の第３の実施形態の予測値をプロットした折れ線グラフである。
【図２４】本発明の第４の実施形態のシミュレーション要求データに含まれる情報を示す図である。
【図２５】本発明の第４の実施形態のイベント効果予測処理を示すフローチャートである。
【発明を実施するための形態】
【００１６】
本発明を実施するための形態について図面を参照して詳細に説明する。図１は、本実施形態の時系列データ分析システム１の一構成例を示す全体図である。時系列データ分析システムは、時系列で集計された値（以下、「時系列値」という）を示す時系列データを分析する。本実施形態では、時系列で集計されたブログの書き込み数を時系列値とする。
【００１７】
図１を参照すると、時系列データ分析システム１は、利用者ＰＣ１０と、分析・可視化サーバ２０と、ブログＤＢ３０と、予測パラメータＤＢ４０とを有する。
【００１８】
分析・可視化サーバ２０は、ブログ抽出・集計部２０１、トレンド成分学習部２０２、周期成分学習部２０３、および分析結果作成部２０４を有する。
【００１９】
ブログＤＢ３０は、ブログデータを収集し、それらのデータをカテゴリに分類して格納する。
【００２０】
利用者ＰＣ１０は、時系列値の分析を要求する時系列分析要求データを受け付ける。時系列分析要求データが入力されると、利用者ＰＣ１０は、その時系列分析要求データを分析・可視化サーバ２０へ送信する。
【００２１】
図２は、時系列分析要求データに含まれる情報を示す図である。同図を参照すると、時系列分析要求データは、「カテゴリ」、「キーワード」、「教師データ期間」、および「予測期間」を示す情報を含む。
【００２２】
「カテゴリ」は、集計対象のブログの属性である。例えば、飲料、電化製品、携帯電話、ビール、有名人など、キーワードのジャンルが指定される。「キーワード」は、ブログにおいて検索される語句である。「教師データ期間」は、ブログが集計される期間である。「予測期間」は、ブログの書き込み数の予測値が算出される期間である。
【００２３】
図３は、ブログ抽出・集計部２０１の一構成例を示すブロック図である。同図を参照すると、ブログ抽出・集計部２０１は、集計部２０１１、および特異点補正部２０１２を有する。
【００２４】
集計部２０１１は、図２に示した時系列分析要求データから、カテゴリ、キーワード、教師データ期間、および予測範囲を取得する。集計部２０１１は、教師データ期間における、そのカテゴリに属するブログデータをブログＤＢ３０に格納されたブログデータから抽出する。集計部２０１１は、それらのブログデータから、キーワードを含むブログの書き込みを検索し、教師データ期間内の各日付において、検索されたブログの書き込み数を集計する。
【００２５】
特異点補正部２０１２は、集計結果から外れ値を検出し、補正する。外れ値は、不定期なイベントによる書き込み数の急増やシステム障害による書き込み数の急減・欠損など、時系列値が突発的に増減したときの値である。外れ値は、特異点とも呼ばれる。
【００２６】
外れ値の抽出において、特異点補正部２０１２は、長さが最小値から最大値までの範囲で可変のバースト期間について、その長さを最小値に設定する。本実施形態では、最小値は１日に設定され、利用者は最大値のみを設定する。
【００２７】
そして、特異点補正部２０１２は、バースト期間の初日を教師データ期間の初日とし、下記の式によりバースト値を算出する。以降、特異点補正部２０１２は、バースト期間の各日付をバースト期間の初日として、各日付に対応するバースト値を算出する。同様の手順で、特異点補正部２０１２は、最大値に達するまで、バースト期間の長さを可変し、各日付に対応するバースト値を算出する。
【００２８】
【数１】
2011081491.tif 000003
【００２９】
【数２】
2011081491.tif 000004
【００３０】
【数３】
2011081491.tif 000005
上記（式１）において、「ＢＵＲＳＴ」はバースト値であり、「α」は、長期間選択調整係数（例えば、２．５）である。「mintime」は、バースト期間内の最初の日付である。「maxtime」は、バースト期間内の最後の日付である。
【００３１】
上記（式３）において、「ｙi」は、「mintime」から数えて「i」番目の日付において集計されたブログの書き込み数である。「ｍ」は、バースト期間の最大値以上で、教師データ期間の長さ以下の値（例えば、７）である。つまり、「ＳＡＢＵＮi」は、「i」番目の日付において集計されたブログの書き込み数「ｙi」と、「i」番目の日付の「ｍ」日前から「i」番目の日付までの間におけるブログの書き込み数の平均値と、の間の差分である。
【００３２】
バースト期間の長さ「Ｗ」と、「ｍ」の値とは、外れ値が適切に検出されるように調整される。例えば、外れ値を１日単位で検出したい場合は、「Ｗ」は０に設定される。また、外れ値が２日以上連続で出現することが多い場合、利用者は「Ｗ」を１以上に設定すればよい。「ｍ」は、「Ｗ」の大きさに応じて適切な値が設定される。
【００３３】
あるバースト期間内に外れ値がある場合、そのバースト期間のバースト値の絶対値は、バースト期間内に外れ値のない場合よりも大きくなる。このため、特異点補正部２０１２は、バースト期間ごとに、バースト値を算出することで、どのバースト期間内に外れ値があるのかを効率的に判断できる。
【００３４】
特異点補正部２０１２は、バースト値が高い順に各バースト期間をソートする。特異点補正部２０１２は、バースト値が高い所定個のバースト期間とバースト値の低い所定個のバースト期間とを、特異点を含むバースト期間として抽出する。例えば、上位、下位１０個のバースト期間が抽出される。
【００３５】
但し、バースト期間の抽出においては、特異点補正部２０１２は、バースト値の絶対値の高い期間を優先してバースト期間を順に抽出する。そして、既に抽出したバースト期間と日付が重複するバースト期間は、抽出対象から除外する。
【００３６】
特異点補正部２０１２は、抽出したバースト期間内の時系列値を外れ値として、それらの外れ値を補正する。例えば、特異点補正部２０１２は、直近データの平均値で補正する。但し、特異点補正部２０１２は、外れ値が最後のデータで、直後のデータがない場合、その外れ値を除外し、以下に説明する時系列分析において、その値を用いない。特異点補正部２０１２は、補正後の集計結果を示す学習データを出力する。
【００３７】
図１に戻り、トレンド成分学習部２０２、周期成分学習部２０３、および分析結果作成部２０４について説明する。
【００３８】
トレンド成分学習部２０２は、特異点補正部２０１２が出力した学習データから、トレンドパラメータを算出する。
【００３９】
トレンドパラメータとは、教師データ期間におけるブログ書き込み数の増加傾向、または減少傾向を示す関数（以下、「トレンド関数」という）のパラメータである。例えば、トレンドパラメータとして、トレンド関数の係数や指数が算出される。
【００４０】
トレンド関数は、例えば、下記の式で表わされる関数である。トレンド関数は、ブログ書き込み数の予測において、増減傾向のみが抽出された成分であるから、トレンド成分とも呼ばれる。
【００４１】
【数４】
2011081491.tif 000006
上記（式４）において、関数「ｙ（ｘ）」は、トレンド関数である。この「ｙ（ｘ）」の値は、教師データ期間における最初の日付から「ｘ」日目におけるブログ書き込み数となる。「α」、「β」は、トレンドパラメータである。
【００４２】
トレンド成分学習部２０２は、例えば、最小ニ乗法を使用することにより、測定データを上記（式４）に近似し、トレンドパラメータを算出する。
【００４３】
具体的には、下記の式に、測定値を代入することにより、トレンドパラメータを算出する。
【００４４】
【数５】
2011081491.tif 000007
【００４５】
【数６】
2011081491.tif 000008
上記（式５）、（式６）において、「ｔ」は教師データ期間内の最初の日からの日数であり、「Σ」は、教師データ期間における総和をとることを示している。「α」
、「β」は、トレンドパラメータである。「ｙ」は、集計されたブログの書き込み数である。
【００４６】
トレンド成分学習部２０２は、算出したトレンドパラメータを適用したトレンド関数を導出し、学習データの各日付の書き込み数から、トレンド関数にその日付に対応する日数（「ｘ」）を代入することで得られた書き込み数（「ｙ（ｘ）」）を減算する。
【００４７】
周期成分学習部２０３は、トレンド関数の算出値を減算した後の学習データをフーリエ解析することにより、周期パラメータを算出する。
【００４８】
周期パラメータとは、教師データ期間におけるブログ書き込み数、すなわち時系列値の周期的な変動を示す関数（以下、「周期関数」という）のパラメータである。例えば、周期パラメータとして、周期関数の係数や指数が算出される。
【００４９】
周期関数は、例えば、下記の式で表わされる関数である。周期関数は、ブログ書き込み数の予測において、周期性のみが抽出された成分であるから、周期成分とも呼ばれる。
【００５０】
【数７】
2011081491.tif 000009
【００５１】
【数８】
2011081491.tif 000010
【００５２】
【数９】
2011081491.tif 000011
【００５３】
【数１０】
2011081491.tif 000012
【００５４】
【数１１】
2011081491.tif 000013
【００５５】
【数１２】
2011081491.tif 000014
【００５６】
【数１３】
2011081491.tif 000015
上記（式７）〜（式１３）において、「Ｙ（ｔ）」は周期関数である。この「Ｙ（ｔ）」の値は、教師データ期間における最初の日付から「ｔ」日目におけるブログ書き込み数となる。「An」、および「Bn」は、１年を周期とする周期パラメータである。「an」、および「bn」は、１週間を周期とする周期パラメータである。「ｎ」は、フーリエ係数である。例えば、周期成分学習部２０３は、１〜５の「ｎ」について、周期パラメータを算出し、それぞれを加算する。
【００５７】
分析結果作成部２０４は、トレンド関数と周期関数とを加算して、ブログ書き込み数の推移を予測するための予測式を生成する。
【００５８】
また、分析結果作成部２０４は、その予測式の各パラメータ、すなわちトレンドパラメータおよび周期パラメータを示す予測式データを予測パラメータＤＢ４０に格納する。
【００５９】
そして、分析結果作成部２０４は、時系列分析要求データの示す予測期間内の各日付のブログ書き込み数を、予測式を使用して算出する。分析結果作成部２０４は、予測期間内で算出した値を予測値として、各予測値を示す分析結果データを利用者ＰＣ１０へ送信する。
【００６０】
また、分析結果作成部２０４は、時系列分析要求データの示すカテゴリ、キーワードが、以前にトレンド成分、周期成分を算出したカテゴリ、キーワードと同じであれば、以前に求めた予測式を用いて予測値を算出する。
【００６１】
図４は、ブログ抽出・集計部２０１が作成する学習データの示す内容をまとめた表である。同図を参照すると、学習データは、「カテゴリ」、および「キーワード」に対応付けて、複数の「日付」および「書き込み件数」を示すデータを含む。「カテゴリ」は、抽出対象のブログデータが属するカテゴリである。「キーワード」は、ブログデータにおいて検索されたキーワードである。「日付」は、集計期間内の各日付である。「書き込み件数」は、キーワードを含むブログが書き込まれた件数である。
【００６２】
例えば、時系列分析要求データの示すカテゴリが「食料品」で、キーワードとして「値上げ」が指定された場合について考える。ブログ抽出・集計部２０１は、「食料品」のカテゴリに属するブログデータについて、指定されたキーワードを検索する。ブログ抽出・集計部２０１は、キーワード「値上げ」を含むブログデータを抽出し、教師データ期間内、各日付のデータ数をそれぞれ集計する。
【００６３】
図５および図６を参照して、学習データから特異点を検出し、補正する方法について説明する。ここでは、バースト期間の長さを１日から７日に可変するものとする。また、上記（式３）における「ｍ」を７に設定する。
【００６４】
図５に示すように、ブログ抽出・集計部２０１は、バースト期間の長さ（バースト範囲）を１から７に可変し、教師データ期間内の各日付をバースト期間の初日としたバースト値を算出する。
【００６５】
図５を参照すると、例えば、バースト範囲が１日で、「mintime」が２００９年８月２５日の場合のバースト値は「０」である。バースト範囲が２日で、「mintime」が２００９年８月２５日の場合、「maxtime」は２００９年８月２６日となり、バースト値は「２２．５」である。
【００６６】
ブログ抽出・集計部２０１は、図６（ａ）のように、各期間において算出したバースト値について降順でソートする。ブログ抽出・集計部２０１は、バースト値の上位から順にバースト期間を抽出する。
【００６７】
ここで、バースト値が同一の値となる期間が複数あった場合、バースト範囲の長い方が優先して抽出される。なお、バースト値の同じ期間が複数あっても、ブログ抽出・集計部２０１がいずれも抽出する構成としてもよい。
【００６８】
図６（ａ）を参照すると、最もバースト値が高いのは、バースト期間の初日を２００９年８月２８日、バースト範囲を１日として算出されたバースト値である。２番目にバースト値が高いのは、バースト期間の初日を２００９年９月１４日、バースト範囲を２日として算出されたバースト値である。３番目にバースト値が高いのは、２００９年９月１４日に対応するバースト値である。まず、これらの３つの期間が、外れ値を含みうるバースト期間として抽出される。
【００６９】
そして、４番目にバースト値が高いのは、バースト期間の初日を２００９年８月２７日、バースト範囲を２日として算出されたバースト値であった。このバースト期間は、既に抽出されたバースト期間内の日付である８月２８日を含む。このため、バースト期間の初日を２００９年８月２７日、バースト範囲を２日とするバースト期間は、抽出対象から除外される。図６（ａ）における斜線部分は、抽出対象から除外されたバースト期間である。これらのバースト期間が除外された結果を図６（ｂ）に示す。
【００７０】
図５における、太い実線で囲まれた数字は、抽出されたバースト期間に対応するバースト値である。斜線で囲まれた数字は、この抽出済みのバースト期間と日付が重複するバースト期間に対応するバースト値である。これらの斜線で囲まれた部分に対応するバースト期間は、抽出対象から除外される。
【００７１】
ブログ抽出・集計部２０１は、図５、図６（ａ）の斜線部分に対応するバースト期間を除いた後、上位１０個のバースト期間を抽出する。次に、最もバースト値が低いバースト値から順に、同様の方法で下位１０個のバースト期間が抽出される。
【００７２】
なお、上位、下位に関わらず、最も絶対値の大きなバースト値に対応するバースト期間から順に抽出される構成としてもよい。
【００７３】
このように、日付が重複するバースト期間を抽出しないことにより、ブログ抽出・集計部２０１は、外れ値を精度よく検出できる。
【００７４】
そして、ブログ抽出・集計部２０１は、抽出したバースト期間内の値を、その前後の平均値に補正する。
【００７５】
図７は、予測式データの示す内容をまとめた図である。同図に示すように、予測式データは、「カテゴリ」、「キーワード」、「教師データ期間」、「トレンドパラメータ」、および「周期パラメータ」を示す情報を含む。
【００７６】
「カテゴリ」は、抽出対象のブログデータが属するカテゴリである。「キーワード」は、ブログデータにおいて検索されたキーワードである。「教師データ期間」は、予測式を算出するために時系列値が集計された期間である。「トレンドパラメータ」は、最小二乗法により算出されたトレンド関数のパラメータである。「周期パラメータ」は、フーリエ解析により算出された周期関数のパラメータである。
【００７７】
続いて、分析・可視化サーバ２０の動作について説明する。図８は、分析・可視化サーバ２０の動作を示すフローチャートである。この動作は、分析・可視化サーバ２０が時系列分析要求データを受信したときに開始する。
【００７８】
ブログ抽出・集計部２０１は、時系列分析要求データに基づいてブログデータを抽出し、集計するブログ抽出・集計処理を実行する（ステップＳ１）。
【００７９】
トレンド成分学習部２０２は、集計結果を上記（４）式に近似することにより、トレンド成分を抽出する（ステップＳ２）。周期成分学習部２０３は、周期成分を抽出するための周期成分抽出処理を実行する（ステップＳ３）。
【００８０】
分析結果作成部２０４は、トレンドパラメータ、周期パラメータを予測パラメータＤＢ４０に格納し、トレンド関数および周期関数を加算して予測式を作成する（ステップＳ４）。分析結果作成部２０４は、予測式から、指定された予測期間における時系列値の予測値を算出し、分析結果データを利用者ＰＣ１０へ出力する（ステップＳ５）。ステップＳ５の後、分析・可視化サーバ２０は動作を終了する。
【００８１】
図９は、ブログ抽出・集計処理を示すフローチャートである。同図を参照すると、ブログ抽出・集計部２０１は、教師データ期間内において、時系列分析要求データで指定されたカテゴリ、キーワードと全部一致または部分一致する単語を含むブログを抽出し、集計する（ステップＳ１１）。
【００８２】
ブログ抽出・集計部２０１は、教師データ期間をバースト期間ごとに分割し、各バースト期間におけるバースト値を算出する。そして、ブログ抽出・集計部２０１は、バースト値が高い所定個のバースト期間を抽出することにより特異点を検出し、補正する（ステップＳ１２）。ステップＳ１２の後、ブログ抽出・集計部２０１は、ブログ抽出・集計処理を終了する。
【００８３】
図１０は、周期成分抽出処理を示すフローチャートである。同図を参照すると、周期成分学習部２０３は、上記（式８）、（式９）を使用して年単位の周期パラメータAn、Bnを算出する（ステップＳ３１）。そして、周期成分学習部２０３は、上記（式１０）、（式１１）を使用して年単位の周期パラメータan、bnを算出する（ステップＳ３２）。ステップＳ３２の後、周期成分学習部２０３は、周期成分抽出処理を終了する。
【００８４】
図１１〜図１５を参照して、分析・可視化サーバ２０の動作結果の一例について説明する。図１１は、分析・可視化サーバ２０が、時系列分析要求データに基づいて、教師データ期間でブログの書き込み数を集計した結果をプロットした折れ線グラフである。同図における、横軸が、教師データ期間内の各集計日であり、縦軸が、集計日に集計されたブログの書き込み数である。
【００８５】
分析・可視化サーバ２０は、教師データ期間を、複数のバースト期間に分割し、上記（式１）を使用して、各バースト期間におけるバースト値を算出する。同図に示すように、Ｂ１において書き込み数が突発的に増加し、Ｂ２において書き込み数が突発的に減少しているので、これらのＢ１，Ｂ２を含むバースト期間において、バースト値の絶対値が高くなる。分析・可視化サーバ２０は、これらバースト期間内の外れ値を補正する。
【００８６】
図１２は特異点を補正した後の集計結果を示す折れ線グラフである。同図における、横軸が、教師データ期間内の各集計日であり、縦軸が、集計日に集計されたブログの書き込み数である。同図に示すように、特異点Ｂ１，Ｂ２は、いずれも補正されている。
【００８７】
図１３は、分析・可視化サーバ２０が、集計結果からトレンド成分を求めた結果を示す図である。同図における、横軸が、教師データ期間内の各集計日であり、縦軸が、集計日に集計されたブログの書き込み数である。集計結果から、最小二乗法により、同図の一点鎖線で表わされる対数関数が導出される。同図に示すように、ブログの書き込み件数の推移には、時間の経過とともに徐々に増加してゆく傾向、すなわちトレンドがある。
【００８８】
図１４は、分析・可視化サーバ２０が集計結果からトレンド成分を除去し、周期成分を求めた結果を示す図である。同図における、横軸が、教師データ期間内の各集計日であり、縦軸が、集計日に集計されたブログの書き込み数である。集計結果から、フーリエ解析より、同図のニ点鎖線で表わされる周期関数が導出される。同図に示すように、ブログの書き込み件数は、年単位で周期的に変化し、週単位でも周期的に変化する。
【００８９】
図１５は、分析・可視化サーバ２０が、予測式を使用して予測値を算出した結果を示す折れ線グラフである。同図における、横軸が、教師データ期間内の各集計日であり、縦軸が、集計日に集計されたブログの書き込み数である。同図の実線の部分は、特異点を補正した集計結果であり、同図の二点鎖線の部分は、予測式で算出した予測期間内の予測値をプロットした結果である。同図に示すように、分析・可視化サーバ２０は、予測式を使用することにより、今後のブログ書き込み数の推移を予測することができる。
【００９０】
分析・可視化サーバ２０は、異なる条件で集計された集計結果や、異なる条件で算出された分析結果を、比較のために同時に表示することもできる。
【００９１】
図１６は、同じカテゴリにおいて、類似する製品「製品Ａ」、「製品Ｂ」をキーワードとして集計された集計結果と算出された分析結果とを示す折れ線グラフである。同図において、実線の部分が、これらのキーワードについての集計結果である。破線の部分が、これらのキーワードについての予測値をプロットしたものである。
【００９２】
利用者は、各製品に関する書き込み数がどのように推移するかを比較し、どのように商品の販売、宣伝などを実行するかを検討することができる。
【００９３】
なお、本実施形態では、分析・可視化サーバ２０は、ブログの書き込み数を時系列値として収集しているが、時系列で収集された値であれば、他の時系列値を収集し、分析することもできる。例えば、時系列値は、ブログ以外のＷｅｂ関係のデータであってもよい。Ｗｅｂ関係のデータは、Ｗｅｂページの記事数や、記事が検索された回数、Ｗｅｂサイトへのアクセスログなどである。
【００９４】
また、時系列値は、株価、為替とすることも可能である。若しくは、分析・可視化サーバ２０は、ＣＰＵ使用率などのシステムログデータなどを時系列値として分析することもできる。
【００９５】
本実施形態では、分析・可視化サーバ２０は、１日ごとにブログ書き込み数を集計しているが、集計する時期は、教師データ期間内の時期であれば、１日に限らない。例えば、分析・可視化サーバ２０は、１時間ごとにブログ書き込み数を集計してもよい。
【００９６】
本実施形態では、分析・可視化サーバ２０は、上記（式１）〜（式３）を使用して外れ値を検出しているが、突発的に増減する値を検出できるのであれば、これらの数式以外の式を使用して外れ値を検出してもよい。
【００９７】
例えば、上記（式３）において、平均値と「ｙi」との間の差分を「ＳＡＢＵＮi」としているが、差分の絶対値を「ＳＡＢＵＮi」として算出してもよい。この場合、分析・可視化サーバ２０は、バースト値の上位、所定個のバースト期間を抽出する。
【００９８】
本実施形態では、分析・可視化サーバ２０は、上記（式４）のように、集計結果を対数関数に近似しているが、対数関数以外の関数に近似してもよいのは勿論である。例えば、分析・可視化サーバ２０は、下記に示すような式を用いて線形近似を行ってもよい。
【００９９】
【数１４】
2011081491.tif 000016
または、分析・可視化サーバ２０は、下記に示すような式を用いて累乗近似を行ってもよい。
【０１００】
【数１５】
2011081491.tif 000017
もしくは、分析・可視化サーバ２０は、下記に示すような式を用いて指数近似を行ってもよい。
【０１０１】
【数１６】
2011081491.tif 000018
上記（式１４）〜（式１６）において、「α」、「β」はトレンドパラメータ、「ｙ」はブログ書き込み数、「ｘ」は集計日である。
【０１０２】
本実施形態では、分析・可視化サーバ２０は、１年を周期とする周期成分と、１週間を周期とする周期成分とを求めているが、これらと異なる周期の周期成分を求めてもよいのは勿論である。例えば、ブログの書き込み数が１時間ごとに集計される場合、分析・可視化サーバ２０は、周期を１日とする周期成分を算出することもできるし、周期を１年でなく、２年や３年としてもよい。
【０１０３】
本実施形態の分析・可視化サーバ２０は、トレンド成分の算出において、最小ニ乗法を使用しているが、集計結果を回帰分析できるのであれば、最小二乗法以外の手法により、トレンド成分を算出してもよい。
【０１０４】
本実施形態の分析・可視化サーバ２０は、フーリエ解析を行うことにより、周期成分を算出しているが、周期的な変動パターンを求めることができるのであれば、フーリエ解析以外の手法により周期成分を算出してよい。例えば、分析・可視化サーバ２０は、ＡＲ（Autoregressive model）モデルを使用することにより、周期成分を求めることが可能である。
【０１０５】
以上説明したように、本実施形態によれば、時系列分析装置は、Ｗｅｂ関係のデータに特有の増減傾向や周期性を表す関数のパラメータを算出するので、時系列値に、周期性のほか、増減傾向があっても、予測値を正確に算出できる。
【０１０６】
（第２の実施形態）
本発明の第２の実施形態について説明する。本実施形態の時系列分析システムは、イベントの影響度を更に算出できる点で、第１の実施形態の時系列分析システムと異なる。なお、本実施形態の時系列分析システムについて、図１〜図１６で説明した構成と同様な構成の詳細な説明を省略し、第１の実施形態と異なる点について詳しく説明する。
【０１０７】
本実施形態の利用者ＰＣ１０は、時系列分析要求データの他、イベントの影響度の算出を要求する影響度分析要求データの入力を受け付ける。イベントの影響度とは、イベントを実行後の、実際のブログの書き込み数が予測値に対して増減した度合を示す値である。例えば、分析・可視化サーバ２０は、実際の書き込み数を予測値で割った値を影響度とする。
【０１０８】
具体的には、イベントの開始から一定期間が経過した時点におけるイベントの影響度を求めたいとき、利用者は、イベントの種類、イベントの開始時期、予測期間、および教師データ期間を利用者ＰＣ１０に入力する。ここで、イベントの開始時期以前の所定期間が教師データ期間として設定される。また、イベントの開始時期以降の所定期間が予測期間として設定される。本実施形態では、教師データ期間の長さは、予測期間の長さ以上に設定されることが望ましい。
【０１０９】
例えば、イベントの開始から１カ月感が経過し、その１カ月間におけるイベントの影響度を算出したい場合、利用者は、イベント開始以前の１年間を教師データ期間、イベント開始時から１カ月間の期間を予測期間として入力する。
【０１１０】
利用者ＰＣ１０は、イベントの種類、イベントの開始時期、予測期間、および教師データ期間を示す情報を含む影響度分析要求データを分析・可視化サーバ２０へ出力する。
【０１１１】
図１７は、影響度分析要求データに含まれる情報を示す図である。同図を参照すると、影響度分析要求データには、「カテゴリ」、「キーワード」、「イベント種類」、「教師データ期間」、および「予測期間」を示す情報を含む。
【０１１２】
「カテゴリ」、「キーワード」は、第１の実施形態と同様である。「教師データ期間」は、イベント開始時期以前の所定期間である。「予測期間」は、イベント開始時期以降の所定期間である。
【０１１３】
影響度分析要求データが入力されたとき、分析・可視化サーバ２０は、第１の実施形態と同様の方法で予測式を求め、その予測式から、予測期間内の各日付の予測値を算出する。
【０１１４】
一方、分析・可視化サーバ２０は、予測期間内におけるブログデータをブログＤＢ３０から読み出し、予測期間における実際のブログ書き込み数を集計して実測値とする。そして、分析・可視化サーバ２０は、予測期間内の実測値の合計を、予測期間内の予測値の合計で除算して得た値をイベントの影響度とする。
【０１１５】
分析・可視化サーバ２０は、予測値、実測値を表示するための分析結果データを利用者ＰＣ１０へ送信する。また、分析・可視化サーバ２０は、イベントの種類、影響度を示す情報を更に含む予測式データを予測パラメータＤＢ４０に格納する。
【０１１６】
図１８は、予測式データの示す内容をまとめた図である。同図に示すように、予測式データは、「カテゴリ」、「キーワード」、「教師データ期間」、「トレンドパラメータ」、および「周期パラメータ」を示す情報に加え、「イベント」および「影響度」を示す情報を含む。
【０１１７】
「イベント」は、分析対象のイベントの種類である。「影響度」は、実測値と予測値とから得られるイベントの影響度である。
【０１１８】
図１９は、本実施形態の分析・可視化サーバ２０の動作を示すフローチャートである。この動作は、分析・可視化サーバ２０が影響度分析要求データを受信したときに開始する。
【０１１９】
分析・可視化サーバ２０は、教師データ期間のブログ書き込み数の集計結果から、予測期間における予測値を算出するための時系列分析処理を実行する（ステップＴ１）。図１９における時系列分析処理は、図８に示した時系列分析処理と同様の処理である。
【０１２０】
分析・可視化サーバ２０は、予測期間内の実測値の合計を、予測期間内の予測値の合計で除算することにより、イベントの影響度を算出する（ステップＴ２）。
【０１２１】
図２０を参照して、イベントの影響度の算出結の一例について説明する。同図は、予測値と実測値とを示す折れ線グラフである。同図における実線は、実際に集計されたブログの書き込み数をプロットしたものである。一点鎖線は、予測値をプロットしたものである。「Ｐ１」の時期はイベントの開始時期である。斜線部分は、予測値と実測値との間の差分を示す部分である。
【０１２２】
図２０に示すように、イベントの開始時期「Ｐ１」から、予測期間内において、予測値に対し、実際のブログ書き込み数が増加している。この増加分がイベントの影響により増加した書き込み数と推測される。
【０１２３】
なお、本実施形態では、分析・可視化サーバ２０は、実測値の合計を予測値の合計で除算した値を影響度としているが、予測値に対する実測値の増減の度合を示す値であれば、他の方法で影響度を算出してもよい。例えば、予測期間内のある日付における実測値と予測値との間の差分自体を影響度としてもよい。
【０１２４】
以上説明したように、本実施形態によれば、分析・可視化サーバ２０は、イベントの効果を定量的に測定することができる。
【０１２５】
（第３の実施形態）
本発明の第３の実施形態について説明する。本実施形態の時系列分析システムは、成長度スコアを更に算出できる点で、第１の実施形態の時系列分析システムと異なる。なお、本実施形態の時系列分析システムについて、図１〜図１６で説明した構成と同様な構成の詳細な説明を省略し、第１の実施形態と異なる点について詳しく説明する。
【０１２６】
本実施形態の利用者ＰＣ１０は、時系列分析要求データの他、成長度スコアの算出を要求する成長度スコア算出要求データの入力を受け付ける。
【０１２７】
成長度スコアとは、ある時期を教師データ期間の終日として算出された予測値に対する、それ以降の時期を教師データ期間の終日として算出された予測値の増減の度合を示す値である。例えば、教師データ期間の終日が遅い方の予測値の合計を、教師データ期間の終日が早い方の予測値の合計で除算した値を成長度スコアとする。
【０１２８】
分析・可視化サーバ２０は、成長度スコア算出要求データに基づいて、成長度スコアを算出する。
【０１２９】
ここでは、最後にブログ書き込み数が集計された日付を「現在」とする。利用者は、「比較対象時期」として、その現在以前の日付を設定する。利用者は、現在を基準として現在以降の所定期間を予測期間に設定し、現在以前の所定期間を教師データ期間に設定する。
【０１３０】
例えば、利用者は、現在より１週間前の日付を比較対象時期に設定し、現在以降の２カ月間を予測期間に設定し、現在以前の１年間を教師データ期間に設定する。
【０１３１】
利用者ＰＣ１０は、比較対象時期、予測期間、および教師データ期間を示す情報を含む成長度スコア算出要求データを分析・可視化サーバ２０へ出力する。
【０１３２】
図２１は、成長度スコア算出要求データに含まれる情報を示す図である。同図を参照すると、影響度分析要求データには、「カテゴリ」、「キーワード」、「教師データ期間」、「予測期間」、および「比較対象時期」を示す情報を含む。
【０１３３】
「カテゴリ」、「キーワード」は、第１の実施形態と同様である。「教師データ期間」、「予測期間」は、現在を基準に設定された教師データ期間、予測期間である。「比較対象時期」は、現在以前の所定の日付である。
【０１３４】
図２２は、本実施形態の分析・可視化サーバ２０の動作を示すフローチャートである。この動作は、分析・可視化サーバ２０が成長度スコア算出要求データを受信したときに開始する。
【０１３５】
分析・可視化サーバ２０は、成長度スコア算出要求データの示す、現在を基準とした「教師データ期間」、「予測期間」を設定する（ステップＴ３）。
【０１３６】
分析・可視化サーバ２０は、教師データ期間のブログ書き込み数の集計結果から、予測期間における予測値を算出するための時系列分析処理を実行する（ステップＴ４）。図２２における時系列分析処理は、図８に示した時系列分析処理と同様の処理である。
【０１３７】
次いで、分析・可視化サーバ２０は、成長度スコア算出要求データの示す「教師データ期間」の長さは変えずに、その終日が「比較対象時期」になるように、「教師データ期間」を設定し直す。
【０１３８】
また、分析・可視化サーバ２０は、成長度スコア算出要求データの示す「予測期間」の長さは変えずに、その初日が「比較対象時期」になるように、「予測期間」を設定し直す。
【０１３９】
つまり、分析・可視化サーバ２０は、比較対象時期を基準とした教師データ期間、予測期間を設定する（ステップＴ５）。
【０１４０】
分析・可視化サーバ２０は、新たに設定された教師データ期間、予測期間に基づいて、時系列分析処理を実行する（ステップＴ６）。
【０１４１】
分析・可視化サーバ２０は、現在を基準として算出された予測期間内の予測値の合計を、比較対象時期を基準として算出された予測期間内の予測値の合計で除算することにより、成長度スコアを算出する（ステップＴ７）。ステップＴ７の後、分析・可視化サーバ２０は、動作を終了する。
【０１４２】
図２３を参照して、成長度スコアの算出結の一例について説明する。同図は、現在、比較対象時期のそれぞれを基準として算出された予測値をプロットした折れ線グラフである。同図における細い実線は、現在を基準とした教師データ期間の集計結果の一部をプロットしたものである。太い実線は、現在を基準とした予測期間内の予測値をプロットしたものである。一点鎖線は、比較対象時期を基準とした予測期間内の予測値をプロットしたものである。比較対象時期は１週間前に設定されたものとする。
【０１４３】
図２３に示すように、１週間前を基準として算出された予測値よりも、現在を基準として算出された予測値は増加している。この予測値の増加率が成長度スコアとして算出される。
【０１４４】
なお、本実施形態では、分析・可視化サーバ２０は、現在を基準として算出された予測期間内の予測値の合計を、比較対象時期を基準として算出された予測期間内の予測値の合計で除算することにより、成長度スコアを算出している。
【０１４５】
しかし、過去の時点の予測値に対する、それ以降における予測値の増減の度合を示す値であれば、他の算出方法で算出した値を成長度スコアとしてもよい。例えば、分析・可視化サーバ２０は、現在を基準として算出された予測期間内の予測値の合計と、比較対象時期を基準として算出された予測期間内の予測値の合計との差分を成長度スコアとしてもよい。
【０１４６】
また、本実施形態では、比較対象時期を１つだけ設定しているが、複数設定してもよいのは勿論である。例えば、分析・可視化サーバ２０は、現在、その１週間前、１カ月前の３つの時点が設定されたとき、１週間前に対する現在の成長度、１カ月前に対する現在の成長度を算出する。
【０１４７】
以上説明したように、本実施形態によれば、過去の予測値と、それ以降の予測値との差を成長度として定量化することができる。
【０１４８】
（第４の実施形態）
本発明の第４の実施形態について説明する。本実施形態の時系列分析システムは、イベントの効果を更に予測できる点で、第２の実施形態の時系列分析システムと異なる。なお、本実施形態の時系列分析システムについて、図１７〜図２０で説明した構成と同様な構成の詳細な説明を省略し、第２の実施形態と異なる点について詳しく説明する。
【０１４９】
本実施形態の利用者ＰＣ１０は、時系列分析要求データの他、イベントの効果の予測を要求するシミュレーション要求データの入力を受け付ける。
【０１５０】
利用者は、カテゴリ、キーワード、イベントの種類、イベントの開始時期、および予測期間を設定する。利用者ＰＣ１０は、これらを示す情報を含むシミュレーション要求データを送信する。
【０１５１】
分析・可視化サーバ２０は、予測パラメータＤＢ４０に格納されたデータの中から、シミュレーション要求データの示すカテゴリと同一カテゴリに属し、類似するイベントを検索する。類似するイベントとは、キーワードの一部または全部が一致し、イベント種類が同じイベントである。分析・可視化サーバ２０は、検索したイベントの予測式の成分を読み出す。また、分析・可視化サーバ２０は、ブログＤＢ３０のブログデータから、検索したイベントのイベント開始時期以降の実測値を取得する。
【０１５２】
分析・可視化サーバ２０は、予測式から算出された予測値と、実測値とを表示するためのデータを利用者ＰＣ１０へ送信する。
【０１５３】
図２４はシミュレーション要求データに含まれる情報を示す図である。同図を参照すると、シミュレーション要求データには、「カテゴリ」、「キーワード」、「イベント種類」、「イベント開始時期」、および「予測期間」を示す情報を含む。「カテゴリ」、「キーワード」、「イベント種類」、「イベント開始時期」、および「予測期間」の内容は、第２の実施形態と同様である。
【０１５４】
図２５は、本実施形態の分析・可視化サーバ２０の実行するイベント効果予測処理を示すフローチャートである。イベント効果予測処理は、分析・可視化サーバ２０がシミュレーション要求データを受信したときに開始する。
【０１５５】
分析・可視化サーバ２０は、シミュレーション要求データの示すカテゴリと同一カテゴリに属し、類似するイベントを検索する。分析・可視化サーバ２０は、検索したイベントの予測式の成分を読み出す。複数のイベントが検索された場合、それらのイベントの成分の平均値を算出する（ステップＴ１１）。
【０１５６】
分析・可視化サーバ２０は、取得した成分から得た予測式で、シミュレーション要求データの示す予測期間における予測値を算出する（ステップＴ１２）。
【０１５７】
分析・可視化サーバ２０は、ブログＤＢ３０から、検索したイベントのブログデータを読み出し、そのイベントの開始時期以降のブログ書き込み数の実測値を集計する。複数のイベントが検索された場合、それらのイベントの実測値の平均値を算出する。分析・可視化サーバ２０は、イベント開始時期から、シミュレーション要求データの示す予測期間と同じ長さの期間が経過するまでの期間について集計する。
【０１５８】
分析・可視化サーバ２０は、予測値と実測値とを表示するためのデータを利用者ＰＣ１０へ送信する（ステップＴ１３）。ステップＴ１３の後、分析・可視化サーバ２０は、イベント効果予測処理を終了する。
【０１５９】
以上説明したように、分析・可視化サーバ２０は、過去の類似イベントの予測値と、実測値とを求めるので、イベントを実行した場合の効果をシミュレートできる。
【符号の説明】
【０１６０】
１    時系列データ分析システム
１０    利用者ＰＣ
２０    分析・可視化サーバ
３０    ブログＤＢ
４０    予測パラメータＤＢ
２０１    ブログ抽出・集計部
２０２    トレンド成分学習部
２０３    周期成分学習部
２０４    分析結果作成部
２０１１    測定部
２０１２    特異点検出・補正部
【図１】
2011081491.tif 000022
【図２】
2011081491.tif 000023
【図３】
2011081491.tif 000024
【図４】
2011081491.tif 000025
【図５】
2011081491.tif 000026
【図６】
2011081491.tif 000027
【図７】
2011081491.tif 000028
【図８】
2011081491.tif 000029
【図９】
2011081491.tif 000030
【図１０】
2011081491.tif 000031
【図１１】
2011081491.tif 000032
【図１２】
2011081491.tif 000033
【図１３】
2011081491.tif 000034
【図１４】
2011081491.tif 000035
【図１５】
2011081491.tif 000036
【図１６】
2011081491.tif 000037
【図１７】
2011081491.tif 000038
【図１８】
2011081491.tif 000039
【図１９】
2011081491.tif 000040
【図２０】
2011081491.tif 000041
【図２１】
2011081491.tif 000042
【図２２】
2011081491.tif 000043
【図２３】
2011081491.tif 000044
【図２４】
2011081491.tif 000045
【図２５】
2011081491.tif 000046