Call:
glm(formula = response ~ diet + room, family = binomial(link = "logit"), 
    data = separation_data)

Coefficients:
                 Estimate Std. Error z value Pr(>|z|)
(Intercept)     2.536e+01  1.193e+05       0        1
dietHigh fibre -5.072e+01  1.377e+05       0        1
roomRoom 2     -5.072e+01  1.822e+05       0        1
roomRoom 3      5.072e+01  1.822e+05       0        1
roomRoom 4      1.186e-11  1.377e+05       0        1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 2.2181e+01  on 15  degrees of freedom
Residual deviance: 2.3267e-10  on 11  degrees of freedom
AIC: 10

Number of Fisher Scoring iterations: 25

fit$converged: TRUE
Fitted range: 2.220446e-16 to 1.000000e+00