度数1のスプラインモデルを異なる結び目でコーディングするにはどうすればよいですか？

Question

私は、複数のノットを持ついくつかの独立変数を持つ多変数スプラインモデルをコーディングしようとしています。スプラインを持つ変数は常に1の次数を持ちます。私はいくつかのコードを持っていますが、私がrでスプライン回帰をしていないので信頼しているかどうかわかりません（独自の「ブラックボックス」ソフトウェアでのみ）。以下はコードです。

私はスプラインで6,000の投稿の多くをチェックしました。私は混乱しているので、非常に多くの異なるコードを参照してください。

誰でもこのコードが私がしたいことをしているか（度= 1 /別の結び目） a）私に教えてください。 b）これを行うには良い方法がありますか？

fit1 <- glm(freq ~ channel + term2 + pay_plan_bucket_2 + 
      state + eff_year + marital_status + 
      vehicle_type + insured_age_bucket + I(pmax(0, insured_age_bucket- 
26)) + I(pmax(0, insured_age_bucket - 70)) + 
      vehicle_length_bucket + I(pmax(0, vehicle_length_bucket - 45)) + 
      veh_age + I(pmax(veh_age -7)) + I(pmax(veh_age - 18)) + 
      rba_bucket + I(pmax(0, rba_bucket - 3500)) + I(pmax(0, rba_bucket - 27000)) + 
      credit_tier_bucket + I(pmax(0, credit_tier_bucket - 3)), 
     family=quasipoisson(link="log"), 
     data=comp_training_set_newpayplan) 

ありがとうございます。

Answer 1

スプラインに対するあなたのブルートフォースのアプローチはおそらく正しいでしょう。スプラインパッケージのbsを使用して出力を確認します。例：bs(credit_tier_bucket, knots=3, degree=1)を式の単一項として使用します。基底がここで行ったのとは異なる形で形成されるので、両方のモデルから予測値を集め、等価であることを確認してスプラインをコード化する2つのアプローチが同等の推定と推論を提供するようにします。