Hi,
I am currently a user of MCSim (
http://www.gnu.org/software/mcsim/),
but would love to have a Python-based tool to carry out MCMC and Monte
Carlo simulations.
We are conducting MCMC simulations of pharmacokinetics (similar to
those described in
http://www.ehponline.org/members/2000/suppl-2/275-282bois/bois-full.html).
To replace MCSim, I am envisioning the use of something like SciPy to
solve the differential equation model coupled with PyMC to carry out
the MCMC analysis and fit the mathematical model to data sets and
compute posterior parameter distributions.
Does this sound appropriate for PyMC?
In general, do you know how PyMC compares to MCSim?
Does PyMC have hierarchichal model capabilities?
Thank you for your help and insights.
Cheers,
brad