предсказания GLM значения

Question

My glm формулу:

fit.full <- glm(y ~ a + b + c + d + e + f, data=df, family=binomial(logit)) 

glm выход:

Estimate Std. Error z value Pr(>|z|) 
a -0.0443 0.0182 -2.43 0.01530 * 
b1 -0.3247 0.0898 -3.62 0.00030 *** 
c2 0.0948 0.0322 2.94 0.00326 ** 

Затем я решил применить glm к уменьшенной модели для сравнения:

fit.reduced <- glm(y ~ a + b1 + c2, data=df, family=binomial(logit)) 

но, прежде код произвел ошибку:

Error in eval(expr, envir, enclos) : object 'b1' not found

Я должен был написать код ниже:

fit.reduced <- glm(y ~ a + b + c, data=df, family=binomial(logit)) 

Есть ли решение ввести только: а, b1 и с2? Кстати, переменные b и c имеют 3 и 5 значений соответственно. Благодарю.

Answer 1

Вывод из вашей первой модели включает в себя как «b1», так и «c2» в качестве меток линий для коэффициентов, но это не означает, что есть столбцы с этими именами. Это означает, что эти переменные являются факторами с уровнями, которые имеют эти имена. Настройки «из коробки» glm производят «контрасты обработки», поэтому отображаются коэффициенты уровней базового фактора с именами переменных, за которыми следует их коэффициент. модель с этими уровнями, так как не было бы ничего, чтобы сравнить их два, если бы вы удалили все случаи с базовыми значениями уровня.

Если вы хотите запустить модель с одинаковыми именами столбцов, но удаление определенных уровней из рассмотрения, затем используйте параметр subset, который поддерживает glm.