StatsModels R-like formula for a general quadratic model

Specific Domains Statistics
6

Based on your solution, the problem can be solved by

a = sum(InteractionTerm(term.((only(x),only(x)))) for x ∈ combinations([:a, :b, :c],1))
b = sum((&)(term.(x)...) for x ∈ combinations([:a, :b, :c],2))
c = sum(term(only(x)) for x ∈ combinations([:a, :b, :c],1))
formula = Term(:y) ~ a + b +c + ConstantTerm(1)

Although an even shorter form would be nice indeed.

1 Like