研究課題
本研究の主眼は,高頻度観測されるジャンプ型確率微分方程式モデルの漸近推測理論の基盤構築である.本年度は主として以下の結果を得た.(1)擬似安定確率微分方程式モデル(SSDE)の一種の擬似最尤推定量(安定近似型推定量,SQMLE)の漸近分布論の基礎構築について,尺度係数関数が定数,かつ擬似安定指数βが未知の場合に,βの√n-収束率を持つ推定量(構成法は知られている)をプラグインした場合の誤差評価を通じ,SQMLEの漸近分布(尺度混合正規分布)が影響を受けないための十分条件を与えた.提案手法では,散過程の場合には理論上不可能なトレンド係数の一致推定も可能である.(2)SSDEの特殊な場合である擬似安定レヴィ過程モデルの場合に局所漸近正規性を導出し,系として特に(1)で筆者が提案しているSQMLEの漸近有効性が得られた.本結果の証明では,尤度関数の解析評価においてポアソン空間上の特殊なマリアヴァン解析(部分積分の公式)を用いたため,確率解析の数理統計への応用についても新たな知見を与えた.(3)エルゴード型ジャンプ型確率微分方程式モデルの係数パラメータ,およびレヴィ測度の汎関数パラメータの二段階推定量を構成し,その漸近分布を導出して近似信頼領域の構成を行った.本結果は,二段階性によるパラメータの次元削減を通じての推定精度安定化へつながるだけでなく,無限次元パラメータであるレヴィ測度のモデル誤特定に対する頑健性も有する.当初の本年度計画には(1),(2)の実施のみ記載していたが,(3)も並行して行い,研究成果が得られた.(1),(3)については現在執筆中で,平成27年度前半には投稿できる見込みである.(2)については2014年11月にALEA journalへ投稿後,採択に好意的な査読レポートを受け取っており,目下改訂中である.
2: おおむね順調に進展している
大規模高頻度従属データ解析に際し,本研究課題の主眼である微小時間安定分布近似に基づいた(非正規)安定型擬似最尤推定方式の基本的原理が固まりつつある.現在,近似安定分布の指数を既知とした場合の擬似最尤推定量の漸近分布を特定し,それが正規型の漸近有効な場合を拡張していることまで解明できており,目下,安定指数が未知の場合の取り扱いについて考察を詰めている.ノイズの対称性や係数の誤特定の可能性など,状況を子細に分類すると先が長そうなので,基本的な設定で一旦論文にまとめて投稿するための準備を進めている状況である.擬似安定レヴィ過程の場合について,非常に一般的な設定で局所漸近正規性を導出できた.例えば,擬似コーシーレヴィ過程(例.Normal inverse Gaussianレヴィ過程)に対してコーシー擬似尤度が漸近最適推定量を与えることが理論的に示されたことになる.レヴィ測度が無限次元パラメータである中で,高頻度データがより大きい情報を与える有限次元パラメータを部分的に推定可能とできたことになり,モデル頑健推定といった方法論構築やSSDEの多段階推定方式の理論構成など,今後多方面への拡張の足掛かりになると期待している.擬似安定過程とは限らないより広いモデル族に対し,正規型擬似尤度とマルチンゲール推定関数型モーメントフィッティングを併用することで,明示的な二段階推定方式を開発し,その漸近分布を示した.一般にこの推定方式は上述の安定型擬似最尤推定量と比べて漸近効率が低いが,計算負荷も低く実装容易であるため,計算機上に比較的容易に構築可能であるという実用上の利点を持つ.
平成27年度半ばまでに,【研究実績の概要】(1)および(3)に記載した研究成果をまとめあげ,投稿する.これらの内容については引き続き,当初の研究計画に記載した方向の改良・拡張を研究する.特に,現在モデルの係数関数形が正しく特定されている状況のみを扱っているが,誤特定された場合には漸近挙動が本質的に変わる可能性があるため,別で精査する.さらに,対応する統計的確率場の裾確率の一様評価も確保し,統計的予測やモデル評価(情報量規準の構成)のための理論基盤の整備をより強固なものとすることを目指す.また(3)については,現在仮定しているサンプリングデザインの緩和や,多段階適合型推定方式への拡張可能性を調べ,理論のみならず,数値計算的な面においても現時点で得られている結果の洗練化を図る.また【研究実績の概要】(2)については,当初の研究計画の通り,非線形状態依存係数関数を持つSSDEの場合において局所漸近混合正規性を導出し,(1)で提案している安定型擬似最尤推定方式の漸近有効性の証明を目指す.おそらくは(1)で筆者が提案しているSQMLEが,この場合にも漸近有効であろうと予想している.さらに,確率微分方程式モデルのオンライン推定方式を開発し,その漸近挙動の解明を目指す.これは大規模従属データからの逐次情報搾取や,推定量のプロットを通じたモデル構造変化のトラッキングに役立つ.
前年度半ばに,今年度の国際会議への参加およびオーガナイザー業務が複数確定したことに加え,本研究課題の一部であるAlexey Kulik氏との共同研究推進にさらなる渡航費もしくは招聘費用が必要となる可能性が高いため,特に渡航費が当初の予定より多めに必要となった.このため研究費の繰越しを行った.すでに今年度渡航先での討論予定も立っており,これらは本研究推進において重要である.
国際研究集会Dynstoch meeting 2015 (Lund, Sweden, 5月末),国際会議International Statistical Institute (ISI) 2015(Rio de Janeiro, Brazil),国際会議the 9th Conference of the Asian Regional Section of the IASC (IASC-ARS 2015)(Singapore)において講演,オーガナイザー業務を行い,本研究に関わる討論を行うことで研究推進を図る.また,本研究課題の一部であるAlexey Kulik氏との共同研究推進のための旅費に充てる.
すべて 2015 2014
すべて 雑誌論文 (4件) (うち査読あり 3件、 オープンアクセス 2件、 謝辞記載あり 2件) 学会発表 (5件) (うち招待講演 2件)
Lecture Notes in Mathematics
巻: 2128 ページ: 179-286
10.1007/978-3-319-12373-8_3
日本統計学会誌和文誌
巻: 44 ページ: 471-495
Journal of Statistical Software
巻: 57 ページ: 1-51
A Mathematical Approach to Research Problems of Science and Technology Mathematics for Industry
巻: 5 ページ: 219-238
10.1007/978-4-431-55060-0_17