stan VS PyMC

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/

I believe the pymc versions were resolved into developing version 4 of pymc. Development at https://github.com/pymc-devs/pymc

It still depends on theano now evolved and renamed

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