Я пытаюсь создать полиномиальную модель логистической регрессии, используя nnet::multinom(). У меня есть 2 независимые переменные (числовые от 0 до 10) и зависимая переменная (фактор с 4 уровнями 1,2,3,4).
Проблема в том, что при проверке мультиколлинеарности функцией vif() из car я получаю следующее предупреждение:
Warning message: No intercept: vifs may not be sensible.
Вот пример кода:
set.seed(150)
example <- data.frame(replicate(2,sample(0:10,1000,rep=TRUE)))
example$dependent <- factor(sample(1:4, 1000, rep=TRUE))
example$dependent <- levels(relevel(example$dependent, ref = "1"))
mymod <- multinom(dependent~., data=example)
vifs <- vif(mymod)
vifs
выход:
X1 X2
5.297042 5.257941
Это дает разумные значения, однако с моими фактическими данными они варьируются от 2,289297e ^ 1 до 5,609178e ^ 15.
Любой совет? Я не понимаю сообщения о том, что перехвата нет, потому что я думал, что перехват выполняется автоматически в модели, если вы его не удалите.
