Could you say more about the efficient sampling through convolutions? Bad interactive latency was a major reason I didn’t spend any time exploring Bayesian approaches for my quick-and-dirty resorter script.
I take it you are not using a conjugate approach for fast exact estimates, nor a Laplacian approximation, nor full slow MCMC, but I’m not sure what sort of speed you actually get from your approach.
The bottleneck for sampling from the posterior is sampling v from the truncated normal distribution, which I’m doing using this implementation of minimax tilting rejection sampling. It’s not exactly instant, but on my laptopt it’s a few CPU-ms/sample up to ~100 comparisons and 50 CPU-ms up to ~200 comparisons. This probably prevents it from being useful for lists with many hundreds of items, but I’ve found it to be fine for interactive use since selecting the next comparison doesn’t require any sampling. The only thing which is prohibitively slow (for me at least) is computing the entropy of the full posterior since it involves so many evaluations of Φm.
Could you say more about the efficient sampling through convolutions? Bad interactive latency was a major reason I didn’t spend any time exploring Bayesian approaches for my quick-and-dirty
resorter
script.I take it you are not using a conjugate approach for fast exact estimates, nor a Laplacian approximation, nor full slow MCMC, but I’m not sure what sort of speed you actually get from your approach.
The bottleneck for sampling from the posterior is sampling v from the truncated normal distribution, which I’m doing using this implementation of minimax tilting rejection sampling. It’s not exactly instant, but on my laptopt it’s a few CPU-ms/sample up to ~100 comparisons and 50 CPU-ms up to ~200 comparisons. This probably prevents it from being useful for lists with many hundreds of items, but I’ve found it to be fine for interactive use since selecting the next comparison doesn’t require any sampling. The only thing which is prohibitively slow (for me at least) is computing the entropy of the full posterior since it involves so many evaluations of Φm.