Hi all,
I upgraded this week to R 4.0.1 after my previous machine quit on me; I had been using 3.6.2. I’m working on a paper using sequential models with continuous parameterization, and when I switched to 4.0.1, something changed in the call between brms and Stan: Stan began to automatically attempt to estimate parameter values for disc. Before the switch, disc was only estimated if I explicitly included it in the model formula. This does not appear to harm anything: Stan attempts to estimate disc and passes all the expected warnings (“chains have not mixed,” “ESS is too low”). However, as my goal is to write accessible code that won’t confuse someone attempting to replicate my results, I have resorted to dropping warnings into my code telling users not to worry about the warnings, the lack of mixing for the mysterious disc on traceplots, etc. I wonder if I’m missing a simple fix, though I haven’t found any discussion of one in the Stan forums or in Buerkner & Vuorre 2018. Is there something I’m missing?
Toy example below. Note that this change with R 4.0 appears to affect both sratio and cratio models.
Thanks in advance!
dat <- data.frame(x = rnorm(500, 0, 1), y = ordered(sample(1:3, 500, replace = T)))
this <- brm(y ~ x, family = cratio, data = dat)
- Operating System: Windows 10
- brms Version: 2.13.0 (R 4.0.1, run with RStudio)