glmmTMB_ordbeta_reprex.R

library(dplyr)
library(glmmTMB)

nobs <- 100
ngroup <- 5
dat <- tibble(
  x = sample(1:3, size = nobs, replace = TRUE, prob = c(0.15, 0.5, 0.35)),
  y = runif(nobs),
  z = case_when(
    x == 1 ~ 0,
    x == 2 ~ y,
    x == 3 ~ 1
  ),
  g = rep(seq_len(ngroup), length.out = nobs)
)
glmmTMB(z ~ (1|g), data = dat, family = ordbeta)