本研究課題では,ノンパラメトリック・ベイズ手法に基づく時系列データ解析手法の開発を目指し, ディリクレ過程混合ARMAモデルによる時系列データのクラスタリング法を提案している.先行研究においては,有限混合ARMAモデルを用いたクラスタリング手法が提案されている.有限混合ARMAモデルに比べ,ディリクレ過程混合ARMAモデルを用いると,クラスタ数をあらかじめ定める必要がないことや,一貫したベイズアプローチであるため,パラメータの不確かさを推定に適切に反映する手法となることなどが期待される. このディリクレ過程混合ARMAモデルのためのサンプリング法を,最も簡単なギブス・サンプリングに基づいて開発した.ARMAモデルのパラメータには共役事前分布が存在しないため,直接的なアプローチでは,メトロポリス・ヘイスティングス(MH)法を基本とするサンプリング法を考えることになる.そして,その場合,複数の時系列データに対するARMAパラメータのサンプリングをMH法により行うことになるが,その事後分布は孤立傾向のある高い多数の峰を持つことが予想されるため,それには困難が伴うだろう.その問題を回避するために本研究では,ギブス・サンプリングを基本とするサンプリングアルゴリズムを開発した. ギブス・サンプリングを基本とする場合には,(1) 時系列データの周辺尤度の評価法と,(2) 単一の時系列データに対するARMAパラメータの事後分布からのサンプリング法が課題となる.本研究では,(1) に対して Annealed Importance Sampling (Neal 2001)が応用可能であること,(2) に対しては,単一時系列データに対するMH法によるシミュレーションをデータ数だけ並列的に行うことで効率のよいサンプリングが実現できることを明らかにした.
|