The standard way (I think) to extract specific parameters from brsm (and other) fits is to use extractor functions like coef, fixef and ranef. This extracts the population and group level parameters. Is there a standard way to extract the family level parameters?
A simple example, where the family level parameter is sigma:
my_model <- brm(hp ~ mpg + (1 | cyl), mtcars)
fixef(my_model)
ranef(my_model)
coef(my_model) # just combines fixef and ranef
sigma(my_model) #This isn't defined for brms
We can use the draws directly to get the answer, e.g.:
mean(as.data.frame(my_model)$sigma)
- Operating System: Linux mint 19.1
- brms Version: 2.7.0