r - Rで切片を削除する方法

Question

切片なしでプロビットモデルを作成する必要があります。では、どうすればRのプロビットモデルから切片を削除できますか？

Answer 1

プロビットモデルをどのように適合させるつもりなのかはわかりませんが、Rの式表記を使用してモデルを記述する場合は、+ 0または- 1式の一部として切片を抑制することができます。

mod <- foo(y ~ 0 + x1 + x2, data = bar)

また

mod <- foo(y ~ x1 + x2 - 1, data = bar)

(もちろん両方とも疑似 R コードを使用 - モデリング関数とデータ/変数を置き換えます。)

これがモデルに適合している場合は、次のglm()ようになります。

mod <- glm(y ~ x1 + x2 - 1, data = bar, family = binomial(link = "probit"))

それを行う必要があります（必要に応じて、データと変数名を再度置き換えてください。）

Answer 2

また、既存の数式オブジェクトがある場合は、次のようfooに切片を削除できます。update

foo <- y ~ x1 + x2
bar <- update(foo, ~ . -1)
# bar == y ~ x1 + x2 - 1