Does PyMC work for ABC SMC?

Questions
,
4

Thank you @junpenglao. I saw this, originally could not get it to work, and just updated to v3.11.4 and it now works. I’m able to get my model running in the ABC SMC framework as described in that article. Unfortunately, it seems to hang after the first stage:

Initializing SMC sampler...
Sampling 2 chains in 1 job
/Users/lollipop/miniconda3/envs/pymc311/lib/python3.9/site-packages/theano/tensor/elemwise.py:826: RuntimeWarning: divide by zero encountered in log1p
  variables = ufunc(*ufunc_args, **ufunc_kwargs)
Stage:   0 Beta: 0.006
INFO:pymc3:Stage:   0 Beta: 0.006

After this first stage, which takes around 2-5 minutes on my MacBook Pro M1 Max, it hasn’t nudged after >20 minutes. I’m new to PyMC and there is very little documentation on ABC with PyMC (as @cluhmann alluded to with the pm.Simulator(), so excuse my initial lack of detail.

Is there a way to get a verbose output from either pm.Simulator or pm.sample_smc? Or, is there a check or a routine I can perform to see what’s going on?

Thanks again!