### Multiple categorical regression

A multi-variable Bayesian categorical regression model using the iris data.

```{r multiple_linear_multilogit_data}
data(iris)
X <- as_data(cbind(1, iris[, 1:4]))
y <- model.matrix(~ Species - 1, iris)
P <- ncol(X)
K <- ncol(y)
```

```{r multiple_linear_multilogit_greta}
beta <- normal(0, 5, dim = c(P, K - 1))
eta <- X %*% beta
prob <- imultilogit(eta)
distribution(y) <- categorical(prob)
```
