Ответ 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"))
должен сделать это (опять же подставляя в свои данные и имена переменных соответственно).