abctools (version 0.2-2)

AS.test: Test for relative approximate sufficiency between two posterior samples.

Description

The function tests to determine adding a (set of) statistics is informative in ABC inference.

Usage

AS.test(grid = 10, x1, x2, supp)

Value

extremea boolean value indicating whether the alternative posterior sample is more informative than the first (i.e. the extra summary statistics add information.

Details

After dividing each posterior sample into a number of bins (specified by grid), the function computes the ratio of the posterior densities. This is seen as a measure of information added (sufficiency) by using the alternative posterior sample instead of the first posterior sample. If the ratio exceeds a particular threshold (a number of standard deviations away from the expected counts in each bin), then the alternative set of summaries is seen as being more informative.

References

Blum, M. G. B, Nunes, M. A., Prangle, D. and Sisson, S. A. (2013) A comparative review of dimension reduction methods in approximate Bayesian computation. Stat. Sci. (to appear). Joyce, P. and P. Marjoram (2008) Approximately sufficient statistics and Bayesian computation. Stat. Appl. Gen. Mol. Biol. 7 Article 26.

See Also

AS.select

Examples

#create two fake posterior samples:

x1<-runif(10000)
x2<-rnorm(10000)

AS.test(x1=x1,x2=x2,supp=range(x2))