Are you aware of a function which can be used to compute contrasts based on brmsfit models
similar to lmerTest::contest function e.g.
I assume I could compute a contrasts just by:
postsample <- posterior_samples(fit1,"^b") as.matrix(postsample) %*% as.matrix(c(1,1,0,-1))
but than I also will need to write all the functions to create visualizations and summaries.
Could you point me in the right direction?
Please also provide the following information in addition to your question:
- Operating System: Windows Linux