julia> using GLM

julia> data = (a = [1, 4, 9], b = [2, 5, 7], c = [3, 6, 11], y = [1, 1, 0])
(a = [1, 4, 9], b = [2, 5, 7], c = [3, 6, 11], y = [1, 1, 0])

julia> X = [data.a data.b data.c]
3×3 Matrix{Int64}:
 1  2   3
 4  5   6
 9  7  11

julia> glm(X, data.y, Bernoulli(), LogitLink());

julia> form = @formula(y ~ a + b + c);

julia> glm(form, data, Bernoulli(), LogitLink());
ERROR: PosDefException: matrix is not positive definite; Cholesky factorization failed.
[...]