Modelling second price auction results

You can try to use a Simulator with sample_smc, it’s the dirtiest way to do inference without a closed form likelihood, when you have a generating process.

Might not scale very well though: pymc.Simulator — PyMC dev documentation