Am I right that the models specified in the linked blog post could be specified and fit using regular old pymc3 (no Theano scan), but with two caveats:
-
You can only generate the smoothed posterior distributions, not the filtered, and
-
The sampling might not be as efficient as when using FFBS