(57)【要約】  本技術は、効率良く、高精度に目的変数の値を推定することができるようにする情報処理装置、情報処理方法、およびプログラムに関する。ログ取得部は、推定対象の目的変数に対応する目的時系列データと、目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得する。モデルパラメータ更新部は、取得された目的時系列データと複数の説明時系列データとを用いて、確率モデルのパラメータを学習する。ログ選択部は、学習により得られた確率モデルのパラメータに基づいて、ログ取得部が取得する説明時系列データに対応する説明変数を選択する。推定部は、選択部の選択結果に基づいてログ取得部が取得した複数の説明時系列データを用いて、目的変数の値を推定する。本技術は、例えば、機器の消費電力を推定する情報処理装置に適用できる。
2013187295.tif 000002
【特許請求の範囲】
【請求項１】
推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得する取得部と、
取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータを学習する学習部と、
学習により得られた前記確率モデルのパラメータに基づいて、前記取得部が取得する前記説明時系列データに対応する前記説明変数を選択する選択部と、
前記選択部の選択結果に基づいて前記取得部が取得した複数の前記説明時系列データを用いて、前記目的変数の値を推定する推定部と
を備える情報処理装置。
【請求項２】
前記学習部は、前記目的変数と前記複数の説明変数との関係を隠れマルコフモデルで学習する
請求項１に記載の情報処理装置。
【請求項３】
前記目的変数が、前記隠れマルコフモデルの隠れ状態と１対１に対応する線形回帰係数と、前記説明変数との線形回帰モデルで表される
請求項２に記載の情報処理装置。
【請求項４】
前記選択部は、前記線形回帰係数が所定の閾値より小さい前記説明変数を、前記取得部が時系列データを取得しない説明変数として選択する
請求項３に記載の情報処理装置。
【請求項５】
情報処理装置が、
推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得し、
取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータを学習し、
学習により得られた前記確率モデルのパラメータに基づいて、取得する前記説明時系列データに対応する前記説明変数を選択し、
選択結果に基づいて取得した複数の前記説明時系列データを用いて、前記目的変数の値を推定する
ステップを含む情報処理方法。
【請求項６】
コンピュータを、
推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得する取得部と、
取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータを学習する学習部と、
学習により得られた前記確率モデルのパラメータに基づいて、前記取得部が取得する前記説明時系列データに対応する前記説明変数を選択する選択部と、
前記選択部の選択結果に基づいて前記取得部が取得した複数の前記説明時系列データを用いて、前記目的変数の値を推定する推定部
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本技術は、情報処理装置、情報処理方法、およびプログラムに関し、特に、効率良く、高精度に目的変数の値を推定することができるようにする情報処理装置、情報処理方法、およびプログラムに関する。
【背景技術】
【０００２】
電子機器の部品であるCPUの負荷率等を説明変数として、複数の説明変数と、その電力係数とからなる線形式で電子機器の消費電力をモデル化し、部品の稼働状態に基づいて、電子機器の消費電力を推定する技術が提案されている（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１０−２２５３３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、特許文献１の技術では、時刻ｔの消費電力を推定する場合、時刻tにおける部品の稼働状態（を示す値）のみが用いられる。したがって、特許文献１の技術では、時刻ｔに至るまでの稼働状況の履歴を考慮した消費電力の推定ができなかった。【０００５】
また、特許文献１の技術では、多重共線性の問題が起こることを回避するために、消費電力の推定に使用する取得データ（説明変数に対応するデータ）の種類を、人が判断し、取捨選択する必要があった。【０００６】
本技術は、このような状況に鑑みてなされたものであり、現在までの稼働状況の履歴を考慮しつつ、取得データを自動で取捨選択して、効率良く、高精度に目的変数の値を推定することができるようにするものである。
【課題を解決するための手段】
【０００７】
本技術の一側面の情報処理装置は、推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得する取得部と、取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータを学習する学習部と、学習により得られた前記確率モデルのパラメータに基づいて、前記取得部が取得する前記説明時系列データに対応する前記説明変数を選択する選択部と、前記選択部の選択結果に基づいて前記取得部が取得した複数の前記説明時系列データを用いて、前記目的変数の値を推定する推定部とを備える。【０００８】
本技術の一側面の情報処理方法は、情報処理装置が、推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得し、取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータを学習し、学習により得られた前記確率モデルのパラメータに基づいて、取得する前記説明時系列データに対応する前記説明変数を選択し、選択結果に基づいて取得した複数の前記説明時系列データを用いて、前記目的変数の値を推定するステップを含む。【０００９】
本技術の一側面のプログラムは、コンピュータを、推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得する取得部と、取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータを学習する学習部と、学習により得られた前記確率モデルのパラメータに基づいて、前記取得部が取得する前記説明時系列データに対応する前記説明変数を選択する選択部と、前記選択部の選択結果に基づいて前記取得部が取得した複数の前記説明時系列データを用いて、前記目的変数の値を推定する推定部として機能させるためのものである。【００１０】
本技術の一側面においては、推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとが取得され、取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータが学習され、学習により得られた前記確率モデルのパラメータに基づいて、取得する前記説明時系列データに対応する前記説明変数が選択され、その選択結果に基づいて取得された複数の前記説明時系列データを用いて、前記目的変数の値が推定される。【００１１】
なお、プログラムは、伝送媒体を介して伝送することにより、又は、記録媒体に記録して、提供することができる。【００１２】
情報処理装置は、独立した装置であっても良いし、１つの装置を構成している内部ブロックであっても良い。
【発明の効果】
【００１３】
本技術の一側面によれば、効率良く、高精度に目的変数の値を推定することができる。
【図面の簡単な説明】
【００１４】
【図１】機器消費電力推定処理の概要について説明する図である。
【図２】本技術が適用された情報処理装置の一実施の形態の構成例を示すブロック図である。
【図３】データ収集処理を説明するフローチャートである。
【図４】モデルパラメータ学習処理を説明するフローチャートである。
【図５】消費電力推定処理を説明するフローチャートである。
【図６】HMMのグラフィカルモデルを示す図である。
【図７】モデルパラメータ更新処理を説明するフローチャートである。
【図８】本技術が適用されたコンピュータの一実施の形態の構成例を示すブロック図である。
【発明を実施するための形態】
【００１５】
［情報処理装置の処理概要］
初めに、図１を参照して、本技術が適用された情報処理装置で実現される機器消費電力推定処理の概要について説明する。
【００１６】
後述する図２に示される情報処理装置１は、機器（電子機器）内部の所定の部品の稼働状態を示す時系列データを取得する。ここで、取得される時系列データは、例えば、図１に示すように、CPUの使用率、メモリ（RAM）のアクセスレート（acess rate）、リムーバブルメディアの書き込み／読み出し回数、などである。また、情報処理装置１は、稼働状態を示す時系列データを取得したときの機器の消費電力の時系列データも同時に取得する。
【００１７】
そして、情報処理装置１は、複数種類の稼働状態と消費電力との関係を所定の学習モデルで予め学習する。情報処理装置１により学習された学習モデルを、以下では、消費電力変動モデルともいう。
【００１８】
情報処理装置１は、消費電力変動モデル（のパラメータ）を学習により決定した後は、学習された消費電力変動モデルを用いて、新たに入力された複数種類の稼働状態を示す時系列データのみに基づいて、現時点の機器の消費電力を推定することが可能となる。そして、情報処理装置１は、例えば、推定結果である現在の消費電力をディスプレイにリアルタイムに表示する。
【００１９】
情報処理装置１の処理は、大きく二つの処理に分かれる。一つは、複数種類の稼働状態と消費電力との関係を所定の学習モデルで学習する学習処理であり、もう一つは、学習処理により得られた学習モデルを用いて、機器の現在の消費電力を推定する消費電力推定処理である。
【００２０】
機器（電子機器）は、例えば、スマートフォン、タブレット端末、などの携帯型の端末や、据え置き型のパーソナルコンピュータなどである。また、機器は、テレビジョン受像機、コンテンツ記録再生装置などでもよい。情報処理装置１は、消費電力を推定する対象の機器の一部として、その機器に含まれていてもよいし、推定対象の機器と異なる装置で構成され、機器と接続されて実行されるものでもよい。また、情報処理装置１は、複数の装置で構成される情報処理システムとして構成することもできる。
【００２１】
［情報処理装置の機能ブロック図］
図２は、情報処理装置１の機能的構成を示すブロック図である。
【００２２】
情報処理装置１は、消費電力計測部１１、消費電力時系列入力部１２、ログ取得部１３、機器制御部１４、ログ時系列入力部１５、時系列履歴保存部１６、モデル学習部１７、消費電力推定部１８、および推定消費電力表示部１９により構成される。
【００２３】
消費電力計測部１１は、例えば、電力計（クランプメータ）、テスタ、オシロスコープなどで構成され、機器の電源ラインに接続し、各時刻における機器の消費電力を計測し、計測結果を消費電力時系列入力部１２に出力する。
【００２４】
消費電力時系列入力部１２は、消費電力計測部１１から供給される各時刻の消費電力値を所定時間蓄積し、消費電力値の時系列データを生成する。生成される消費電力値の時系列データ（以下、消費電力時系列データともいう。）は、取得したときの時刻と消費電力値との組が、所定の期間だけ集められたデータとなる。
【００２５】
ログ取得部１３は、機器内部の所定部品の稼働状態を示すデータをログ情報として取得する。ログ取得部１３は、複数種類のログ情報を同時に取得し、ログ時系列入力部１５に出力する。ログ取得部１３が取得するログ情報の種類としては、例えば、CPU使用率、GPU使用率、Wifi通信量、移動体通信回線の通信量（３G通信量）、ディスプレイ輝度、起動中の各アプリケーションのリストとCPU使用比率の対データ、などが挙げられるが、これらに限定されるわけではない。
【００２６】
機器制御部１４は、消費電力変動モデルの学習処理において、実際に想定される様々な状態における消費電力を消費電力変動モデルとして学習するため、様々な状態を作り出す機器の制御を行う。例えば、機器制御部１４は、ゲーム、表計算処理ソフトなど、複数種類のアプリケーションを同時に起動させて処理を実行させたり、データ通信を実行および停止させたりなどして、想定される稼働状態の組み合わせを機器に実行させる。
【００２７】
ログ時系列入力部１５は、学習処理において、ログ取得部１３から供給される各時刻の稼働状態を示すログ情報を所定時間蓄積し、その結果得られるログ時系列データを、時系列履歴保存部１６に出力する。
【００２８】
また、ログ時系列入力部１５は、消費電力推定処理において、ログ取得部１３から供給される各時刻のログ情報を所定時間蓄積し、その結果得られるログ時系列データを、消費電力推定部１８に出力する。
【００２９】
なお、ログ情報の種類が、例えば、起動中の各アプリケーションのリストとCPU使用比率の対データである場合には、起動中の各アプリケーションのリストが、機器制御部１４からログ時系列入力部１５に供給され、CPU使用比率が、ログ取得部１３からログ時系列入力部１５に供給される。
【００３０】
ログ時系列入力部１５は、異常値を除去する処理などのデータ加工を必要に応じて実行する。異常値を除去する処理としては、例えば、 “On-line outlier detection and data cleaning” ,Computers and Chemical Engineering 28 (2004) ,1635-1647, by Liu et al で提案されている処理を採用することができる。
【００３１】
時系列履歴保存部１６は、消費電力時系列入力部１２から供給される消費電力時系列データと、ログ時系列入力部１５から供給されるログ時系列データを、保存（記憶）する。時系列履歴保存部１６に保存された消費電力時系列データとログ時系列データは、モデル学習部１７が消費電力変動モデルを学習（更新）する際に利用される。
【００３２】
モデル学習部１７は、モデルパラメータ更新部２１、モデルパラメータ記憶部２２、およびログ選択部２３により構成される。
【００３３】
モデルパラメータ更新部２１は、時系列履歴保存部１６に保存された消費電力時系列データとログ時系列データを用いて、消費電力変動モデルを学習し、その結果得られる消費電力変動モデルのパラメータを、モデルパラメータ記憶部２２に記憶させる。以下では、消費電力変動モデルのパラメータを、単に、モデルパラメータともいう。
【００３４】
また、モデルパラメータ更新部２１は、新たな消費電力時系列データとログ時系列データが時系列履歴保存部１６に保存された場合、その新たな時系列データを用いてモデルパラメータ記憶部２２に保持されている消費電力変動モデルのパラメータを更新する。
【００３５】
モデルパラメータ更新部２１は、消費電力変動モデルを表す確率モデルとして、機器の動作状態を隠れ状態Sとして表す隠れマルコフモデル（HMM:Hidden Markov Model）と関連ベクタマシン（RVM:Relevance Vector Machine）とを混合した確率モデルであるHMM+RVMを採用する。HMM+RVMの詳細については後述する。
【００３６】
モデルパラメータ記憶部２２は、モデルパラメータ更新部２１により更新（学習）された消費電力変動モデルのパラメータを記憶する。モデルパラメータ記憶部２２に記憶されている消費電力変動モデルのパラメータは、消費電力推定部１８に供給される。
【００３７】
ログ選択部２３は、ログ取得部１３が取得している複数種類のログ情報のうち、不要なログ情報（の種類）を選択制御する。より具体的には、ログ選択部２３は、モデルパラメータ記憶部２２に記憶されている消費電力変動モデルのパラメータ（の値）に基づいて、不要なログ情報を判定する。そして、ログ選択部２３は、判定結果に基づいて、不要と判定されたログ情報を取得しないようにログ取得部１３を制御する。
【００３８】
消費電力推定部１８は、学習処理により得られた消費電力変動モデルのパラメータを、モデルパラメータ記憶部２２から取得する。そして、消費電力推定部１８は、消費電力推定処理において、ログ時系列入力部１５から供給される、現在時刻から一定時間前までのログ時系列データを、学習した消費電力変動モデルに入力し、現在時刻の消費電力値を推定する。推定結果の消費電力値は、推定消費電力表示部１９に供給される。
【００３９】
推定消費電力表示部１９は、消費電力推定部１８から供給される現在時刻の消費電力値を、所定の方法で表示する。例えば、推定消費電力表示部１９は、現在時刻の消費電力値をデジタル表示したり、一定時間前から現在時刻までの消費電力値の推移をグラフ化して表示する。
【００４０】
情報処理装置１は、以上のように構成されている。
【００４１】
［データ収集処理のフローチャート］
図３のフローチャートを参照して、情報処理装置１の学習処理の一部であり、モデルパラメータを計算するためのデータを収集するデータ収集処理について説明する。
【００４２】
初めに、ステップＳ１において、消費電力計測部１１が、機器の消費電力の計測を開始する。ステップＳ１の処理後、消費電力が一定時間間隔で計測され、計測結果が、消費電力時系列入力部１２に、順次、出力される。
【００４３】
ステップＳ２において、機器制御部１４が、複数種類のアプリケーションを起動し、実行させる。
【００４４】
ステップＳ３において、ログ取得部１３が、複数種類のログ情報の取得を開始する。ステップＳ３の処理後、複数種類のログ情報が一定時間間隔で取得され、ログ時系列入力部１５に、順次、出力される。
【００４５】
ステップＳ１乃至Ｓ３の各処理は、任意の順番で実行することができる。また、ステップＳ１乃至Ｓ３の各処理を同時に実行してもよい。
【００４６】
ステップＳ４において、消費電力時系列入力部１２が、消費電力計測部１１から供給された各時刻の消費電力値を所定時間蓄積し、消費電力時系列データを生成する。消費電力時系列入力部１２は、生成した消費電力時系列データを、時系列履歴保存部１６に供給する。
【００４７】
ステップＳ５において、ログ時系列入力部１５が、ログ取得部１３から供給される各時刻のログ情報を所定時間蓄積し、ログ時系列データを生成する。ログ時系列入力部１５は、生成したログ時系列データを、時系列履歴保存部１６に供給する。
【００４８】
ステップＳ６において、時系列履歴保存部１６は、消費電力時系列入力部１２から供給された消費電力時系列データと、ログ時系列入力部１５から供給されたログ時系列データを保存する。
【００４９】
以上のステップＳ１乃至Ｓ６の処理が実行されると、ステップＳ２で機器制御部１４が制御した所定の一つの動作条件での学習用データ（消費電力時系列データとログ時系列データの組）が、時系列履歴保存部１６に保存される。
【００５０】
情報処理装置１は、動作条件を異なるように様々に変更して、上述のデータ収集処理を所定回数繰り返すことで、想定される様々な稼働状態における学習用データを蓄積する。換言すれば、情報処理装置１は、ステップＳ２の処理を様々に変更して、ステップＳ１乃至Ｓ６の処理を所定回数繰り返すことで、想定される様々な稼働状態における学習用データを、時系列履歴保存部１６に保存する。
【００５１】
［モデルパラメータ学習処理のフローチャート］
次に、図４のフローチャートを参照して、情報処理装置１の学習処理の一部であり、データ収集処理により収集された学習用データを用いてモデルパラメータを求めるモデルパラメータ学習処理について説明する。
【００５２】
初めに、ステップＳ２１において、モデルパラメータ更新部２１が、現在のモデルパラメータを、モデルパラメータ記憶部２２から取得する。モデルパラメータ更新部２１が、消費電力変動モデルの学習を初めて行う場合、モデルパラメータ更新部２１には、モデルパラメータの初期値が記憶されている。
【００５３】
ステップＳ２２において、モデルパラメータ更新部２１が、時系列履歴保存部１６に保存されている消費電力時系列データとログ時系列データを取得する。
【００５４】
ステップＳ２３において、モデルパラメータ更新部２１が、モデルパラメータ記憶部２２から取得した現在のモデルパラメータを初期値として、時系列履歴保存部１６から取得した新たな消費電力時系列データとログ時系列データを用いて、モデルパラメータを更新する。
【００５５】
ステップＳ２４において、モデルパラメータ更新部２１が、更新後のモデルパラメータを、モデルパラメータ記憶部２２に供給し、記憶させる。モデルパラメータ記憶部２２は、モデルパラメータ更新部２１から供給された更新後のモデルパラメータを、現在のモデルパラメータに上書きして記憶する。
【００５６】
ステップＳ２５において、ログ選択部２３が、モデルパラメータ記憶部２２に記憶されている更新後のモデルパラメータに基づいて、不要なログ情報を判定する。そして、ログ選択部２３が、判定結果に基づいて、不要と判定されたログ情報を取得しないようにログ取得部１３を制御する。このログ選択部２３の選択制御は、ログ取得部１３がログ情報の取得処理（ステップＳ３の処理）を次に実行するときから反映される。
【００５７】
以上のようにして、時系列履歴保存部１６に保存された新たな消費電力時系列データとログ時系列データを用いたモデルパラメータの学習（更新）が実行される。
【００５８】
［消費電力推定処理のフローチャート］
次に、図５のフローチャートを参照して、学習したモデルパラメータを用いて、現在の稼働状態における消費電力を推定する消費電力推定処理について説明する。
【００５９】
初めに、ステップＳ４１において、消費電力推定部１８が、学習処理により得られたモデルパラメータを、モデルパラメータ記憶部２２から取得する。
【００６０】
ステップＳ４２において、ログ取得部１３が、現在時刻における複数種類のログ情報を取得し、ログ時系列入力部１５に出力する。ステップＳ４２の処理では、ログ選択部２３により選択制御されたログ情報の種類についてのみ取得される。
【００６１】
ステップＳ４３において、ログ時系列入力部１５が、ログ取得部１３から供給される現在時刻のログ情報を一時保存するとともに、現在時刻から一定時間前までのログ時系列データを、消費電力推定部１８に供給する。現在時刻のログ情報がログ取得部１３から供給されることにより、保存しておく必要がなくなった古いログ情報は削除される。
【００６２】
ステップＳ４４において、消費電力推定部１８が、学習した消費電力変動モデルを用いた消費電力推定処理を実行する。すなわち、消費電力推定部１８は、ログ時系列入力部１５からのログ時系列データを消費電力変動モデルに入力し、現在時刻の消費電力値を推定（算出）する。推定結果の消費電力値は、推定消費電力表示部１９に供給される。
【００６３】
ステップＳ４５において、推定消費電力表示部１９は、消費電力推定部１８から供給された現在時刻の消費電力値（の推定値）を、所定の方法で表示して、処理を終了する。
【００６４】
以上のステップＳ４１乃至Ｓ４５の処理が、例えば、ログ取得部１３で新たなログ情報が取得されるごとに実行される。
【００６５】
［HMM+RVMの詳細説明］
次に、本実施の形態において、機器の消費電力変動を学習する学習モデルとして採用しているHMM+RVMの詳細について説明する。
【００６６】
まず、一般的なHMMの確率モデルを示す。図６は、HMMのグラフィカルモデルを示している。
【００６７】
観測される消費電力時系列データをY=｛Y1，Y2，Y3，・・・，Yt，・・・，YT｝、ログ時系列データをX=｛X1，X2，X3，・・・，Xt，・・・，XT｝、背後に想定される機器の隠れ状態に対応する隠れ変数の時系列データをS=｛S1，S2，S3，・・・，St，・・・ST｝とすると、隠れ変数Stと観測データXt，Ytの同時確率は、次式（１）で与えられる。
【数１】
2013187295.tif 000003
式（１）において、Ytは、消費電力計測部１１の時刻ｔにおける消費電力の計測値ytであり、１次元のデータである。Xtは、ログ取得部１３で取得される時刻ｔにおける複数種類（Dx個）のログ情報xt1,xt2,・・・,xtDxであり、Dx次元のベクトルを表す。
【００６８】
式（１）中の、P(S1)は初期確率、P(St｜St-1)は隠れ状態St-1から隠れ状態Stへの状態遷移確率、P(Yt｜St)およびP(Xt｜St)は観測確率を表す。観測確率P(Yt｜St)およびP(Xt｜St)は、それぞれ、以下の式（２）および式（３）で計算される。
【数２】
2013187295.tif 000004
ここで、式（２）のρStは隠れ状態Stでの消費電力（出力Y）の平均値を表し、βは、出力Yにかかるガウスノイズの大きさ（分散）を表す。出力Yにかかるガウスノイズの大きさβは、隠れ状態に依存しないものとするが、依存するように定義することも容易にできる。同様に、式（３）のμStは隠れ状態Stでのログ情報（入力X）の平均値を表し、ΣStは、入力Xの分散を表す。また、Tは、転置を表す。
【００６９】
以上の一般的なHMMの確率モデルに対して、本実施の形態において消費電力変動モデルとして採用するHMM+RVMの確率モデルは、以下の式（４）で表される。
【数３】
2013187295.tif 000005
【００７０】
式（４）のｗStは、隠れ状態Stにおける入力Xと出力Yの線形回帰係数を表し、P(ｗS)は、線形回帰係数ｗSの事前確率分布を表す。線形回帰係数ｗSの事前確率分布P(ｗS)は、式（５）のとおり、平均ν0、分散α-1（αの逆行列）のガウス分布で仮定される。
【数４】
2013187295.tif 000006
なお、平均ν0は、０に設定されるとともに、分散α-1は対角行列とする。
【００７１】
また、式（４）の観測確率P(Xt,Yt｜St,ｗSt)は、式（６）のように観測確率P(Yt｜Xt,ｗSt)と観測確率P(Xt｜St)の積で表すことができ、観測確率P(Yt｜Xt,ｗSt)および観測確率P(Xt｜St)は、それぞれ、式（７）および式（８）で表される。
【数５】
2013187295.tif 000007
【００７２】
式（７）は、出力Ytが、隠れ状態Stの線形回帰係数ｗStを用いた入力Xtの線形回帰モデルで表されることを意味し、式（８）は、入力Xtが、平均μSt、分散ΣStのガウス分布で表されることを意味している。したがって、隠れ状態Stは、HMMのような、機器の隠れた状態を表すというよりも、消費電力値（出力Yt）とログ情報（入力Xt）との関係（確率的な関係性）を表す線形回帰モデルの変数（隠れ変数）であると言うことができる。また、出力Ytは、線形回帰モデルの目的変数であり、入力Xtとしての複数種類のログ情報が、線形回帰モデルの説明変数に相当する。
【００７３】
図４を参照して説明したモデルパラメータ学習処理のステップＳ２３におけるモデルパラメータの更新では、モデルパラメータ更新部２１が、これらのパラメータ｛ｗ，ν0，α，β，μ，Σ，P(S|S')，P(S1)｝を更新する。モデルパラメータの更新は、HMMで使用される反復型のアルゴリズムであるEMアルゴリズムと同様の要領で実行される。
【００７４】
［アルゴリズム更新処理の詳細フロー］
図７のフローチャートを参照して、図４のステップＳ２３として実行される、モデルパラメータ更新処理の詳細について説明する。
【００７５】
初めに、ステップＳ６１において、モデルパラメータ更新部２１は、モデルパラメータの初期値を設定する。モデルパラメータの初期値は、例えば、乱数を用いて決定された値など、所定の方法で設定される。
【００７６】
ステップＳ６２において、モデルパラメータ更新部２１は、各隠れ状態Ｓに付随した線形回帰係数ｗSを更新する。線形回帰係数ｗSは、点推定ではなく、分布で推定して更新する。すなわち、各隠れ状態Ｓに付随した線形回帰係数ｗSの分布をガウス分布ｑ（ｗS）とすると、線形回帰係数ｗSのガウス分布ｑ（ｗS）は、次式（９）で計算される。
【数６】
2013187295.tif 000008
ここで、＜・＞q(S)は、隠れ状態Ｓについての期待値を表す。式（９）のP(X,Y,S｜ｗ)は、ｗ＝｛ｗS1,ｗS2,ｗS3,・・・,ｗSt,・・・,ｗST｝として、式（１０）で表される。
【数７】
2013187295.tif 000009
【００７７】
より具体的には、ガウス分布ｑ（ｗS）の平均をλS、分散をτSとすると、平均λSと分散τSは、以下の式（１１）および式（１２）で更新される。
【数８】
2013187295.tif 000010
ここで、ｑt(S)は、時刻ｔにおいて隠れ状態Sにある確率を表し、後述する式（２３）で表される。
【００７８】
ステップＳ６３において、モデルパラメータ更新部２１は、次式（１３）により、出力Yのガウスノイズの大きさβを更新する。
【数９】
2013187295.tif 000011
【００７９】
ステップＳ６４において、モデルパラメータ更新部２１は、次式（１４）により、ある時刻において隠れ状態S'にいるとき、次の時刻で隠れ状態Sに遷移する確率である遷移確率P(S|S')を更新する。
【数１０】
2013187295.tif 000012
ここで、ｑt(S',S)は、時刻ｔとｔ＋１において、それぞれ隠れ状態S'とSに存在する確率を表す。
【００８０】
ステップＳ６５において、モデルパラメータ更新部２１は、次式（１５）乃至式（１７）により、初期確率分布P(S1)、隠れ状態Sの入力Xの平均μSおよび分散ΣSを更新する。
【数１１】
2013187295.tif 000013
【００８１】
ステップＳ６６において、モデルパラメータ更新部２１は、次式（１８）で表される、隠れ状態Sの確率ｑ（S）を計算する。
【数１２】
2013187295.tif 000014
【００８２】
具体的には、モデルパラメータ更新部２１は、時刻ｔにおいて状態Sにある確率ｑt(S)を、以下の手順で計算する。
【００８３】
まず、モデルパラメータ更新部２１は、状態Stのフォワード尤度α（St）とバックワード尤度β（St）を、次式（１９）および式（２０）により計算する。
【数１３】
2013187295.tif 000015
ここで、ｐ(Xt,Yt｜St)は、
【数１４】
2013187295.tif 000016
で計算される。
【００８４】
次に、モデルパラメータ更新部２１は、時刻ｔとｔ＋１において、それぞれ隠れ状態SとS'に存在する確率ｑt(S,S')を、次式（２２）により計算する。
【数１５】
2013187295.tif 000017
そして、得られた確率ｑt(S,S')を用いて、モデルパラメータ更新部２１は、式（２３）により、時刻ｔにおいて状態Sにある確率ｑt(S)を計算する。
【数１６】
2013187295.tif 000018
【００８５】
ステップＳ６７において、モデルパラメータ更新部２１は、線形回帰係数ｗSの事前確率分布P(ｗS)のパラメータ、すなわち、平均ν0および分散α-1を更新する。
【００８６】
平均ν0および分散α-1は、概念的には、以下の条件を満たすものとなる。
【数１７】
2013187295.tif 000019
ここで、N（ｗS；ν0，α-1）は、確率変数ｗSが平均ν0および分散α-1の正規分布に従うことを表し、KL（ｑ（ｗS）||N（ｗS；ν0，α-1））は、ｑ（ｗS）とN（ｗS；ν0，α-1）のカルバックライブラーダイバージェンスを表す。また、Argminは、全ての隠れ状態SについてのKL（ｑ（ｗS）||N（ｗS；ν0，α-1））の和（Σ）を最小にする変数（ν0および分散α）を意味する。
【００８７】
事前確率分布P(ｗS)の平均ν0および分散α-1は、具体的には、以下の式で計算することができる。
【数１８】
2013187295.tif 000020
【００８８】
ステップＳ６８において、モデルパラメータ更新部２１は、モデルパラメータの収束条件を満たしているかを判定する。例えば、モデルパラメータ更新部２１は、ステップＳ６２乃至６８の処理の繰り返し回数が予め設定した所定の回数に到達した場合、または、モデルパラメータの更新による状態尤度の変化量が所定値以内である場合に、モデルパラメータの収束条件を満たしていると判定する。
【００８９】
ステップＳ６８で、モデルパラメータの収束条件をまだ満たしていないと判定された場合、処理はステップＳ６２に戻り、ステップＳ６２乃至Ｓ６８の処理が繰り返される。
【００９０】
一方、ステップＳ６８で、モデルパラメータの収束条件を満たしたと判定された場合、モデルパラメータ更新部２１は、モデルパラメータ更新処理を終了する。
【００９１】
なお、更新するモデルパラメータの計算順序は、必ずしも上述したステップＳ６２乃至Ｓ６７の順番で行う必要はなく、任意の順番で行うことができる。
【００９２】
以上のモデルパラメータ更新処理により、モデルパラメータが更新されると、ステップＳ６７で計算された、線形回帰係数ｗSの事前確率分布P(ｗS)の分散α-1のうち、多くのα-1kが無限大となる。分散α-1kが無限大であれば、平均ν0,kは０に設定されているので、すべての線形回帰係数ｗSの第ｋ成分が強制的に０になる。これは、入力Xの第ｋ成分の重要性が低いことを表し、第ｋ成分の入力Xを使用しなくてもよいことを意味する。
【００９３】
そこで、図４を参照して説明したモデルパラメータ学習処理のステップＳ２５で実行される、不要なログ情報の判定処理では、線形回帰係数ｗSの第ｋ成分が所定の閾値（例えば、０．０１など）よりも小さいか否かを、ログ選択部２３が判定する。そして、所定の閾値より小さい線形回帰係数ｗSの成分に対応するログ情報が、不要なログ情報と判定され、ログ選択部２３は、不要と判定したログ情報を、次回以降使用しないように、ログ取得部１３を選択制御する。
【００９４】
なお、取得するログ情報の種類を新たに追加した場合には、追加したログ情報を加えて、学習処理、すなわち、図３のデータ収集処理および図４のモデルパラメータ学習処理を再び実行する必要がある。この場合、モデルパラメータの初期値には、モデルパラメータ記憶部２２に記憶されている現時点のモデルパラメータを用いることができる。
【００９５】
[消費電力推定処理の詳細]
次に、以上のようにして学習されたモデルパラメータを用いて、図５のステップＳ４４で実行される、消費電力推定処理の詳細について説明する。
【００９６】
消費電力推定処理においては、ログ時系列入力部１５から、ログ時系列データのみが取得される。ここで取得された消費電力推定処理時のログ時系列データを｛Xd1，Xd2，Xd3，・・・，Xdt，・・・，XdT｝で表すと、消費電力推定部１８は、以下の式（２７）を満たす隠れ状態S*tを求める。
【数１９】
2013187295.tif 000021
ここで、時刻ｔにおける隠れ状態Stとログ情報Xdtの同時確率分布P(｛St，Xdt｝)は、    
【数２０】
2013187295.tif 000022
で表される。
【００９７】
式（２７）は、取得されたログ時系列データ｛Xd1，Xd2，Xd3，・・・，Xdt，・・・，XdT｝が観測される尤度を、最も大にする状態遷移の過程（最尤状態系列）における時刻ｔの隠れ状態StをS*tとして求めるものである。最尤状態系列は、ビタビアルゴリズムを用いて求めることができる。
【００９８】
ビタビアルゴリズム、および、EMアルゴリズム（Baum-Welchアルゴリズム）の詳細については、例えば、“パターン認識と機械学習（下）”，Ｃ．Ｍ．ビショップ著，（英語原書：“Pattern Recognition and Machine Learning (Information Science and Statistics) ”，Christopher M. BishopSpringer, New York, 2006.）のP.347,P.333に記載されている。
【００９９】
消費電力推定部１８は、式（２７）を満たす隠れ状態S*tを求めた後、以下の式（２９）により、消費電力の推定値Y*tを求める（推定する）。
【数２１】
2013187295.tif 000023
したがって、消費電力の推定値Y*tは、隠れ状態S*tの線形回帰係数ｗSの平均λStと、入力Xtの内積により、求められる。
【０１００】
消費電力推定部１８は、図５のステップＳ４４において、以上のようにして、現在時刻の消費電力値Y*tを推定する。
【０１０１】
[HMMへの変形例]
これまで説明したHMM+RVMのアルゴリズムは、HMMに、RVMを適用したものということができる。したがって、上述したHMM+RVMのモデルパラメータを所定の条件に設定すると、通常のHMMにもなる。そこで、消費電力変動モデルとして、通常のHMMを適用し、推定を行う場合について以下で説明する。
【０１０２】
初めに、上述したHMM+RVMからHMMへの変形について説明する。
【０１０３】
HMM+RVMで入力Xとしたログ時系列データXの第（Dx+１）次元目に新しい成分を加えたものを、HMMの入力X~と考える。すなわち、HMMの時刻ｔにおける入力X~tは、
【数２２】
2013187295.tif 000024
とする。
【０１０４】
また、線形回帰係数ｗSの事前確率分布P(ｗS)の分散α-1（αの逆行列）を、以下のように、第（Dx+１）成分を無限大、それ以外の成分を０に設定する。これにより、線形回帰係数ｗSの事前確率分布P(ｗS)のパラメータが固定化される。
【数２３】
2013187295.tif 000025
【０１０５】
この場合、上述した式（４）で表現されていたHMM+RVMの確率モデルは、以下の式（３２）のように表すことができる。
【数２４】
2013187295.tif 000026
【０１０６】
ここで、式（３２）の観測確率P(Xt,Yt｜St)は、
【数２５】
2013187295.tif 000027
と変形できるから、結局、式（３２）は、式（１）として示した通常のHMMの確率モデルとなる。なお、式（２）の観測確率P(Yt｜St)のρStは、線形回帰係数ｗSの第（Dx+１）成分に相当する。
【０１０７】
消費電力変動モデルとして、通常のHMMを採用した場合には、上述したように、線形回帰係数ｗSの事前確率分布P(ｗS)のパラメータは固定化されるので、図７のステップＳ６８の処理は省略される。また、事前確率分布P(ｗS)のパラメータの固定化により、図４のモデルパラメータ学習処理におけるステップＳ２５の処理も省略される。
【０１０８】
したがって、消費電力変動モデルとして通常のHMMを採用した場合には、ログ選択部２３による不要なログ情報の種類の選択制御はできなくなる。換言すれば、HMM+RVMでは、多種類のログ情報を入力Xとして与えた場合でも、情報処理装置１が、電力推定に必要なログ情報のみを（自動で）選択することができる。これにより、人が、消費電力の推定に使用する取得データの種類を判断し、取捨選択する必要はなく、人の負担を軽減することができる。また、一回の学習処理の終了後、次からのデータ収集処理、モデルパラメータ学習処理、および消費電力推定処理においては、不要なログ情報は取得しなくても済むので、取得するログの情報量および計算量も削減され、処理時間も短縮させることができる。すなわち、本技術のHMM+RVMを用いた確率モデルによれば、推定に有用なログ情報のみを使用して、効率の良く推定することができる。
【０１０９】
消費電力変動モデルとして、通常のHMMを採用した場合の消費電力推定処理では、線形回帰係数ｗSの事前確率分布P(ｗS)の分散α-1（αの逆行列）が、式（３１）で構成されることから、消費電力の推定値Y*tを求める式は、HMM+RVMの式（２９）から、以下の式（３４）のように簡略化される。
【数２６】
2013187295.tif 000028
【０１１０】
消費電力変動モデルとして通常のHMMを採用した場合には、モデルパラメータの学習処理において、消費電力時系列データを重視した学習を行うことができる。式（３５）は、学習モデルがHMMである場合の隠れ状態Sの確率分布ｑ（S）の式である。
【数２７】
2013187295.tif 000029
ここで、ωxは、観測されるログ時系列データに対する重み係数であり、ωyは、消費電力時系列データに対する重み係数である。式（３５）において、消費電力時系列データの重み係数ωyをログ時系列データの重み係数ωxより大きく設定することにより、消費電力時系列データを重視した学習を行うことができる。
【０１１１】
上述した実施の形態においては、時刻ｔの観測データXtには、その時刻ｔに得られたログ情報をそのまま用いたが、必要に応じて、ログ情報に対して所定のデータ加工した後の値を、時刻ｔの観測データXtとして扱うようにすることができる。
【０１１２】
例えば、時刻ｔと、それからΔｔ時間前の時刻ｔ-ΔｔまでのBt個のログ情報を、時刻ｔの観測データXtとしてもよい。この場合、観測データXtは、Dx行Bt列のベクトルで構成される。
【０１１３】
また例えば、時刻ｔにおける複数種類（Dx個）のログ情報xt1,xt2,・・・,xtDxのうち、所定の種類のログ情報xti(ｉ＝1，2，・・・，Dxのいずれか)については、関数f(x)＝log(１＋x)を用いて変換した値log(１＋xti)を、入力Xの構成要素としてもよい。
【０１１４】
本実施の形態の学習モデルには、HMMをベースとした学習モデルを採用している。HMMを採用することにより、現在時刻の稼働状態だけでなく、現在時刻の稼働状態になるまでの過去の稼働状態の履歴をも考慮した推定が可能となるので、特許文献１のような線形式で推定する場合より、高精度に消費電力を推定することができる。例えば、ある一定期間CPUに負荷がかかり続けると、それによりCPUの温度が上がったり、ファンが回転するなどして消費電力が増大することがある。本実施の形態の学習モデルによれば、一定期間CPUが高負荷状態であることを履歴として学習した推定結果を出力することができる。
【０１１５】
上述した実施の形態では、本技術の学習モデルであるHMM+RVMを、消費電力時系列データとログ時系列データを観測時系列データとして、消費電力を推定する処理に適用した例について説明した。
【０１１６】
しかし、本技術のHMM+RVMを用いた推定処理は、消費電力以外の推定にも適用することができる。他の適用例について簡単に説明する。
【０１１７】
例えば、本技術は、ロボット等の物体の姿勢を推定する姿勢推定処理に適用することができる。この場合、ロボット等の物体に取り付けられた複数の加速度センサから得られる時系列センサデータと、物体の姿勢を示す位置の時系列データである位置時系列データを、それぞれ、学習処理における入力Xおよび出力Yとする。そして、推定処理では、時系列センサデータを用いて、現在の物体の姿勢を推定することができる。
【０１１８】
本技術を適用した姿勢を推定する学習モデルの学習処理によれば、ログ選択部２３が、不要な加速度センサのセンサデータを除外することができる。現在時刻のセンサデータのみを用いた線形式では、現在の物体の姿勢を判定することは困難であるが、センサデータの時系列データを用いて、加速度の値の累積（積分）を利用することにより、物体の姿勢がどのように変化したかを推定することが可能となる。
【０１１９】
また、本技術は、例えば、ビデオコンテンツの特徴量系列から、人の感じる「うるささ」を推定する処理に適用することができる。この場合、ビデオコンテンツの特徴量（例えば、音量、画像の特徴量など）の時系列データと、人の感じる「うるささ」の時系列データを、それぞれ、学習処理における入力Xおよび出力Yとする。
【０１２０】
人の感じる「うるささ」は、音の種類や前後の文脈に依存するので、単純な音量とは異なる。また、「うるささ」という指標は、ビデオコンテンツの文脈（それまでのビデオの内容とその音声）に依存する。例えば、盛り上がりシーンでは、音量が「３」から「４」に変化したとしても、人はうるさいとは感じないが、音量「１」の静かなシーンで突然音量が「４」になったとすると、人はうるさいと感じる場合がある。したがって、現在時刻の特徴量だけでは、正確な「うるささ」の推定が困難であり、ビデオコンテンツの特徴量の時系列データを入力Xとすることで、人の感じる「うるささ」を、高精度に推定することが可能となる。
【０１２１】
そして、本技術を適用した「うるささ」を推定する学習モデルの学習処理によれば、ログ選択部２３が、取得しているビデオコンテンツの特徴量のうち、不要な特徴量を判定し、使用しないように制御することができる。
【０１２２】
さらに、本技術は、例えば、ユーザのスマートフォンの操作状況を示す操作時系列データから、そのユーザのその日のテレビの視聴時間を推定する処理に適用することができる。この場合、時間インデックスｔを日にち（日単位）とし、ユーザのスマートフォンの操作状況を示す操作時系列データと、そのユーザのテレビの視聴時間の時系列データを、それぞれ、学習処理における入力Xおよび出力Yとする。そして、推定処理では、ある日のマートフォンの操作状況を示す操作時系列データを用いて、その日のテレビの視聴時間を推定することができる。
【０１２３】
人間の行動パターンには、「ちょうど一週間前にテレビを沢山見ていた場合には、その日も沢山見る可能性が高い」など、履歴依存性がある。したがって、一日のユーザの行動状況だけで推定するよりも、数日間、数週間といった所定期間の時系列データを用いることで、高精度に推定することができる。なお、推定する項目は、「テレビの視聴時間」以外に、例えば、「車に乗っている時間」や、「フェイスブック（登録商標）等のSNS(Social Networking Service)にログインしている時間」などでもよい。
【０１２４】
上述した一連の処理は、ハードウエアにより実行することもできるし、ソフトウエアにより実行することもできる。一連の処理をソフトウエアにより実行する場合には、そのソフトウエアを構成するプログラムが、コンピュータにインストールされる。ここで、コンピュータには、専用のハードウエアに組み込まれているコンピュータや、各種のプログラムをインストールすることで、各種の機能を実行することが可能な、例えば汎用のパーソナルコンピュータなどが含まれる。
【０１２５】
図８は、上述した一連の処理をプログラムにより実行するコンピュータのハードウエアの構成例を示すブロック図である。
【０１２６】
コンピュータにおいて、CPU（Central Processing Unit）１０１，ROM（Read Only Memory）１０２，RAM（Random Access Memory）１０３は、バス１０４により相互に接続されている。
【０１２７】
バス１０４には、さらに、入出力インタフェース１０５が接続されている。入出力インタフェース１０５には、入力部１０６、出力部１０７、記憶部１０８、通信部１０９、及びドライブ１１０が接続されている。
【０１２８】
入力部１０６は、キーボード、マウス、マイクロホンなどよりなる。出力部１０７は、ディスプレイ、スピーカなどよりなる。記憶部１０８は、ハードディスクや不揮発性のメモリなどよりなる。通信部１０９は、ネットワークインタフェースなどよりなる。ドライブ１１０は、磁気ディスク、光ディスク、光磁気ディスク、或いは半導体メモリなどのリムーバブル記録媒体１１１を駆動する。
【０１２９】
以上のように構成されるコンピュータでは、CPU１０１が、例えば、記憶部１０８に記憶されているプログラムを、入出力インタフェース１０５及びバス１０４を介して、RAM１０３にロードして実行することにより、上述した一連の処理が行われる。
【０１３０】
コンピュータでは、プログラムは、リムーバブル記録媒体１１１をドライブ１１０に装着することにより、入出力インタフェース１０５を介して、記憶部１０８にインストールすることができる。また、プログラムは、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線または無線の伝送媒体を介して、通信部１０９で受信し、記憶部１０８にインストールすることができる。その他、プログラムは、ROM１０２や記憶部１０８に、あらかじめインストールしておくことができる。
【０１３１】
なお、コンピュータが実行するプログラムは、本明細書で説明する順序に沿って時系列に処理が行われるプログラムであっても良いし、並列に、あるいは呼び出しが行われたとき等の必要なタイミングで処理が行われるプログラムであっても良い。
【０１３２】
なお、本明細書において、フローチャートに記述されたステップは、記載された順序に沿って時系列的に行われる場合はもちろん、必ずしも時系列的に処理されなくとも、並列に、あるいは呼び出しが行われたとき等の必要なタイミングで実行されてもよい。
【０１３３】
なお、本明細書において、システムとは、複数の構成要素（装置、モジュール（部品）等）の集合を意味し、すべての構成要素が同一筐体中にあるか否かは問わない。したがって、別個の筐体に収納され、ネットワークを介して接続されている複数の装置、及び、１つの筐体の中に複数のモジュールが収納されている１つの装置は、いずれも、システムである。
【０１３４】
本技術の実施の形態は、上述した実施の形態に限定されるものではなく、本技術の要旨を逸脱しない範囲において種々の変更が可能である。
【０１３５】
例えば、上述した複数の実施の形態の全てまたは一部を組み合わせた形態を採用することができる。
【０１３６】
例えば、本技術は、１つの機能をネットワークを介して複数の装置で分担、共同して処理するクラウドコンピューティングの構成をとることができる。
【０１３７】
また、上述のフローチャートで説明した各ステップは、１つの装置で実行する他、複数の装置で分担して実行することができる。
【０１３８】
さらに、１つのステップに複数の処理が含まれる場合には、その１つのステップに含まれる複数の処理は、１つの装置で実行する他、複数の装置で分担して実行することができる。
【０１３９】
なお、本技術は以下のような構成も取ることができる。
（１）
推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得する取得部と、
取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータを学習する学習部と、
学習により得られた前記確率モデルのパラメータに基づいて、前記取得部が取得する前記説明時系列データに対応する前記説明変数を選択する選択部と、
前記選択部の選択結果に基づいて前記取得部が取得した複数の前記説明時系列データを用いて、前記目的変数の値を推定する推定部と
を備える情報処理装置。
（２）
前記学習部は、前記目的変数と前記複数の説明変数との関係を隠れマルコフモデルで学習する
前記（１）に記載の情報処理装置。
（３）
前記目的変数が、前記隠れマルコフモデルの隠れ状態と１対１に対応する線形回帰係数と、前記説明変数との線形回帰モデルで表される
前記（２）に記載の情報処理装置。
（４）
前記選択部は、前記線形回帰係数が所定の閾値より小さい前記説明変数を、前記取得部が時系列データを取得しない説明変数として選択する
前記（３）に記載の情報処理装置。
（５）
情報処理装置が、
推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得し、
取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータを学習し、
学習により得られた前記確率モデルのパラメータに基づいて、取得する前記説明時系列データに対応する前記説明変数を選択し、
選択結果に基づいて取得した複数の前記説明時系列データを用いて、前記目的変数の値を推定する
ステップを含む情報処理方法。
（６）
コンピュータを、
推定対象の目的変数に対応する時系列データである目的時系列データと、前記目的変数を説明する複数の説明変数に対応する時系列データである複数の説明時系列データとを取得する取得部と、
取得された前記目的時系列データと前記複数の説明時系列データとを用いて、確率モデルのパラメータを学習する学習部と、
学習により得られた前記確率モデルのパラメータに基づいて、前記取得部が取得する前記説明時系列データに対応する前記説明変数を選択する選択部と、
前記選択部の選択結果に基づいて前記取得部が取得した複数の前記説明時系列データを用いて、前記目的変数の値を推定する推定部
として機能させるためのプログラム。
【符号の説明】
【０１４０】
１  情報処理装置，  １１  消費電力計測部，  １２  消費電力時系列入力部，  １３  ログ取得部，  １５  ログ時系列入力部，  １６  時系列履歴保存部，  １７  モデル学習部，  １８  消費電力推定部，  １９  推定消費電力表示部，  ２１  モデルパラメータ更新部，  ２２  モデルパラメータ記憶部，  ２３  ログ選択部
【図１】
2013187295.tif 000030
【図２】
2013187295.tif 000031
【図３】
2013187295.tif 000032
【図４】
2013187295.tif 000033
【図５】
2013187295.tif 000034
【図６】
2013187295.tif 000035
【図７】
2013187295.tif 000036
【図８】
2013187295.tif 000037
【国際調査報告】
2013187295.tif 000038
2013187295.tif 000039
2013187295.tif 000040
2013187295.tif 000041