# StatsModels R-like formula for a general quadratic model

**URL:** <https://discourse.julialang.org/t/statsmodels-r-like-formula-for-a-general-quadratic-model/64441>\
**Category:** Statistics\
**Created:** [July 10, 2021, 11:29pm UTC](https://discourse.julialang.org/t/statsmodels-r-like-formula-for-a-general-quadratic-model/64441 "2021-07-10T23:29:14Z")\
**Posts on this page:** 1\
**Showing post:** 2

<div class="post-metadata">

**Author:** ![xiaodai](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/xiaodai/32/15937_2.png) [@xiaodai](https://discourse.julialang.org/u/xiaodai)\
**Post date:** [July 11, 2021, 1:54am UTC](https://discourse.julialang.org/t/statsmodels-r-like-formula-for-a-general-quadratic-model/64441/2 "2021-07-11T01:54:11Z")

</div>

Doesn’t seem to say in the manual [Modeling tabular data · StatsModels.jl](https://juliastats.org/StatsModels.jl/stable/formula/#Modeling-tabular-data)

But I was curious.

I think you can do it via constructing terms directly.

```julia
using Base.Iterators
function allinteractions(terms::Vector{Symbol})
    tmp = sum(Term.(terms))
    tmp + sum([t1&t2 for(t1, t2) in Iterators.product(tmp, tmp)])
end

formula = Term(:y)~allinteractions([:a, :b, :c])
using GLM
lm(formula, df)

```

---

_[View the full topic](https://discourse.julialang.org/t/statsmodels-r-like-formula-for-a-general-quadratic-model/64441)._