I will give PyStan a try over the holidays. Also, many thanks for all of the great comments here (statistical modeling is not in my toolbox yet, and the discussion here is grounding).