That is not exactly what I want to model.
Anyway, if I use Bambi, the model formula is quite cumbersome, as I need to write the model as a string, without using matrix algebra, somehow like this:
model = bmb.Model("y ~ var_1*covar_1 + var_2*cov_2 + ... + var_n*covar_k", data, family="gaussian")
Or is there also a possibility to use dot products, like
model = bmb.Model(y = np.dot(X, Z), data, family=“gaussian”)
?