stan VS probability

Alternatives: stan and edward

remove.packages(c("StanHeaders", "rstan")) install.packages("rstan", repos = c("https://mc-stan.org/r-packages/", getOption("repos")))

My approach to problems like this is to write down the proposed model mathematically first, in extreme detail. I find hierarchical form to be the easiest way to break it down piece by piece. Once I have the maths then I turn it into a Stan model. Last step is to use the Stan output to answer the research questions.

Here is a good breakdown of recommendations from Andrew Gelman.

For instance my first choice in these cases is always a Bayesian inference tool like Stan. In my experience as someone who’s more of a programmer than mathematician/statistician, Bayesian tools like this make it much easier to not accidentally fool yourself with assumptions, and they can be pretty good at catching statistical mistakes.

I tend to be most impressed by tools and libraries. The stuff that has most impressed me in my time in ML is stuff like pytorch and Stan, tools that allow expression of a wide variety of statistical (and ML, DL models, if you believe there's a distinction) models and inference from those models. These are the things that have had the largest effect in my own work, not in the sense of just using these tools, but learning from their design and emulating what makes them successful.

oh its certainly used in practice. you should look into frameworks like Stan[1] and pyro[2]. i think bayesian models are seen as more explainable so they will be used in industries that value that sort of thing

[1] https://mc-stan.org/

TensorFlow-Probability

If you are interested in probabilistic programming, causal modeling and bayesian graphical modeling, I recommend checking out Tensorflow Probability (https://www.tensorflow.org/probability).

As others have said, R for academia, Python for industry. However, i'd also throw Stan into the mix, along with other PPL frameworks like Tensorflow Probability and Pyro. The latter two will require you to learn Python first, though.

This tutorial explains what is probabilistic programming & provides a review of 5 frameworks (PPLs) using an example taken from Chapter 4 of Statistical Rethinking by Dr. Richard McElreath. Frameworks (PPLs) reviewed are - Stan (https://mc-stan.org/) PyMC3 (https://docs.pymc.io/) Tensorflow Probability (https://www.tensorflow.org/probability) Pyro/NumPyro (https://pyro.ai/) Turing.jl (https://turing.ml/stable/) I also provide the basic review of a great library called arviz (https://arviz-devs.github.io/arviz/), which can be used for all the above-mentioned PPLs to do Exploratory Data Analysis of Bayesian Models. Here is the link to the notebook in which I have implemented the example model using the above Frameworks/PPLs https://colab.research.google.com/drive/1zgR2b0j2waGi1ppnIe1rw7emkbBXtMqF?usp=sharing

You implemented much of your model yourself. I have come across this GItHub issue recently where people reported a large slowdown when they were running their self-built models eagerly. It seems to be somewhat similar to yours. One answer (this one: https://github.com/tensorflow/probability/issues/629#issuecomment-551875968) seems to allow people to achieve speedups.