Multivariate GLMs about turingglm.jl HOT 2 OPEN

Does StatsModels.jl and MixedModels.jl formulas support multivariate stuff in the right hand side of ~?

from turingglm.jl.

I think that this would be a great idea, but there may be some practical issues. I don't believe that MixedModels.jl supports multivariate (not multiple) linear models. Even MultivariateStats.jl doesn't seem to support it.

If I enter a formula with multiple responses variables, it appears to work, at least minimally:

julia> @formula(y + z ~ 1 + x)
FormulaTerm
Response:
  y(unknown)
  z(unknown)
Predictors:
  1
  x(unknown)

But, it looks like there have been some questions about easily implementing multivariate models in Turing more generally; e.g., https://discourse.julialang.org/t/singular-exception-with-lkjcholesky/85020 (there are other discussion too, some of which are linked therein). People seem to be having issues with using a prior commonly used for the covariance matrix in multivariate regression models.

from turingglm.jl.