While I’m on a microblogging spree, I’ve been meaning to link to this informative comparison of pymc, emcee, and pystan: http://jakevdp.github.io/blog/2014/06/14/frequentism-and-bayesianism-4-bayesian-in-python/
Tag Archives: pymc
MCMC in Python: Another thing that turns out to be hard
Here is an interesting StackOverflow question, about a model for data distributed as the sum of two uniforms with unknown support. I was surprised how hard it was for me.
I think the future of probabilistic programming should be to make a model for this easy to code.
Comments Off on MCMC in Python: Another thing that turns out to be hard
Filed under statistics
MCMC in Python: Never… no… always check for convergence
I’ve had no teaching responsibilities over the last quarter, and I must miss it. I’ve found myself responding to PyMC questions on StackOverflow more than ever before. It is an interesting window into what is hard in Bayesian computation. Checking (and achieving) MCMC convergence is one thing that is hard. Here are some questions and answers that include it:
http://stackoverflow.com/questions/24294203/difference-between-bugs-model-and-pymc/24347102#24347102
http://stackoverflow.com/questions/24242660/pymc3-multiple-observed-values/24271760#24271760
http://stackoverflow.com/questions/24402834/fitting-power-law-function-with-pymc/24413323#24413323
Comments Off on MCMC in Python: Never… no… always check for convergence
Filed under statistics
MCMC in Python: sampling in parallel with PyMC
Comments Off on MCMC in Python: sampling in parallel with PyMC
Filed under software engineering
MCMC in Python: Estimating failure rates from observed data
A question and answer on CrossValidated, which make me reflect on the danger of knowing enough statistics to be dangerous.
Comments Off on MCMC in Python: Estimating failure rates from observed data
Filed under statistics
MCMC in Python: How to make a custom sampler in PyMC
The PyMC documentation is a little slim on the topic of defining a custom sampler, and I had to figure it out for some DisMod work over the years. Here is a minimal example of how I did it, in answer to a CrossValidated question.
Comments Off on MCMC in Python: How to make a custom sampler in PyMC
Filed under MCMC
MCMC in Python: How to set a custom prior with joint distribution on two parameters in PyMC
Question and answer on Stackoverflow. Motivated by question and answer on CrossValidated about modeling incidence rates.
Comments Off on MCMC in Python: How to set a custom prior with joint distribution on two parameters in PyMC
Filed under Uncategorized
