回帰の定義自然3次スプライン


16

私はHastieらによる「統計学習データマイニング、推論、および予測の要素」という本からスプラインについて学んでいます。145ページで、Natural 3次スプラインが境界の結び目を越えて線形であることを発見しました。スプラインにはノット、があり、本のそのようなスプラインについては次のとおりです。Kξ1、ξ2、。。。ξKここに画像の説明を入力してください

質問1: 4自由度はどのように解放されますか?私はこの部分を取得しません。

質問2:の定義において場合次いで。この式で著者がやろうとしていることは何ですか?これは、スプラインが境界ノットを超えて線形であることを確認するのにどのように役立ちますか?dk(バツ)k=KdK(バツ)=00

回答:


17
  1. まず、通常の3次スプラインを考えてみましょう。ノットの各ペアの間は立方体で、境界ノットの外側は立方体です。最初のキュービック(最初の境界ノットの左側)の4dfから開始し、各ノットは1つの新しいパラメーターを追加します(キュービックスプラインと導関数の連続性と2次導関数は3つの制約を追加し、1つの自由パラメーターを残します)。のパラメータノット。KK+4K

    自然な3次スプラインは両端で線形です。これにより、3次および2次部分が0に制限され、それぞれdfが1ずつ減少します。つまり、曲線の両端で2 dfになり、が減少します。KK+4K

    ノンパラメトリック曲線の推定に自由度の合計数(など)を費やすことができると決めたと想像してください。自然なスプラインを課すことは、通常の3次スプラインよりも4少ない自由度(同じノット数)を使用するため、これらのパラメーターを使用すると、境界ノット間の曲線をモデル化するために、さらに4ノット(したがって4パラメーター)を持つことができます。Ppp

  2. の定義はあることに注意してください(すべてに基底関数があるため)。したがって、そのリストの最後の基底関数、です。したがって、定義に必要な最高のは、です。(つまり、使用しないので、が何をする可能性があるかをする必要はありません。)、K = 1 、2 、。。。、K − 2 K N K = d K − 2 − d K − 1 k d k k = K − 1 d KNk+2k=1,2,...,K−2KNK=dK−2−dK−1kdkk=K−1dK


4

ノット例では、「これにより4つの自由度(両方の境界領域にそれぞれ2つの制約)が解放されます」というアサーションの詳細を説明します。関連する間隔は、および(したがって、間隔とノットがあります)。ξ 1、ξ 2 ] - ∞ 、ξ 1 [ ] ξ 1、ξ 2 [ ] ξ 2、+ ∞ [ | 私| = 3 | 私| − 1 = 22ξ1、ξ2]−∞、ξ1[]ξ1、ξ2[]ξ2、+∞[|私|=3|私|−1=2

(一般的な)3次スプラインの場合

規則性の制約がない場合、方程式があります。4|私|=12

1(ξ 1 ≤ X < ξ 2)。1(ξ 1 ≤ X < ξ 2)X 。1

1(バツ<ξ1)  ;  1(バツ<ξ1)バツ  ;  1(バツ<ξ1)バツ2  ;  1(バツ<ξ1)バツ3  ;
1(ξ 2 ≤ X )。1(ξ 2 ≤ X )X 。1(ξ 2 ≤ X )X 2。1(ξ 2 ≤ X )X 3
1(ξ1≤バツ<ξ2)  ;  1(ξ1≤バツ<ξ2)バツ  ;  1(ξ1≤バツ<ξ2)バツ2  ;  1(ξ1≤バツ<ξ2)バツ3  ;
1(ξ2≤バツ)  ;  1(ξ2≤バツ)バツ  ;  1(ξ2≤バツ)バツ2  ;  1(ξ2≤バツ)バツ3。

制約を追加することにより(キュービックスプラインは規則性を仮定します)、を追加する必要があります線形係数の制約。Crr=2(r+1)×(|私|−1)=3×(|私|−1)=6

最終的にの自由度になります。12−6=6

自然な3次スプラインの場合

「自然な3次スプラインは、追加の制約を追加します。つまり、その関数は境界ノットを超えて線形になります。」

規則性の制約がない場合、方程式があります(つの多項式と3次多項式が関係するため、両方の境界領域で方程式を削除しました)。4|私|−4=12−442

1(バツ<ξ1)  ;  1(バツ<ξ1)バツ  ;  
1(ξ1≤バツ<ξ2)  ;  1(ξ1≤バツ<ξ2)バツ  ;  1(ξ1≤バツ<ξ2)バツ2  ;  1(ξ1≤バツ<ξ2)バツ3  ;
1(ξ2≤バツ)  ;  1(ξ2≤バツ)バツ。

制約は以前と同じなので、線形係数に制約を追加する必要があります。3×(|私|−1)=6

最終的にはの自由度になります。8−6=2

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.