Back

Overinflation and overconcentration: why Cauchy perturbation kernels are the right choice for ABC-SMC

Sturrock, M.; Shahrezaei, V.

2026-07-09 systems biology
10.64898/2026.06.24.734205 bioRxiv
Show abstract

Approximate Bayesian computation sequential Monte Carlo (ABC-SMC) propagates its particles with a perturbation kernel, and with the standard Normal kernel it degrades sharply as the parameter dimension grows, a failure usually attributed to dimension itself. We show instead that it is governed by the quality of the summary statistics, with dimension entering only through a separate and milder mechanism, and that the two must act together for the Normal kernel to break. The first ingredient is covariance overinflation: the kernel covariance, estimated from the particle cloud, overshoots the true posterior covariance by a factor set by information loss in the summary statistics. We derive this overscaling factor in closed form for a Gaussian model with sufficient statistics and show that it stays modest at any dimension, shrinking toward its baseline value as the tolerance tightens; the extreme values seen in practice (of order 103) are a signature of insufficient summaries, not of dimension. The second ingredient is perturbation overconcentration: the normalised Normal step size concentrates around one as the dimension grows, so every proposal overshoots by the same factor. Either ingredient alone is harmless; only their combination breaks the Normal kernel. A Cauchy kernel (multivariate t with one degree of freedom) removes the concentration, keeping a positive acceptance rate under arbitrary overscaling at a bounded worst-case cost of 1.87x in expected squared jump distance. In a Metropolis-Hastings framework we derive closed-form acceptance rates for both kernels that illustrate the advantage of the Cauchy kernel in this limit. A series of full ABC-SMC computational experiments on five problems at d = 12, including a hierarchical gene-expression model, show the Cauchy reducing the sliced Wasserstein distance to the reference posterior by factors of up to 50 with the same simulation budget. Since the summary statistics are commonly insufficient for the models that require ABC, overinflation is structural and the Cauchy perturbation kernel is the right default for problems in higher dimensions.

Matching journals

The top 11 journals account for 50% of the predicted probability mass.

1
PLOS Computational Biology
1863 papers in training set
Top 2%
12.8%
2
Bioinformatics
1204 papers in training set
Top 3%
8.0%
3
PLOS ONE
5266 papers in training set
Top 31%
4.9%
4
The Journal of Chemical Physics
56 papers in training set
Top 0.1%
4.4%
5
BMC Bioinformatics
457 papers in training set
Top 2%
3.5%
6
Nature Communications
5641 papers in training set
Top 35%
3.3%
7
eLife
5828 papers in training set
Top 33%
3.3%
8
Journal of Mathematical Biology
40 papers in training set
Top 0.2%
3.3%
9
Cell Systems
201 papers in training set
Top 1%
3.3%
10
Journal of The Royal Society Interface
235 papers in training set
Top 2%
2.8%
11
Journal of Theoretical Biology
162 papers in training set
Top 1%
2.4%
50% of probability mass above
12
Proceedings of the National Academy of Sciences
2444 papers in training set
Top 21%
2.4%
13
GENETICS
483 papers in training set
Top 2%
2.1%
14
The Annals of Applied Statistics
19 papers in training set
Top 0.1%
2.1%
15
Biophysical Journal
631 papers in training set
Top 3%
2.1%
16
Scientific Reports
3612 papers in training set
Top 50%
1.9%
17
Systematic Biology
144 papers in training set
Top 0.5%
1.9%
18
Physical Review E
112 papers in training set
Top 0.8%
1.8%
19
Bulletin of Mathematical Biology
92 papers in training set
Top 0.8%
1.8%
20
Nature Methods
385 papers in training set
Top 4%
1.7%
21
Nature Genetics
286 papers in training set
Top 3%
1.7%
22
Patterns
78 papers in training set
Top 2%
1.1%
23
Biostatistics
24 papers in training set
Top 0.2%
1.1%
24
Statistics in Medicine
40 papers in training set
Top 0.4%
1.1%
25
Imaging Neuroscience
282 papers in training set
Top 4%
1.0%
26
Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
12 papers in training set
Top 0.1%
1.0%
27
Genome Research
468 papers in training set
Top 6%
0.9%
28
Genome Biology
637 papers in training set
Top 8%
0.9%
29
Frontiers in Systems Biology
10 papers in training set
Top 0.1%
0.9%
30
Communications Biology
993 papers in training set
Top 29%
0.9%