Input:
x <- data.frame("Question" = c(0,0,0,0,1,1,1,1,0),
"Age" = 1:9, "Color" = c("R","B","G", 'G','B','R','B','B','G'))
model <- glm(Question ~., data = x, family = binomial('logit'))
summary(model)$coefficients
Output:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -55.323075 156049.02 -3.545237e-04 0.9997171
Age 15.774626 39420.88 4.001592e-04 0.9996807
ColorG -111.445840 263204.29 -4.234195e-04 0.9996622
ColorR 0.112095 47453169.66 2.362223e-09 1.0000000
How can I exclude ColorR instead of ColorB from the coefficients in the model?
