web-dev-qa-db-ja.com

Rでインターセプトを削除する方法

インターセプトなしでプロビットモデルを作成する必要があります。それでは、Rのプロビットモデルから切片を削除するにはどうすればよいですか?

16
Kazo

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"))
_

それを行う必要があります(再度、データ名と変数名を適切に置き換えます。)

34
Gavin Simpson

また、既存の数式オブジェクトfooがある場合、次のようにupdateを使用してインターセプトを削除できます。

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