jakeret / abcpmc

Approximate Bayesian Computation Population Monte Carlo
GNU General Public License v3.0
40 stars 19 forks source link

Update distance function #14

Open NiallMac opened 7 years ago

NiallMac commented 7 years ago

Is it possible to update the distance metric between iterations? My motivation for this is the case where one has multiple summary statistics entering the distance metric, but it is not obvious how to weight them. A reasonable choice would involve the covariance matrix of the summary statistics, but this is not known in advance. I'd like to use e.g. the sample covariance of the summary statistics in the previous iteration to weight them in the distance metric.

jakeret commented 7 years ago

I never thought about this but in principle I think this should be possible. Once an interation is computed:

for pool in abcpmc_sampler.sample(prior, eps):

you could set you new distance function:

sampler.dist = new_dist_with_covariance_matrix