-
Notifications
You must be signed in to change notification settings - Fork 4
Open
Description
Would it be possible to add support for ordinal models fit with the bayesian BRMS package:
# List required packages
Pkgs <- c("tidyverse", "ordinal", "brms", "rstan")
# Load packages ====
lapply(Pkgs, require, c = T)
# Load data
data <- ordinal::wine
# Model formula
fmla <- bf(Rating ~ temp + contact +temp:contact + 1|judge) )
# Compile and fit model
mod <- brm(fmla,
data = data,
family = cumulative("probit"),
inits = 0,
iter = 2000,
warmup = 1000,
chains = 2,
cores = 2,
sample_prior = TRUE,
save_all_pars = TRUE)
)
Metadata
Metadata
Assignees
Labels
No labels