(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-11-05
(45)【発行日】2024-11-13
(54)【発明の名称】データ算出装置、データ算出方法およびプログラム
(51)【国際特許分類】
G06N 20/00 20190101AFI20241106BHJP
【ＦＩ】
G06N20/00
【請求項の数】
6
(21)【出願番号】P 2022550307
(86)(22)【出願日】2020-09-18
(86)【国際出願番号】 JP2020035546
(87)【国際公開番号】W WO2022059189
(87)【国際公開日】2022-03-24
【審査請求日】2023-02-22
(73)【特許権者】
【識別番号】000004237
【氏名又は名称】日本電気株式会社
(74)【代理人】
【識別番号】100149548
【弁理士】
【氏名又は名称】松沼 泰史
(74)【代理人】
【識別番号】100181135
【弁理士】
【氏名又は名称】橋本 隆史
(72)【発明者】
【氏名】木佐森 慶一
(72)【発明者】
【氏名】清川 裕
【審査官】渡辺 順哉
(56)【参考文献】
【文献】国際公開第２０２０／１０５４６８（ＷＯ，Ａ１）
【文献】国際公開第２０１６／０３１１７４（ＷＯ，Ａ１）
【文献】吉本潤一郎 ほか，変分法的ベイズ推定法に基づく正規化ガウス関数ネットワークと階層的モデル選択法，計測自動制御学会論文集，日本，社団法人計測自動制御学会，2003年05月31日，Vol.39, No.5，pp.503-512，ISSN 0453-4654
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０６Ｎ ３／００－９９／００
(57)【特許請求の範囲】
【請求項1】
対象の部分を模擬して、前記対象を模擬する上位モデルのパラメータ値を出力する下位モデルについて、前記上位モデルのパラメータ値の事後分布を示すデータに基づく、前記下位モデルのパラメータ値の事後分布を示すデータを算出する下位事後分布データ算出手段と、前記下位モデルのパラメータ値の事後分布を示すデータに基づくパラメータ値を前記下位モデルに適用してシミュレーションを実行する下位モデルシミュレーション実行手段と、
前記下位モデルを用いたシミュレーションの結果に基づく前記上位モデルのパラメータ値のサンプルデータを前記上位モデルに適用してシミュレーションを実行する上位モデルシミュレーション実行手段と、
前記上位モデルシミュレーション実行手段による前記シミュレーションの結果として得られる前記上位モデルの出力のサンプルデータと、前記上位モデルの出力に対するターゲットデータとの類似度に基づいて、前記上位モデルの出力に対するターゲットデータに基づく、前記上位モデルのパラメータ値の事後分布を示すデータを算出する上位事後分布データ算出手段と、
を備えるデータ算出装置。
【請求項2】
前記上位事後分布データ算出手段は、前記上位モデルのパラメータ値のサンプルデータのカーネル平均に基づいてカーネルハーディングを実行し、前記上位モデルのパラメータ値の事後分布を示すデータを算出する、請求項１に記載のデータ算出装置。
【請求項3】
前記下位事後分布データ算出手段は、前記下位モデルのパラメータ値のサンプルデータのカーネル平均に基づいてカーネルハーディングを実行し、前記下位モデルのパラメータ値の事後分布を示すデータを算出する、請求項１または請求項２に記載のデータ算出装置。
【請求項4】
前記下位事後分布データ算出手段が算出する前記下位モデルのパラメータ値の事後分布を示すデータに基づいて、前記対象を制御する制御手段をさらに備える、請求項１から３の何れか一項に記載のデータ算出装置。
【請求項5】
コンピュータが、対象の部分を模擬して、前記対象を模擬する上位モデルのパラメータ値を出力する下位モデルについて、前記上位モデルのパラメータ値の事後分布を示すデータに基づく、前記下位モデルのパラメータ値の事後分布を示すデータを算出し、
前記下位モデルのパラメータ値の事後分布を示すデータに基づくパラメータ値を前記下位モデルに適用してシミュレーションを実行し、
前記下位モデルを用いたシミュレーションの結果に基づく前記上位モデルのパラメータ値のサンプルデータを前記上位モデルに適用してシミュレーションを実行し、
前記上位モデルを用いたシミュレーションの結果として得られる前記上位モデルの出力のサンプルデータと、前記上位モデルの出力に対するターゲットデータとの類似度に基づいて、前記上位モデルの出力に対するターゲットデータに基づく、前記上位モデルのパラメータ値の事後分布を示すデータを算出する、
データ算出方法。
【請求項6】
コンピュータに、対象の部分を模擬して、前記対象を模擬する上位モデルのパラメータ値を出力する下位モデルについて、前記上位モデルのパラメータ値の事後分布を示すデータに基づく、前記下位モデルのパラメータ値の事後分布を示すデータを算出することと、
前記下位モデルのパラメータ値の事後分布を示すデータに基づくパラメータ値を前記下位モデルに適用してシミュレーションを実行することと、
前記下位モデルを用いたシミュレーションの結果に基づく前記上位モデルのパラメータ値のサンプルデータを前記上位モデルに適用してシミュレーションを実行することと、
前記上位モデルを用いたシミュレーションの結果として得られる前記上位モデルの出力のサンプルデータと、前記上位モデルの出力に対するターゲットデータとの類似度に基づいて、前記上位モデルの出力に対するターゲットデータに基づく、前記上位モデルのパラメータ値の事後分布を示すデータを算出することと、
を実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、データ算出装置、データ算出方法およびプログラムに関する。
【背景技術】
【0002】
シミュレータの出力相当するターゲットデータに基づいて、そのシミュレータがターゲットデータを出力するようなパラメータ値を推定する技術が提案されている。例えば、特許文献１には、カーネルＡＢＣ（Kernel ABC; Kernel Approximate Bayesian Computation）およびカーネルハーディング（Kernel Herding）を用いて、ターゲットデータに応じたパラメータ値の事後分布を算出することが記載されている。
【先行技術文献】
【特許文献】
【0003】
【文献】国際公開第２０２０－１０５４６８号
【発明の概要】
【発明が解決しようとする課題】
【0004】
シミュレーションの対象の規模が大きい場合など、値を推定したいパラメータの個数が多い場合、パラメータ値の推定に時間がかかることが考えられる。特に、値を推定したいパラメータの個数の増加に応じて、パラメータ値の推定に要する時間が指数的に増加することが考えられる。このため、値を推定したいパラメータの個数が多い場合、現実的な時間内にパラメータの推定値を得られないことが考えられる。
【0005】
本発明の目的の１つは、上述の課題を解決することのできるデータ算出装置、データ算出方法およびプログラムを提供することである。
【課題を解決するための手段】
【0006】
本発明の第１の態様によれば、データ算出装置は、対象の部分を模擬して、前記対象を模擬する上位モデルのパラメータ値を出力する下位モデルについて、前記上位モデルのパラメータ値の事後分布を示すデータに基づく、前記下位モデルのパラメータ値の事後分布を示すデータを算出する下位事後分布データ算出手段と、前記下位モデルのパラメータ値の事後分布を示すデータに基づくパラメータ値を前記下位モデルに適用してシミュレーションを実行する下位モデルシミュレーション実行手段と、前記下位モデルを用いたシミュレーションの結果に基づく前記上位モデルのパラメータ値のサンプルデータを前記上位モデルに適用してシミュレーションを実行する上位モデルシミュレーション実行手段と、前記上位モデルシミュレーション実行手段による前記シミュレーションの結果として得られる前記上位モデルの出力のサンプルデータと、前記上位モデルの出力に対するターゲットデータとの類似度に基づいて、前記上位モデルの出力に対するターゲットデータに基づく、前記上位モデルのパラメータ値の事後分布を示すデータを算出する上位事後分布データ算出手段と、を備える
【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】
以下、本発明の実施形態を説明するが、以下の実施形態は請求の範囲にかかる発明を限定するものではない。また、実施形態の中で説明されている特徴の組み合わせの全てが発明の解決手段に必須であるとは限らない。
図１は、実施形態に係るデータ算出装置の機能構成を示す概略ブロック図である。図１に示す構成で、データ算出装置１００は、通信部１１０と、表示部１２０と、操作入力部１３０と、記憶部１７０と、制御部１８０とを備える。制御部１８０は、初期設定部１８１と、ループ処理部１８２と、上位シミュレータ部１８３と、上位事後分布データ算出部１８４と、下位シミュレータ部１８５と、下位事後分布データ算出部１８６とを備える。
また、図１には、シミュレーションモデルによる模擬（シミュレーション、Simulation）の対象である対象９００が示されている。対象９００のシミュレーションモデルを、対象９００のモデル、あるいは単にモデルとも称する。
図１は、実施形態に係るデータ算出装置の機能構成を示す概略ブロック図である。図１に示す構成で、データ算出装置１００は、通信部１１０と、表示部１２０と、操作入力部１３０と、記憶部１７０と、制御部１８０とを備える。制御部１８０は、初期設定部１８１と、ループ処理部１８２と、上位シミュレータ部１８３と、上位事後分布データ算出部１８４と、下位シミュレータ部１８５と、下位事後分布データ算出部１８６とを備える。
また、図１には、シミュレーションモデルによる模擬（シミュレーション、Simulation）の対象である対象９００が示されている。対象９００のシミュレーションモデルを、対象９００のモデル、あるいは単にモデルとも称する。
【0013】
データ算出装置１００は、モデルの出力に対するターゲットデータを用いて、モデルの出力値がターゲットデータの値に近付くように、モデルのパラメータ値を推定する。データ算出装置１００が、例えばワークステーション（Workstation）またはパソコン（Personal Computer；ＰＣ）等のコンピュータを用いて構成されていてもよい。
【0014】
データ算出装置１００が、１つの装置として構成されていてもよいし、複数の装置の組み合わせにて構成されていてもよい。例えば、上位シミュレータ部１８３、および、下位シミュレータ部１８５が、モデルごとに別々のシミュレータ装置として構成されていてもよい。
ターゲットデータの値をターゲット値とも称する。モデルの出力値がターゲット値に近付くように、モデルのパラメータ値を推定する処理をデータ同化とも称する。データ同化におけるパラメータ値の推定をパラメータ値の調整とも称する。
ターゲットデータの値をターゲット値とも称する。モデルの出力値がターゲット値に近付くように、モデルのパラメータ値を推定する処理をデータ同化とも称する。データ同化におけるパラメータ値の推定をパラメータ値の調整とも称する。
【0015】
ターゲットデータは、対象９００における実測データであってもよい。この場合、データ算出装置１００がパラメータ値を推定することで、対象９００を高精度に模擬できるようにモデルを調整することができる。
また、推定結果のパラメータ値を参照することで、対象９００の状態を推定することができる。例えば、対象９００が工場であり、データ算出装置１００による推定対象のパラメータに、工場の生産ラインごとの生産量である場合、データ算出装置１００が推定するパラメータ値を、実測データから推定される、工場の生産ラインごとの生産量として用いることができる。
また、推定結果のパラメータ値を参照することで、対象９００の状態を推定することができる。例えば、対象９００が工場であり、データ算出装置１００による推定対象のパラメータに、工場の生産ラインごとの生産量である場合、データ算出装置１００が推定するパラメータ値を、実測データから推定される、工場の生産ラインごとの生産量として用いることができる。
【0016】
あるいは、ターゲットデータは、対象９００における目標値を示すデータであってもよい。この場合、データ算出装置１００による推定結果のパラメータを参照することで、目標を達成するための方策を作成することができる。さらには、データ算出装置１００によれば、作成した方策に基づいて工場内の各装置を制御するなど、対象９００を制御することができる。
【0017】
例えば、対象９００が工場であり、データ算出装置１００による推定対象のパラメータに、工場の生産ラインごとの生産量である場合、データ算出装置１００が推定するパラメータ値を、ターゲットデータに示される目標値を達成するための、工場の生産ラインごとの生産量として用いることができる。
データ算出装置１００によれば、工場の生産ラインごとに推定した生産量に基づいて生産ラインごとに機器を制御するなど、推定結果に基づく制御を行うことができる。
データ算出装置１００によれば、工場の生産ラインごとに推定した生産量に基づいて生産ラインごとに機器を制御するなど、推定結果に基づく制御を行うことができる。
【0018】
データ算出装置１００は、対象９００を模擬するモデルとして階層的なモデルを用いる。例えば、データ算出装置１００が、対象９００の全体を模擬する上位モデルと、対象９００の一部を模擬する１つ以上の下位モデルとを含むモデルを用いるようにしてもよい。
【0019】
ここでいう「階層」は、上下関係を定めることができることである。例えば上記のように、全体と部分との関係に基づいて、対象９００の全体を模擬するモデルを上位とし、対象９００の部分を模擬するモデルを下位とすることができる。この上下関係は固定的であってもよいし、可変であってもよい。例えば、会社を模擬するモデルで、経費の流れに関しては経理部のモデルが上位となるが、資材の流れに関しては購買部のモデルが上位になるというように、データ算出の目的に応じて上下関係が異なってもよい。
【0020】
さらに、データ算出装置１００は、下位モデルが上位モデルのパラメータ値を出力する、階層的なモデルを用いる。これにより、データ算出装置１００は、モデルの出力のターゲットデータに基づいてモデルのパラメータ値を推定する処理を、モデルの階層構造に従って繰り返す。
【0021】
例えば、データ算出装置１００は、上位モデルの出力のターゲットデータに基づいて上位モデルのパラメータ値を推定する。そして、データ算出装置１００は、推定した上位モデルのパラメータ値を下位モデルの出力のターゲットデータとして用いて、下位モデルのパラメータ値を推定する。
【0022】
このように、データ算出装置１００がモデルの階層構造に応じてモデルのパラメータ値の推定を階層的に行うことで、対象９００を１つのモデルのみで模擬する場合よりも、個々のモデルのパラメータの個数を少なくすることができる。例えば、上位モデルのパラメータが、複数の生産ラインを含む部門ごとの生産量をパラメータとして用いることで、個々の生産ラインの生産量をパラメータとして用いる場合よりも、上位モデルのパラメータの個数が少なくて済む。
【0023】
ここで、モデルのパラメータの個数が多い場合、パラメータ値の推定に時間がかかることが考えられる。特に、パラメータの個数の増加に応じて、パラメータ値の推定に要する時間が指数的に増加することが考えられる。このため、パラメータの個数が多い場合、現実的な時間内にパラメータの推定値を得られないことが考えられる。
これに対して、データ算出装置１００では、上記のように個々のモデルのパラメータの個数が比較的少ない階層的なモデルを用いることで、パラメータの個数が多い１つのモデルを用いる場合よりも計算時間が短くて済むと期待される。
これに対して、データ算出装置１００では、上記のように個々のモデルのパラメータの個数が比較的少ない階層的なモデルを用いることで、パラメータの個数が多い１つのモデルを用いる場合よりも計算時間が短くて済むと期待される。
【0024】
対象９００は、データ算出装置１００による模擬の対象である。階層的なシミュレーションモデルを用いて模擬可能ないろいろなものを対象９００とすることができる。
例えば、対象９００は工場であってもよい。そして、データ算出装置１００が、工場全体を模擬する上位のモデルと、工場の一部を模擬する１つ以上の下位のモデルとを含む階層的なモデルを用いるようにしてもよい。
あるいは、対象９００が、生産拠点、物流設備、および、販売拠点を含むサプライチェーンであってもよい。そして、データ算出装置１００が、サプライチェーン全体を模擬する上位のモデルと、生産拠点、物流設備、販売拠点のそれぞれを模擬する下位のモデルとを含む階層的なモデルを用いるようにしてもよい。
例えば、対象９００は工場であってもよい。そして、データ算出装置１００が、工場全体を模擬する上位のモデルと、工場の一部を模擬する１つ以上の下位のモデルとを含む階層的なモデルを用いるようにしてもよい。
あるいは、対象９００が、生産拠点、物流設備、および、販売拠点を含むサプライチェーンであってもよい。そして、データ算出装置１００が、サプライチェーン全体を模擬する上位のモデルと、生産拠点、物流設備、販売拠点のそれぞれを模擬する下位のモデルとを含む階層的なモデルを用いるようにしてもよい。
【0025】
あるいは、対象９００が、ある地域であってもよい。そして、データ算出装置１００が、地域全体について感染の広がりの状況を模擬する上位のモデルと、その地域の一部または交通機関での感染の状況を模擬する下位のモデルとを含む階層的なモデルを用いるようにしてもよい。
対象９００は、実在していてもよいし、仮想的なものであってもよい。
対象９００は、実在していてもよいし、仮想的なものであってもよい。
【0026】
図２は、データ算出装置１００が用いる対象９００のモデルの例を示す図である。図２の例で、モデル２００は、１つの上位モデル２１０と、Ｓ個の下位モデル２２０－１から２２０－Ｓまでとを含む。Ｓは、Ｓ≧１の整数であり、下位モデル２２０の個数を示す。下位モデル２２０－１から２２０－Ｓまでを総称して下位モデル２２０とも表記する。
【0027】
上位モデル２１０への入力をベクトルＸで示し、上位モデル２１０のパラメータをベクトルＺで示し、上位モデル２１０の出力をベクトルＹで示す。ベクトルＸ、Ｙ、Ｚは何れも実数ベクトルである。ベクトルＸ、Ｙの何れも次元数は１次元以上であればよい。ベクトルＺの次元数はＳ次元以上であればよい。
ベクトルＸを入力ベクトルＸとも称する。ベクトルＹを出力ベクトルＹとも称する。ベクトルＺをパラメータベクトルＺとも称する。
ベクトルＸを入力ベクトルＸとも称する。ベクトルＹを出力ベクトルＹとも称する。ベクトルＺをパラメータベクトルＺとも称する。
【0028】
入力ベクトルＸの要素に、整数値または二値などの離散値をとる要素が含まれていてもよい。同様に、出力ベクトルＹの要素に、整数値または二値などの離散値をとる要素が含まれていてもよい。パラメータベクトルＺの要素に、整数値または二値などの離散値をとる要素が含まれていてもよい。
上位モデル２１０をｆで示し、ｆ（Ｘ，Ｚ）＝Ｙと表記する。あるいは、入力ベクトルＸの表記を省略してｆ（Ｚ）＝Ｙとも表記する。
上位モデル２１０をｆで示し、ｆ（Ｘ，Ｚ）＝Ｙと表記する。あるいは、入力ベクトルＸの表記を省略してｆ（Ｚ）＝Ｙとも表記する。
【0029】
下位モデル２２０－ｓへの入力をベクトルＸｓで示し、下位モデル２２０－ｓのパラメータをベクトルθｓで示し、下位モデル２２０－ｓの出力をベクトルＺｓで示す。ｓは、１≦ｓ≦Ｓの整数である。
ベクトルＸｓ、Ｚｓ、θｓは何れも実数ベクトルである。Ｘｓ、ベクトルＺｓ、θｓの何れも次元数は１次元以上であればよい。ベクトルＸｓの次元数は、０次元以上であればよい。
ベクトルＸｓ、Ｚｓ、θｓは何れも実数ベクトルである。Ｘｓ、ベクトルＺｓ、θｓの何れも次元数は１次元以上であればよい。ベクトルＸｓの次元数は、０次元以上であればよい。
【0030】
入力ベクトルＸｓの要素に、整数値または二値などの離散値をとる要素が含まれていてもよい。同様に、出力ベクトルＺｓの要素に、整数値または二値などの離散値をとる要素が含まれていてもよい。パラメータベクトルθｓの要素に、整数値または二値などの離散値をとる要素が含まれていてもよい。
下位モデル２２０をｇｓで示し、ｇｓ（Ｘｓ，θｓ）＝Ｚｓと表記する。あるいは、入力ベクトルＸｓの表記を省略し、ｇｓ（θｓ）＝Ｚｓとも表記する。
下位モデル２２０をｇｓで示し、ｇｓ（Ｘｓ，θｓ）＝Ｚｓと表記する。あるいは、入力ベクトルＸｓの表記を省略し、ｇｓ（θｓ）＝Ｚｓとも表記する。
【0031】
データ算出装置１００は、入力ベクトルＸの値と出力ベクトルＹの値との組み合わせを示すターゲットデータを取得し、上位モデル２１０の出力値がターゲット値に近付くようにパラメータベクトルＺの値を調整する。
データ算出装置１００による調整前のパラメータベクトルＺの確率分布を事前分布πＺと表記し、事前分布πＺによるパラメータベクトルＺの値をサンプルＺ＜ＰＲＩ＞と表記する。ベクトルの確率分布は、ベクトルの要素ごとの確率分布の組み合わせで示されてもよい。
データ算出装置１００による調整前のパラメータベクトルＺの確率分布を事前分布πＺと表記し、事前分布πＺによるパラメータベクトルＺの値をサンプルＺ＜ＰＲＩ＞と表記する。ベクトルの確率分布は、ベクトルの要素ごとの確率分布の組み合わせで示されてもよい。
【0032】
事前分布πＺは、ユーザが対象９００の情報に基づいて設定するようにしてもよい。あるいは、事前分布πＺとして一様分布など所定の分布を用いるようにしてもよい。
データ算出装置１００による調整後のパラメータベクトルＺの確率分布を事後分布ｐ（Ｚ｜Ｙ）と表記し、事後分布ｐ（Ｚ｜Ｙ）によるパラメータベクトルＺの値をサンプルＺ＜ＰＯＳＴ＞と表記する。ここでの事後分布ｐ（Ｚ｜Ｙ）は、出力ベクトルＹのターゲット値が与えられたときの、パラメータベクトルＺの値の事後分布を意味する。
データ算出装置１００による調整後のパラメータベクトルＺの確率分布を事後分布ｐ（Ｚ｜Ｙ）と表記し、事後分布ｐ（Ｚ｜Ｙ）によるパラメータベクトルＺの値をサンプルＺ＜ＰＯＳＴ＞と表記する。ここでの事後分布ｐ（Ｚ｜Ｙ）は、出力ベクトルＹのターゲット値が与えられたときの、パラメータベクトルＺの値の事後分布を意味する。
【0033】
また、データ算出装置１００は、上位モデル２１０への入力値（入力ベクトルＸの値）に基づいて下位モデル２２０－ｓへの入力値（入力ベクトルＸｓの値）を決定する。例えば、入力ベクトルＸｓのベクトル空間が、入力ベクトルＸのベクトル空間の部分空間になっていてもよい。すなわち、入力ベクトルＸｓが入力ベクトルＸの要素のうち一部の要素を用いて構成されていてもよい。この場合、データ算出装置１００は、入力ベクトルＸと入力ベクトルＸｓとの関係に応じて、入力ベクトルＸの要素の値を入力ベクトルＸｓの要素の値に設定する。
あるいは、入力ベクトルＸｓの要素の値に、入力ベクトルＸの値に基づく計算値が含まれていてもよい。
あるいは、入力ベクトルＸｓの要素の値に、入力ベクトルＸの値に基づく計算値が含まれていてもよい。
【0034】
また、データ算出装置１００は、サンプルＺ＜ＰＯＳＴ＞に基づいて下位モデル２２０－ｓの出力のターゲット値を決定する。例えば、出力ベクトルＺｓのベクトル空間が、パラメータベクトルＺのベクトル空間の部分空間になっていてもよい。すなわち、出力ベクトルＺｓがパラメータベクトルＺの要素のうち一部の要素を用いて構成されていてもよい。この場合、データ算出装置１００は、パラメータベクトルＺと出力ベクトルＺｓとの関係に応じて、サンプルＺ＜ＰＯＳＴ＞の要素の値を下位モデル２２０－ｓの出力のターゲットデータにおける出力ベクトルＺｓの要素の値に設定する。
【0035】
パラメータベクトルＺの要素のうち、出力ベクトルＺｓに紐付けられない要素があってもよい。すなわち、上位モデル２１０のパラメータとして、上位モデル２１０での計算に用いられ、下位モデル２２０には値が出力されないパラメータがあってもよい。
データ算出装置１００がサンプルＺ＜ＰＯＳＴ＞に基づいて決定する下位モデル２２０－ｓの出力のターゲット値をサンプルＺ＜ＰＯＳＴ＞ ｓと表記する。
データ算出装置１００がサンプルＺ＜ＰＯＳＴ＞に基づいて決定する下位モデル２２０－ｓの出力のターゲット値をサンプルＺ＜ＰＯＳＴ＞ ｓと表記する。
【0036】
このように、データ算出装置１００が、上位モデル２１０のパラメータ値を下位モデル２２０の出力値のターゲットデータに設定することで、データ同化を階層的に行うことができる。具体的には、データ算出装置１００が、上位モデル２１０でのパラメータ値の調整後のサンプルＺ＜ＰＯＳＴ＞に基づいて下位モデル２２０の出力のターゲット値（サンプルＺ＜ＰＯＳＴ＞
ｓ）を決定することで、上位モデル２１０での調整結果を下位モデル２２０でのパラメータ値の調整に反映させることができる。
【0037】
パラメータベクトルＺの次元数が、パラメータベクトルθ１からθＳまでの次元数の合計よりも小さいことで、対象９００を１つのモデルのみで模擬する場合よりも、個々のモデルのパラメータの個数を少なくすることができる。この点で、データ算出装置１００によれば比較的短時間でデータ同化を行うことができ、パラメータベクトルθ１からθＳまでの値を推定できると期待される。特に、データ算出装置１００によれば、モデルのパラメータの個数が多いことで現実的な時間内にパラメータの推定値を得られないことを回避できると期待される。
【0038】
また、データ算出装置１００は、下位モデル２２０のそれぞれについて、出力ベクトルＺｓの値がサンプルＺ＜ＰＯＳＴ＞
ｓに近付くようにパラメータベクトルθｓの値を調整するデータ同化を行う。
データ算出装置１００による調整前のパラメータベクトルθｓの確率分布を事前分布πθｓと表記し、事前分布πθｓによるパラメータベクトルθｓの値をサンプルθ＜ＰＲＩ＞ ｓと表記する。
データ算出装置１００による調整前のパラメータベクトルθｓの確率分布を事前分布πθｓと表記し、事前分布πθｓによるパラメータベクトルθｓの値をサンプルθ＜ＰＲＩ＞ ｓと表記する。
【0039】
事前分布πθｓは、ユーザが対象９００の情報に基づいて設定するようにしてもよい。あるいは、事前分布πθｓとして一様分布など所定の分布を用いるようにしてもよい。
データ算出装置１００による調整後のパラメータベクトルＺの確率分布を事後分布ｐ（θｓ｜Ｚｓ）と表記し、事後分布ｐ（θｓ｜Ｚｓ）によるパラメータベクトルＺの値をサンプルθ＜ＰＯＳＴ＞ ｓと表記する。
データ算出装置１００による調整後のパラメータベクトルＺの確率分布を事後分布ｐ（θｓ｜Ｚｓ）と表記し、事後分布ｐ（θｓ｜Ｚｓ）によるパラメータベクトルＺの値をサンプルθ＜ＰＯＳＴ＞ ｓと表記する。
【0040】
また、データ算出装置１００は、サンプルθ＜ＰＯＳＴ＞
ｓを下位モデル２２０－ｓのパラメータに適用してシミュレーションを実行し、入力ベクトルＸｓに対する出力ベクトルＺｓを算出する。パラメータのサンプルをモデルのパラメータに適用するとは、サンプル値をパラメータ値として設定することである。サンプルをモデルのパラメータに適用することを、サンプルをモデルに適用するとも称する。サンプルθ＜ＰＯＳＴ＞
ｓを適用したシミュレーションで算出される出力ベクトルＺｓの値をサンプルＺ＜ＰＲＥＤ＞
ｓと表記する。
【0041】
そして、データ算出装置１００は、サンプルＺ＜ＰＲＥＤ＞
１からサンプルＺ＜ＰＲＥＤ＞
Ｓまでの値に基づいて、パラメータベクトルＺの値のサンプルを生成する。このサンプルは、データ算出装置１００が、上位モデル２１０を用いたデータ同化を繰り返し行う際に、下位モデル２２０でのデータ同化の結果を反映したサンプルとして用いられる。
【0042】
パラメータベクトルＺの何れの要素も、サンプルＺ＜ＰＲＥＤ＞
１からサンプルＺ＜ＰＲＥＤ＞
Ｓまでの要素のうち高々１つの要素と紐付けられている場合、この関係に応じて、データ算出装置１００が、サンプルＺ＜ＰＲＥＤ＞
１からサンプルＺ＜ＰＲＥＤ＞
Ｓまで要素の値を、生成するサンプルの要素の値として用いるようにしてもよい。
【0043】
あるいは、パラメータベクトルＺの要素のうち、サンプルＺ＜ＰＲＥＤ＞
１からサンプルＺ＜ＰＲＥＤ＞
Ｓまでの複数の要素と紐付けられる要素があってもよい。この場合、データ算出装置１００が、サンプルＺ＜ＰＲＥＤ＞
１からサンプルＺ＜ＰＲＥＤ＞
Ｓまでの複数の要素の値の平均をとるなど所定の方法で、生成するサンプルの要素の値を算出するようにしてもよい。
【0044】
このように、データ算出装置１００が、サンプルＺ＜ＰＲＥＤ＞
１からサンプルＺ＜ＰＲＥＤ＞
Ｓまでの値に基づいて、パラメータベクトルＺのサンプルを生成することで、階層的なデータ同化を繰り返すことができる。
具体的には、データ算出装置１００は、生成したサンプルの値を上位モデル２１０のパラメータベクトルＺに適用してシミュレーションを実行する。これにより、データ算出装置１００は、生成したパラメータのサンプルに基づく出力ベクトルＹのサンプルを取得する。そして、データ算出装置１００は、得られた出力ベクトルＹのサンプルとターゲット値との関係に基づいて、上位モデル２１０の出力値がターゲット値に近付くように、パラメータベクトルＺの値を推定する。データ算出装置１００は、推定したパラメータベクトルＺの値をサンプルＺ＜ＰＯＳＴ＞として用いて、上述したのと同様に、下位モデル２２０におけるデータ同化を行い、データ同化の結果に基づいてパラメータベクトルＺのサンプルを生成する。データ算出装置１００は、所定の終了条件が成立するまで、階層的なデータ同化を繰り返す。
具体的には、データ算出装置１００は、生成したサンプルの値を上位モデル２１０のパラメータベクトルＺに適用してシミュレーションを実行する。これにより、データ算出装置１００は、生成したパラメータのサンプルに基づく出力ベクトルＹのサンプルを取得する。そして、データ算出装置１００は、得られた出力ベクトルＹのサンプルとターゲット値との関係に基づいて、上位モデル２１０の出力値がターゲット値に近付くように、パラメータベクトルＺの値を推定する。データ算出装置１００は、推定したパラメータベクトルＺの値をサンプルＺ＜ＰＯＳＴ＞として用いて、上述したのと同様に、下位モデル２２０におけるデータ同化を行い、データ同化の結果に基づいてパラメータベクトルＺのサンプルを生成する。データ算出装置１００は、所定の終了条件が成立するまで、階層的なデータ同化を繰り返す。
【0045】
このように、データ算出装置１００が階層的なデータ同化を繰り返すことで、パラメータ値を高精度に推定できると期待される。特に、パラメータベクトルＺの確率分布、および、パラメータベクトルθ１からθＳまでのそれぞれの確率分布が、何れもターゲット値に応じたあるパラメータ値に収束していくことが期待される。ここでいう確率分布があるパラメータ値に収束するとは、その値における確率が大きくなることである。
【0046】
例えば、データ同化によって、出力ベクトルＹの値がターゲット値に近い値となるパラメータＺの値の確率が比較的大きくなることで、データ同化を再度行う際に、出力ベクトルＹのサンプリング値としてターゲット値に近いものが多くなる。これにより、再度のデータ同化によって、出力ベクトルＹの値がターゲット値に近い値となるパラメータＺの値の確率がさらに大きくなると期待される。
パラメータベクトルθ１からθＳまでについても、データ同化を繰り返すことで、ターゲット値に応じたパラメータ値の確率が大きくなっていくと期待される。
パラメータベクトルθ１からθＳまでについても、データ同化を繰り返すことで、ターゲット値に応じたパラメータ値の確率が大きくなっていくと期待される。
【0047】
以下では、データ算出装置１００が、モデルの出力値をターゲット値に近付けるようにパラメータ値を推定する処理として、カーネルＡＢＣおよびカーネルハーディングを用いる場合を例に説明する。
ただし、データ算出装置１００が、モデルの出力値をターゲット値に近付けるようにパラメータ値を推定する方法は、特定の方法に限定されない。例えば、データ算出装置１００が、モデルの出力値をターゲット値に近付けるようにパラメータ値を推定する処理として、マルコフチェーンモンテカルロ（Markov Chain Monte Carlo；ＭＣＭＣ）、または、ＡＢＣ（Approximate Bayesian Computation）を用いるようにしてもよい。
ただし、データ算出装置１００が、モデルの出力値をターゲット値に近付けるようにパラメータ値を推定する方法は、特定の方法に限定されない。例えば、データ算出装置１００が、モデルの出力値をターゲット値に近付けるようにパラメータ値を推定する処理として、マルコフチェーンモンテカルロ（Markov Chain Monte Carlo；ＭＣＭＣ）、または、ＡＢＣ（Approximate Bayesian Computation）を用いるようにしてもよい。
【0048】
通信部１１０は、他の装置と通信を行う。例えば、データ算出装置１００がターゲットデータとしてセンサによる実測データを取得する場合、通信部１１０が、対象９００に設置されたセンサから、そのセンサによる実測データを受信するようにしてもよい。
【0049】
表示部１２０は、例えば液晶パネルまたはＬＥＤ（Light Emitting Diode、発光ダイオード）パネル等の表示画面を備え、各種画像を表示する。例えば、表示部１２０が、データ算出装置１００によるパラメータ値の推定結果を表示するようにしてもよい。
操作入力部１３０は、例えばキーボードおよびマウス等の入力デバイスを備え、ユーザ操作を受け付ける。例えば、操作入力部１３０が、ターゲットデータとして対象９００における目標値を入力するユーザ操作を受け付けるようにしてもよい。
操作入力部１３０は、例えばキーボードおよびマウス等の入力デバイスを備え、ユーザ操作を受け付ける。例えば、操作入力部１３０が、ターゲットデータとして対象９００における目標値を入力するユーザ操作を受け付けるようにしてもよい。
【0050】
記憶部１７０は、各種データを記憶する。例えば、記憶部１７０は、対象９００のモデルとして上記のように階層化されたモデルを記憶する。記憶部１７０は、データ算出装置１００が備える記憶デバイスを用いて構成される。
制御部１８０は、データ算出装置１００の各部を制御して各種処理を行う。制御部１８０の機能は、データ算出装置１００が備えるＣＰＵ（Central Processing Unit、中央処理装置）が記憶部１７０からプログラムを読み出して実行することで実行されてもよい。
制御部１８０は、データ算出装置１００の各部を制御して各種処理を行う。制御部１８０の機能は、データ算出装置１００が備えるＣＰＵ（Central Processing Unit、中央処理装置）が記憶部１７０からプログラムを読み出して実行することで実行されてもよい。
【0051】
初期設定部１８１は、階層的なデータ同化のための初期設定を行う。
具体的には、初期設定部１８１は、階層的なデータ同化を繰り返すためのループカウンタとして用いられる変数ｌに１を代入する。
また、初期設定部１８１は、パラメータベクトルＺの事前分布πＺに基づくサンプリングを行う。すなわち、事前分布πＺに基づいてサンプルＺ＜ＰＲＩ＞を算出する。初期設定部１８１は、サンプルＺ＜ＰＲＩ＞をＭ個算出する。Ｍは、Ｍ≧１の整数であり、サンプルの個数を示す。初期設定部１８１のサンプリングによって得られるＭ個のサンプルＺ＜ＰＲＩ＞（ｌ）を、｛Ｚ＜ＰＲＩ＞（１） ｍ｝ｍ＝１ Ｍと表記する。（）内の値は、変数ｌの値である。すなわち、（）内の値は、階層的なデータ同化の何回目の実行であるかを示す。「Ｚ＜ＰＲＩ＞（１） ｍ」の下付の「ｍ」はＭ個のサンプルＺ＜ＰＲＩ＞（１）を識別する識別番号を示す。
具体的には、初期設定部１８１は、階層的なデータ同化を繰り返すためのループカウンタとして用いられる変数ｌに１を代入する。
また、初期設定部１８１は、パラメータベクトルＺの事前分布πＺに基づくサンプリングを行う。すなわち、事前分布πＺに基づいてサンプルＺ＜ＰＲＩ＞を算出する。初期設定部１８１は、サンプルＺ＜ＰＲＩ＞をＭ個算出する。Ｍは、Ｍ≧１の整数であり、サンプルの個数を示す。初期設定部１８１のサンプリングによって得られるＭ個のサンプルＺ＜ＰＲＩ＞（ｌ）を、｛Ｚ＜ＰＲＩ＞（１） ｍ｝ｍ＝１ Ｍと表記する。（）内の値は、変数ｌの値である。すなわち、（）内の値は、階層的なデータ同化の何回目の実行であるかを示す。「Ｚ＜ＰＲＩ＞（１） ｍ」の下付の「ｍ」はＭ個のサンプルＺ＜ＰＲＩ＞（１）を識別する識別番号を示す。
【0052】
ループ処理部１８２は、階層的なデータ同化を繰り返すループの処理を行う。具体的には、ループ処理部１８２は、階層的なデータ同化を繰り返し回数に応じて、変数ｌの値をカウントアップする。また、ループ処理部１８２は、ループの終了条件が成立しているか否かを判定する。
ここでのループの終了条件として、いろいろな条件を用いることができる。例えば、ループの終了条件が、変数ｌの値が所定の閾値以上であることであってもよい。あるいは、ループの終了条件が、上位モデル２１０の出力値が所定の条件以上にターゲット値に近いと評価されることであってもよい。
ここでのループの終了条件として、いろいろな条件を用いることができる。例えば、ループの終了条件が、変数ｌの値が所定の閾値以上であることであってもよい。あるいは、ループの終了条件が、上位モデル２１０の出力値が所定の条件以上にターゲット値に近いと評価されることであってもよい。
【0053】
上位シミュレータ部１８３は、上位モデル２１０を用いたシミュレーションを実行する。
特に、上位シミュレータ部１８３は、下位モデル２２０を用いたシミュレーションの結果として得られるサンプルＺ＜ＰＲＥＤ＞ １からＺ＜ＰＲＥＤ＞ Ｓまでに基づく更新後のパラメータベクトルＺの値を上位モデル２１０に適用してシミュレーションを実行する。
上位シミュレータ部１８３は、上位シミュレータ部１８３上位モデルシミュレーション実行手段の例に該当する。
特に、上位シミュレータ部１８３は、下位モデル２２０を用いたシミュレーションの結果として得られるサンプルＺ＜ＰＲＥＤ＞ １からＺ＜ＰＲＥＤ＞ Ｓまでに基づく更新後のパラメータベクトルＺの値を上位モデル２１０に適用してシミュレーションを実行する。
上位シミュレータ部１８３は、上位シミュレータ部１８３上位モデルシミュレーション実行手段の例に該当する。
【0054】
サンプルＺ＜ＰＲＥＤ＞（ｌ）
１からＺ＜ＰＲＥＤ＞（ｌ）
Ｓまでに基づく更新後のパラメータベクトルＺの値をサンプルＺ＜ＰＲＩ＞（ｌ＋１）と表記する。上述したように、（）内の値は、変数ｌの値である。
上位シミュレータ部１８３は、階層的なデータ同化のｌ回目の実行にて得られるサンプルＺ＜ＰＲＥＤ＞（ｌ） １からＺ＜ＰＲＥＤ＞（ｌ） Ｓに基づいて、階層的なデータ同化のｌ＋１回目のためのサンプルＺ＜ＰＲＩ＞（ｌ＋１）を取得する。データ算出装置１００について上述したように、上位シミュレータ部１８３が、サンプルＺ＜ＰＲＥＤ＞（ｌ） ｓの要素の値を、サンプルＺ＜ＰＲＩ＞（ｌ＋１）の要素の値として用いるようにしてもよい。あるいは、データ算出装置１００について上述したように、上位シミュレータ部１８３が、サンプルＺ＜ＰＲＥＤ＞ １からサンプルＺ＜ＰＲＥＤ＞ Ｓまでの複数の要素の値の平均をとるなど、所定の方法でサンプルＺ＜ＰＲＥＤ＞（ｌ＋１） ｓの要素の値を算出するようにしてもよい。
上位シミュレータ部１８３は、階層的なデータ同化のｌ回目の実行にて得られるサンプルＺ＜ＰＲＥＤ＞（ｌ） １からＺ＜ＰＲＥＤ＞（ｌ） Ｓに基づいて、階層的なデータ同化のｌ＋１回目のためのサンプルＺ＜ＰＲＩ＞（ｌ＋１）を取得する。データ算出装置１００について上述したように、上位シミュレータ部１８３が、サンプルＺ＜ＰＲＥＤ＞（ｌ） ｓの要素の値を、サンプルＺ＜ＰＲＩ＞（ｌ＋１）の要素の値として用いるようにしてもよい。あるいは、データ算出装置１００について上述したように、上位シミュレータ部１８３が、サンプルＺ＜ＰＲＥＤ＞ １からサンプルＺ＜ＰＲＥＤ＞ Ｓまでの複数の要素の値の平均をとるなど、所定の方法でサンプルＺ＜ＰＲＥＤ＞（ｌ＋１） ｓの要素の値を算出するようにしてもよい。
【0055】
上位シミュレータ部１８３は、サンプルＺ＜ＰＲＩ＞（ｌ＋１）を上位モデル２１０に適用してシミュレーションを実行し、サンプルＹ＜ＰＲＩ＞（ｌ＋１）を算出する。
上位シミュレータ部１８３は、サンプルＺ＜ＰＲＥＤ＞（ｌ） １からＺ＜ＰＲＥＤ＞（ｌ） Ｓまでの組み合わせをＭ個取得して、サンプルＺ＜ＰＲＩ＞（ｌ＋１）をＭ個算出し、サンプルＹ＜ＰＲＩ＞（ｌ＋１）をＭ個算出する。上位シミュレータ部１８３が算出するＭ個のサンプルＺ＜ＰＲＩ＞（ｌ）を、｛Ｚ＜ＰＲＩ＞（ｌ） ｍ｝ｍ＝１ Ｍと表記する。上位シミュレータ部１８３が算出するＭ個のサンプルＹ＜ＰＲＩ＞（ｌ）を、｛Ｙ＜ＰＲＩ＞（ｌ） ｍ｝ｍ＝１ Ｍと表記する。「Ｚ＜ＰＲＩ＞（ｌ） ｍ」の下付の「ｍ」、および、「Ｙ＜ＰＲＩ＞（ｌ） ｍ」の下付の「ｍ」は、いずれも、Ｍ個のサンプルを識別する識別番号を示す。
上位シミュレータ部１８３は、サンプルＺ＜ＰＲＥＤ＞（ｌ） １からＺ＜ＰＲＥＤ＞（ｌ） Ｓまでの組み合わせをＭ個取得して、サンプルＺ＜ＰＲＩ＞（ｌ＋１）をＭ個算出し、サンプルＹ＜ＰＲＩ＞（ｌ＋１）をＭ個算出する。上位シミュレータ部１８３が算出するＭ個のサンプルＺ＜ＰＲＩ＞（ｌ）を、｛Ｚ＜ＰＲＩ＞（ｌ） ｍ｝ｍ＝１ Ｍと表記する。上位シミュレータ部１８３が算出するＭ個のサンプルＹ＜ＰＲＩ＞（ｌ）を、｛Ｙ＜ＰＲＩ＞（ｌ） ｍ｝ｍ＝１ Ｍと表記する。「Ｚ＜ＰＲＩ＞（ｌ） ｍ」の下付の「ｍ」、および、「Ｙ＜ＰＲＩ＞（ｌ） ｍ」の下付の「ｍ」は、いずれも、Ｍ個のサンプルを識別する識別番号を示す。
【0056】
上位事後分布データ算出部１８４は、上位モデル２１０の出力ベクトルＹに対するターゲットデータに基づく、上位モデル２１０のパラメータベクトルＺの値の事後分布を示すデータを算出する。上位事後分布データ算出部１８４は、上位事後分布データ算出手段の例に該当する。
具体的には、上位事後分布データ算出部１８４は、出力ベクトルＹの値がターゲット値に近付くようにパラメータベクトルＺの値を推定するデータ同化を、カーネルＡＢＣおよびカーネルハーディングを用いて行う。上位事後分布データ算出部１８４は、このデータ同化で、Ｍ個のサンプル｛Ｚ＜ＰＲＩ＞（ｌ） ｍ｝ｍ＝１ Ｍに応じて、サンプルＺ＜ＰＯＳＴ＞（ｌ）をＭ個算出する。上述したように、（）内の値は、変数ｌの値である。
具体的には、上位事後分布データ算出部１８４は、出力ベクトルＹの値がターゲット値に近付くようにパラメータベクトルＺの値を推定するデータ同化を、カーネルＡＢＣおよびカーネルハーディングを用いて行う。上位事後分布データ算出部１８４は、このデータ同化で、Ｍ個のサンプル｛Ｚ＜ＰＲＩ＞（ｌ） ｍ｝ｍ＝１ Ｍに応じて、サンプルＺ＜ＰＯＳＴ＞（ｌ）をＭ個算出する。上述したように、（）内の値は、変数ｌの値である。
【0057】
上位事後分布データ算出部１８４が算出するＭ個のサンプルＺ＜ＰＯＳＴ＞（ｌ）を、｛Ｚ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍと表記する。「Ｚ＜ＰＯＳＴ＞（ｌ）
ｍ」の下付の「ｍ」は、いずれも、Ｍ個のサンプル「Ｚ＜ＰＯＳＴ＞（ｌ）」を識別する識別番号を示す。Ｍ個のサンプル｛Ｚ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍは、ターゲットデータに基づく、上位モデル２１０のパラメータベクトルＺの値の事後分布ｐ（Ｚ｜Ｙ）を示すデータの例に該当する。
【0058】
上位事後分布データ算出部１８４は、カーネルＡＢＣの実行にて、ターゲット値と、｛Ｙ＜ＰＲＩ＞（ｌ）
ｍ｝ｍ＝１
Ｍの各々とを再生核ヒルベルト空間（Reproducing Kernel Hilbert Space；ＲＫＨＳ）における点に射影する。そして、上位事後分布データ算出部１８４は、ターゲット値が射影された点と、サンプルＹ＜ＰＲＩ＞（ｌ）
ｍが射影された点との再生核ヒルベルト空間における距離に基づいて、カーネル平均を算出する。上位事後分布データ算出部１８４は、算出したカーネル平均にカーネルハーディングを適用して、｛Ｚ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍを算出する。
再生核ヒルベルト空間における距離は、サンプルＹ＜ＰＲＩ＞（ｌ） ｍとターゲット値との類似度を示す指標の例に該当する。サンプルＹ＜ＰＲＩ＞（ｌ） ｍとターゲット値との類似度が高いほど、再生核ヒルベルト空間における距離は小さくなる。
再生核ヒルベルト空間における距離は、サンプルＹ＜ＰＲＩ＞（ｌ） ｍとターゲット値との類似度を示す指標の例に該当する。サンプルＹ＜ＰＲＩ＞（ｌ） ｍとターゲット値との類似度が高いほど、再生核ヒルベルト空間における距離は小さくなる。
【0059】
下位シミュレータ部１８５は、下位モデル２２０を用いたシミュレーションを実行する。
特に、下位シミュレータ部１８５は、データ算出装置１００について上述したように、下位モデル２２０－ｓにおけるデータ同化で得られるサンプルθ＜ＰＯＳＴ＞ ｓを下位モデル２２０－ｓに適用してシミュレーションを実行し、サンプルＺ＜ＰＲＥＤ＞ ｓを算出する。下位モデル２２０－ｓにおけるデータ同化は、上述した、下位モデル２２０－ｓの出力ベクトルＺｓの値がターゲット値（サンプルＺ＜ＰＯＳＴ＞ ｓ）に近付くように下位モデル２２０のパラメータベクトルθｓの値を推定するデータ同化である。
特に、下位シミュレータ部１８５は、データ算出装置１００について上述したように、下位モデル２２０－ｓにおけるデータ同化で得られるサンプルθ＜ＰＯＳＴ＞ ｓを下位モデル２２０－ｓに適用してシミュレーションを実行し、サンプルＺ＜ＰＲＥＤ＞ ｓを算出する。下位モデル２２０－ｓにおけるデータ同化は、上述した、下位モデル２２０－ｓの出力ベクトルＺｓの値がターゲット値（サンプルＺ＜ＰＯＳＴ＞ ｓ）に近付くように下位モデル２２０のパラメータベクトルθｓの値を推定するデータ同化である。
【0060】
下位シミュレータ部１８５は、下位事後分布データ算出部１８６がデータ同化で算出するＭ個のサンプルθ＜ＰＯＳＴ＞
ｓを取得して、Ｍ個のサンプルＺ＜ＰＲＥＤ＞
ｓを算出する。下位シミュレータ部１８５は、下位モデルシミュレーション実行手段の例に該当する。
下位事後分布データ算出部１８６が算出するＭ個のθ＜ＰＯＳＴ＞ ｓを、｛θ＜ＰＯＳＴ＞（ｌ） ｍ｝ｍ＝１ Ｍと表記する。上述したように、（）内の値は、変数ｌの値である。「θ＜ＰＯＳＴ＞（ｌ） ｍ」の下付の「ｍ」は、Ｍ個のサンプルθ＜ＰＯＳＴ＞（ｌ）を識別する識別番号を示す。
下位シミュレータ部１８５が算出するＭ個のＺ＜ＰＲＥＤ＞ ｓを、｛Ｚ＜ＰＲＥＤ＞（ｌ） ｓ，ｍ｝ｍ＝１ Ｍと表記する。上述したように、（）内の値は、変数ｌの値である。「Ｚ＜ＰＲＥＤ＞（ｌ） ｓ，ｍ」の下付の「ｍ」は、Ｍ個のサンプルＺ＜ＰＲＥＤ＞（ｌ） ｓを識別する識別番号を示す。
下位事後分布データ算出部１８６が算出するＭ個のθ＜ＰＯＳＴ＞ ｓを、｛θ＜ＰＯＳＴ＞（ｌ） ｍ｝ｍ＝１ Ｍと表記する。上述したように、（）内の値は、変数ｌの値である。「θ＜ＰＯＳＴ＞（ｌ） ｍ」の下付の「ｍ」は、Ｍ個のサンプルθ＜ＰＯＳＴ＞（ｌ）を識別する識別番号を示す。
下位シミュレータ部１８５が算出するＭ個のＺ＜ＰＲＥＤ＞ ｓを、｛Ｚ＜ＰＲＥＤ＞（ｌ） ｓ，ｍ｝ｍ＝１ Ｍと表記する。上述したように、（）内の値は、変数ｌの値である。「Ｚ＜ＰＲＥＤ＞（ｌ） ｓ，ｍ」の下付の「ｍ」は、Ｍ個のサンプルＺ＜ＰＲＥＤ＞（ｌ） ｓを識別する識別番号を示す。
【0061】
下位事後分布データ算出部１８６は、上位モデル２１０のパラメータベクトルＺの事後分布ｐ（Ｚ｜Ｙ）を示すデータに基づく、下位モデル２２０のパラメータ値の事後分布を示すデータを算出する。下位事後分布データ算出部１８６は、下位事後分布データ算出手段の例に該当する。
具体的には、下位事後分布データ算出部１８６は、下位モデル２２０－ｓの出力値が、ターゲット値であるサンプルＺ＜ＰＯＳＴ＞（ｌ）に近付くようにパラメータθｓの値を推定するデータ同化を、カーネルＡＢＣおよびカーネルハーディングを用いて行う。
下位事後分布データ算出部１８６は、データ同化にて、上述したようにＭ個のサンプル｛θ＜ＰＯＳＴ＞（ｌ） ｍ｝ｍ＝１ Ｍを算出する。
具体的には、下位事後分布データ算出部１８６は、下位モデル２２０－ｓの出力値が、ターゲット値であるサンプルＺ＜ＰＯＳＴ＞（ｌ）に近付くようにパラメータθｓの値を推定するデータ同化を、カーネルＡＢＣおよびカーネルハーディングを用いて行う。
下位事後分布データ算出部１８６は、データ同化にて、上述したようにＭ個のサンプル｛θ＜ＰＯＳＴ＞（ｌ） ｍ｝ｍ＝１ Ｍを算出する。
【0062】
制御部１８０が、通信部１１０を介して対象９００に制御信号を送信するなど、対象９００を制御するようにしてもよい。例えば、制御部１８０が、モデル２００のパラメータの推定結果に基づいて、あるいはモデル２００を用いたシミュレーションの実行結果に基づいて、対象９００の各装置など対象９００の各部を制御するようにしてもよい。さらに例えば、制御部１８０が、下位事後分布データ算出部１８６が算出する下位モデル２２０のパラメータ値の事後分布を示すデータ（サンプル｛θ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍ）に基づいて、対象９００を制御するようにしてもよい。
制御部１８０は、制御手段の例に該当する。
制御部１８０は、制御手段の例に該当する。
【0063】
図３は、データ算出装置１００が階層的なデータ同化を行う処理手順の例を示すフローチャートである。
（ステップＳ１０１）
初期設定部１８１は、ループカウンタとして用いられる変数ｌの値を初期設定する。具体的には、初期設定部１８１は、変数ｌの値を１に設定する。変数ｌの値は、ステップＳ１０３からＳ１１０までのループの繰り返し回数を示す。
ステップＳ１０１の後、処理がステップＳ１０２へ遷移する。
（ステップＳ１０１）
初期設定部１８１は、ループカウンタとして用いられる変数ｌの値を初期設定する。具体的には、初期設定部１８１は、変数ｌの値を１に設定する。変数ｌの値は、ステップＳ１０３からＳ１１０までのループの繰り返し回数を示す。
ステップＳ１０１の後、処理がステップＳ１０２へ遷移する。
【0064】
（ステップＳ１０２）
初期設定部１８１は、パラメータベクトルＺの事前分布πＺからのＭ個のサンプル｛Ｚ＜ＰＲＩ＞（１） ｍ｝ｍ＝１ Ｍを算出する。
ステップＳ１０２の後、処理がステップＳ１０３へ遷移する。
初期設定部１８１は、パラメータベクトルＺの事前分布πＺからのＭ個のサンプル｛Ｚ＜ＰＲＩ＞（１） ｍ｝ｍ＝１ Ｍを算出する。
ステップＳ１０２の後、処理がステップＳ１０３へ遷移する。
【0065】
（ステップＳ１０３）
ループ処理部１８２は、ループの終了条件が成立しているか否かを判定する。
終了条件が成立していないとループ処理部１８２が判定した場合（ステップＳ１０３：ＮＯ）、処理がステップＳ１０４へ遷移する。
一方、終了条件が成立しているとループ処理部１８２が判定した場合（ステップＳ１０３：ＹＥＳ）、データ算出装置１００は、図３の処理を終了する。
ループ処理部１８２は、ループの終了条件が成立しているか否かを判定する。
終了条件が成立していないとループ処理部１８２が判定した場合（ステップＳ１０３：ＮＯ）、処理がステップＳ１０４へ遷移する。
一方、終了条件が成立しているとループ処理部１８２が判定した場合（ステップＳ１０３：ＹＥＳ）、データ算出装置１００は、図３の処理を終了する。
【0066】
（ステップＳ１０４）
上位シミュレータ部１８３は、サンプルＺ＜ＰＲＩ＞（ｌ） ｍ（ｍ＝１，・・・，Ｍ）を上位モデル２１０のパラメータに設定してシミュレーションを実行し、出力のサンプルＹ＜ＰＲＩ＞（ｌ） ｍを算出する。Ｙ＜ＰＲＩ＞（ｌ） ｍは、式（１）のように示される。
上位シミュレータ部１８３は、サンプルＺ＜ＰＲＩ＞（ｌ） ｍ（ｍ＝１，・・・，Ｍ）を上位モデル２１０のパラメータに設定してシミュレーションを実行し、出力のサンプルＹ＜ＰＲＩ＞（ｌ） ｍを算出する。Ｙ＜ＰＲＩ＞（ｌ） ｍは、式（１）のように示される。
【0067】
【数1】
【0068】
ステップＳ１０４の後、処理がステップＳ１０５へ遷移する。
【0069】
（ステップＳ１０５）
上位事後分布データ算出部１８４は、カーネル平均μ＾（ｌ） Ｚ｜Ｙを算出する。カーネル平均μ＾（ｌ） Ｚ｜Ｙは、式（２）のように示される。
上位事後分布データ算出部１８４は、カーネル平均μ＾（ｌ） Ｚ｜Ｙを算出する。カーネル平均μ＾（ｌ） Ｚ｜Ｙは、式（２）のように示される。
【0070】
【数2】
【0071】
ｋｚは、カーネルを示す。
括弧内の「・」は、再生核ヒルベルト空間における関数の変数が特定のものに限定されないことを示す。
重みｗｍは、式（３）のように示される。
括弧内の「・」は、再生核ヒルベルト空間における関数の変数が特定のものに限定されないことを示す。
重みｗｍは、式（３）のように示される。
【0072】
【数3】
【0073】
「Ｉ」は単位行列を表す。
δは正則化定数（Regularization Constant）を示す。δ＞０である。
Ｇｙはグラム行列（Gram Matrix）を示し、式（４）のように表される。
δは正則化定数（Regularization Constant）を示す。δ＞０である。
Ｇｙはグラム行列（Gram Matrix）を示し、式（４）のように表される。
【0074】
【数4】
【0075】
ｋｙ（Ｙ）はｙのカーネルを示し、式（５）のように表される。
【0076】
【数5】
【0077】
上付のＴは、行列またはベクトルの転置を示す。
ステップＳ１０５の後、処理がステップＳ１０６へ遷移する。
ステップＳ１０５の後、処理がステップＳ１０６へ遷移する。
【0078】
（ステップＳ１０６）
上位事後分布データ算出部１８４は、カーネル平均μ＾（ｌ） Ｚ｜ＹによるパラメータＺのサンプル｛Ｚ＜ＰＯＳＴ＞（ｌ） ｍ｝ｍ＝１ Ｍを算出する。
カーネル平均μ＾（ｌ） Ｚ｜Ｙは、出力ベクトルＹとしてターゲットデータが与えられた場合の事後分布ｐ（Ｚ｜Ｙ）と一対一に紐付けられる。カーネルハーディング（Kernel Herding）の処理を行うことで、事後分布によるサンプルを得られる。
例えば、上位事後分布データ算出部１８４は、式（６）に基づいてＺ＜ＰＯＳＴ＞（ｌ） ｍを算出する。
上位事後分布データ算出部１８４は、カーネル平均μ＾（ｌ） Ｚ｜ＹによるパラメータＺのサンプル｛Ｚ＜ＰＯＳＴ＞（ｌ） ｍ｝ｍ＝１ Ｍを算出する。
カーネル平均μ＾（ｌ） Ｚ｜Ｙは、出力ベクトルＹとしてターゲットデータが与えられた場合の事後分布ｐ（Ｚ｜Ｙ）と一対一に紐付けられる。カーネルハーディング（Kernel Herding）の処理を行うことで、事後分布によるサンプルを得られる。
例えば、上位事後分布データ算出部１８４は、式（６）に基づいてＺ＜ＰＯＳＴ＞（ｌ） ｍを算出する。
【0079】
【数6】
【0080】
「ａｒｇｍａｘＺｈｊ（Ｚ）」は、「ｈｊ（Ｚ）」が最大となる「Ｚ」を出力する関数である。
ｈｊは、式（７）のように再帰的に示される。
ｈｊは、式（７）のように再帰的に示される。
【0081】
【数7】
【0082】
式（７）の「Ｈ」は、再生核ヒルベルト空間を示す。
ｈ０は、式（８）のように示される。
ｈ０は、式（８）のように示される。
【0083】
【数8】
【0084】
ステップＳ１０６の後、処理がステップＳ１０７へ遷移する。
【0085】
（ステップＳ１０７）
上位事後分布データ算出部１８４は、Ｍ個のサンプルＺ＜ＰＯＳＴ＞（ｌ） ｍ（ｍ＝１，・・・，Ｍ）を、サンプルＺ＜ＰＯＳＴ＞（ｌ）に集約する。例えば、上位事後分布データ算出部１８４は、式（９）に示されるように、Ｍ個のサンプルＺ＜ＰＯＳＴ＞（ｌ） ｍ（ｍ＝１，・・・，Ｍ）を、サンプルＺ＜ＰＯＳＴ＞（ｌ）に集約する。
上位事後分布データ算出部１８４は、Ｍ個のサンプルＺ＜ＰＯＳＴ＞（ｌ） ｍ（ｍ＝１，・・・，Ｍ）を、サンプルＺ＜ＰＯＳＴ＞（ｌ）に集約する。例えば、上位事後分布データ算出部１８４は、式（９）に示されるように、Ｍ個のサンプルＺ＜ＰＯＳＴ＞（ｌ） ｍ（ｍ＝１，・・・，Ｍ）を、サンプルＺ＜ＰＯＳＴ＞（ｌ）に集約する。
【0086】
【数9】
【0087】
ステップＳ１０７の後、処理がステップＳ１０８へ遷移する。
【0088】
（ステップＳ１０８）
データ算出装置１００は、下位モデル２２０ごとの処理を実行する。上述したように、データ算出装置１００は、下位モデル２２０－ｓの処理でサンプルＺ＜ＰＲＥＤ＞（ｌ） ｓを算出する。
ステップＳ１０８の後、処理がステップＳ１０９に遷移する。
データ算出装置１００は、下位モデル２２０ごとの処理を実行する。上述したように、データ算出装置１００は、下位モデル２２０－ｓの処理でサンプルＺ＜ＰＲＥＤ＞（ｌ） ｓを算出する。
ステップＳ１０８の後、処理がステップＳ１０９に遷移する。
【0089】
（ステップＳ１０９）
上位シミュレータ部１８３は、下位モデル２２０ごとの処理結果を上位モデル２１０の処理に受け渡す。具体的には、データ算出装置１００について上述したように、上位シミュレータ部１８３は、下位モデル２２０ごとに得られたサンプルＺ＜ＰＲＥＤ＞（ｌ） ｓの値に基づいて、サンプルＺ＜ＰＲＩ＞（ｌ＋１）の値を設定する。例えば、上位シミュレータ部１８３は、式（１０）に示されるように、サンプルＺ＜ＰＲＩ＞（ｌ＋１）の値を設定する。
上位シミュレータ部１８３は、下位モデル２２０ごとの処理結果を上位モデル２１０の処理に受け渡す。具体的には、データ算出装置１００について上述したように、上位シミュレータ部１８３は、下位モデル２２０ごとに得られたサンプルＺ＜ＰＲＥＤ＞（ｌ） ｓの値に基づいて、サンプルＺ＜ＰＲＩ＞（ｌ＋１）の値を設定する。例えば、上位シミュレータ部１８３は、式（１０）に示されるように、サンプルＺ＜ＰＲＩ＞（ｌ＋１）の値を設定する。
【0090】
【数10】
【0091】
ステップＳ１０９の後、処理がステップＳ１１０に遷移する。
【0092】
（ステップＳ１１０）
ループ処理部１８２は、ループカウンタとして用いられる変数ｌの値をカウントアップする。すなわち、ループ処理部１８２は、変数ｌの値を１増加させる。
ステップＳ１１０の後、処理がステップＳ１０２へ遷移する。
ループ処理部１８２は、ループカウンタとして用いられる変数ｌの値をカウントアップする。すなわち、ループ処理部１８２は、変数ｌの値を１増加させる。
ステップＳ１１０の後、処理がステップＳ１０２へ遷移する。
【0093】
図４は、データ算出装置１００が下位モデル２２０ごとに行う処理の手順の例を示すフローチャートである。データ算出装置１００は、図３のステップＳ１０８で図４の処理を行う。
（ステップＳ２０１）
上位事後分布データ算出部１８４は、上位モデル２１０のパラメータベクトルＺの値を下位モデル２２０に受け渡す。具体的には、データ算出装置１００について上述したように、上位事後分布データ算出部１８４は、サンプルＺ＜ＰＯＳＴ＞（ｌ）に基づいてサンプルＺ’＜ＰＯＳＴ＞（ｌ） ｓを取得する。例えば、上位事後分布データ算出部１８４は、式（１１）に基づいてＺ’＜ＰＯＳＴ＞（ｌ） ｓを取得する。
（ステップＳ２０１）
上位事後分布データ算出部１８４は、上位モデル２１０のパラメータベクトルＺの値を下位モデル２２０に受け渡す。具体的には、データ算出装置１００について上述したように、上位事後分布データ算出部１８４は、サンプルＺ＜ＰＯＳＴ＞（ｌ）に基づいてサンプルＺ’＜ＰＯＳＴ＞（ｌ） ｓを取得する。例えば、上位事後分布データ算出部１８４は、式（１１）に基づいてＺ’＜ＰＯＳＴ＞（ｌ） ｓを取得する。
【0094】
【数11】
【0095】
Ｚ＜ＰＯＳＴ＞（ｌ）［ｓ，：］は、行列Ｚ＜ＰＯＳＴ＞（ｌ）のｓ行目の行ベクトルを示す。
式（１１）では、Ｚ＜ＰＯＳＴ＞（ｌ） ｍと区別するために、Ｚに’を付している。
Ｚ’＜ＰＯＳＴ＞（ｌ） ｓは、式（１２）のように示される。
式（１１）では、Ｚ＜ＰＯＳＴ＞（ｌ） ｍと区別するために、Ｚに’を付している。
Ｚ’＜ＰＯＳＴ＞（ｌ） ｓは、式（１２）のように示される。
【0096】
【数12】
【0097】
Ｚ’＜ＰＯＳＴ＞（ｌ）
ｓ、ｍは、式（１３）のように示される。
【0098】
【数13】
【0099】
Ｚ＜ＰＯＳＴ＞（ｌ）［ｓ，ｍ］は、行列Ｚ＜ＰＯＳＴ＞（ｌ）のｓ行目かつｍ列目の要素を示す。
式（１３）の例で、Ｚ＜ＰＯＳＴ＞（ｌ）は、Ｚ’＜ＰＯＳＴ＞（ｌ） ｓ、ｍを要素とする行列として構成されている。上位事後分布データ算出部１８４は、式（１１）に示されるように、Ｚ＜ＰＯＳＴ＞（ｌ）のｓ行目を読み出すことで、下位モデル２２０－ｓのターゲットデータとして用いるためのＭ個のサンプルＺ’＜ＰＯＳＴ＞（ｌ） ｓ、１，・・・，Ｚ’＜ＰＯＳＴ＞（ｌ） ｓ、Ｍを取得する。
ステップＳ２０１の後、処理がステップ２０２へ遷移する。
式（１３）の例で、Ｚ＜ＰＯＳＴ＞（ｌ）は、Ｚ’＜ＰＯＳＴ＞（ｌ） ｓ、ｍを要素とする行列として構成されている。上位事後分布データ算出部１８４は、式（１１）に示されるように、Ｚ＜ＰＯＳＴ＞（ｌ）のｓ行目を読み出すことで、下位モデル２２０－ｓのターゲットデータとして用いるためのＭ個のサンプルＺ’＜ＰＯＳＴ＞（ｌ） ｓ、１，・・・，Ｚ’＜ＰＯＳＴ＞（ｌ） ｓ、Ｍを取得する。
ステップＳ２０１の後、処理がステップ２０２へ遷移する。
【0100】
（ステップＳ２０２）
下位シミュレータ部１８５は、事前分布πθｓからのＭ個のサンプル｛θ＜ＰＲＩ＞（ｌ） ｓ，ｍ｝ｍ＝１ Ｍを算出する。
ステップＳ２０２の後、処理がステップＳ２０３へ遷移する。
下位シミュレータ部１８５は、事前分布πθｓからのＭ個のサンプル｛θ＜ＰＲＩ＞（ｌ） ｓ，ｍ｝ｍ＝１ Ｍを算出する。
ステップＳ２０２の後、処理がステップＳ２０３へ遷移する。
【0101】
（ステップＳ２０３）
下位シミュレータ部１８５は、サンプルθ＜ＰＲＩ＞（ｌ） ｓ，ｍを下位モデル２２０－ｓに適用してシミュレーションを実行し、サンプルＺ＜ＰＲＩ＞（ｌ） ｓ，ｍを算出する。サンプルＺ＜ＰＲＩ＞（ｌ） ｓ，ｍは、式（１４）のように示される。
下位シミュレータ部１８５は、サンプルθ＜ＰＲＩ＞（ｌ） ｓ，ｍを下位モデル２２０－ｓに適用してシミュレーションを実行し、サンプルＺ＜ＰＲＩ＞（ｌ） ｓ，ｍを算出する。サンプルＺ＜ＰＲＩ＞（ｌ） ｓ，ｍは、式（１４）のように示される。
【0102】
【数14】
【0103】
ステップＳ２０３の後、処理がステップＳ２０４へ遷移する。
【0104】
（ステップＳ２０４）
下位事後分布データ算出部１８６は、カーネル平均μ＾（ｌ） θ｜Ｚ，ｓを算出する。カーネル平均μ＾（ｌ） θ｜Ｚ，ｓは、式（１５）のように示される。
下位事後分布データ算出部１８６は、カーネル平均μ＾（ｌ） θ｜Ｚ，ｓを算出する。カーネル平均μ＾（ｌ） θ｜Ｚ，ｓは、式（１５）のように示される。
【0105】
【数15】
【0106】
ｋθは、カーネルを示す。
重みｗｓ，ｍは、式（１６）のように示される。
重みｗｓ，ｍは、式（１６）のように示される。
【0107】
【数16】
【0108】
Ｇｚ，ｓはグラム行列を示し、式（１７）のように表される。
【0109】
【数17】
【0110】
Ｋθ，ｓ（Ｚ）はＺのカーネルを示し、式（１８）のように表される。
【0111】
【数18】
【0112】
ステップＳ２０４の後、処理がステップＳ２０５へ遷移する。
【0113】
（ステップＳ２０５）
下位事後分布データ算出部１８６は、カーネル平均μ＾（ｌ） θ｜Ｚ，ｓによるパラメータθｓのサンプル｛θ＜ＰＯＳＴ＞（ｌ） ｓ，ｍ｝ｍ＝１ Ｍを算出する。
例えば、下位事後分布データ算出部１８６は、式（１９）に基づいてθ＜ＰＯＳＴ＞（ｌ） ｓ，ｍを算出する。
下位事後分布データ算出部１８６は、カーネル平均μ＾（ｌ） θ｜Ｚ，ｓによるパラメータθｓのサンプル｛θ＜ＰＯＳＴ＞（ｌ） ｓ，ｍ｝ｍ＝１ Ｍを算出する。
例えば、下位事後分布データ算出部１８６は、式（１９）に基づいてθ＜ＰＯＳＴ＞（ｌ） ｓ，ｍを算出する。
【0114】
【数19】
【0115】
ｈｓ，ｊは、式（２０）のように再帰的に示される。
【0116】
【数20】
【0117】
ｈｓ，０は、式（２１）のように示される。
【0118】
【数21】
【0119】
ステップＳ２０５の後、処理がステップＳ２０６へ遷移する。
【0120】
（ステップＳ２０６）
下位シミュレータ部１８５は、サンプルθ＜ＰＯＳＴ＞（ｌ） s,ｍを下位モデル２２０－ｓのパラメータに設定してシミュレーションを実行し、サンプルＺ＜ＰＲＥＤ＞（ｌ） ｓ，ｍを算出する。Ｚ＜ＰＲＥＤ＞（ｌ） ｓ，ｍは、式（２２）のように示される。
下位シミュレータ部１８５は、サンプルθ＜ＰＯＳＴ＞（ｌ） s,ｍを下位モデル２２０－ｓのパラメータに設定してシミュレーションを実行し、サンプルＺ＜ＰＲＥＤ＞（ｌ） ｓ，ｍを算出する。Ｚ＜ＰＲＥＤ＞（ｌ） ｓ，ｍは、式（２２）のように示される。
【0121】
【数22】
【0122】
ステップＳ２０６の後、処理がステップＳ２０７へ遷移する。
【0123】
（ステップＳ２０７）
下位事後分布データ算出部１８６は、Ｍ個のサンプルＺ＜ＰＲＥＤ＞（ｌ） ｓ，ｍ（ｍ＝１，・・・，Ｍ）を、サンプルＺ＜ＰＲＥＤ＞（ｌ） ｓに集約する。
ステップＳ２０７の後、データ算出装置１００は、図４の処理を終了する。
下位事後分布データ算出部１８６は、Ｍ個のサンプルＺ＜ＰＲＥＤ＞（ｌ） ｓ，ｍ（ｍ＝１，・・・，Ｍ）を、サンプルＺ＜ＰＲＥＤ＞（ｌ） ｓに集約する。
ステップＳ２０７の後、データ算出装置１００は、図４の処理を終了する。
【0124】
データ算出装置１００が行う処理について、例を用いてさらに説明する。
図５は、対象９００の例に該当する生産システム９１０の構成を示す図である。図５に示す構成で、生産システム９１０は、第一サブシステム９１１－１と、第二サブシステム９１１－２とを備える。第一サブシステム９１１－１と第二サブシステム９１１－２とを総称して、サブシステム９１１と表記する。
図５は、対象９００の例に該当する生産システム９１０の構成を示す図である。図５に示す構成で、生産システム９１０は、第一サブシステム９１１－１と、第二サブシステム９１１－２とを備える。第一サブシステム９１１－１と第二サブシステム９１１－２とを総称して、サブシステム９１１と表記する。
【0125】
第一サブシステム９１１－１は、第一生産ライン９１２－１、第二生産ライン９１２－２、第三生産ライン９１２－３、第四生産ライン９１２－４、および、第五生産ライン９１２－５を備える。第二サブシステム９１１－２は、第六生産ライン９１２－６、第七生産ライン９１２－７、第八生産ライン９１２－８、および、第九生産ライン９１２－９を備える。第一生産ライン９１２－１から第九生産ライン９１２－９までを総称して、生産ライン９１２と表記する。
【0126】
図６は、生産システム９１０を模擬するモデルの例を示す図である。図６の例では、サブシステム９１１ごとに下位モデル２２０が設けられており、下位モデル２２０の個数は２つとなっている。
上位モデル２１０では、２つの下位モデルに応じて、２つのパラメータｚ１およびｚ２が設けられている。Ｚ＝（ｚ１，ｚ２）と表される。
上位モデル２１０では、２つの下位モデルに応じて、２つのパラメータｚ１およびｚ２が設けられている。Ｚ＝（ｚ１，ｚ２）と表される。
【0127】
下位モデル２２０－１では、第一サブシステム９１１－１が備える第一生産ライン９１２－１から第五生産ライン９１２－５までの５つの生産ライン９１２に応じて、５つのパラメータθ１，１からθ１，５までが設けられている。θ１＝（θ１，１，・・・，θ１，５）と表される。
【0128】
下位モデル２２０－２では、第二サブシステム９１１－２が備える第六生産ライン９１２－６から第九生産ライン９１２－９までの４つの生産ライン９１２に応じて、４つのパラメータθ２，１からθ２，４までが設けられている。θ２＝（θ２，１，・・・，θ２，４）と表される。
【0129】
モデルの出力およびパラメータとして、製品または部品など各部での生産単位１つあたりの生産にかかる時間を用いる。
パラメータｚ１、ｚ１は、それぞれ第一サブシステム９１１－１、第二サブシステム９１１－２で、生産単位１つあたりの生産にかかる時間を示す。
パラメータθ１，１からθ１，５までは、第一生産ライン９１２－１から第五生産ライン９１２－５までで、生産単位１つあたりの生産にかかる時間を示す。パラメータθ２，１からθ２，４までは、第六生産ライン９１２－６から第九生産ライン９１２－９までで、生産単位１つあたりの生産にかかる時間を示す。
パラメータｚ１、ｚ１は、それぞれ第一サブシステム９１１－１、第二サブシステム９１１－２で、生産単位１つあたりの生産にかかる時間を示す。
パラメータθ１，１からθ１，５までは、第一生産ライン９１２－１から第五生産ライン９１２－５までで、生産単位１つあたりの生産にかかる時間を示す。パラメータθ２，１からθ２，４までは、第六生産ライン９１２－６から第九生産ライン９１２－９までで、生産単位１つあたりの生産にかかる時間を示す。
【0130】
ターゲットデータに目標時間を設定する場合、このターゲット時間、および、上位モデル２１０から下位モデル２２０への出力（サンプル｛Ｚ＜ＰＯＳＴ＞（ｌ）ｍ｝ｍ＝１
Ｍ）は、タクトタイム（Takt Time）のような、時間についての要求を示すと考えられる。一方、下位モデル２２０から上位モデル２１０への出力（サンプル｛Ｚ＜ＰＲＥＤ＞（ｌ）
ｓ，ｍ｝ｍ＝１
Ｍ）は、処理に必要な時間を示すと考えられる。
この場合、上位モデル２１０と下位モデル２２０とがデータをやり取りすることで、要求される時間と必要な時間との両方に対応できる解を探索すると見ることができる。
この場合、上位モデル２１０と下位モデル２２０とがデータをやり取りすることで、要求される時間と必要な時間との両方に対応できる解を探索すると見ることができる。
【0131】
図７は、事前分布に基づく上位モデル２１０のパラメータ値のサンプルの例を示す図である。図７の横軸は、パラメータｚ１の値を示す。縦軸は、パラメータｚ２の値を示す。
領域Ａ１１に、パラメータｚ１の事前分布が示されている。領域Ａ１２に、パラメータｚ２の事前分布が示されている。領域Ａ１３の点の各々は、（ｚ１，ｚ２）のサンプルの例を示している。領域Ａ１３に示されるサンプルは、｛Ｚ＜ＰＲＩ＞（１） ｍ｝ｍ＝１ １００と表される。
図７は、図３のステップＳ１０２で初期設定部１８１が算出する、事前分布に基づくパラメータＺのサンプルの例を示している。
領域Ａ１１に、パラメータｚ１の事前分布が示されている。領域Ａ１２に、パラメータｚ２の事前分布が示されている。領域Ａ１３の点の各々は、（ｚ１，ｚ２）のサンプルの例を示している。領域Ａ１３に示されるサンプルは、｛Ｚ＜ＰＲＩ＞（１） ｍ｝ｍ＝１ １００と表される。
図７は、図３のステップＳ１０２で初期設定部１８１が算出する、事前分布に基づくパラメータＺのサンプルの例を示している。
【0132】
図８は、上位モデル２１０の出力値のサンプルの第１例を示す図である。図８の横軸は、上位モデル２１０の出力値を示す。縦軸は度数を示す。
図８は、図７に示されるパラメータベクトルＺの値ごとに上位モデル２１０を用いたシミュレーションを実行して得られる、上位モデル２１０の出力値を度数分布にて示している。図３のステップＳ１０４で、上位シミュレータ部１８３が、このシミュレーションを行う。
上位事後分布データ算出部１８４は、図５のステップＳ１０５およびＳ１０６の処理で、図８に示されるモデル出力が、実際値である「５」に近付くような、パラメータＺの値を推定する。
図８は、図７に示されるパラメータベクトルＺの値ごとに上位モデル２１０を用いたシミュレーションを実行して得られる、上位モデル２１０の出力値を度数分布にて示している。図３のステップＳ１０４で、上位シミュレータ部１８３が、このシミュレーションを行う。
上位事後分布データ算出部１８４は、図５のステップＳ１０５およびＳ１０６の処理で、図８に示されるモデル出力が、実際値である「５」に近付くような、パラメータＺの値を推定する。
【0133】
図９は、事後分布に基づく上位モデル２１０のパラメータ値のサンプルの例を示す図である。図９の横軸は、パラメータｚ１の値を示す。縦軸は、パラメータｚ２の値を示す。
領域Ａ２１に、パラメータｚ１の事後分布が示されている。領域Ａ２２に、パラメータｚ２の事後分布が示されている。領域Ａ２３の点の各々は、（ｚ１，ｚ２）のサンプルの例を示している。領域Ａ２３に示されるサンプルは、｛Ｚ＜ＰＯＳＴ＞（１） ｍ｝ｍ＝１ １００と表される。
領域Ａ２１に、パラメータｚ１の事後分布が示されている。領域Ａ２２に、パラメータｚ２の事後分布が示されている。領域Ａ２３の点の各々は、（ｚ１，ｚ２）のサンプルの例を示している。領域Ａ２３に示されるサンプルは、｛Ｚ＜ＰＯＳＴ＞（１） ｍ｝ｍ＝１ １００と表される。
【0134】
図９は、図８に示されるモデル出力および実際値を用いて、上位事後分布データ算出部１８４が推定するパラメータ値の例を示している。図７に示されるパラメータＺの事前分布のサンプルの場合、領域Ａ１３全体にプロットされている。これに対し、図９に示されるパラメータＺの事後分布のサンプルの場合、ｚ１の値が１０から１５程度、かつ、ｚ２の値が１０から１５程度の領域に多くプロットされている。
上位事後分布データ算出部１８４によるデータ同化によって、パラメータＺの値が収束していると考えられる。
上位事後分布データ算出部１８４によるデータ同化によって、パラメータＺの値が収束していると考えられる。
【0135】
図１０は、事前分布に基づく下位モデル２２０のパラメータ値のサンプルの例を示す図である。
図１０は、下位モデル２２０－１のパラメータθ１，１からθ１，５まで、および下位モデル２２０－２のパラメータθ２，１からθ２，４までのそれぞれについて、事前分布に基づくサンプルを、横軸にパラメータ値をとり縦軸に度数をとった度数分布にて示している。
下位事後分布データ算出部１８６は、図４のステップＳ２０２で、図１０に例示される事前分布に基づいて下位モデル２２０のパラメータ値をサンプリングする。
図１０は、下位モデル２２０－１のパラメータθ１，１からθ１，５まで、および下位モデル２２０－２のパラメータθ２，１からθ２，４までのそれぞれについて、事前分布に基づくサンプルを、横軸にパラメータ値をとり縦軸に度数をとった度数分布にて示している。
下位事後分布データ算出部１８６は、図４のステップＳ２０２で、図１０に例示される事前分布に基づいて下位モデル２２０のパラメータ値をサンプリングする。
【0136】
図１１は、下位モデル２２０の出力値のサンプルの例を示す図である。図１１の横軸は、下位モデル２２０の出力値を示す。縦軸は確率を示す。
図１１は、図１０に示されるパラメータ値を下位モデル２２０に適用してシミュレーションを実行して得られる下位モデル２２０の出力値を確率分布にて示している。図４のステップＳ２０３で、下位シミュレータ部１８５が、このシミュレーションを行う。
図１１は、図１０に示されるパラメータ値を下位モデル２２０に適用してシミュレーションを実行して得られる下位モデル２２０の出力値を確率分布にて示している。図４のステップＳ２０３で、下位シミュレータ部１８５が、このシミュレーションを行う。
【0137】
図１２は、事後分布に基づく下位モデル２２０のパラメータ値のサンプルの例を示す図である。
図１２は、下位モデル２２０－１のパラメータθ１，１からθ１，５まで、および下位モデル２２０－２のパラメータθ２，１からθ２，４までのそれぞれについて、事後分布に基づくサンプルを、横軸にパラメータ値をとり縦軸に度数をとった度数分布にて示している。
図１２は、下位モデル２２０－１のパラメータθ１，１からθ１，５まで、および下位モデル２２０－２のパラメータθ２，１からθ２，４までのそれぞれについて、事後分布に基づくサンプルを、横軸にパラメータ値をとり縦軸に度数をとった度数分布にて示している。
【0138】
図４のステップＳ２０４およびＳ２０５で、下位事後分布データ算出部１８６がデータ同化を行って、事後分布に基づく下位モデル２２０のパラメータ値のサンプルを算出する。
図１０と図１２とを比較すると、図１０よりも図１２のほうが、パラメータ値が偏在している。下位事後分布データ算出部１８６によるデータ同化によって、パラメータθ１およびθ２の値が収束していると考えられる。
図１０と図１２とを比較すると、図１０よりも図１２のほうが、パラメータ値が偏在している。下位事後分布データ算出部１８６によるデータ同化によって、パラメータθ１およびθ２の値が収束していると考えられる。
【0139】
図１３は、下位モデル２２０の出力値のサンプルの例を示す図である。図１３の横軸は、パラメータｚ１の値を示す。縦軸は、パラメータｚ２の値を示す。
図１３は、図１２に示されるパラメータ値を下位モデル２２０に適用してシミュレーションを実行した場合の、下位モデル２２０の出力値を示している。図４のステップＳ２０６で、下位シミュレータ部１８５が、このシミュレーションを行う。
図１３は、図１２に示されるパラメータ値を下位モデル２２０に適用してシミュレーションを実行した場合の、下位モデル２２０の出力値を示している。図４のステップＳ２０６で、下位シミュレータ部１８５が、このシミュレーションを行う。
【0140】
領域Ａ３３の点の各々は、（ｚ１，ｚ２）のサンプルの例を示している。領域Ａ３３に示されるサンプルは、｛Ｚ＜ＰＲＥＤ＞（１）
ｓ，ｍ｝ｍ＝１
１００と表される。
領域Ａ３１に、領域Ａ３３に示されるサンプルに基づくパラメータｚ１の分布が示されている。領域Ａ３２に、領域Ａ３３に示されるサンプルに基づくパラメータｚ２の分布が示されている。
領域Ａ３１に、領域Ａ３３に示されるサンプルに基づくパラメータｚ１の分布が示されている。領域Ａ３２に、領域Ａ３３に示されるサンプルに基づくパラメータｚ２の分布が示されている。
【0141】
図１３では、図９の場合よりもさらにパラメータ値が偏在しており、ｚ１の値が６から１０程度、かつ、ｚ２の値が７から１５程度の領域に多くプロットされている。
下位事後分布データ算出部１８６によるデータ同化の結果が反映されることで、パラメータＺの値がさらに収束していると考えられる。
下位事後分布データ算出部１８６によるデータ同化の結果が反映されることで、パラメータＺの値がさらに収束していると考えられる。
【0142】
図１４は、上位モデル２１０の出力値のサンプルの第２例を示す図である。図１４の横軸は、上位モデル２１０の出力値を示す。縦軸は度数を示す。
図１４は、図１３に示されるパラメータ値を上位モデル２１０に適用してシミュレーションを実行した場合の、上位モデル２１０の出力の例を度数分布で示している。図３のステップＳ１０３からＳ１１０のループの２回目の実行におけるステップＳ１０４で、上位シミュレータ部１８３が、このシミュレーションを行う。
図７の場合よりも図１３の場合のほうがパラメータＺの値が収束していることに応じて、図８の場合よりも図１４の場合のほうが、上位モデル２１０の出力が収束している。
図１４は、図１３に示されるパラメータ値を上位モデル２１０に適用してシミュレーションを実行した場合の、上位モデル２１０の出力の例を度数分布で示している。図３のステップＳ１０３からＳ１１０のループの２回目の実行におけるステップＳ１０４で、上位シミュレータ部１８３が、このシミュレーションを行う。
図７の場合よりも図１３の場合のほうがパラメータＺの値が収束していることに応じて、図８の場合よりも図１４の場合のほうが、上位モデル２１０の出力が収束している。
【0143】
図１５は、上位モデル２１０のパラメータ値のサンプルの第３例を示す図である。図１５の横軸は、パラメータｚ１の値を示す。縦軸は、パラメータｚ２の値を示す
図１５は、図１４に示される上位モデル２１０の出力値に基づいて、上位モデル２１０の出力値が実際に近付くようにパラメータ値を推定するデータ同化を行った場合の、上位モデル２１０のパラメータの値を示している。図３のステップＳ１０３からＳ１１０のループの２回目の実行における、ステップＳ１０４からＳ１０５で、上位事後分布データ算出部１８４が、このデータ同化を行う。
図１５は、図１４に示される上位モデル２１０の出力値に基づいて、上位モデル２１０の出力値が実際に近付くようにパラメータ値を推定するデータ同化を行った場合の、上位モデル２１０のパラメータの値を示している。図３のステップＳ１０３からＳ１１０のループの２回目の実行における、ステップＳ１０４からＳ１０５で、上位事後分布データ算出部１８４が、このデータ同化を行う。
【0144】
領域Ａ４３の点の各々は、（ｚ１，ｚ２）のサンプルの例を示している。領域Ａ３３に示されるサンプルは、｛Ｚ＜ＰＲＥＤ＞（１）
ｓ，ｍ｝ｍ＝１
１００と表される。
領域Ａ４１に、領域Ａ４３に示されるサンプルに基づくパラメータｚ１の分布が示されている。領域Ａ４２に、領域Ａ３３に示されるサンプルに基づくパラメータｚ２の分布が示されている。
領域Ａ４１に、領域Ａ４３に示されるサンプルに基づくパラメータｚ１の分布が示されている。領域Ａ４２に、領域Ａ３３に示されるサンプルに基づくパラメータｚ２の分布が示されている。
【0145】
図１５では、図１３の場合よりもさらにパラメータ値が偏在しており、ｚ１の値が９程度、かつ、ｚ２の値が１２程度の領域に点が集中している。
データ算出装置１００が適切なパラメータＺの値の推定に成功したと考えられる。
データ算出装置１００が適切なパラメータＺの値の推定に成功したと考えられる。
【0146】
図１６は、データ算出装置１００による下位モデル２２０のパラメータ値の推定結果の例を示す図である。
図１６は、図１５に示される上位モデル２１０のパラメータ値を下位モデル２２０の出力のターゲット値として、下位モデル２２０の出力値がターゲット値に近付くように下位モデル２２０のパラメータ値を推定するデータ同化を行った場合の、下位モデル２２０のパラメータ値の例を、横軸にパラメータ値をとり縦軸に度数をとった度数分布にて示している。
図１６は、図１５に示される上位モデル２１０のパラメータ値を下位モデル２２０の出力のターゲット値として、下位モデル２２０の出力値がターゲット値に近付くように下位モデル２２０のパラメータ値を推定するデータ同化を行った場合の、下位モデル２２０のパラメータ値の例を、横軸にパラメータ値をとり縦軸に度数をとった度数分布にて示している。
【0147】
図４の処理の２回目の実行におけるステップＳ２０４およびＳ２０５で、下位事後分布データ算出部１８６が、このデータ同化を行う。
図１６では、図１２の場合よりも、パラメータ値がさらに偏在している。データ算出装置１００が、適切なパラメータθ１およびθ２の値の推定に成功したと考えられる。
図１６では、図１２の場合よりも、パラメータ値がさらに偏在している。データ算出装置１００が、適切なパラメータθ１およびθ２の値の推定に成功したと考えられる。
【0148】
モデル２００が、３層以上のモデルとして構成されていてもよい。この場合、最上位層のモデルを上位モデルと称し、最下位層のモデルを下位モデルと称し、それ以外の層のモデルを中間モデルと称する。
モデルの層数が３層以上の場合も、上位モデルに関する処理は、図３を参照して説明したのと同様である。上位モデルと直接データの入出力を行うモデルが中間モデルである場合は、図３の説明で「下位モデル」を中間モデルと読み替える。
モデルの層数が３層以上の場合も、上位モデルに関する処理は、図３を参照して説明したのと同様である。上位モデルと直接データの入出力を行うモデルが中間モデルである場合は、図３の説明で「下位モデル」を中間モデルと読み替える。
【0149】
モデルの層数が３層以上の場合も、下位モデルに関する処理は、図４を参照して説明したのと同様である。下位モデルと直接データの入出力を行うモデルが中間モデルである場合は、図４の説明で「上位モデル」を中間モデルと読み替える。
【0150】
【0151】
図１７は、モデルの層数が３層以上である場合の中間層におけるデータの入出力の例を示す図である。図１７の例で、中間モデル２３２は、上位層モデル２３１および下位層モデル２３３とデータの入出力を行う。上位層モデル２３１は、中間モデル２３２のすぐ上の層のモデルである。上位層モデル２３１は、上位モデルであってもよいし、中間モデルであってもよい。
【0152】
下位層モデル２３３は、中間モデル２３２のすぐ下の層のモデルである。下位層モデル２３３は、下位モデルであってもよいし、中間モデルであってもよい。
１つの上位層モデル２３１とデータの入出力を行う中間モデルの個数は、１つ以上であればよい。１つの中間モデル２３２とデータの入出力を行う下位層モデルの個数は、１つ以上であればよい。
モデルの層数は一律でなくてもよい。例えば、上位モデルのすぐ下の層のモデルとして、中間モデルと下位モデルとの両方があってもよい。中間モデルのすぐ下の層のモデルとして、中間モデルと下位モデルとの両方があってもよい。
１つの上位層モデル２３１とデータの入出力を行う中間モデルの個数は、１つ以上であればよい。１つの中間モデル２３２とデータの入出力を行う下位層モデルの個数は、１つ以上であればよい。
モデルの層数は一律でなくてもよい。例えば、上位モデルのすぐ下の層のモデルとして、中間モデルと下位モデルとの両方があってもよい。中間モデルのすぐ下の層のモデルとして、中間モデルと下位モデルとの両方があってもよい。
【0153】
入力ベクトルについては、データ算出装置１００が、中間モデル２３２への入力となる入力ベクトルＸ（ｕ）、および、下位層モデル２３３への入力となる入力ベクトルＸ（ｄ）を設定する。これは、上述した、データ算出装置１００が、下位モデル２２０－ｓの入力ベクトルＸｓを設定する処理を階層的に行うことに相当する。
【0154】
モデルパラメータについてのサンプルについては、上位層モデル２３１から中間モデル２３２へ出力するサンプルを、Ｚ’＜ＰＯＳＴ＞
（ｕ）と表記する。中間モデル２３２から下位層モデル２３３へ出力するサンプルを、Ｚ’＜ＰＯＳＴ＞
（ｄ）と表記する。
Ｚ’＜ＰＯＳＴ＞ （ｕ）およびＺ’＜ＰＯＳＴ＞ （ｄ）は、図４の処理におけるＺ’＜ＰＯＳＴ＞（ｌ） ｓ、１，・・・，Ｚ’＜ＰＯＳＴ＞（ｌ） ｓ、Ｍに相当する。Ｚ’＜ＰＯＳＴ＞ （ｕ）は、中間モデル２３２におけるデータ同化でターゲットデータとして用いられる。Ｚ’＜ＰＯＳＴ＞ （ｕ）は、下位層モデル２３３におけるデータ同化でターゲットデータとして用いられる。
Ｚ’＜ＰＯＳＴ＞ （ｕ）およびＺ’＜ＰＯＳＴ＞ （ｄ）は、図４の処理におけるＺ’＜ＰＯＳＴ＞（ｌ） ｓ、１，・・・，Ｚ’＜ＰＯＳＴ＞（ｌ） ｓ、Ｍに相当する。Ｚ’＜ＰＯＳＴ＞ （ｕ）は、中間モデル２３２におけるデータ同化でターゲットデータとして用いられる。Ｚ’＜ＰＯＳＴ＞ （ｕ）は、下位層モデル２３３におけるデータ同化でターゲットデータとして用いられる。
【0155】
下位層モデル２３３から中間モデル２３２へ出力するサンプルを、Ｚ＜ＰＲＥＤ＞
（ｄ）と表記する。中間モデル２３２から上位層モデル２３１へ出力するサンプルを、Ｚ＜ＰＲＥＤ＞
（ｕ）と表記する。
Ｚ＜ＰＲＥＤ＞ （ｄ）およびＺ＜ＰＲＥＤ＞ （ｕ）は、図３の処理および図４の処理におけるＺ＜ＰＲＥＤ＞（ｌ） ｓに相当する。Ｚ＜ＰＲＥＤ＞ （ｄ）は、中間モデル２３２における、事後分布に基づくサンプルＺ＜ＰＯＳＴ＞ （ｄ）を更新するために用いられる。ここでいう事後分布に基づくサンプルは、データ同化によって得られるサンプルであり、モデルの出力値のターゲットデータが与えられたときの、事後分布に基づくパラメータのサンプルに該当する。
Ｚ＜ＰＲＥＤ＞ （ｄ）およびＺ＜ＰＲＥＤ＞ （ｕ）は、図３の処理および図４の処理におけるＺ＜ＰＲＥＤ＞（ｌ） ｓに相当する。Ｚ＜ＰＲＥＤ＞ （ｄ）は、中間モデル２３２における、事後分布に基づくサンプルＺ＜ＰＯＳＴ＞ （ｄ）を更新するために用いられる。ここでいう事後分布に基づくサンプルは、データ同化によって得られるサンプルであり、モデルの出力値のターゲットデータが与えられたときの、事後分布に基づくパラメータのサンプルに該当する。
【0156】
上位層モデル２３１が中間モデルである場合、Ｚ＜ＰＲＥＤ＞
（ｕ）は、上位層モデル２３１における、事後分布に基づくサンプルＺ＜ＰＯＳＴ＞
（ｄ）を更新するために用いられる。一方、上位層モデル２３１が上位モデルである場合、Ｚ＜ＰＲＥＤ＞
（ｕ）は、上位層モデル２３１の、次回のループ実行における、事前分布に基づくサンプルを生成するために用いられる。
【0157】
データ算出装置１００は、図４のステップＳ２０１からＳ２０５に相当する処理によって、中間モデル２３２の出力Ｚｕの値がターゲットデータＺ＜ＰＯＳＴ＞
（ｕ）の値に近付くように、パラメータＺ（ｄ）の値を推定するデータ同化を行う。これにより、データ算出装置１００は、サンプルＺ＜ＰＯＳＴ＞
（ｄ）を算出する。
【0158】
そして、データ算出装置１００は、図３のステップＳ１０７からＳ１０９に相当する処理によって、下位層モデル２３３でのデータ同化の結果が反映されたサンプルＺ＜ＰＲＥＤ＞
（ｄ）を算出する。
図３のステップＳ１０９では、データ算出装置１００は、Ｚ（ＰＲＥＤ）（ｌ） ｓを用いてＺ＜ＰＲＩ＞（ｌ＋１）を算出する。これに対し、中間モデル２３２に対する処理では、データ算出装置１００は、Ｚ＜ＰＲＥＤ＞ （ｄ）を用いて事後分布に基づくサンプルＺ＜ＰＯＳＴ＞ （ｄ）を更新する。
図３のステップＳ１０９では、データ算出装置１００は、Ｚ（ＰＲＥＤ）（ｌ） ｓを用いてＺ＜ＰＲＩ＞（ｌ＋１）を算出する。これに対し、中間モデル２３２に対する処理では、データ算出装置１００は、Ｚ＜ＰＲＥＤ＞ （ｄ）を用いて事後分布に基づくサンプルＺ＜ＰＯＳＴ＞ （ｄ）を更新する。
【0159】
そして、データ算出装置１００は、図４のステップＳ２０６からＳ２０７に相当する処理によって、更新されたサンプルＺ＜ＰＯＳＴ＞
（ｄ）を中間モデル２３２に適用してシミュレーションを実行し、Ｚ＜ＰＲＥＤ＞
（ｕ）を算出する。
このように、中間モデル２３２に関する処理では、中間モデル２３２におけるデータ同化の後、下位層モデル２３３におけるデータ同化の結果をサンプルＺ＜ＰＯＳＴ＞ （ｄ）に反映させてから、Ｚ＜ＰＲＥＤ＞ （ｕ）を算出する点で、図４の処理と異なる。それ以外の点では、中間モデルに関する処理は、図４の処理と同様である。
このように、中間モデル２３２に関する処理では、中間モデル２３２におけるデータ同化の後、下位層モデル２３３におけるデータ同化の結果をサンプルＺ＜ＰＯＳＴ＞ （ｄ）に反映させてから、Ｚ＜ＰＲＥＤ＞ （ｕ）を算出する点で、図４の処理と異なる。それ以外の点では、中間モデルに関する処理は、図４の処理と同様である。
【0160】
以上のように、上位事後分布データ算出部１８４は、上位モデル２１０の出力に対するターゲットデータに基づく、上位モデル２１０のパラメータ値の事後分布を示すデータ｛Ｚ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍを算出する。下位事後分布データ算出部１８６は、下位モデル２２０について、上位モデル２１０のパラメータ値の事後分布を示すデータ｛Ｚ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍに基づく、下位モデル２２０のパラメータ値の事後分布を示すデータ｛θ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍを算出する。
【0161】
下位事後分布データ算出部１８６が、上位モデル２１０のパラメータ値を用いて下位モデル２２０のパラメータ値の事後分布を示すデータを算出することで、データ同化を階層的に行うことができる。具体的には、下位事後分布データ算出部１８６が、上位モデル２１０でのパラメータ値の調整後のサンプルＺ＜ＰＯＳＴ＞に基づいて決定される、下位モデル２２０の出力のターゲット値（サンプルＺ＜ＰＯＳＴ＞
ｓ）を用いることで、上位モデル２１０でのパラメータ値の調整結果を下位モデル２２０でのパラメータ値の調整に反映させることができる。
【0162】
パラメータベクトルＺの次元数が、パラメータベクトルθ１からθＳまでの次元数の合計よりも小さいことで、対象９００を１つのモデルのみで模擬する場合よりも、個々のモデルのパラメータの個数を少なくすることができる。この点で、データ算出装置１００によれば比較的短時間でデータ同化を行うことができ、パラメータベクトルθ１からθＳまでの値を推定できると期待される。特に、データ算出装置１００によれば、モデルのパラメータの個数が多いことで現実的な時間内にパラメータの推定値を得られないことを回避できると期待される。
【0163】
また、下位モデル２２０同士間では直接的なデータのやり取りは無く、下位モデル２２０の各々は並列処理可能である。データ算出装置１００によれば、この点でも比較的短時間でデータ同化を行えると期待される。
【0164】
また、下位シミュレータ部１８５は、下位モデル２２０のパラメータ値の事後分布を示すデータ｛θ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍに基づくパラメータ値θ＜ＰＯＳＴ＞を下位モデル２２０に適用してシミュレーションを実行する。上位シミュレータ部１８３は、下位モデル２２０を用いたシミュレーションの結果｛Ｚ＜ＰＲＥＤ＞（ｌ）
ｓ，ｍ｝ｍ＝１
Ｍに基づく上位モデル２１０のパラメータ値のサンプルデータ｛Ｚ＜ＰＲＩ＞（ｌ＋１）
ｍ｝ｍ＝１
Ｍを上位モデル２１０に適用してシミュレーションを実行する。上位事後分布データ算出部１８４は、上位モデル２１０を用いたシミュレーションの結果として得られる上位モデル２１０の出力のサンプルデータ｛Ｙ＜ＰＲＩ＞（ｌ＋１）
ｍ｝ｍ＝１
Ｍとターゲットデータとの類似度に基づいて、上位モデル２１０の出力に対するターゲットデータに基づく、上位モデル２１０のパラメータ値の事後分布を示すデータ｛Ｚ＜ＰＯＳＴ＞（ｌ＋１）
ｍ｝ｍ＝１
Ｍを算出する。
【0165】
このように、データ算出装置１００では、下位モデル２２０におけるデータ同化の結果を上位モデル２１０のパラメータ値に反映させることができ、これにより、階層的なデータ同化を繰り返すことができる。データ算出装置１００が階層的なデータ同化を繰り返すことで、パラメータ値を高精度に推定できると期待される。特に、パラメータベクトルＺの確率分布、および、パラメータベクトルθ１からθＳまでのそれぞれの確率分布が、何れもターゲット値に応じたあるパラメータ値に収束していくことが期待される。上述したように、確率分布があるパラメータ値に収束するとは、その値における確率が大きくなることである。
【0166】
例えば、データ同化によって、出力ベクトルＹの値がターゲット値に近い値となるパラメータＺの値の確率が比較的大きくなることで、データ同化を再度行う際に、出力ベクトルＹのサンプリング値としてターゲット値に近いものが多くなる。これにより、再度のデータ同化によって、出力ベクトルＹの値がターゲット値に近い値となるパラメータＺの値の確率がさらに大きくなると期待される。
パラメータベクトルθ１からθＳまでについても、データ同化を繰り返すことで、ターゲット値に応じたパラメータ値の確率が大きくなっていくと期待される。
パラメータベクトルθ１からθＳまでについても、データ同化を繰り返すことで、ターゲット値に応じたパラメータ値の確率が大きくなっていくと期待される。
【0167】
また、上位シミュレータ部１８３によるシミュレーションの実行で、上位モデル２１０の出力｛Ｙ＜ＰＲＩ＞（ｌ＋１）
ｍ｝ｍ＝１
Ｍが、下位モデル２２０でのデータ同化の結果に応じて特定の値の近傍に集中することが考えられる。上位事後分布データ算出部１８４が、この出力｛Ｙ＜ＰＲＩ＞（ｌ＋１）
ｍ｝ｍ＝１
Ｍを用いてデータ同化を行うことで、上位モデル２１０の出力の特定の値の近傍について集中的にデータ同化の解となるパラメータ値を探索する。データ算出装置１００によれば、この点で、単一階層のモデルを用いてパラメータ空間を均一的に探索する場合よりも効率的に、パラメータ値を推定できると期待される。
【0168】
また、上位事後分布データ算出部１８４は、上位モデル２１０のパラメータ値のサンプルデータのカーネル平均に基づいてカーネルハーディングを実行し、上位モデル２１０のパラメータ値の事後分布を示すデータ｛Ｚ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍを算出する。
【0169】
このように、上位事後分布データ算出部１８４がカーネル平均に基づいて、上位モデル２１０のパラメータ値の事後分布を示すデータを算出することで、与えられるＭ個の事前分布のサンプルに対して、同数であるＭ個の事後分布のサンプルを算出することができる。データ算出装置１００によれば、この点で、パラメータ値の推定を効率的に行うことができる。
【0170】
また、下位事後分布データ算出部１８６は、下位モデル２２０のパラメータ値のサンプルデータのカーネル平均に基づいてカーネルハーディングを実行し、下位モデルのパラメータ値の事後分布を示すデータ｛θ＜ＰＯＳＴ＞（ｌ）
ｍ｝ｍ＝１
Ｍを算出する。
このように、下位事後分布データ算出部１８６がカーネル平均に基づいて、下位モデル２２０のパラメータ値の事後分布を示すデータを算出することで、与えられるＭ個の事前分布のサンプルに対して、同数であるＭ個の事後分布のサンプルを算出することができる。データ算出装置１００によれば、この点で、パラメータ値の推定を効率的に行うことができる。
このように、下位事後分布データ算出部１８６がカーネル平均に基づいて、下位モデル２２０のパラメータ値の事後分布を示すデータを算出することで、与えられるＭ個の事前分布のサンプルに対して、同数であるＭ個の事後分布のサンプルを算出することができる。データ算出装置１００によれば、この点で、パラメータ値の推定を効率的に行うことができる。
【0171】
また、制御部１８０は、下位事後分布データ算出部１８６が算出する下位モデル２２０のパラメータ値の事後分布を示すデータに基づいて、対象９００を制御する。
上記のように、データ算出装置１００が階層的なモデルを用いることで、個々のモデルのパラメータの個数が比較的少ない。これにより、対象９００が比較的複雑な場合でも、下位事後分布データ算出部１８６は、下位モデル２２０のパラメータ値を推定することができる。制御部１８０は、下位事後分布データ算出部によるパラメータ値の推定結果を参照することで、対象９００の状態を把握して対象９００を比較的高精度に制御することができる。
上記のように、データ算出装置１００が階層的なモデルを用いることで、個々のモデルのパラメータの個数が比較的少ない。これにより、対象９００が比較的複雑な場合でも、下位事後分布データ算出部１８６は、下位モデル２２０のパラメータ値を推定することができる。制御部１８０は、下位事後分布データ算出部によるパラメータ値の推定結果を参照することで、対象９００の状態を把握して対象９００を比較的高精度に制御することができる。
【0172】
図１８は、実施形態に係るデータ算出装置の構成例を示す図である。図１８に示す構成で、データ算出装置６１０は、上位事後分布データ算出部６１１と、下位事後分布データ算出部６１２とを備える。
かかる構成で、上位事後分布データ算出部６１１は、対象を模擬する上位モデルの出力に対するターゲットデータに基づく、上位モデルのパラメータ値の事後分布を示すデータを算出する。下位事後分布データ算出部６１２は、対象の部分を模擬して上位モデルのパラメータ値を出力する下位モデルについて、上位モデルのパラメータ値の事後分布を示すデータに基づく、下位モデルのパラメータ値の事後分布を示すデータを算出する。
上位事後分布データ算出部６１１は、上位事後分布データ算出手段の例に該当する。下位事後分布データ算出部６１２は、下位事後分布データ算出手段の例に該当する。
かかる構成で、上位事後分布データ算出部６１１は、対象を模擬する上位モデルの出力に対するターゲットデータに基づく、上位モデルのパラメータ値の事後分布を示すデータを算出する。下位事後分布データ算出部６１２は、対象の部分を模擬して上位モデルのパラメータ値を出力する下位モデルについて、上位モデルのパラメータ値の事後分布を示すデータに基づく、下位モデルのパラメータ値の事後分布を示すデータを算出する。
上位事後分布データ算出部６１１は、上位事後分布データ算出手段の例に該当する。下位事後分布データ算出部６１２は、下位事後分布データ算出手段の例に該当する。
【0173】
下位事後分布データ算出部６１２が、上位モデルのパラメータ値を用いて下位モデルのパラメータ値の事後分布を示すデータを算出することで、データ同化を階層的に行うことができる。具体的には、下位事後分布データ算出部６１２が、上位モデルでのデータ同化後のパラメータ値のサンプルに基づいて決定される、下位モデルの出力のターゲット値を用いることで、上位モデルでのデータ同化の結果を下位モデルでのデータ同化に反映させることができる。
【0174】
上位モデルのパラメータの個数が、下位モデルのパラメータの個数の合計よりも少ないことで、対象を１つのモデルのみで模擬する場合よりも、個々のモデルのパラメータの個数を少なくすることができる。この点で、データ算出装置６１０によれば比較的短時間でデータ同化を行うことができ、パラメータの値を推定できると期待される。特に、データ算出装置６１０によれば、モデルのパラメータの個数が多いことで現実的な時間内にパラメータの推定値を得られないことを回避できると期待される。
【0175】
また、下位モデル同士間では直接的なデータのやり取りは無く、下位モデルの各々は並列処理可能である。データ算出装置６１０によれば、この点でも比較的短時間でデータ同化を行えると期待される。
【0176】
上位事後分布データ算出部６１１は、例えば、図１に示される上位事後分布データ算出部１８４等の機能を用いて実現することができる。下位事後分布データ算出部６１２は、例えば、図１に示される下位事後分布データ算出部１８６等の機能を用いて実現することができる。
【0177】
図１９は、実施形態に係るデータ算出装置の構成のもう１つの例を示す図である。図１９に示す構成で、データ算出装置６２０は、下位モデルシミュレーション実行部６２１と、上位モデルシミュレーション実行部６２２とを備える。
かかる構成で、下位モデルシミュレーション実行部６２１は、対象を模擬する上位モデルのパラメータ値を出力する下位モデルに、下位モデルのパラメータ値の分布を示すデータを適用してシミュレーションを実行し、上位モデルのパラメータ値の分布を示すデータを算出する。上位モデルシミュレーション実行部６２２は、上位モデルに、算出したパラメータ値の分布を示すデータを適用して対象のシミュレーションを実行する。
かかる構成で、下位モデルシミュレーション実行部６２１は、対象を模擬する上位モデルのパラメータ値を出力する下位モデルに、下位モデルのパラメータ値の分布を示すデータを適用してシミュレーションを実行し、上位モデルのパラメータ値の分布を示すデータを算出する。上位モデルシミュレーション実行部６２２は、上位モデルに、算出したパラメータ値の分布を示すデータを適用して対象のシミュレーションを実行する。
【0178】
下位モデルシミュレーション実行部６２１は、下位モデルシミュレーション実行手段の例に該当する。上位モデルシミュレーション実行部６２２は、上位モデルシミュレーション実行手段の例に該当する。
データ算出装置６２０によれば、下位モデルを用いたパラメータ値の推定結果を上位モデルに反映させることができ、この点で、対象のシミュレーションを高精度に実行することができる。
データ算出装置６２０によれば、下位モデルを用いたパラメータ値の推定結果を上位モデルに反映させることができ、この点で、対象のシミュレーションを高精度に実行することができる。
【0179】
下位モデルシミュレーション実行部６２１は、例えば、図１に示される下位シミュレータ部１８５等の機能を用いて実現することができる。上位モデルシミュレーション実行部６２２は、例えば、図１に示される上位シミュレータ部１８３等の機能を用いて実現することができる。
【0180】
図２０は、実施形態に係るデータ算出方法における処理手順の例を示す図である。図２０に示すデータ算出方法は、上位モデルについて事後分布を示すデータを算出する工程（ステップＳ６１１）と、下位モデルについて事後分布を示すデータを算出する工程（ステップＳ６１２）とを含む。
【0181】
上位モデルについて事後分布を示すデータを算出する工程（ステップＳ６１１では、対象を模擬する上位モデルの出力に対するターゲットデータに基づく、上位モデルのパラメータ値の事後分布を示すデータを算出する。
下位モデルについて事後分布を示すデータを算出する工程（ステップＳ６１２）では、対象の部分を模擬して上位モデルのパラメータ値を出力する下位モデルについて、上位モデルのパラメータ値の事後分布を示すデータに基づく、下位モデルのパラメータ値の事後分布を示すデータを算出する。
下位モデルについて事後分布を示すデータを算出する工程（ステップＳ６１２）では、対象の部分を模擬して上位モデルのパラメータ値を出力する下位モデルについて、上位モデルのパラメータ値の事後分布を示すデータに基づく、下位モデルのパラメータ値の事後分布を示すデータを算出する。
【0182】
図２０に示すデータ算出方法によれば、上位モデルのパラメータ値を用いて下位モデルのパラメータ値の事後分布を示すデータを算出することで、データ同化を階層的に行うことができる。具体的には、上位モデルでのデータ同化後のパラメータ値のサンプルに基づいて決定される、下位モデルの出力のターゲット値を用いることで、上位モデルでのデータ同化の結果を下位モデルでのデータ同化に反映させることができる。
【0183】
上位モデルのパラメータの個数が、下位モデルのパラメータの個数の合計よりも少ないことで、対象を１つのモデルのみで模擬する場合よりも、個々のモデルのパラメータの個数を少なくすることができる。図２０に示すデータ算出方法によれば、この点で、比較的短時間でデータ同化を行うことができ、パラメータの値を推定できると期待される。特に、図２０に示すデータ算出方法によれば、モデルのパラメータの個数が多いことで現実的な時間内にパラメータの推定値を得られないことを回避できると期待される。
【0184】
また、下位モデル同士間では直接的なデータのやり取りは無く、下位モデルの各々は並列処理可能である。図２０に示すデータ算出方法によれば、この点でも比較的短時間でデータ同化を行えると期待される。
【0185】
図２１は、少なくとも１つの実施形態に係るコンピュータの構成を示す概略ブロック図である。
図２１に示す構成で、コンピュータ７００は、ＣＰＵ７１０と、主記憶装置７２０と、補助記憶装置７３０と、インタフェース７４０と、不揮発性記録媒体７５０とを備える。
図２１に示す構成で、コンピュータ７００は、ＣＰＵ７１０と、主記憶装置７２０と、補助記憶装置７３０と、インタフェース７４０と、不揮発性記録媒体７５０とを備える。
【0186】
上記のデータ算出装置１００およびデータ算出装置６１０のうち何れか１つ以上またはその一部が、コンピュータ７００に実装されてもよい。その場合、上述した各処理部の動作は、プログラムの形式で補助記憶装置７３０に記憶されている。ＣＰＵ７１０は、プログラムを補助記憶装置７３０から読み出して主記憶装置７２０に展開し、当該プログラムに従って上記処理を実行する。また、ＣＰＵ７１０は、プログラムに従って、上述した各記憶部に対応する記憶領域を主記憶装置７２０に確保する。各装置と他の装置との通信は、インタフェース７４０が通信機能を有し、ＣＰＵ７１０の制御に従って通信を行うことで実行される。また、インタフェース７４０は、不揮発性記録媒体７５０用のポートを有し、不揮発性記録媒体７５０からの情報の読出、および、不揮発性記録媒体７５０への情報の書込を行う。
【0187】
データ算出装置１００がコンピュータ７００に実装される場合、制御部１８０およびその各部の動作は、プログラムの形式で補助記憶装置７３０に記憶されている。ＣＰＵ７１０は、プログラムを補助記憶装置７３０から読み出して主記憶装置７２０に展開し、当該プログラムに従って上記処理を実行する。
【0188】
また、ＣＰＵ７１０は、プログラムに従って、記憶部１７０に対応する記憶領域を主記憶装置７２０に確保する。
通信部１１０による他の装置との通信は、インタフェース７４０が通信機能を有し、ＣＰＵ７１０の制御に従って動作することで実行される。
表示部１２０による表示は、インタフェース７４０が表示装置を有し、ＣＰＵ７１０の制御に従って各種画像を表示することで実行される。
操作入力部１３０によるユーザ操作の受け付けは、インタフェース７４０が例えばキーボードおよびマウスなどの入力デバイスを有してユーザ操作を受け付け、受け付けたユーザ操作を示す情報をＣＰＵ７１０へ出力することで実行される。
通信部１１０による他の装置との通信は、インタフェース７４０が通信機能を有し、ＣＰＵ７１０の制御に従って動作することで実行される。
表示部１２０による表示は、インタフェース７４０が表示装置を有し、ＣＰＵ７１０の制御に従って各種画像を表示することで実行される。
操作入力部１３０によるユーザ操作の受け付けは、インタフェース７４０が例えばキーボードおよびマウスなどの入力デバイスを有してユーザ操作を受け付け、受け付けたユーザ操作を示す情報をＣＰＵ７１０へ出力することで実行される。
【0189】
データ算出装置６１０がコンピュータ７００に実装される場合、上位事後分布データ算出部６１１、および、下位事後分布データ算出部６１２の動作は、プログラムの形式で補助記憶装置７３０に記憶されている。ＣＰＵ７１０は、プログラムを補助記憶装置７３０から読み出して主記憶装置７２０に展開し、当該プログラムに従って上記処理を実行する。
【0190】
また、ＣＰＵ７１０は、プログラムに従って、データ算出装置６１０が行う処理のための記憶領域を主記憶装置７２０に確保する。
データ算出装置６１０と他の装置との通信は、インタフェース７４０が通信機能を有し、ＣＰＵ７１０の制御に従って動作することで実行される。
データ算出装置６１０とユーザとのインタラクションは、インタフェース７４０が入力デバイスおよび出力デバイスを有し、ＣＰＵ７１０の制御に従って出力デバイスにて情報をユーザに提示し、入力デバイスにてユーザ操作を受け付けることで実行される。
データ算出装置６１０と他の装置との通信は、インタフェース７４０が通信機能を有し、ＣＰＵ７１０の制御に従って動作することで実行される。
データ算出装置６１０とユーザとのインタラクションは、インタフェース７４０が入力デバイスおよび出力デバイスを有し、ＣＰＵ７１０の制御に従って出力デバイスにて情報をユーザに提示し、入力デバイスにてユーザ操作を受け付けることで実行される。
【0191】
上述したプログラムのうち何れか１つ以上が不揮発性記録媒体７５０に記録されていてもよい。この場合、インタフェース７４０が不揮発性記録媒体７５０からプログラムを読み出すようにしてもよい。そして、ＣＰＵ７１０が、インタフェース７４０が読み出したプログラムを直接実行するか、あるいは、主記憶装置７２０または補助記憶装置７３０に一旦保存して実行するようにしてもよい。
【0192】
なお、データ算出装置１００、および、データ算出装置６１０が行う処理の全部または一部を実行するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各部の処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ（Read Only Memory）、ＣＤ－ＲＯＭ（Compact Disc Read Only Memory）等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ（Read Only Memory）、ＣＤ－ＲＯＭ（Compact Disc Read Only Memory）等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
【0193】
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【産業上の利用可能性】
【0194】
本発明の実施形態は、データ算出装置、データ算出方法および記録媒体に適用してもよい。
【符号の説明】
【0195】
１００、６１０ データ算出装置
１１０ 通信部
１２０ 表示部
１３０ 操作入力部
１７０ 記憶部
１８０ 制御部
１８１ 初期設定部
１８２ ループ処理部
１８３ 上位シミュレータ部
１８４、６１１ 上位事後分布データ算出部
１８５ 下位シミュレータ部
１８６、６１２ 下位事後分布データ算出部
１１０ 通信部
１２０ 表示部
１３０ 操作入力部
１７０ 記憶部
１８０ 制御部
１８１ 初期設定部
１８２ ループ処理部
１８３ 上位シミュレータ部
１８４、６１１ 上位事後分布データ算出部
１８５ 下位シミュレータ部
１８６、６１２ 下位事後分布データ算出部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】