I was searching for this answer and I'm really suprised that haven't found it. I just want to peform three level logistic regression in R.

Let's define some artificial data:

```
set.seed(42)
y <- sample(0:2, 100, replace = T)
x <- rnorm(100)

```

My variable `y` is containing three numbers - 0, 1 and 2. So I thought that the simplest way would be just to use:

```
glm(y ~ x, family = binomial("logit"))

```

However I got information that `y` should be in interval [0,1]. Do you know how I can perform this regression ?

Please notice - I know that it's not so straightforward to perform multilevel logistic regression, there are several techniques how to do so e.g. One vs all. But as I was seeking for it, I wasn't able to find any.