I found this very old demo (added as regressions-demo-older.py). Might be useful to incorporate into pystata-demo