Brad Reisfeld breisfeld

Export to Joplin

The main code was lightly adapted from here. Consult these instructions, but use the code from this project.

Some notes:

I hereby claim:

To claim this, I am signing this object:

	import pymc

	# load the model
	import pk_model as mod

	# fit the model with mcmc
	mc = pymc.MCMC(mod)
	mc.use_step_method(pymc.AdaptiveMetropolis, [mod.kcp, mod.kpc, mod.ke])
	mc.sample(iter=2000, burn=1000, thin=20, verbose=1)