Brms model crashes when using specific variable

#1
  • Operating System: Windows 10
  • brms Version: 2.7.0

Here’s a sample of the dataset I’m working with and the offending model:

structure(list(imm = c(2, 0, 0, 6, 3, 0, 42, 0, 0, 0), week = c(4, 
5, 2, 11, 11, 3, 7, 7, 10, 4), rep = structure(c(2L, 3L, 2L, 
2L, 3L, 2L, 2L, 1L, 1L, 2L), .Label = c("1", "2", "3"), class = "factor"), 
    trt = structure(c(2L, 1L, 4L, 5L, 4L, 5L, 5L, 1L, 1L, 3L), .Label = c("1", 
    "2", "3", "4", "5"), class = "factor"), bindenom = c(10, 
    8, 11, 8, 9, 9, 11, 11, 6, 10), nlive = c(6, 0, 3, 0, 0, 
    7, 6, 0, 0, 6), plantid = c(15, 21, 32, 49, 11, 49, 51, 36, 
    34, 4)), row.names = c(NA, -10L), class = c("tbl_df", "tbl", 
"data.frame")) -> hwf

priors <- set_prior("normal(0,5)", class = "b")

brm(bf(nlive|trials(bindenom) ~ trt + week + rep,
       zi ~ trt), data = hwf, family = zero_inflated_binomial(), chains = 1)

The model will run fine if I either omit the zi term OR if I add more variables to the predictor. Using trt alone, however, crashes the RStudio session.

#2

On my windows 10 laptop with brms 2.8.0+, R 3.6, and the latest stable Rtools (3.5) it compiles and fits just fine. I don’t know what causes this problem on your machine but you may want to update brms, R, and/or Rtools.

1 Like
#3

Updated R to 3.6 and brms to 2.8.0, no issues now. Thanks!