Fitting joint distributions using custom distribution as part of mixture

Since in this toy example the two dimension of y are independent, I suggest you to separate the 2 and write 2 mixtures, each with observed y[0, :] and y[1, :], and then try combining them.