インターセプトなしでプロビットモデルを作成する必要があります。それでは、Rのプロビットモデルから切片を削除するにはどうすればよいですか?
howとは言わないが、プロビットモデルに適合させるつもりであるが、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"))
_
それを行う必要があります(再度、データ名と変数名を適切に置き換えます。)
また、既存の数式オブジェクトfoo
がある場合、次のようにupdate
を使用してインターセプトを削除できます。
foo <- y ~ x1 + x2
bar <- update(foo, ~ . -1)
# bar == y ~ x1 + x2 - 1