I’ve started playing around with brms and Likert scale data and I’ve stumbled upon something that I need someone with a bigger brain to sanity check for me. I’ve read @paul.buerkner’s tutorial paper (which I can highly recommend) but I couldn’t find a solution to my problem.

If I understand things correctly I should model Likert scale (e.g., 1…5) data as:

`Response | resp_cat(5) ~ ...`

(see what I did with the … @paul.buerkner ;)

But what happens if, among all subjects, no one answered, e.g., 3, and all answers are either 1-2 or 4-5? Or for that matter (a bit more common) no one has provided the answer 1 or 5? We’re not talking about missing data here, i.e., it’s missing but it should be missing…

`resp_cat()`

requires a number, i.e., the number of categories, but does it figure out by itself that, e.g., 3 is missing, or does that not matter? Intuitively, I feel that it’s a big difference if nobody has answered 1 or 3, and that it affects(?) the cumulative prob.

In `formula-helpers.R`

the code for `resp_cat()`

only checks things like, e.g., `is.numeric(x)`

, while I would expect that I should be forced to list each category explicitly, i.e., `c(seq(1,5))`

.

Any insight would be nice! Thanks :)