[論文レビュー] Local Asymptotics of P-splines

Luo Xiao, Yingxing Li+2 more|arXiv (Cornell University)|Jan 3, 2012

Advanced Numerical Analysis Techniques参考文献 18被引用数 12

ひとこと要約

本稿は、一般の順序のBスプラインおよび差分罰則を用いたPスプラインの局所漸近理論を確立し、罰則の順序 $m$ のみに依存する同等核を持つ核推定器としてPスプラインが漸近的に振る舞うことを示している。罰則の順序 $m$ に依存するが、スプライン次数 $p$ には依存しない。漸近的バイアスおよび分散の明示的表現を導出し、十分に速い速度で発散する場合、 Knot の数の発散レートが漸近的分布に影響しないことを確認した。また、平滑化パラメータ $ \lambda$ の最適収束レートを特定した。結果は、$p$-$m$ パrameter空間全体に一般化され、Pスプライン推論の理論的基盤を提供する。

ABSTRACT

This report studies local asymptotics of P-splines with $p$th degree B-splines and a $m$th order difference penalty. Earlier work with $p$ and $m$ restricted is extended to the general case. Asymptotically, penalized splines are kernel estimators with equivalent kernels depending on $m$, but not on $p$. A central limit theorem provides simple expressions for the asymptotic mean and variance. Provided it is fast enough, the divergence rate of the number of knots does not affect the asymptotic distribution. The optimal convergence rate of the penalty parameter is given.

研究の動機と目的

• スプライン次数 $p$ および罰則順序 $m$ の制限されたケースを超えて、一般の $p$-$m$ 設定におけるPスプラインの漸近理論を拡張すること。
• Pスプラインが漸近的に、スプライン次数 $p$ に依存せず、罰則順序 $m$ のみに依存する同等核を持つ核推定器として振る舞うことを確立すること。
• Pスプライン推定量の漸近的平均および分散の明示的表現を、内部点および境界点において導出すること。
• Knot の数 $K$ の発散レートが十分に速い限り、漸近的分布に影響しないことを確認すること。
• 漸近的正規性を保証し、平均二乗誤差を最小化するための平滑化パラメータ $\lambda$ の最適収束レートを特定すること。

提案手法

• Pスプラインスムーザーの漸近的挙動を分析するために同等核法を用い、スムーザ行列 $\mathbf{\Lambda}$ を重みなし成分と罰則付き成分の和として表現する。
• Pスプラインの同等核 $g(x - x_i)/h_n$ を導出し、Bスプライン基底関数の性質と漸近展開を用いて、これがスプライン次数 $p$ に依存せず、罰則順序 $m$ のみに依存することを示す。
• 中心極限定理を適用して $\hat{\mu}(x)$ の漸近的分布を導出し、$\mu^{(m)}(x)$ および $\sigma^2(x)$ を用いた漸近的バイアスおよび分散の明示的表現を得る。
• 複素解析およびフーリエ変換技術を用いて、罰則作用素のグリーン関数を含む積分を評価し、特に単位根および留数積分を用いる。
• 境界付近でのスムーザの挙動を、Bスプラインの対称性および台の性質を用いて、エッジ効果を含む核近似に拡張して分析する。
• 不規則に配置された設計点を扱うために、それらを正則なグリッドにマッピングするバインニング手順を導入し、ややきめの設計条件のもとで漸近的性質を保持することを示す。

実験結果

リサーチクエスチョン

• RQ1 Pスプライン推定量の漸近的バイアスおよび分散は、スプライン次数 $p$ および罰則順序 $m$ にどのように依存するか？
• RQ2 Knot の数 $K$ の増加レートがPスプライン推定量の漸近的分布に影響を与えるか？
• RQ3 Pスプラインの同等核を明示的に特徴づけられるか？また、それは $p$ に依存するのか、それとも $m$ のみに依存するのか？
• RQ4 漸近的正規性および最小平均二乗誤差を達成するための平滑化パラメータ $\lambda$ の最適収束レートは何か？
• RQ5 境界点におけるPスプラインの漸近的性質は、内部点とどのように異なるか？

主な発見

• Pスプラインの同等核は、スプライン次数 $p$ に依存せず、罰則順序 $m$ のみに依存する。これは、$p$ が局所的な滑らかさに影響するが、漸近的核の形状に影響しないことを確認する。
• 漸近的バイアスは $\lambda h_n^{2m}$ に比例し、漸近的分散は $\sigma^2(x)/(n h_n)$ に比例する。それぞれ $\mu^{(m)}(x)$ および $\sigma^2(x)$ に明示的に依存する。
• 平滑化パラメータの最適収束レートは $\lambda \sim n^{-2m/(2m+1)}$ であり、正規性条件のもとで平均二乗誤差を最小化する。
• Knot の数 $K$ の発散レートが $K \to \infty$ かつ $K h_n \to \infty$ を満たす限り、漸近的分布に影響しない。
• 境界付近では、有効バンド幅が小さくなるため、漸近的分散は内部点よりも大きくなるが、バイアスのオーダーは同じである。
• 不規則に配置されたデータに対しても、バインニング手順により漸近的分布が保持され、ややきめの設計条件のもとで非等間隔設計への理論の拡張が可能になる。