Multivariate Polynomial Regression

The tutorial is great, but here's where I'm at a loss: if I have multiple independent variables that I want to consider in the regression, how can I do that in a way that's similar to sklearns PolynomialFeatures: https://scikit-learn.org/stable/modules/generated/sklearn.preprocessing.PolynomialFeatures.html#sklearn.preprocessing.PolynomialFeatures
I've been using MLJLinearModels.jl where applicable, but would it be better to just sub-in the SKLearnMLJ interface?
If there's a way to do it that easily builds off of the tutorial and I've just missed that, pointing me in the right direction is also appreciated.

Inizia leggendo il tutorial di regressione esistente e confronta gli approcci MLJLinearModels e SKLearnMLJ citati nell’issue. Determina se il workflow richiesto appartiene al tutorial e quale percorso Julia supportato dovrebbe essere documentato. Il lavoro è completato quando è disponibile un esempio riproducibile di regressione polinomiale multivariata approvato da un maintainer oppure una spiegazione chiara dell’approccio supportato.