Using multiple index variables in Turing model

Specific Domains Statistics
1

Hello everyone,

I would like to use two index variables in a model, but I don’t know how to do it… I would like to do something like this:

@model function m2(A, G, D)
    α ~ MvNormal([[0 0]
                   [0 0]], 1)
    p = α[G, D]
    A ~ arraydist(LazyArray(@~ BernoulliLogit.(p)))
end

But this obviously gives an error. Is there a way to do this?