如何删除R中的拦截

Question

我需要在没有拦截的情况下创建一个probit模型.那么,如何从R中的probit模型中删除截距？

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