rstan VS brms

I have a Colab notebook with an R kernel that I'm using to share with students for remote lessons in statistics. This notebook relies on "rstanarm", which is pretty massive with the number of dependencies - it takes ~50minutes to install into a fresh Colab session with install.packages(). It seems the issue is that many of the dependencies of this package need to be compiled from source, which takes a long time on Linux distributions like Colab.

You need to install the preview version of rstan: https://github.com/stan-dev/rstan/wiki/Configuring-C---Toolchain-for-Windows

If you're on Windows, there are extra hoop to jump through, I'm afraid https://github.com/stan-dev/rstan/wiki/

And here are the instructions on how to build RStan from source: https://github.com/stan-dev/rstan/wiki

I ran the R code from here

In R, the rstanarm package should do you well. You'll need to install rstan and make sure you have a C++ complier set up as well (instructions here: https://github.com/stan-dev/rstan/wiki/RStan-Getting-Started)

[2] https://paul-buerkner.github.io/brms/

Maybe something like the Logistic Normal ? (e.g. see this issue from brms). If that fits what you are looking for, you can use brms to generate the Stan code for you (brms::make_stan_code()) and work from that.

Okay so first off, I recommend that you read [this](https://link.springer.com/article/10.3758/s13423-016-1221-4) article about "The Bayesian New Statistics", which highlights estimation rather than hypothesis testing from a Bayesian perspective (see Fig. 1, second row, second column). Instead of a t-test, then, we can *estimate the difference* between two groups/variables. If you want to go deeper than JASP etc, I recommend that you use [brms](https://paul-buerkner.github.io/brms/), or, if you want to go even deeper, [Stan](https://mc-stan.org/) (brms is a front-end to Stan).

I think the brms package (https://github.com/paul-buerkner/brms) or the blavaan package (http://ecmerkle.github.io/blavaan/) have support for SEM. I've never done it myself, so I unfortunately can't give you any direction for that in particular. However, I have used stan in multi-level meta-analysis regression (combining multiple CRISPRa experiments to find determinants of CRISPRa activity, see https://github.com/timydaley/CRISPRa-sgRNA-determinants/blob/master/metaAnalysis/NeuronAndSelfRenewalMetaMixtureRegression.Rmd) and had some success.

brms baby

Could also be done using brms and the gr term. See this for the motivation behind this syntax.

Anyway, I found I can include weights into the brm function by using gr(RE, by = var) to deal with the heterogeneous variance and it should automatically assume that each observation within a group is correlated according to the brms reference manual.

This is what actually should work- I must be declaring my variables incorrectly. The issue I'm having is that what you refer to as lin , I tried calling a few things, from b to LinPred (which worked in the link here: brms issue 47). When I've tried doing this, I receive errors that say "The following variables are missing from the dataset....[insert variable used to symbolize linear part of the model)". But I believe you're code is on the right path for what needs to be done- I'll try altering my syntax to be sure it resembles yours let you know if it works.