The number of terms should not change the type of coeffnames,
it should always be a tuple of tuples
1 RHS term
df = (y=1:10, a = 1:10, b=1:10);
f = @formula(y ~ a);
f = apply_schema(f, schema(f, df));
gives:
julia> coefnames(f)
("y", "a")
2 RHS terms
df = (y=1:10, a = 1:10, b=1:10);
f = @formula(y ~ a + b);
f = apply_schema(f, schema(f, df));
gives
julia> coefnames(f)
("y", ["a", "b"])