In a similar manner the rejection sampler can be setup and run using
options.mcmc.adaptive_rejection=1; % automatically adjust the normalizing likelihood options.mcmc.nite=100000; options=sippi_rejection(data,prior,forward,options);
This site is hosted by sourceforge.net |