[論文レビュー] Exploring recursion for optimal estimators under cascade rotation

Jan Kowalski, Jacek Wesołowski|arXiv (Cornell University)|Dec 2, 2014

Statistical and numerical algorithms参考文献 14被引用数 3

ひとこと要約

本稿では、欠けを含むリターン型の回転パターンを有する繰り返し調査における母平均の最適線形推定量の解析的再帰公式を提示する。長年の課題であったこの問題を、代数的作用素アプローチを用いて解決し、再帰の次数が回転パターンにおける最大のギャップのサイズに1を加えたものに等しいことを示した。問題は、チェビシェフ多項式に関連する多項式の根を特定し、行列のランクを解析することに還元され、再帰の次数は回転パターンにおける最大ギャップのサイズに1を加えたものに等しいことが示された。

ABSTRACT

We are concerned with optimal linear estimation of means on subsequent occasions under sample rotation where evolution of samples in time is designed through a cascade pattern. It has been known since the seminal paper of Patterson (1950) that when the units are not allowed to return to the sample after leaving it for certain period (there are no gaps in the rotation pattern), one step recursion for optimal estimator holds. However, in some important real surveys, e.g. Current Population Survey in the US or Labour Force Survey in many countries in Europe, units return to the sample after being absent in the sample for several occasions (there are gaps in rotation patterns). In such situations difficulty of the question of the form of the recurrence for optimal estimator increases drastically. This issue has not been resolved yet. Instead alternative sub-optimal approaches were developed, as K-composite estimation (see e.g. Hansen et al. (1955)), AK-composite estimation (see e.g. Gurney and Daly (1965) or time series approach (see e.g. Binder and Hidiroglou (1988)). In the present paper we overcome this long-standing difficulty, that is, we present analytical recursion formulas for the optimal linear estimator of the mean for schemes with gaps in rotation patterns. It is achieved under some technical conditions: ASSUMPTION I and ASSUMPTION II (numerical experiments suggest that these assumptions might be universally satisfied). To attain the goal we develop an algebraic operator approach which allows to reduce the problem of recursion for the optimal linear estimator to two issues: (1) localization of roots (possibly complex) of a polynomial Q_p defined in terms of the rotation pattern (Q_p happens to be conveniently expressed through Chebyshev polynomials of the first kind), (2) rank of a matrix S defined in terms of the rotation pattern and the roots of the polynomial Q_p.

研究の動機と目的

• 回転パターンにギャップを含む繰り返し調査における最良線形不偏推定量（BLUE）の再帰的公式を導出するという、長年の未解決問題を解決すること。
• 元々非再入単位にのみ有効であったパターソン（1950）の1ステップ再帰を、欠けの後再入する単位を含む場合に拡張すること。
• 推定再帰問題を多項式の根の局所化と行列ランク解析に還元する一般化された代数的フレームワークを構築すること。
• 再帰の次数が回転パターンにおける最大ギャップのサイズによって決定されることを確立すること。
• 技術的仮定の下で再帰係数の正確な公式を提供し、数値的証拠からこれらの仮定が実世界の調査設計において広く満たされている可能性を示すこと。

提案手法

• 最適線形推定量を過去の観測値と過去の推定量の再帰的線形結合として定式化し、パターソンの1ステップ再帰を一般化する。
• 時間発展をモデル化するため、シフト作用素と反転作用素を用いた代数的作用素フレームワークを導入する。
• チェビシェフ多項式の第一種を用いて表現可能な多項式 $ Q_p(x) $ を定義し、その根が再帰の構造を決定することを示す。
• 再帰係数を求める問題を2つの核心的タスクに還元する：(1) $ Q_p $ の根の特定、(2) 回転パターンとこれらの根から導かれる行列 $ S $ のランクの計算。
• 対角行列および相関行列に関連する行列の逆行列とトレースの恒等式を適用し、推定量の分散の式を導出する。
• 作用素方程式から導かれる制約を用いて、得られた推定量が不偏かつ最適であることを検証する。

実験結果

リサーチクエスチョン

• RQ1 ギャップを含む回転スキーム（単位が欠けた後再入する）において、最良線形不偏推定量（BLUE）の再帰的公式を導出できるか？
• RQ2 このようなスキームにおける再帰の次数は何かが決定するか—特に、回転パターンにおける最大ギャップのサイズに関連しているか？
• RQ3 再帰係数は回転パターンおよび相関構造を用いて解析的にどのように表現できるか？
• RQ4 代数的作用素アプローチが、BLUEの有効かつ最適な再帰を導くための条件は何か？
• RQ5 回転パターンから導かれる多項式 $ Q_p $ を、チェビシェフ多項式と系統的に関連づけることで、根の局所化を簡略化できるか？

主な発見

• 最適線形推定量の再帰の次数は、回転パターンにおける最大ギャップのサイズに1を加えたものに正確に等しい。
• 再帰係数は、第一種チェビシェフ多項式で表現可能な多項式 $ Q_p $ の根から導かれる公式を用いて正確に計算可能である。
• 最適性は2つの技術的仮定（仮定IおよびII）のもとで保証され、数値実験からこれらが実世界の調査設計において広く満たされている可能性が示唆されている。
• 最適推定量の分散は、$ Q_p $ の根と行列 $ S $ を含む方程式系の解から得られる $ heta_{0,0} $ に等しいことが示され、推定量の最小分散が確認された。
• 本手法により、パターソン（1950）の1ステップ再帰がギャップを含む場合にまで一般化され、従来は得られなかった完全な解析的解が得られた。
• 実世界の回転スキームの例を通じて、理論的枠組みが現在の人口調査および労働力調査などの調査に実際に適用可能であることが確認された。