- Operating System: Microsoft Windows [Version 10.0.17134.345]
- brms Version: 2.5.0
I want to fit a binomial brms model containing the weight option. I have followed the instructions discussed here but when I try using the code and data below I get the error:
Error: $ operator is invalid for atomic vectors
dt = read.table(header = TRUE, text = "
n r r/n group treat c2 c1 w
62 3 0.048387097 1 0 0.1438 1.941115288 1.941115288
96 1 0.010416667 1 0 0.237 1.186583128 1.186583128
17 0 0 0 0 0.2774 1.159882668 3.159882668
41 2 0.048780488 1 0 0.2774 1.159882668 3.159882668
212 170 0.801886792 0 0 0.2093 1.133397521 1.133397521
143 21 0.146853147 1 1 0.1206 1.128993008 1.128993008
143 0 0 1 1 0.1707 1.128993008 2.128993008
143 33 0.230769231 0 1 0.0699 1.128993008 1.128993008
73 62 1.260273973 0 1 0.1351 1.121927228 1.121927228
73 17 0.232876712 0 1 0.1206 1.121927228 1.121927228")
mod ← brm(r|trials(n)|weights(w) ~ treat*c2+(1|group), data=dt, family=binomial(link=logit))
I would appreciate any help to solve this error.