tjhladish / AbcSmc

Sequential Monte Carlo Approximate Bayesian Computation with Partial Least Squares parameter estimator
GNU General Public License v3.0
15 stars 7 forks source link

read-last-set-only filtering mode #50

Open pearsonca opened 1 year ago

pearsonca commented 1 year ago

Currently, every pass with AbcSmc must ingest the entire history of complete jobs to process the latest set => propose next set. This is because the weights are computed iteratively: the latest set weights actually depend on the entire history of weights. We don't store those weights, nor the parameters for noising draws.

We should continue to store the history, but if we also start storing weights and doubled variances, then we can simplify the code that deals with filtering to just take the last complete step. This could allow us to substantially trim down the AbcSmc object itself, getting rid of all the vectors-of-stuff and the assorted look-up-by-set accessors.

At the moment, we have no fitting-algorithmic reason to use all the data, though we do look at it for diagnostic reasons. We don't generally need an AbcSmc object, however, to do the diagnostic steps.